A method and related equipment for calculating induced polarization parameters based on interference signals in working areas

Abstract

The present invention provides a method and related equipment for calculating induced polarization parameters based on a work area interference signal, comprising: collecting interference time series data of an interference signal when no spread spectrum signal is sent and spread spectrum time series data fed back when a spread spectrum signal is sent on the ground of a target work area, wherein the spread spectrum signal comprises multiple main frequencies and multiple combined frequencies, and the combined frequency is formed by fusing multiple main frequencies; calculating the spectrum energy distribution of the interference time series data to obtain an interference equivalent spectrum; for each main frequency among the multiple main frequencies, obtaining a weighting coefficient of the main frequency according to the interference equivalent spectrum and updating the number of main frequencies in each combined frequency; obtaining a frequency value of each combined frequency according to the weighting coefficient of each main frequency; obtaining an equivalent spectrum value of each combined frequency according to the spectrum energy distribution of the spread spectrum time series data; calculating the induced polarization parameters of the spread spectrum induced polarization signal according to the frequency value and the equivalent spectrum value; solving the problem of excessive error in induced polarization parameters measured during spread spectrum induced polarization exploration.

Description

Technical Field

[0001] The present invention relates to the technical field of geophysical electrical prospecting, and in particular to a method for calculating induced polarization parameters based on a working area interference signal and related equipment. Background Art

[0002] The spread spectrum signal is a pseudo-random m-sequence signal, which appears as a pseudo-random step waveform containing only high and low levels in the time domain. In the frequency domain, the energy of the entire frequency domain is concentrated on multiple spread spectrum IP signal main frequency points, and the energy of other non-spread spectrum IP signal main frequency points is basically 0. The amplitudes of all spread spectrum IP signal main frequency points are basically the same, and the frequency points are evenly spaced on the linear frequency coordinate. If the order of the spread spectrum wave is m, then the number of spread spectrum IP signal main frequencies contained in the spread spectrum wave is 2 m-1 , the main frequency value of the first spread spectrum IP signal is the fundamental frequency value, and the frequency interval value between the main frequencies of adjacent spread spectrum IP signals is equal to the fundamental frequency value.

[0003] Spread-spectrum IP (SIPP) is an exploration method that uses spread-spectrum signals to detect differences in resistivity and polarizability of underground targets. This method offers advantages such as high frequency resolution, multiple measurement parameters, and strong anti-coupling capabilities. During field operations, to improve anti-interference capabilities, the main frequencies of the spread-spectrum IP signals are typically combined. A fixed number of adjacent main frequencies are averaged to form a single frequency, ultimately yielding three to five combined frequencies. The frequency and spectrum values ​​of these combined frequencies are then used to calculate IP parameters to detect the target.

[0004] The existing spread spectrum IP exploration acquisition and processing method does not take into account the interference conditions of the work area. The main frequency of the spread spectrum IP signal in all work areas is processed in the same way, so that IP parameters with small errors can be obtained in some weak interference areas, but large errors can be obtained in some strong interference areas. Summary of the Invention

[0005] The present invention provides a method and related equipment for calculating induced polarization parameters based on interference signals in a working area, the purpose of which is to solve the problem of large errors in measured induced polarization parameters during spread spectrum induced polarization exploration.

[0006] In order to achieve the above object, the present invention provides a method for calculating induced polarization parameters based on a working area interference signal, comprising:

[0007] Step 1: On the ground of the target work area, collect interference time series data of the interference signal when no spread spectrum signal is transmitted and spread spectrum time series data fed back when the spread spectrum signal is transmitted. The spread spectrum signal includes multiple main frequencies and multiple combined frequencies. The combined frequency is formed by fusing multiple main frequencies.

[0008] Step 2: Calculate the spectrum energy distribution of the interference time series data to obtain the interference equivalent spectrum;

[0009] Step 3: for each main frequency among the multiple main frequencies, obtain a weighting coefficient of the main frequency according to the interference equivalent spectrum;

[0010] Step 4: for each of the multiple combined frequencies, update the number of main frequencies in the combined frequency according to the interference equivalent spectrum;

