A method, device, medium and equipment for obtaining induced polarization parameters using spread spectrum signals

Abstract

The present invention provides a method, device, medium and equipment for obtaining induced polarization parameters using spread spectrum signals, including: sending a spread spectrum signal to a detection target and obtaining time series data of the spread spectrum signal; obtaining the spectrum value of each main frequency in the spread spectrum signal according to the time series data, and obtaining the normalized complex value of each main frequency based on the spectrum value of each main frequency; obtaining the frequency value of each combined frequency in the spread spectrum signal according to the frequency value of each main frequency; calculating the normalized complex value of each combined frequency according to the normalized complex value of each main frequency; calculating the induced polarization parameters of all combined frequencies according to the frequency values ​​of all combined frequencies and the normalized complex values ​​of all combined frequencies, the induced polarization parameters including apparent resistivity, dispersion rate and relative phase; realizing the function of calculating the induced polarization parameters of the spread spectrum signal, solving the problem of distortion of induced polarization parameters measured by conventional induced polarization methods due to low frequency domain resolution, weak anti-interference ability and large coupling influence.

Description

Technical Field

[0001] The present invention relates to the technical field of geophysical electrical prospecting, and in particular to a method, device, medium and equipment for obtaining induced polarization parameters by using spread spectrum signals. Background Art

[0002] Induced polarization (IP) is a geophysical exploration method that solves geological problems by studying the differences in resistivity and polarizability of rocks and minerals under the influence of electric fields in the time or frequency domain. Compared to other electrical and electromagnetic methods, IP can more easily distinguish mineral-induced anomalies in solid mineral exploration through polarization characteristics, such as anomalies caused by dense massive lead-zinc ores and disseminated metal sulfides. Conventional frequency-domain IP generally uses square wave, dual-frequency, or pseudo-random multi-frequency wave signals. The differences between the different transmitted signals are mainly reflected in the spectral energy distribution, spectral density, and calculation parameters. Among them, the energy of the square wave signal is mainly distributed on the odd harmonics, and the energy decreases rapidly with the increase of the harmonic number. Generally, during induced polarization acquisition, only the data of the first 2 to 3 frequency points can be used to calculate the apparent resistivity, phase and dispersion rate; the dual-frequency wave contains 2 main frequencies, and the interval between the two frequencies is large, and 2 apparent resistances, 2 phases and 1 dispersion rate parameters can be calculated; the multi-frequency wave generally contains 3 to 7 frequencies, and the interval between the frequencies is large, and the apparent resistivity, phase and dispersion rate parameters can be calculated.

[0003] The above methods mainly have the following defects:

[0004] (1) The number of transmitted signal frequencies is small and the frequency intervals are large, which makes it impossible to accurately depict the variation patterns of the IP parameters of rocks and minerals on the spectrum;

[0005] (2) The energy of the main frequency of the transmitted signal is uneven, and the signal-to-noise ratio of some frequencies is very low;

[0006] (3) The frequency interval is large, the energy is dispersed, and reliable data cannot be obtained in areas with strong interference;

[0007] (4) The calculated dispersion rate and phase parameters are greatly affected by coupling, resulting in distortion of the measurement data. Summary of the Invention

[0008] The present invention provides a method, device, medium and equipment for obtaining induced polarization parameters using spread spectrum signals. The purpose is to obtain induced polarization parameters using spread spectrum signals to solve the problem of distortion of induced polarization parameters measured by conventional induced polarization methods due to low frequency domain resolution, weak anti-interference ability and large coupling influence.

[0009] In order to achieve the above object, the present invention provides a method for obtaining IP parameters using a spread spectrum signal, comprising:

[0010] Step 1: Send a spread spectrum signal to the detection target and obtain time series data of the spread spectrum signal, the time series data including voltage time series data and current time series data;

[0011] Step 2: Obtain the spectrum value of each main frequency in the spread spectrum signal according to the time series data, and obtain the normalized complex value of each main frequency based on the spectrum value of each main frequency;

[0012] Step 3, obtaining the frequency value of each combination frequency in the spread spectrum signal according to the frequency value of each main frequency;