[0011] Step 5: Obtain the frequency value of each combined frequency according to the weighting coefficients of the multiple main frequencies;

[0012] Step 6: Obtain the equivalent spectrum value of each combination frequency according to the spectrum energy distribution of the spread spectrum time series data, and calculate the IP parameters of the spread spectrum IP signal according to the frequency value and the equivalent spectrum value. The IP parameters include apparent resistivity, relative phase and dispersion rate.

[0013] More specifically, step 1 includes:

[0014] Select N sampling points on the ground within the target work area, and collect interference time series data T of the interference signal when no spread spectrum signal is sent at each sampling point;

[0015] The spread spectrum timing data fed back when the spread spectrum signal is sent is collected at the designed points in the target work area.

[0016] More specifically, step 2 includes:

[0017] For the interference time series data T corresponding to each sampling point, perform Fourier transform on the interference time series data T and normalize it to obtain S:

[0018] S=DFT(T,Nf)×2 / Nf

[0019] Where Nf is the number of sampling points of Fourier transform;

[0020] Calculate the average value of all sampling points at each main frequency point to obtain the interference equivalent spectrum S d .

[0021] More specifically, step 3 includes:

[0022] For each of the multiple main frequencies, obtain the main frequency in the interference equivalent spectrum S d The amplitude value A on

[0023] The weighting coefficient w of the main frequency is calculated according to the amplitude value A. The calculation formula is:

[0024]

[0025] More specifically, step 4 includes:

[0026] For each of the multiple combined frequencies, obtain the number N of main frequencies in the combined frequency according to the interference equivalent spectrum. j ;

[0027] For each of the multiple main frequencies in the combined frequency, calculate the average value AV of the interference equivalent spectrum amplitude of the main frequency j ;

[0028] Comparison of the average value of the interference equivalent spectrum amplitude AV j and the preset threshold Y a size;

[0029] When the average value of the interference equivalent spectrum amplitude AV j Greater than the preset threshold Y a When the number of main frequencies N j Update to N j =3×N c / 2;

[0030] When the average value of the interference equivalent spectrum amplitude AV j Less than the preset threshold Y a When the number of main frequencies N j Update to N j =N c .

[0031] More specifically, step 5 includes:

[0032] According to the number N of main frequencies in each combination frequency j And the weighted coefficient w of each main frequency, calculate the frequency value F of each combination frequency j , the calculation formula is as follows:

[0033]

[0034] Among them, j is the index value of the combination frequency, N j is the number of main frequencies in the j+1th combination frequency, w i is the weighting coefficient of the ith main frequency, f i is the frequency value of the main frequency, P j is the starting index of the main frequency in the j+1th combination frequency, and its calculation formula is as follows:

[0035]

[0036] where N z N is the number of combined frequencies required in the end. c N is the default number of combined frequencies and fused main frequencies. c =2 m-1 / N z , m is the order of the spread spectrum signal.

[0037] More specifically, step 6 includes:

[0038] According to the spectrum energy distribution of the spread spectrum time series data, the number N of main frequencies in each combination frequency j And the weighting coefficient w of each main frequency, calculate the equivalent spectrum value C of each combination frequency j , the calculation formula is as follows:

[0039]

[0040] Among them, P j is the starting index of the main frequency in the j+1th combination frequency, w i is the weighting coefficient of the ith main frequency, X i is the spectrum value of the main frequency of the spread spectrum IP signal.

[0041] The present invention also provides a device for calculating induced polarization parameters based on a working area interference signal, comprising:

[0042] The measurement module is used to collect interference time series data of the interference signal when the spread spectrum signal is not transmitted and the spread spectrum time series data fed back when the spread spectrum signal is transmitted on the ground of the target work area. The spread spectrum signal includes multiple main frequencies and multiple combined frequencies. The combined frequency is formed by fusing multiple main frequencies.

[0043] The first calculation module is used to calculate the spectrum energy distribution of the interference time series data to obtain the interference equivalent spectrum;

[0044] A first acquisition module is configured to acquire, for each main frequency among the multiple main frequencies, a weighting coefficient of the main frequency according to an interference equivalent spectrum;