[0013] Step 4, calculating the normalized complex value of each combination frequency based on the normalized complex value of each main frequency;

[0014] Step 5: Calculate the IP parameters of all the combined frequencies according to the frequency values ​​of all the combined frequencies and the normalized complex values ​​of all the combined frequencies. The IP parameters include apparent resistivity, dispersion rate, and relative phase.

[0015] More specifically, step 2 includes:

[0016] Step 21, obtaining the fundamental frequency f0 and order N of the spread spectrum signal according to the time series data;

[0017] Step 22: Calculate the number of main frequencies Nf=2 according to the order N of the spread spectrum signal. N-1 ;

[0018] Step 23: Calculate the frequency corresponding to each main frequency as f according to the basic frequency f0 of the spread spectrum signal. i =f0*i, i is the main frequency number from low frequency to high frequency;

[0019] Step 24, perform Fourier transform on the time series data to obtain a transform result, and read the kth value in the transform result as the spectrum value Z of each main frequency, where k is the frequency point index of each main frequency.

[0020] Furthermore, the normalized complex value of each main frequency is obtained based on the spectrum value of each main frequency, including:

[0021] The spectrum value U of each main frequency in the voltage time series data is obtained by using the calculation method of the spectrum value of each main frequency in the spread spectrum signal;

[0022] The spectrum value I of each main frequency in the current time series data is obtained by using the calculation method of the spectrum value of each main frequency in the spread spectrum signal;

[0023] According to the spectrum value U of each main frequency in the voltage time series data and the spectrum value I of each main frequency in the current time series data, the normalized complex value of each main frequency is calculated as X=U / I.

[0024] Furthermore, the frequency value F of each combination frequency j The calculation formula is:

[0025]

[0026] Among them, p is the starting index of the starting main frequency, the value of p is (j-1)*L, the value range of L is [1,n], n is the number of combined frequencies, w i is the frequency weighted value, w i The value of f is 1 / n, i The frequency of the main frequency.

[0027] Furthermore, the normalized complex value C of each combination frequency j The calculation formula is:

[0028]

[0029] Among them, p is the starting index of the starting main frequency, the value of p is (j-1)*L, the value range of L is [1,n], n is the number of combined frequencies, w i is the frequency weighted value, w i The value of X is 1 / n, i is the normalized complex value of the ith main frequency.

[0030] Furthermore, the resistivity ρ j The calculation formula is:

[0031] ρ j =k*|C j |

[0032] Where k is the device coefficient calculated from the receiving point location and the power supply point location, C j is the normalized complex value of the j-th combination frequency;

[0033] Dispersion rate fs j The calculation formula is:

[0034] fs j =(|C j |-|C h |) / |C j |

[0035] or

[0036] fs j =(|C j |-|C h |) / |C h |

[0037] Among them, C j is the normalized complex value of the jth combination frequency, C h is the h-th combination frequency normalized complex value, h = j + 1;

[0038] Relative phase cp j The calculation formula is:

[0039]

[0040] k=F h / F j

[0041]

[0042] in, is the normalized complex-valued phase of the jth combined frequency, is the phase of the hth frequency-normalized complex value, F h is the frequency value of the hth combination frequency, F j is the frequency value of the jth combination frequency.

[0043] The present invention also provides a device for obtaining IP parameters using a spread spectrum signal, comprising:

[0044] A first acquisition module is used to send a spread spectrum signal to the detection target and acquire time series data of the spread spectrum signal, where the time series data includes voltage time series data and current time series data;

[0045] A second acquisition module is used to obtain the spectrum value of each main frequency in the spread spectrum signal according to the time series data, and obtain the normalized complex value of each main frequency based on the spectrum value of each main frequency;

[0046] A first calculation module is used to obtain the frequency value of each combination frequency in the spread spectrum signal according to the frequency value of each main frequency;

[0047] A second calculation module, configured to calculate a normalized complex value of each combination frequency based on the normalized complex value of each main frequency;