[0045] An updating module, configured to update, for each of the plurality of combined frequencies, the number of main frequencies in the combined frequency according to an interference equivalent spectrum;

[0046] A second acquisition module is used to obtain the frequency value of each combined frequency according to the weighting coefficients of the multiple main frequencies;

[0047] The second calculation module is used to obtain the equivalent spectrum value of each combination frequency according to the spectrum energy distribution of the spread spectrum time series data, and calculate the IP parameters of the spread spectrum IP signal according to the frequency value and the equivalent spectrum value. The IP parameters include apparent resistivity, relative phase and dispersion rate.

[0048] The present invention also provides a computer-readable storage medium for storing a computer program, and the computer program is executed to implement a method for calculating induced polarization parameters based on a work area interference signal.

[0049] The present invention also provides a device for calculating induced polarization parameters based on a work area interference signal, which is used to implement a method for calculating induced polarization parameters based on a work area interference signal, comprising:

[0050] memory and processor;

[0051] Memory is used to store computer programs;

[0052] The processor is used to execute the computer program stored in the memory.

[0053] The above solution of the present invention has the following beneficial effects:

[0054] The present invention collects interference time series data of an interference signal when no spread spectrum signal is sent and spread spectrum time series data fed back when a spread spectrum signal is sent on the ground of a target work area, wherein the spread spectrum signal includes multiple main frequencies and multiple combined frequencies, and the combined frequency is formed by fusing multiple main frequencies; calculates the spectrum energy distribution of the interference time series data to obtain an interference equivalent spectrum, thereby pre-estimating the reliability of the measurement data; obtains a weighting coefficient of the main frequency based on the interference equivalent spectrum and updates the number of main frequencies in the combined frequency, thereby enhancing the signal-to-noise ratio of the combined frequency; obtains a frequency value and an equivalent spectrum value of the combined frequency based on the weighting coefficient of the main frequency and the spectrum energy distribution of the spread spectrum time series data, and calculates the induced polarization parameters of the spread spectrum induced polarization signal based on the frequency value and the equivalent spectrum value; compared with the prior art, the weighting coefficient of each main frequency is calculated based on the amplitude distribution characteristics of the work area interference at the main frequency point of the spread spectrum signal, thereby minimizing the influence of the interference on the final result; when the interference intensity of the main frequencies in the combined frequency is large, the signal-to-noise ratio of the combined frequency is improved by increasing the number of updated main frequencies, thereby reducing the error of the induced polarization parameters.

[0055] Other beneficial effects of the present invention will be described in detail in the subsequent specific implementation section. BRIEF DESCRIPTION OF THE DRAWINGS

[0056] Figure 1 Schematic diagram of a flow chart of an embodiment of the present invention;

[0057] Figure 2 Schematic diagram of measuring background field interference signals in a work area according to an embodiment of the present invention. DETAILED DESCRIPTION

[0058] To make the technical problems, technical solutions, and advantages to be solved by the present invention more clear, the following is a detailed description with reference to the accompanying drawings and specific embodiments. It is obvious that the embodiments described are only some of the embodiments of the present invention, not all of them. All other embodiments derived by persons of ordinary skill in the art based on the embodiments of the present invention without inventive effort are within the scope of protection of the present invention.

[0059] In the description of the present invention, it should be noted that the terms "center," "upper," "lower," "left," "right," "vertical," "horizontal," "inner," and "outer," etc., indicating orientations or positional relationships, are based on the orientations or positional relationships shown in the accompanying drawings and are intended solely to facilitate and simplify the description of the present invention. They are not intended to indicate or imply that the devices or components referred to must have, be constructed, or operate in a specific orientation, and therefore should not be construed as limitations on the present invention. Furthermore, the terms "first," "second," and "third" are used for descriptive purposes only and should not be construed as indicating or implying relative importance.