[0048] The third calculation module is used to calculate the induced polarization parameters of all the combined frequencies according to the frequency values ​​and normalized complex values ​​of all the combined frequencies, wherein the induced polarization parameters include apparent resistivity, dispersion rate and relative phase.

[0049] The present invention also provides a computer-readable storage medium for storing a computer program, and the computer program is executed to implement the above-mentioned method of obtaining induced polarization parameters using a spread spectrum signal.

[0050] The present invention also provides a device for obtaining IP parameters using a spread spectrum signal, which is used to implement the above-mentioned calculation method, including:

[0051] memory and processor;

[0052] Memory is used to store computer programs;

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

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

[0055] The present invention sends a spread spectrum signal to the detection target and obtains the time series data of the spread spectrum signal, obtains the spectrum value of each main frequency in the spread spectrum signal according to the time series data, and obtains the normalized complex value of each main frequency based on the spectrum value of each main frequency; obtains the frequency value of each combination frequency in the spread spectrum signal according to the frequency value of each main frequency; calculates the normalized complex value of each combination frequency according to the normalized complex value of each main frequency; calculates the induced polarization parameters of all combination frequencies according to the frequency values ​​and normalized complex values ​​of all combination frequencies, and the induced polarization parameters include apparent resistivity, dispersion rate and relative phase; realizes the function of obtaining induced polarization parameters using spread spectrum signals, and solves the problem of distortion of induced polarization parameters measured by conventional induced polarization methods due to low frequency domain resolution, weak anti-interference ability and large coupling influence; at the same time, according to exploration needs, the frequency density and number of combined frequencies can also be customized to achieve fine characterization of the induced polarization spectrum characteristics of rocks and minerals, enhance the signal energy of a single main frequency point, and improve anti-interference ability.

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

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

[0058] Figure 2 is a timing diagram of a spread spectrum signal in an embodiment of the present invention;

[0059] Figure 3 Schematic diagram of decimation and combination of four combined frequencies calculated from a 5th-order spread spectrum signal according to an embodiment of the present invention;

[0060] Figure 4 Schematic diagram of decimation combination of five combination frequencies calculated from a five-order spread spectrum signal according to an embodiment of the present invention;

[0061] Figure 5 It is a pseudo-section diagram of a cross section calculated and drawn using an embodiment of the present invention in field measured data. DETAILED DESCRIPTION

[0062] 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.

[0063] 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.

[0064] 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.

[0065] 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.

[0066] In view of the existing problems, the present invention provides a method, device, medium and equipment for obtaining induced polarization parameters by using spread spectrum signals.

[0067] like Figure 1 As shown, an embodiment of the present invention provides a method for obtaining IP parameters using a spread spectrum signal, comprising:

[0068] Step 1: Send a spread spectrum signal to a detection target and obtain time series data of the spread spectrum signal, wherein the time series data includes voltage time series data and current time series data;

[0069] Step 2: obtaining a spectrum value of each main frequency in the spread spectrum signal according to the time series data, and obtaining a normalized complex value of each main frequency based on the spectrum value of each main frequency;

[0070] Step 3, obtaining the frequency value of each combined frequency in the spread spectrum signal according to the frequency value of each main frequency;

[0071] Step 4, calculating the normalized complex value of each combination frequency according to the normalized complex value of each main frequency;

[0072] Step 5: Calculate the IP parameters of all the combined frequencies according to the frequency values ​​of all the combined frequencies and the normalized complex values ​​of all the combined frequencies. The IP parameters include apparent resistivity, dispersion rate and relative phase.

[0073] Specifically, step 2 includes:

[0074] Step 21, obtaining the fundamental frequency f0 and order N of the spread spectrum signal according to the time series data;

[0075] Step 22: Calculate the number of main frequencies Nf=2 according to the order N of the spread spectrum signal. N-1 ;

[0076] Step 23, based on the basic frequency f0 of the spread spectrum signal, calculate the frequency corresponding to each main frequency as f i =f0*i, i is the main frequency number from low frequency to high frequency;

[0077] Step 23: Perform Fourier transform on the time series data to obtain a transform result, and read the kth value in the transform result as the spectrum value Z of each main frequency, where k is the frequency point index of each main frequency.