[0060] In the description of the present invention, it should be noted that, unless otherwise expressly specified or limited, the terms "mounted," "connected," and "connected" should be understood broadly. For example, they may refer to a locking connection, a detachable connection, or an integral connection; they may refer to a mechanical connection or an electrical connection; they may refer to a direct connection or an indirect connection through an intermediate medium; and they may refer to internal communication between two components. Those skilled in the art will understand the specific meanings of the above terms in the present invention based on the specific circumstances.

[0061] In addition, the technical features involved in the different embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0062] In view of the existing problems, the present invention provides a method for calculating induced polarization parameters based on a working area interference signal and related equipment.

[0063] like Figure 1 As shown, an embodiment of the present invention provides a method for calculating induced polarization parameters based on a work area interference signal, comprising:

[0064] Step 1: On the ground of the target work area, collect interference time series data of the interference signal when no spread spectrum signal is transmitted and spread spectrum time series data fed back when the spread spectrum signal is transmitted. The spread spectrum signal includes multiple main frequencies and multiple combined frequencies. The combined frequency is formed by fusing multiple main frequencies.

[0065] Step 2: Calculate the spectrum energy distribution of the interference time series data to obtain the interference equivalent spectrum;

[0066] Step 3: for each main frequency among the multiple main frequencies, obtain a weighting coefficient of the main frequency according to the interference equivalent spectrum;

[0067] Step 4: for each of the multiple combined frequencies, update the number of main frequencies in the combined frequency according to the interference equivalent spectrum;

[0068] Step 5: Obtain the frequency value of each combined frequency according to the weighting coefficients of the multiple main frequencies;

[0069] Step 6: Obtain the equivalent spectrum value of each combination frequency according to the spectrum energy distribution of the spread spectrum time series data, and calculate the IP parameters of the spread spectrum IP signal according to the frequency value and the equivalent spectrum value. The IP parameters include apparent resistivity, relative phase and dispersion rate.

[0070] Specifically, step 1 includes:

[0071] Select N sampling points in the target work area, and collect a section of interference time series data T of the interference signal when no spread spectrum signal is sent at each sampling point;

[0072] The spread spectrum timing data fed back when the spread spectrum signal is sent is collected at the designed points in the target work area.

[0073] like Figure 2 As shown in the figure, the outer frame is the construction area, and the black dots are the sampling points for collecting background field interference signals. Figure 1 Starting from the northwest edge, a 1km x 1km grid is established. Interference signal collection points are located at grid intersections within the construction area (excluding the edge). Designed work area points are generally based on exploration specifications or actual needs. They are typically distributed throughout the construction area at regular intervals, with typical point and line spacing ranging from 10m to 100m and 20m to 200m.

[0074] In the embodiment of the present invention, for the sake of convenience, it is assumed that the spread spectrum signal designed in the target work area is a 5th order spread spectrum signal with a base frequency of 1 / 16 Hz, the sampling rate of the acquisition device is 15 Hz, and the final combination frequency to be obtained is 4. After the sampling point is determined, the interference time series data T is generally collected at the sampling point for about 10 minutes. i , T i is a time series array, i is the i-th sampling point.

[0075] Specifically, in this embodiment of the present invention, step 2 includes:

[0076] According to the designed 5th-order spread spectrum signal with a fundamental frequency of 1 / 16 Hz, the interference time series data T corresponding to each sampling point is Fourier transformed and normalized to obtain S;

[0077] S=DFT(T,Nf)×2 / Nf

[0078] Among them, Nf is the number of sampling points of Fourier transform. Nf is generally the maximum number of full-cycle sampling points that can be taken within the sampling time. The maximum number of spread spectrum cycles contained in 10 minutes is T. n =37, a total of 592 sampling points, that is, Nf = 592; then calculate the average value of all sampling points at each main frequency point to obtain the interference equivalent spectrum S of the entire work area d , interference equivalent spectrum S dis an array of complex numbers of size 592.

[0079] Specifically, in the embodiment of the present invention, the interference equivalent spectrum S of the working area is obtained. d Then, the weighting coefficient of the main frequency in the spread spectrum signal is calculated by the following steps:

[0080] For each of the multiple main frequencies, obtain the main frequency in the interference equivalent spectrum S d The amplitude value A on i ; The amplitude value corresponding to the i+1th main frequency is A i , which is in the interference equivalent spectrum S d The index in is (i+1)×T n , i starts counting from 0;

[0081] According to the amplitude value A i Calculate the weighting coefficient w of the main frequency i , the calculation formula is:

[0082]

[0083] w i is the weighting coefficient corresponding to the i+1th main frequency, where i starts counting from 0.

[0084] Specifically, step 4, updating the data of the main frequency in each combined frequency, is intended to enable the combined frequency to obtain a better signal-to-noise ratio, including:

[0085] For each of the multiple combined frequencies, obtain the number N of main frequencies in the combined frequency according to the interference equivalent spectrum. j ; Generally, each combination frequency needs to integrate the default number of main frequencies N j Set to N c =2 m-1 / N z , where m is the order of the spread spectrum signal, N z The number of combined frequencies required in the end. According to the parameters assumed in the embodiment of the present invention, N c =2 5 -1 / 4=4.

[0086] For each of the multiple main frequencies in the combined frequency, calculate the average value AV of the interference equivalent spectrum amplitude of the main frequency j ; The starting index of the j+1th combined frequency fusion main frequency is N c ×j, according to the parameters assumed in the embodiment of the present invention, the average amplitude of the interference equivalent spectrum of the first main frequency is

[0087] Comparison of the average value of the interference equivalent spectrum amplitude AV j and the preset threshold Y asize;

[0088] When the average value of the interference equivalent spectrum amplitude AV j Greater than the preset threshold Y a When the number of main frequencies N j Update to N j =3×N c / 2;

[0089] When the average value of the interference equivalent spectrum amplitude AV j Less than the preset threshold Y a When the number of main frequencies N j Update to N j =N c , generally the background field is greater than 0.5mV, which can be considered as a large interference, so the Y a Set to 5×10 -4 V.

[0090] Specifically, calculate the weighted coefficient w of each main frequency i and the number of main frequencies N in each combination frequency j After that, the number of main frequencies in each combination frequency N can be j And the weighted coefficient w of each main frequency, calculate the frequency value of each combination frequency, the frequency value F of the j+1th combination frequency j The calculation formula is as follows:

[0091]

[0092] Among them, j is the index value of the combination frequency, N j is the number of main frequencies in the j+1th combination frequency, w i is the weighting coefficient of the ith main frequency, f i is the frequency value of the main frequency, P j is the starting index of the main frequency in the j+1th combination frequency, and its calculation formula is as follows:

[0093]

[0094] where N z N is the number of combined frequencies required in the end. c N is the default number of combined frequencies and fused main frequencies. c =2 m-1 / N z , m is the order of the spread spectrum signal.

[0095] Specifically, step 6 includes:

[0096] According to the spectrum energy distribution of the spread spectrum time series data, the number N of main frequencies in each combination frequency jAnd the weighting coefficient w of each main frequency, calculate the equivalent spectrum value C of each combination frequency j , the calculation formula is as follows:

[0097]

[0098] Among them, P j is the starting index of the main frequency in the j+1th combination frequency, w i is the weighting coefficient of the ith main frequency, X i is the spectrum value of the main frequency of the spread spectrum IP signal.

[0099] Finally, based on the spectral energy distribution of the spread spectrum time series data, the equivalent spectrum value of each combination frequency is obtained, and the induced polarization parameters of the spread spectrum induced polarization signal are calculated based on the frequency value and the equivalent spectrum value. The induced polarization parameters include apparent resistivity, relative phase and dispersion rate. The calculation method of apparent resistivity, relative phase and dispersion rate is a well-known technology in the field, so it will not be repeated in the embodiments of the present invention.

[0100] The embodiment of the present invention collects interference time series data of an interference signal when no spread spectrum signal is transmitted and spread spectrum time series data fed back when a spread spectrum signal is transmitted on the ground of a target work area. The spread spectrum signal includes multiple main frequencies and multiple combined frequencies, and the combined frequencies are formed by fusing multiple main frequencies. The spectral energy distribution of the interference time series data is calculated to obtain an interference equivalent spectrum, thereby pre-estimating the reliability of the measurement data. Based on the interference equivalent spectrum, a weighting coefficient of each main frequency is obtained and the number of main frequencies in each combined frequency is updated, thereby enhancing the signal-to-noise ratio of the combined frequency. Based on the weighting coefficient of the main frequency and the spectral energy distribution of the spread spectrum time series data, the frequency value and equivalent spectrum value of the combined frequency are obtained, and the induced polarization parameters of the spread spectrum induced polarization signal are calculated based on the frequency value and equivalent spectrum value. Compared with the prior art, the weighting coefficient of each main frequency is calculated based on the amplitude distribution characteristics of the work area interference at the main frequency point of the spread spectrum signal, thereby minimizing the impact of the interference on the final result. When the interference intensity of the main frequencies in the combined frequency is large, the signal-to-noise ratio of the combined frequency is improved by increasing the number of updated main frequencies, thereby reducing the error of the induced polarization parameters.

[0101] An embodiment of the present invention further provides a device for calculating induced polarization parameters based on a work area interference signal, comprising:

[0102] The measurement module is used to collect interference time series data of the interference signal when the spread spectrum signal is not transmitted and the spread spectrum time series data fed back when the spread spectrum signal is transmitted on the ground of the target work area. The spread spectrum signal includes multiple main frequencies and multiple combined frequencies. The combined frequency is formed by fusing multiple main frequencies.

[0103] The first calculation module is used to calculate the spectrum energy distribution of the interference time series data to obtain the interference equivalent spectrum;

[0104] A first acquisition module is configured to acquire, for each main frequency among the multiple main frequencies, a weighting coefficient of the main frequency according to an interference equivalent spectrum;

[0105] An updating module, configured to update, for each of the plurality of combined frequencies, the number of main frequencies in the combined frequency according to an interference equivalent spectrum;

[0106] A second acquisition module is used to obtain the frequency value of each combined frequency according to the weighting coefficients of the multiple main frequencies;

[0107] The second calculation module is used to obtain the equivalent spectrum value of each combination frequency based on the spectrum energy distribution of the spread spectrum time series data, and calculate the IP parameters of the spread spectrum IP signal based on the frequency value and the equivalent spectrum value. The IP parameters include apparent resistivity, relative phase and dispersion rate.

[0108] It should be noted that the information interaction, execution process, etc. between the above-mentioned devices / units are based on the same concept as the method embodiments of the embodiments of the present invention. Their specific functions and technical effects can be found in the method embodiment section and will not be repeated here.

[0109] Those skilled in the art can clearly understand that, for the convenience and brevity of description, only the division of the above-mentioned functional units and modules is used as an example for illustration. In actual applications, the above-mentioned functions can be distributed and completed by different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiment can be integrated into one processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit. The above-mentioned integrated unit can be implemented in the form of hardware or in the form of software functional units. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing each other, and are not used to limit the scope of protection of the embodiments of the present invention. The specific working process of the units and modules in the above-mentioned system can refer to the corresponding process in the aforementioned method embodiment, and will not be repeated here.

[0110] An embodiment of the present invention further provides a computer-readable storage medium for storing a computer program, and a method for calculating induced polarization parameters based on a work area interference signal is implemented by executing the computer program.

[0111] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the embodiments of the present invention implement all or part of the processes in the above-mentioned embodiment method, which can be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by a processor, it can implement the steps of the above-mentioned various method embodiments. Among them, the computer program includes computer program code, which can be in source code form, object code form, executable file or some intermediate form. The computer-readable medium can at least include: any entity or device capable of carrying computer program code to a construction device / terminal device, a recording medium, a computer memory, a read-only memory (ROM), a random access memory (RAM), an electric carrier signal, a telecommunication signal and a software distribution medium. For example, a USB flash drive, a mobile hard disk, a magnetic disk or an optical disk. In some jurisdictions, according to legislation and patent practice, computer-readable media cannot be electric carrier signals and telecommunication signals.

[0112] The present invention also provides a device for calculating induced polarization parameters based on a work area interference signal, which is used to implement a method for calculating induced polarization parameters based on a work area interference signal, comprising:

[0113] memory and processor;

[0114] Memory is used to store computer programs;

[0115] The processor is used to execute the computer program stored in the memory.

[0116] It should be noted that the terminal device may be a mobile phone, tablet computer, laptop computer, ultra-mobile personal computer (UMPC), netbook, personal digital assistant (PDA), and other terminal devices. For example, the terminal device may be a station (ST, STAION) in a WLAN, a cellular phone, a cordless phone, a Session Initiation Protocol (SIP) phone, a wireless local loop (WLL) station, a personal digital assistant (PDA), a handheld device with wireless communication capabilities, a computing device or other processing device connected to a wireless modem, a computer, a laptop computer, a handheld communication device, a handheld computing device, a satellite wireless device, etc. The embodiments of the present invention do not impose any restrictions on the specific type of the terminal device.

[0117] The processor may be a central processing unit (CPU), other general-purpose processors, a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor.

[0118] In some embodiments, the memory may be an internal storage unit of the terminal device, such as a hard disk or memory of the terminal device. In other embodiments, the memory may also be an external storage device of the terminal device, such as a plug-in hard disk equipped on the terminal device, a smart memory card (SMC, Smart Media Card), a secure digital (SD, Secure Digital) card, a flash card, etc. Furthermore, the memory may include both an internal storage unit of the terminal device and an external storage device. The memory is used to store an operating system, an application program, a boot loader (BootLoader), data, and other programs, such as the program code of the computer program. The memory may also be used to temporarily store data that has been output or is to be output.

[0119] It should be noted that the information interaction, execution process, etc. between the above-mentioned devices / units are based on the same concept as the method embodiments of the embodiments of the present invention. Their specific functions and technical effects can be found in the method embodiment section and will not be repeated here.

[0120] The above is a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principles of the present invention. These improvements and modifications should also be regarded as within the scope of protection of the present invention.

Claims

1. A method for calculating IP parameters based on a work area interference signal, characterized in that: include: Step 1: On the ground of the target work area, collect interference time series data of the interference signal when no spread spectrum signal is transmitted and spread spectrum time series data fed back when the spread spectrum signal is transmitted, wherein the spread spectrum signal includes multiple main frequencies and multiple combined frequencies, and the combined frequency is formed by fusing multiple main frequencies; Step 2: Calculate the spectrum energy distribution of the interference time series data to obtain an interference equivalent spectrum; Step 3, for each of the plurality of main frequencies, obtaining a weighting coefficient of the main frequency according to the interference equivalent spectrum, includes: For each of the plurality of main frequencies, obtain the interference equivalent spectrum S of the main frequency d The amplitude value A on According to the amplitude value A, the weighting coefficient w of the main frequency is calculated using the following formula: Step 4: for each of the plurality of combined frequencies, updating the number of main frequencies in the combined frequency according to the interference equivalent spectrum; Step 5, obtaining the frequency value of each combined frequency according to the weighting coefficients of the plurality of main frequencies; Step 6: Obtain the equivalent spectrum value of each combination frequency according to the spectrum energy distribution of the spread spectrum time series data, and calculate the IP parameters of the spread spectrum IP signal according to the frequency value and the equivalent spectrum value. The IP parameters include apparent resistivity, relative phase and dispersion rate.

2. The method for calculating induced polarization parameters based on a working area interference signal according to claim 1, characterized in that: The step 1 comprises: Selecting N sampling points on the ground within the target work area, and collecting interference time series data T of the interference signal when no spread spectrum signal is transmitted at each sampling point; The spread spectrum timing data fed back when the spread spectrum signal is sent is collected at the designed points within the target work area.

3. The method for calculating induced polarization parameters based on work area interference signals according to claim 2, characterized in that: The step 2 includes: For the interference time series data T corresponding to each sampling point, the interference time series data T is Fourier transformed and normalized to obtain S: S=DFT(T,Nf)×2 / Nf Where Nf is the number of sampling points of Fourier transform; Calculate the average value of all sampling points at each main frequency point to obtain the interference equivalent spectrum S d .