[0078] The embodiment of the present invention takes a 5th order spread spectrum signal with a fundamental frequency of 0.0625 Hz as a specific embodiment to further illustrate the present invention.

[0079] Specifically, acquiring time series data including voltage time series data and current time series data of the spread spectrum signal;

[0080] According to the time series data, the basic frequency is obtained as f0=0.0625 Hz, the order of the spread spectrum IP signal is 5, and the time series sampling rate is 64 Hz.

[0081] Specifically, according to the formula Nf=2 N-1 , the number of main frequencies is calculated to be Nf=16.

[0082] According to the formula f i =f0*i, i is the main frequency number from low frequency to high frequency, calculate the main frequency f0~f 15 They are 0.0625Hz, 0.125Hz, 0.1875Hz, 0.25Hz, 0.3125Hz, 0.375Hz, 0.4375Hz, 0.5Hz, 0.5625Hz, 0.625Hz, 0.6875Hz, 0.75Hz, 0.8125Hz, 0.875Hz, 0.9375Hz, and 1Hz respectively.

[0083] Perform two cycles of Fourier transform on the voltage time series data and the current time series data to obtain the voltage Fourier transform result U and the current Fourier transform result I, respectively. Both U and I are arrays containing 2048 complex values, and two cycles means 2048 sample points.

[0084] The complex results of the voltage time series data and the current time series data are obtained, which are 16 data respectively, and the corresponding indexes in U and I are 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 0, and 32.

[0085] Specifically, the normalized complex value of each main frequency is obtained based on the spectrum value of each main frequency, including:

[0086] Obtaining a spectrum value U of each main frequency in the voltage time series data by using a calculation method for the spectrum value of each main frequency in the spread spectrum signal;

[0087] Obtaining the spectrum value I of each main frequency in the current time series data by using a calculation method for the spectrum value of each main frequency in the spread spectrum signal;

[0088] According to the spectrum value U of each main frequency in the voltage time series data and the spectrum value I of each main frequency in the current time series data, the normalized complex value of each main frequency is calculated as X=U / I; according to the calculation formula X i =U i / I i , = 1 to 16, calculate the normalized complex value X of 16 main frequencies respectively i .

[0089] Specifically, calculate the frequency value F of the four combined frequencies j , the calculation formula is as follows:

[0090]

[0091] Among them, p is the starting index of the starting main frequency, the value of p is (j-1)*L, the value range of L is [1,n], n is the number of combined frequencies, w i is the frequency weighted value, w i The value of f is 1 / n, i The larger the L value is, the finer the spectrum characteristics of the rock and ore will be described by the IP parameters. The n value is determined according to the interference situation. When the interference is large, the n value needs to be increased. When the interference at the corresponding main frequency point is large, the w value of the main frequency can be reduced. i , and increase the w of other main frequencies i , making

[0092] In the embodiment of the present invention, n=4, p=4*j, w i =0.25, then the four combined frequencies are F0=0.15625Hz, F1=0.8125Hz, F2=1.96875Hz, and F3=3.625Hz.

[0093] Specifically, the normalized complex value C of each combination frequency jThe calculation formula is:

[0094]

[0095] Among them, p is the starting index of the starting main frequency, the value of p is (j-1)*L, the value range of L is [1,n], n is the number of combined frequencies, w i is the frequency weighted value, w i The value of L is 1 / n. The larger the value of L, the finer the description of the spectrum characteristics of the rock and ore by the IP parameters; the value of n is determined according to the interference situation. When the interference is large, the value of n needs to be increased; when the interference of the corresponding main frequency point is large, the frequency weighting value w of the main frequency can be reduced. i , and increase the frequency weighting value w of other main frequencies i , making X i is the normalized complex value of the ith main frequency.

[0096] In the embodiment of the present invention, n=4, p=4*j, w i = 0.2, and obtain the normalized complex values ​​of the four combined frequencies. For subsequent calculations, assume that the calculated normalized complex values ​​C j The amplitudes are 0.98, 0.97, 0.96, and 0.95 respectively. The calculated normalized complex value C j The phase is 1.0mrad, 6.0mrad, 10.0mrad, and 15.0mrad.

[0097] Specifically, the apparent resistivity ρ of the four combined frequencies is calculated j The calculation formula is:

[0098] ρ j =k*|C j |

[0099] Where k is the device coefficient calculated from the receiving point location and the power supply point location, C j is the normalized complex value of the jth combination frequency; assuming k = 100, the apparent resistivity ρ j The values ​​are 98, 97, 96, and 95 respectively.

[0100] Calculate the dispersion rate fs of the first three combined frequencies j The calculation formula is:

[0101] fs j =(|C j |-|C h |) / |C j |

[0102] or

[0103] fsj =(|C j |-|C h |) / |C h |

[0104] Among them, C j is the normalized complex value of the jth combination frequency, C h is the normalized complex value of the hth combination frequency, h = j + 1. According to the aforementioned assumptions, the calculated results of the dispersion rate are 1.02%, 1.03%, and 1.04%;

[0105] like Figure 2 The following is a timing diagram of a 5th-order spread-spectrum signal with a fundamental frequency of 1 / 16 Hz, measured according to an embodiment of the present invention. The horizontal axis represents the sample coordinates, and the vertical axis represents the signal amplitude (in mV). The signal was acquired at a sampling frequency of 15 Hz, so a single cycle of the spread-spectrum wave consists of 240 sampling points. The figure shows that the spread-spectrum signal profile consists of only two levels, high and low, with some jitter in the levels representing induced signals and noise.

[0106] like Figure 3 and Figure 4 As shown in the figure, the graph includes three types of data: X, Y, and Z. X is the Fourier transform result, Y is the extracted primary frequency data, and Z is the combined frequency data. Each column in the figure represents a data point, and the column length indicates the signal energy intensity. The corresponding Fourier transform period in both figures is 2. The X data contains 16 primary frequency signals, separated by a single frequency point. Y is the 16 primary frequency signals extracted from the X data. Figure 3 and Figure 4 Both use four main frequencies for combination, but the starting frequencies of the two combinations are different, resulting in different numbers of final combination frequencies. However, the energy of the final combination frequencies is about four times that of the original signal. The main frequency distribution of different spread spectrum signals at a base frequency of 1 / 16 Hz is shown in Table 1 below:

[0107] Table 1

[0108]

[0109]

[0110] As can be seen from Table 1 above, the 5th order spread spectrum wave contains two frequency distributions of calculation parameters. The calculation parameters of combination A are (n=4, p=2.5+j, w i =0.25, j is the combined frequency index starting from 0), the calculation parameters of combination B are (n=8, p=8*j, w i=0.125, j is the combined frequency index starting from 0). It can be seen from the figure that the frequency interval of the square wave is 0.124Hz; the frequency interval of the dual-frequency wave is 0.75Hz; the minimum frequency interval of the 5th-order multi-frequency wave is 0.0625Hz, and then it continues to increase with the increase of frequency, and the maximum is 0.5Hz; the frequency interval of the spread spectrum wave combination A is 0.0625Hz; the frequency interval of the spread spectrum wave combination B is 0.5Hz. Therefore, the method proposed in the present invention can obtain the densest frequency distribution as shown in combination A, and can also obtain a sparse frequency distribution as shown in combination B. Combination A can be used to calculate the induced polarization parameters in exploration areas where the differences in the spectral characteristics of rocks and minerals are small, and combination B can be used to calculate the induced polarization parameters in areas where the differences in the spectral characteristics of rocks and minerals are large. The comparison of the main frequency energy distribution of different induced polarization signals at the fundamental frequency of 1 / 16Hz is shown in Table 2 below:

[0111] Table 2

[0112] Signal Type <![CDATA[F0 energy]]> <![CDATA[F1 energy]]> <![CDATA[F2 energy]]> <![CDATA[F3 energy]]> <![CDATA[F4 energy]]> Square Wave 1 0.33 0.2 0.14 0.11 Dual-frequency wave 1 1 —— —— —— 5th-order multi-frequency wave 1 1 1 1 1 5th-order spread spectrum wave_Combination A 4 4 4 4 4 5th-order spread spectrum wave_combination B 8 8 —— —— ——