4. The method for calculating IP parameters based on work area interference signals according to claim 3, characterized in that: The step 4 comprises: For each of the plurality of combined frequencies, the number N of main frequencies in the combined frequency is obtained according to the interference equivalent spectrum. j ; For each of the multiple main frequencies in the combined frequency, calculate the interference equivalent spectrum amplitude average value AV of the main frequency j ; Compare the average value AV of the interference equivalent spectrum amplitude j and the preset threshold Y a size; When the average value of the interference equivalent spectrum amplitude AV j Greater than the preset threshold Y a When the number of the main frequency N j Update to N j =3×N c / 2; When the average value of the interference equivalent spectrum amplitude AV j Less than the preset threshold Y a When the number of the main frequency N j Update to N j =N c , N c The default number of main frequencies for the combined frequency fusion.

5. The method for calculating IP parameters based on work area interference signals according to claim 4, characterized in that: The step 5 comprises: According to the number N of main frequencies in each of the combined frequencies j And the weighted coefficient w of each main frequency, calculate the frequency value F of each combination frequency j , the calculation formula is as follows: Among them, j is the index value of the combination frequency, N j is the number of main frequencies in the j+1th combination frequency, w i is the weighting coefficient of the ith main frequency, f i is the frequency value of the main frequency, P j is the starting index of the main frequency in the j+1th combination frequency, and its calculation formula is as follows: where N z N is the number of combined frequencies required in the end. c N is the default number of combined frequencies and fused main frequencies. c =2 m-1 / N z , m is the order of the spread spectrum signal.

6. The method for calculating induced polarization parameters based on a working area interference signal according to claim 5, characterized in that: The step 6 comprises: According to the spectrum energy distribution of the spread spectrum time series data, the number N of the main frequencies in each of the combined frequencies j And the weighting coefficient w of each main frequency, calculate the equivalent spectrum value C of each combination frequency j , the calculation formula is as follows: Among them, P j is the starting index of the main frequency in the j+1th combination frequency, w i is the weighting coefficient of the ith main frequency, X i is the spectrum value of the i-th main frequency in the combined frequency.

7. A device for calculating induced polarization parameters based on interference signals in a working area, characterized in that: include: An acquisition module is configured to collect, on the ground of the target work area, interference time series data of an interference signal when no spread spectrum signal is transmitted and spread spectrum time series data fed back when a spread spectrum signal is transmitted, wherein the spread spectrum signal includes multiple main frequencies and multiple combined frequencies, and the combined frequency is formed by fusing the multiple main frequencies; A first calculation module is used to calculate the spectrum energy distribution of the interference time series data to obtain an interference equivalent spectrum; The first acquisition module is configured to acquire, for each of the plurality of main frequencies, a weighting coefficient of the main frequency according to the interference equivalent spectrum, including: For each of the plurality of main frequencies, obtain the interference equivalent spectrum S of the main frequency d The amplitude value A on According to the amplitude value A, the weighting coefficient w of the main frequency is calculated using the following formula: An updating module, configured to update, for each of the plurality of combined frequencies, the number of main frequencies in the combined frequency according to the interference equivalent spectrum; A second acquisition module is used to obtain the frequency value of each combined frequency according to the weighting coefficients of the plurality of main frequencies; The second calculation module is used to obtain the equivalent spectrum value of each combination frequency according to the spectrum energy distribution of the spread spectrum time series data, and calculate the induced polarization parameters of the spread spectrum induced polarization signal according to the frequency value and the equivalent spectrum value, wherein the induced polarization parameters include apparent resistivity, relative phase and dispersion rate.

8. A computer-readable storage medium for storing a computer program, characterized in that: By executing the computer program, the method for calculating induced polarization parameters based on a working area interference signal as described in any one of claims 1 to 6 is implemented.

9. A device for calculating induced polarization parameters based on a work area interference signal, used to implement the method for calculating induced polarization parameters based on a work area interference signal according to any one of claims 1 to 6, characterized in that: include: memory and processor; The memory is used to store computer programs; The processor is configured to execute the computer program stored in the memory.

Images