[0113] Table 2 shows the theoretical energy distribution of different IP signals corresponding to the frequencies in Table 1 (the actual values ​​will vary slightly due to different signal modulation methods). The values ​​shown in the table are normalized energy values ​​(main frequency energy divided by fundamental frequency energy). The calculation parameters for spread spectrum wave combinations A and B are the same as those in Table 1. The table shows that the main frequency energy of the spread spectrum waves calculated using this scheme is greatly enhanced, and the energy distribution is uniform. In exploration areas with less interference, the parameters of combination A can be used to calculate IP parameters; in areas with greater interference, combination B can be used to calculate IP parameters.

[0114] In induced polarization (IP) exploration, in order to accurately measure the target object's IP phase parameter φ0, it is usually necessary to ensure strict synchronization between the transmitter and receiver. However, in actual construction, due to reasons such as hardware design or GPS synchronization accuracy, there will be an original phase difference φ1 between the transmitter and receiver (φ1 is proportional to the frequency F). Therefore, the actually measured IP phase parameter φ2 = φ0 + φ1, that is, directly measuring the phase will produce an error of φ1.

[0115] Calculate the relative phase cp of the first three combined frequencies j The calculation formula is:

[0116]

[0117] k=F h / F j

[0118]

[0119] in, is the relative phase of the normalized complex value of the jth combination frequency, is the phase of the normalized complex value of the j+1th combination frequency, F h is the frequency value of the j+1th combination frequency, F j The jth frequency combination is calculated as 2.125 mrad, 0.5 mrad, and 3.125 mrad. This relative phase calculation formula can remove some IP coupling, greatly reducing IP parameter distortion, while ensuring that the correct relative phase value can still be calculated when the phase of adjacent frequency points is reversed.

[0120] For the sake of convenience, assume The range is but Taking the phase error φ1 between the transmitter and the receiver into consideration, h =φ h0 +φ h1 and φ j =φ j0 +φ j1 Substitute the relative phase cp j The calculation formula is, then cp j =(φ h0 -k*φ j0 +φ h1 -k*φ j1 ) / (1-k); Since φ1 is proportional to the frequency F, we have φ h1 =F h / F j *φ j1 =k*φ j1 , bring in cp j In the formula, cp j =(φ h0 -k*φ j0 ) / (1-k), that is, the relative phase cp j It is only related to the polarization parameters of the target, eliminating the influence of synchronization errors between the transmitter and receiver. Similarly, the linear term of the induced polarization coupling effect is also proportional to the frequency, and the linear term generated by it can also be eliminated by the calculation formula of this relative phase, so this parameter can reduce the induced polarization coupling effect.

[0121] The pseudo-section diagram drawn using field measured data according to the embodiment of the present invention is as follows: Figure 5As shown, the figure contains the dispersion rate (FS) and relative phase (CP) parameters calculated from the fifth-order spread spectrum wave collected at the same time. The ordinate in the figure is the relative coordinate of the tripole sounding power supply point A divided by 5, and the abscissa is the relative coordinate of the measurement point. The negative values ​​of the dispersion rate and relative phase in the figure are caused by electromagnetic coupling. The larger the ratio of the minimum negative value to the maximum positive value, the stronger the coupling induction. The maximum positive value of the dispersion rate parameter in the figure is approximately 6, the minimum negative value is approximately -44, and the absolute value ratio of negative to positive values ​​is approximately 7.33; the maximum positive value of the relative phase parameter is approximately 47, the minimum negative value is approximately -160, and the negative to positive ratio is approximately 3.4. Therefore, from the comparison of the calculated values, it can be concluded that the electromagnetic coupling induction intensity of the relative phase parameter is much smaller than that of the dispersion rate parameter. The figure also shows that the range of electromagnetic coupling induction of the relative phase parameter is also much smaller than that of the dispersion rate parameter. The range of abnormal morphology displayed by the relative phase parameter is also clearer than that of the dispersion rate parameter. These test results show that the relative phase calculated using the scheme of the present invention can remove some induced polarization coupling, thereby reducing the distortion of the induced polarization parameters.

[0122] The embodiment of the present invention sends a spread spectrum signal to the detection target and obtains the time series data of the spread spectrum signal, obtains the spectrum value of each main frequency in the spread spectrum signal according to the time series data, and obtains the normalized complex value of each main frequency based on the spectrum value of each main frequency; obtains the frequency value of each combination frequency in the spread spectrum signal according to the frequency value of each main frequency; calculates the normalized complex value of each combination frequency according to the normalized complex value of each main frequency; calculates the induced polarization parameters of all combination frequencies according to the frequency values ​​and normalized complex values ​​of all combination frequencies, and the induced polarization parameters include apparent resistivity, dispersion rate and relative phase; realizes the function of obtaining induced polarization parameters using spread spectrum signals, solves the problem of distortion of induced polarization parameters measured by conventional induced polarization methods due to low frequency domain resolution, weak anti-interference ability and large coupling influence; at the same time, according to exploration needs, the frequency density and number of combined frequencies can also be customized to achieve fine characterization of the induced polarization spectrum characteristics of rocks and minerals, enhance the signal energy of a single main frequency point, and improve anti-interference ability.

[0123] An embodiment of the present invention further provides a device for obtaining IP parameters using a spread spectrum signal, comprising:

[0124] A first acquisition module is configured to send a spread spectrum signal to a detection target and acquire time series data of the spread spectrum signal, wherein the time series data includes voltage time series data and current time series data;

[0125] a second acquisition module, configured to acquire a spectrum value of each main frequency in the spread spectrum signal according to the time series data, and obtain a normalized complex value of each main frequency based on the spectrum value of each main frequency;

[0126] A first calculation module is used to obtain the frequency value of each combined frequency in the spread spectrum signal according to the frequency value of each main frequency;

[0127] A second calculation module is used to calculate the normalized complex value of each combination frequency according to the normalized complex value of each main frequency;

[0128] The third calculation module is used to calculate the induced polarization parameters of all the combined frequencies according to the frequency values ​​of all the combined frequencies and the normalized complex values ​​of all the combined frequencies, wherein the induced polarization parameters include apparent resistivity, dispersion rate and relative phase.

[0129] 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.

[0130] 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.

[0131] An embodiment of the present invention further provides a computer-readable storage medium for storing a computer program, and the computer program is executed to implement the above-mentioned method of obtaining induced polarization parameters using a spread spectrum signal.

[0132] 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.

[0133] An embodiment of the present invention further provides a device for obtaining IP parameters using a spread spectrum signal, which is used to implement the above-mentioned calculation method, including:

[0134] memory and processor;

[0135] Memory is used to store computer programs;

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

[0137] 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.

[0138] 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.

[0139] 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.

[0140] 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 present invention. Their specific functions and technical effects can be found in the method embodiment section and will not be described in detail here. 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 considered within the scope of protection of the present invention.

Claims

1. A method for obtaining IP parameters using a spread spectrum signal, characterized in that: include: Step 1: Send a spread spectrum signal to a detection target and obtain time series data of the spread spectrum signal, wherein the time series data includes voltage time series data and current time series data; Step 2: obtaining a spectrum value of each main frequency in the spread spectrum signal according to the time series data, and obtaining a normalized complex value of each main frequency based on the spectrum value of each main frequency; Step 3, obtaining the frequency value of each combined frequency in the spread spectrum signal according to the frequency value of each main frequency; Step 4, calculating the normalized complex value of each combination frequency based on the normalized complex value of each main frequency; Step 5: Calculate the IP parameters of all the combined frequencies according to the frequency values ​​of all the combined frequencies and the normalized complex values ​​of all the combined frequencies. The IP parameters include apparent resistivity, dispersion rate and relative phase.

2. The method for obtaining IP parameters using a spread spectrum signal according to claim 1, wherein: The step 2 includes: Step 21, obtaining the fundamental frequency f0 and order N of the spread spectrum signal according to the time series data; Step 22: Calculate the number of main frequencies Nf=2 according to the order N of the spread spectrum signal. N-1 ; Step 23, based on the basic frequency f0 of the spread spectrum signal, calculate the frequency corresponding to each main frequency as f i =f0*i, i is the main frequency number from low frequency to high frequency; Step 24: Perform Fourier transform on the time series data to obtain a transform result, and read the kth value in the transform result as the spectrum value Z of each main frequency, where k is the frequency point index of each main frequency.

3. The method for obtaining IP parameters using a spread spectrum signal according to claim 2, wherein: Based on the spectrum value of each main frequency, the normalized complex value of each main frequency is obtained, including: Obtaining a spectrum value U of each main frequency in the voltage time series data by using a calculation method for the spectrum value of each main frequency in the spread spectrum signal; Obtaining the spectrum value I of each main frequency in the current time series data by using a calculation method for the spectrum value of each main frequency in the spread spectrum signal; According to the spectrum value U of each main frequency in the voltage time series data and the spectrum value I of each main frequency in the current time series data, the normalized complex value of each main frequency is calculated as X=U / I.

4. The method for obtaining IP parameters using a spread spectrum signal according to claim 2, wherein: The frequency value F of the combined frequency j The calculation formula is: Among them, p is the starting index of the starting main frequency, the value of p is (j-1)*L, the value range of L is [1,n], n is the number of combined frequencies, w i is the frequency weighted value, w i The value of f is 1 / n, i The frequency of the main frequency.

5. The method for obtaining IP parameters using a spread spectrum signal according to claim 3, wherein: The normalized complex value C of the combined frequency j The calculation formula is: Among them, p is the starting index of the starting main frequency, the value of p is (j-1)*L, the value range of L is [1,n], n is the number of combined frequencies, w i is the frequency weighted value, w i The value of X is 1 / n, i is the normalized complex value of the ith main frequency.

6. The method for obtaining IP parameters using a spread spectrum signal according to claim 5, wherein: The apparent resistivity ρ j The calculation formula is: r j =k*|C j | Where k is the device coefficient calculated from the receiving point location and the power supply point location, C j is the normalized complex value of the j-th combination frequency; The dispersion rate fs j The calculation formula is: fs j =(|C j |-|C h |) / |C j | or fs j =(|C j |-|C h |) / |C h | Among them, C j is the normalized complex value of the jth combination frequency, C h is the h-th combination frequency normalized complex value, h = j + 1; The relative phase cp j The calculation formula is: k=F h / F j in, is the normalized complex-valued phase of the jth combined frequency, is the phase of the hth frequency-normalized complex value, F h is the frequency value of the hth combination frequency, F j is the frequency value of the jth combination frequency.

7. A device for obtaining IP parameters using a spread spectrum signal, characterized in that: include: A first acquisition module is configured to send a spread spectrum signal to a detection target and acquire time series data of the spread spectrum signal, wherein the time series data includes voltage time series data and current time series data; a second acquisition module, configured to acquire a spectrum value of each main frequency in the spread spectrum signal according to the time series data, and obtain a normalized complex value of each main frequency based on the spectrum value of each main frequency; A first calculation module is used to obtain the frequency value of each combined frequency in the spread spectrum signal according to the frequency value of each main frequency; A second calculation module is used to calculate the normalized complex value of each combination frequency according to the normalized complex value of each main frequency; The third calculation module is used to calculate the induced polarization parameters of all the combined frequencies according to the frequency values ​​of all the combined frequencies and the normalized complex values ​​of all the combined frequencies, wherein the induced polarization parameters include apparent resistivity, dispersion rate and relative phase.

8. A computer-readable storage medium for storing a computer program, characterized in that: The computer program is executed to implement the method for obtaining induced polarization parameters using a spread spectrum signal as described in any one of claims 1 to 6.

9. A device for obtaining IP parameters using a spread spectrum signal, used to implement the method for obtaining IP parameters using a spread spectrum 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