Articles | Volume 28, issue 2
Nonlin. Processes Geophys., 28, 247–256, 2021
https://doi.org/10.5194/npg-28-247-2021
Nonlin. Processes Geophys., 28, 247–256, 2021
https://doi.org/10.5194/npg-28-247-2021

Research article 19 May 2021

Research article | 19 May 2021

An enhanced correlation identification algorithm and its application on spread spectrum induced polarization data

An enhanced correlation identification algorithm and its application on spread spectrum induced polarization data
Siming He1,2, Jian Guan3, Xiu Ji1,4, Hang Xu1, and Yi Wang2 Siming He et al.
  • 1School of Electrical and Information Engineering, Changchun Institute of Technology, Changchun 130000, China
  • 2College of Instrumentation and Electrical Engineering, Jilin University, Changchun 130000, China
  • 3College of Electronic Science and Engineering, Jilin University, Changchun 130000, China
  • 4National Local Joint Engineering Research Center for Smart Distribution Grid Measurement and Control with Safety Operation Technology, Changchun Institute of Technology, Changchun 130000, China

Correspondence: Yi Wang (wangyijlu@jlu.edu.cn) and Xiu Ji (jixiu523@163.com)

Abstract

In spread spectrum induced polarization (SSIP) data processing, attenuation of background noise from the observed data is the essential step that improves the signal-to-noise ratio (SNR) of SSIP data. The time-domain spectral induced polarization based on pseudorandom sequence (TSIP) algorithm has been proposed to improve the SNR of these data. However, signal processing in background noise is still a challenging problem. We propose an enhanced correlation identification (ECI) algorithm to attenuate the background noise. In this algorithm, the cross-correlation matching method is helpful for the extraction of useful components of the raw SSIP data and suppression of background noise. Then the frequency-domain IP (FDIP) method is used for extracting the frequency response of the observation system. Experiments on both synthetic and real SSIP data show that the ECI algorithm will not only suppress the background noise but also better preserve the valid information of the raw SSIP data to display the actual location and shape of adjacent high-resistivity anomalies, which can improve subsequent steps in SSIP data processing and imaging.

1 Introduction

Induced polarization (IP) technology operated in both the time domain and the frequency domain is useful in exploration for groundwater mapping, mineral exploration, and other environmental studies (Revil et al., 2012, 2019; Høyer et al., 2018). Since the phenomenon of IP in the time domain was first discovered by Liu et al. (2017b), there has been consistent efforts to explore its utilization in various research efforts. In 1959, the frequency-domain IP (FDIP) approach was proposed by Collett et al. (1959) and Seigel (1959), which became a classic, widely used mapping technique. For example, the first variable-frequency approach was proposed by Wait (1959), then the spectrum approach of the complex resistivity was developed by Zonge and Wynn (1975), and the dual-frequency IP approach was presented and developed by He (1993) and Han et al. (2013). Recently, spread spectrum induced polarization (SSIP) is a popular branch of FDIP which uses pseudorandom current pulses of opposite polarity as an excitation source (Chen et al., 2007; Xi et al., 2013, 2014; He et al., 2015). According to the intrinsic broadband characteristics of the source itself, the spectral response of an observation system can be simultaneously calculated at multiple frequencies in electrical exploration (Liu et al., 2019). Thus, this SSIP technology has been gaining attention and application in electrical prospecting (Xi et al., 2014; Lu et al., 2019; Wang and He, 2020).

In field detection experiments, it is still a major problem that IP data are often contaminated with background noise. The background noise can be mainly categorized into two types: the Gaussian noise and the impulsive interference with different percentage of outliers (Liu et al., 2016; Kimiaefar et al., 2018; Li et al., 2019). If the background noise is not effectively reduced, the remnant noise can affect the calculation of complex resistivity and may mislead subsequent interpretations of the subsurface structure.

The field of FDIP denoising has achieved quite good results through the constant research of experts and scholars. There have been many algorithms that can be used to suppress the FDIP random noise (Mo et al., 2017), such as smooth filter (Guo, 2017), Mean stack (Liu, 2015), digital filter (Meng et al., 2015), and robust stacking (Liu et al., 2016). The smooth can effectively attenuate Gaussian noise, but the impulsive interference with intense energy leaves the effectiveness of this algorithm limited. Therefore, an effective attenuating algorithm for background noise is still a challenging task for traditional noise suppression algorithms (Neelamani et al., 2008; Liu et al., 2017a). SSIP method also faces the same issue (Liu and Chen, 2016; Liu et al., 2017b).

Recently, the new algorithm based on a circular cross-correlation method, time-domain spectral induced polarization based on pseudorandom sequence (TSIP) algorithm, has also been used to suppress the SSIP noise (Li et al. 2013; Zhang et al., 2020). Due to its effective denoising ability, the identification method has gained more attention and development. However, the TSIP algorithm is restricted because the excitation signal is sensitive to the random noise. For this problem, we propose an enhanced correlation identification (ECI) algorithm for reducing the noise in SSIP data. The ECI algorithm obtains cross-correlations between the transmitter output signal, the excitation signal, and the response signal. The performance of the ECI algorithm is demonstrated on both synthetic and field SSIP data. Experimental results show that the ECI algorithm can effectively control the root mean square of noise (NRMS) increase, enhance its denoising performance in background noise and improve the valid signal preservation to display the actual location and shape of high-resistivity anomalies with higher resolution.

2 Theory

2.1 The TSIP theoretical model

Figure 1 shows a traditional diagram of the electrical resistivity survey. The transmitter output signal uT(t) is poured from electrode A to electrode B, the excitation signal i(t) flows from electrode A to electrode B, and the response signal u(t) between the electrodes M and N is measured. To simultaneously obtain the spectral response of subsurface at various frequencies, pseudorandom sequence based the excitation signal i(t) is considered. Thus, the spectral response of subsurface be retrieved by the TSIP algorithm, and its spectral response be expressed as (Li et al., 2013):

(1) H e ( ω ) = P u i ( ω ) P i i ( ω ) P S ( ω ) ,

where Pui(ω) is the cross-power spectral density of u(t) and i(t), Pii(ω) the auto-power spectral density of i(t), and PS(ω) is the impulse spectral response of the observing system.

https://npg.copernicus.org/articles/28/247/2021/npg-28-247-2021-f01

Figure 1(a) The observation model of the four-electrode measurement. (b) Its equivalent diagram.

Download

Given this observation mode using low-power signals, the magnetotelluric system is a time-invariant system and let us suppose thatHS(ω) is 1. Equation (1) can further be expressed as

(2) H e ( ω ) = P u i ( ω ) P i i ( ω ) = fft R u i ( τ ) fft R i i ( τ ) = U ( ω ) I ( ω ) ,

where fft[.] denotes fast Fourier transform (FFT), Rui(τ) is the cross-correlation function of u(t) and i(t), Rii(τ) is the autocorrelation function of i(t), U(ω) and I(ω) depict the geometric factor defined by the frequency spectrum of u(t)and the frequency spectrum of i(t) respectively, and τdenotes time delay.

https://npg.copernicus.org/articles/28/247/2021/npg-28-247-2021-f02

Figure 2Schematic diagram using the TSIP algorithm.

Download

In the practical field environment, this observation mode is contaminated by the background noise, as shown in Fig. 2. The output of the sensors Ak (k=1,2,3) can be expressed as

(3)y1=uT(t)+n1(t),(4)y2=u(t)+n2(t),(5)y3=i(t)+n3(t),

where nk(t) is the background noise.

Therefore, according to Eq. (2), the formula of the TSIP algorithm is given as

(6) H e ( ω ) = P y 2 y 3 ( ω ) P y 3 y 3 ( ω ) = fft R y 2 y 3 ( τ ) fft R y 3 y 3 ( τ ) = fft R u i ( τ ) + R u n 2 ( τ ) + R i n 1 ( τ ) fft R i i ( τ ) + R i n 1 ( τ ) + R n 1 n 1 ( τ ) fft R u i ( τ ) fft R i i ( τ ) + R n 3 n 3 ( τ ) .

Equation (6) demonstrates that the TSIP algorithm has a weak denoising effect when n3(t) is the massive intense noise. In other words, the TSIP algorithm depends on the energy intensity of n3(t) present in i(t).

https://npg.copernicus.org/articles/28/247/2021/npg-28-247-2021-f03

Figure 3The schematic diagram of the ECI denoising model.

Download

2.2 The ECI theoretical model

That the denoising ability of the TSIP algorithm is limited is caused by that i(t) is sensitive to n3(t). To solve this problem, the ECI algorithm is proposed in Fig. 3 and its derivation process is as follows.

Firstly, let us suppose that the telluric system is a time-invariant system under low-power signals. For three sensor output signals, their cross-correlation functions are the periodic correlation functions of time τ. When the length of the correlation window NT is specified, 0.0125 s in this experiment. The cross-correlation functions can be expressed as follows:

(7)Ry1y2(τ)=Ey1(t)y2(t-τ)=RuTu(τ)N+Rn1n2(τ),(8)Ry1y3(τ)=Ey1(t)y3(t-τ)=RuTi(τ)N+Rn1n3(τ),

where Rn1n2(τ) and Rn1n3(τ) are the cross-correlations of, n2(t) and n3(t) respectively, and τ is time delay that lies in the range of −NT to NT.

Figure 4 shows the schematic diagram of the ZW-CMDSII instrument (Zhang et al., 2014; He et al., 2014). As is known from the figure, we can conclude that uT(t) is mainly disturbed by the floor noise energy of the instrument, and i(t) and u(t) are mainly contaminated by environmental noise. The floor noise is relatively very low, while environmental noise possesses a much higher energy level. Thus we assume that n1(t)≈0, and can conclude that zero correlation between n1(t) and n2(t), n3(t), Rn1n2(τ)0 and Rn1n3(τ)0.

https://npg.copernicus.org/articles/28/247/2021/npg-28-247-2021-f04

Figure 4Schematic diagram of the instrument.

Download

Based on the above analyses, we can further obtain:

(9)Ry1y2(τ)RuTu(τ)N,(10)Ry1y3(τ)RuTi(τ)N.

Then the cross-power spectrum of Eqs. (9) and (10) can be written as following

(11)Py1y2(ω)PuTu(ω),(12)Py1y3(ω)PuTi(ω).

Finally, according to Eqs. (2) and (11), Eq. (12) can be expressed as following

(13) H e ( ω ) = U ( ω ) I ( ω ) = U ( ω ) U T ( ω ) I ( ω ) U T ( ω ) = P u T u ( ω ) P u T i ( ω ) P y 1 y 2 ( ω ) P y 1 y 3 ( ω ) = P y 1 y 2 ( ω ) P y 1 y 3 ( ω ) e - j φ y 1 y 2 ( ω ) - φ y 1 y 3 ( ω )

where φy1y2(ω) and φy1y3(ω) denotes the difference between y1(t), y2(t) and y3(t).

So, Eq. (13) is the formula of the ECI algorithm. The derivation process of this formula clearly describes that the ECI algorithm can effectively suppress the background noise and be independent on the degree of n3(t) present in i(t).

https://npg.copernicus.org/articles/28/247/2021/npg-28-247-2021-f05

Figure 5(a) Experimental schematic; (b) experimental setup.

Download

3 Experiment on synthetic SSIP data record

We test the ECI algorithm for attenuating background noise of SSIP data sets in comparison with the FDIP algorithm and the TSIP algorithm. For the comparison, the signal-to-noise ratio (SNR), root mean square of noise (NRMS), and relative error (ε) are the objective parameters to judge the performance of denoising, which are calculated as follows:

(14)SNR=10log10i=1My(i)-μy2i=1Mn(i)-μn2,(15)NRMS=i=1M[n(i)]2M,(16)ε=100×ρ1-ρ0ρ0,(17)i(t)=ui(t)/RS,

where μy and μn denote the mean values of the useful signal and the noise separately. y(i) and n(i) are the useful signal and the noise separately, M is the length, ρ0 denotes the complex resistivity calculated without noise, and ρ1 is the complex resistivity calculated with the noise added to ρ0. RS is the value of the sampling resistor (RS=1Ω), and ui(t) is the voltage at the sampling resistor.

To validate the effectiveness of the ECI system, we performed a resistance–capacitance experiment, as shown in Fig. 5. The circuit parameters are chosen to be RA=30.3Ω/5 W, RMN=30.1Ω/5 W, RB=30Ω/5 W, and CMN=470µF. We recorded the applied voltage uT(t), the injected current i(t), and the measured potential signal u(t) as the raw signals. These signals are a three-order spread spectrum pseudorandom sequence at the clock cycle of 0.0125 s, as shown in Fig. 6a–c and Table 1.

https://npg.copernicus.org/articles/28/247/2021/npg-28-247-2021-f06

Figure 6The time waves of (a) the applied voltage uT(t), (b) the measured potential signal u(t), (c) the voltage ui(t) at the sampling resistor, (d) noise (t), and (e) the synthetic signal (ui(t)+noise(t)).

Download

Table 1Amplitude and phase values of complex resistivity obtained with Fig. 6a–c.

Download Print Version | Download XLSX

https://npg.copernicus.org/articles/28/247/2021/npg-28-247-2021-f07

Figure 7Amplitude and phase of complex resistivity values at (a1, b1) 80 Hz, and (a2, b2) 160 Hz, (a3, b3) 320 Hz using the three methods.

Download

Since our experiment is in a stable environment, we consider the system to be linear time-invariant, and the noise from the current and voltage measurements are linearly superpositioned (Pelton and Sill, 1983; Vinegar and Waxman, 1984; Garrouch and Sharma, 1998). Therefore, it is actually equivalent whenever the noise is added to the injected current i(t), the measured potential signal u(t), or the applied voltage uT(t). Therefore, the injected current i(t) is only polluted by the synthetic background noise, including Gaussian and impulsive, as shown in Fig. 6d and e. Thirdly, the complex resistivity of the main frequency is considered and discussed because the main energy of the pseudorandom signal is concentrated on the main frequency (He, 2017). Finally, for detailed comparisons between the ECI algorithm and the others, we add the synthetic Gaussian and impulsive noises to the response signal i(t), respectively.

We use synthetic Gaussian noise with the deviation and mean values of 0.1 and 1.1 as a standard template. The excitation signal i(t) is polluted by synthetic different energy levels of the Gaussian noise. Figure 7 shows that the denoised results are obtained and compared at the three main frequencies when the NRMS ranges from 0.12 to 0.25. The figure shows that as the NRMS increases, the complex resistivity information obtained by each algorithm decreases. However, the amplitude spectrum after ECI processing has the slowest-falling speed, and the phase spectrum has the slowest-falling speed at 80 Hz.

Previous literature has shown that if the percentages of outliers in impulsive noise exceed 50 %, the traditional denoising algorithm will be limited (Liu and Chen, 2016; Liu et al., 2017a). Thus, synthetic impulsive noise is added to the excitation signal i(t) in 10 % steps. Their standard deviations (SDs) and skewness values (SKs) are shown in Fig. 8. As depicted in Fig. 9, the three algorithms have a certain degree of denoising performance versus the different percentages of the synthetic outliers against the raw data. The figure shows that with the discrete points of impulse noise growing, the NRMS is different. The amplitude spectrum and phase spectrum of complex resistivity obtained by each algorithm fluctuate. The amplitude spectrum after ECI processing remained the slowest-falling speed. Although the noise reduction performance of the phase spectrum processed by ECI does not stand out, the overall change of the amplitude spectrum after ECI processing is still slow, especially when the discrete point is more than 60 %.

https://npg.copernicus.org/articles/28/247/2021/npg-28-247-2021-f08

Figure 8The standard deviations (SDs) and skewness values (SKs) of synthetic impulsive noise.

Download

https://npg.copernicus.org/articles/28/247/2021/npg-28-247-2021-f09

Figure 9Complex resistivity values at (a1, b1) 80 Hz, (a2, b2) 160 Hz, and (a3, b3) 320 Hz using the three methods.

Download

4 Experiment on real SSIP data record

To further verify the performance of the ECI algorithm, the Wenner array, which is the traditionally applied system in the field, was selected for performing laboratory tests, as shown in Figs. 10 and 11. SSIP data was acquired with a high-density meter and 20 electrodes at 1 m spacing. A Wenner acquisition sequence was adopted with 55 potential measurements expressed using the green points. The figure shows an example of two high-resistance cavities. The two cavities were presented by the letters A and B, and their calibers were about 1.8 m × 2 m. The two cavities are buried by loess. The loess is measured to have an electronic resistivity of 36 Ωm. The measured excitation signal had a range between 0.04 and 0.19 A approximately. The transmitter output signal is a three-order sequence with 80 Hz frequency, and its voltage is about ±11.8 V. The sampling frequency is 625 kHz. The excitation and response data of 40 periods were recorded at each point.

https://npg.copernicus.org/articles/28/247/2021/npg-28-247-2021-f10

Figure 10Diagram of the field test.

Download

https://npg.copernicus.org/articles/28/247/2021/npg-28-247-2021-f11

Figure 11The schematic of the two high-resistance cavities.

Download

https://npg.copernicus.org/articles/28/247/2021/npg-28-247-2021-f12

Figure 12Inverted resistivity sections of the two high-resistivity anomalies (A and B) at 80 Hz with using (a) the FDIP method, (b) the TSIP algorithm, and (c) the ECI algorithm.

Download

Figure 12 demonstrates the experimental SSIP data processed by the three algorithms, inverted with Res2DInv (Arifin et al., 2019). It can be observed that the location and shape of two abnormal bodies are distinguished only in the ECI algorithm while recognized as one whole in the other algorithms. We believe the reason that ECI has higher detection precision is due to its higher denoising ability.

https://npg.copernicus.org/articles/28/247/2021/npg-28-247-2021-f13

Figure 13Standard deviation (SD) values of the ECI algorithm and the others compared to the data dots from 18 to 50 at 80 Hz.

Download

To verify the reason for the improved detecting precision, the SDs of data points are calculated from 18 to 50 (Fig. 11), as shown in Fig. 13. This figure shows that the 33 SD in ECI processing the SSIP data is the lowest at all points. The average SD values in ECI processing of the SSIP data are 7 % and 3.8 % lower than the FDIP and TSIP methods, respectively. Also, the maximum value of SDs with the ECI method is 5 % and 1.4 % lower than the others, and the minimum value is 8 % and 10 % lower, respectively.

Meanwhile, amplitude–frequency |ρ(f)| and phase–frequency φ(f) characteristics of complex resistivity are calculated by the three algorithms (one period) in survey point no. 38 in Fig. 11.

https://npg.copernicus.org/articles/28/247/2021/npg-28-247-2021-f14

Figure 14Complex resistivity spectrum calculated by the three algorithms (one period) in survey point no. 38.

Download

For example, Fig. 14a1 and a2 show that the amplitude and phase of the complex resistivity spectrum for this point at 80 Hz processed by FDIP are 39.7 Ωm and 0.0881 rad, the amplitude and phase are 40.9 Ωm and 6.12 rad when at 160 Hz, and the amplitude and phase are 38.7 Ωm and 0.253 rad when at 320 Hz. As depicted in Fig. 14, the complex resistivity processed by the ECI shows a linear trend with the three main frequencies. Also, the SD of the amplitude–frequency |ρ(f)| characteristic is 0.10 and 0.49 lower than the others, and the SD of the phase–frequency φ(f) is 3.56 and 0.03 lower. Therefore, we believe that the ECI algorithm has an advantage in suppressing background noise, which benefits the subsequent steps in SSIP data processing and imaging.

5 Discussion

The simulation results indicate that the ECI algorithm has very good performance in noise reduction and robustness. Along with the increase of the Gaussian noise level, we found that the ECI algorithm can, to some extend, overcome the shortcomings of the TSIP algorithm has, i.e., being susceptible to the noise of the current. This result coincided with Eqs. (6) and (13), which provides a novel approach for correlated identification noise reduction. In the impulsive noise experiment, we found that the ECI algorithm still has good noise reduction when the discrete point is more than 60 %, which compensates for the disadvantage of the traditional denoising algorithm. Moreover, these simulation results also reveal that the ECI algorithm should have high robustness.

The standard deviation analysis of the real data indicates that the ECI algorithm improves the accuracy and robustness of the collected data, which are compatible with the simulation analyses. This consistency shows that the ECI algorithm can obtain the location and shape of two abnormal bodies by improving the SNR of SSIP data, which can increase the resolution of inversion results.

6 Conclusions

We propose the ECI algorithm that effectively attenuates the background noise in SSIP data and improves the complex resistivity spectrum. This algorithm uses the correlation function to neutralize the influence of the background noise in the SSIP data, and the spectrum complex resistivity can be calculated at multiple frequencies by the formula of the complex resistivity. Simulation results show that the ECI algorithm can effectively attenuate the background noise and improve the SNR. Subsequently, the practicability of the ECI algorithm is further verified by a field test. The results demonstrate that the SD of the SSIP data is improved, which benefits the accuracy and stability of the collected data. There is a good agreement between the complex resistivity and the geological target body with high resistance, which indicates that the ECI algorithm can help to improve the quality of interpretation and inversion in the survey area. For the amplitude spectrum, the ECI algorithm can more effectively suppress the background noise, including the Gaussian random and impulsive noises. Still, its effect is very limited for the phase spectrum. Therefore, a denoising algorithm based on pseudorandom sequence correlation identification is still left open for more investigation.

Code availability

The code is a collection of routines in MATLAB (MathWorks) and is available upon request to the author (e-mail: hsmfly1982@163.com).

Data availability

All the SSIP data are collected by the ZW-CMDSII and are available upon request to the author (e-mail: hsmfly1982@163.com).

Author contributions

SH and YW designed the study, performed the research, analyzed data, and wrote the paper. JG contributed to language polishing and response. XJ and HX contributed to refining the ideas, carrying out additional analyses, and finalizing this paper.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

We are grateful for the help of Jun Wang, Shi Zhu, Hui Wang, and Jinshi Cui. We thank the editors and the reviewers for the constructive comments that helped to improve this article.

Financial support

This research has been supported by the Key Technology Projects of Science and Technology Department of Jilin Province Scientific (grant no. 20190303015SF), a research project of Jilin Provincial Department of Education (grant nos. JJKH20210692KJ and JJKH20211053KJ), and the Fundamental Research and Theme Funds for Changchun Institute of Technology, China.

Review statement

This paper was edited by Richard Gloaguen and reviewed by three anonymous referees.

References

Arifin, M. H., Kayode, J. S., Izwan, M. K., Zaid, H. A. H., and Hussin, H.: Data for the potential gold mineralization mapping with the applications of Electrical Resistivity Imaging and Induced Polarization geophysical surveys, Data in Brief, 22, 830–835, https://doi.org/10.1016/j.dib.2018.12.086, 2019. 

Chen, R. J., Luo, W. B., and He, J. S.: High precision multi-frequency multi-function receiver for electrical exploration, 2007 8th International Conference on Electronic Measurement and Instruments (ICEMI'07), 16–18 August 2007, Xian, China, IEEE, Expanded Abstracts, 599–602, https://doi.org/10.1109/icemi.2007.4350521, 2007. 

Collett, L. S., Brant, A. A., Bell, W. E., Ruddock, K. A., Seigel, H. O., and Wait, J. R.: Laboratory investigation of overvoltage, Overvoltage research and geophysical applications, International series of monographs on earth sciences, Pergamon, New York, USA, 50–70, https://doi.org/10.1016/b978-0-08-009272-0.50009-1, 1959. 

Garrouch, A. A. and Sharma M. M.: Dielectric dispersion of partially saturated porous media in the frequency range 10 Hz to 10 MHz, The Log Analyst, 39, 48–53, 1998. 

Guo, H.: Study of key technology and data fusion of multi-probe penetration based on gas hydrate exploration, PhD Thesis, China University of Geosciences, Wuhan, China, 2017. 

Han, S. L., Zhang, S. G., Liu, J. X., Hu, J., and Zhang, W. S.: Integrated interpretation of dual frequency induced polarization measurement based on wavelet analysis and metal factor methods, T. Nonferr. Metal. Soc., 23, 1465–1471, https://doi.org/10.1016/S1003-6326(13)62618-7, 2013. 

He, G.: Wide area electromagnetic method and pseudo random signal method, Higher Education Press, Beijing, China,, 2017. 

He, G., Wang, J., Zhang, B. Y., Li, M., and Ma, C.: Design of High-density Electrical Method Data Acquisition System, Instrument Technique and Sensor, 8, 18–19, https://doi.org/10.3969/j.issn.1002-1841.2014.08.007, 2014. 

He, J. H., Yang, Y., Li, D. Q., and Weng, J. B.: Wide field electromagnetic sounding methods, in: Symposium on the Application of Geophysics to Engineering and Environmental Problems (SAGEEP 2015), 22–26 March 2015, Texas, USA, EEGS, Expanded Abstracts, 325–329, https://doi.org/10.4133/sageep.28-047, 2015. 

He, J. S.: Dual-frequency IP, T. Nonferr. Metal. Soc., 3, 1–10, 1993. 

Høyer, A. S., Klint, K. E. S., Fiandaca, G., Maurya, P. K., Christiansen, A. V., Balbarini, N., Bjerg, P. L., Hansen, T. B., and Møller, I.: Development of a high-resolution 3D geological model for landfill leachate risk assessment, Eng. Geol., 249, 45–59, https://doi.org/10.1016/j.enggeo.2018.12.015, 2018. 

Kimiaefar, R., Siahkoohi, S. H., Hajian, A., and Kalhor, A.: Random noise attenuation by Wiener-ANFIS filtering, J. Appl. Geophys., 159, 453–459, https://doi.org/10.1016/j.jappgeo.2018.05.017, 2018. 

Li, G., Liu, X., Tang, J., Li, J., Ren, Z., and Chen, C.: De-noising low-frequency magnetotelluric data using mathematical morphology filtering and sparse representation, J. Appl. Geophys., 172, 103919, https://doi.org/10.1016/j.jappgeo.2019.103919, 2019. 

Li, M., Wei, W., Luo, W., and Xu, Q: Time-domain spectral induced polarization based on pseudo-random sequence, Pure Appl. Geophys., 170, 2257–2262, https://doi.org/10.1007/s00024-012-0624-z, 2013. 

Liu, N.: Preprocessing and Research of denosing methods for marine controlled source electromangnetic data, MSc Thesis, Jilin University, Jilin, China, 2015. 

Liu, W. Q. and Chen, R. J.: Coherence analysis for multi-frequency induced polarization signal processingin strong interference environment, The Chinese Journal of Nonferrous Metals, 26, 655–665, https://doi.org/10.19476/j.ysxb.1004.0609.2016.03.022, 2016 (in Chinese). 

Liu, W. Q., Chen, R. J., Cai, H. Z., and Luo, W. B.: Robust statistical methods for impulse noise suppressing of spread spectrum induced polarization data, with application to a mine site, Gansu province, China, J. Appl. Geophys., 135, 397–407, https://doi.org/10.1016/j.jappgeo.2016.04.020, 2016. 

Liu, W. Q., Chen, R. J., Cai, H. Z., Luo, W. B., and Revil, André: Correlation analysis for spread spectrum induced polarization signal processing in electromagnetically noisy environments, Geophysics, 82, E243–E256, https://doi.org/10.1190/geo2016-0109.1, 2017a. 

Liu, W. Q., Wang, J. L., and Lin, P. R.: Signal processing approaches to obtain complex resistivity and phase at multiple frequencies for the electrical exploration method, B. Geofis. Teor. Appl., 58, 103–114, https://doi.org/10.4430/bgta0190, 2017b. 

Liu, W. Q., Lü, Q. T., Chen, R. J., Lin, P. R., Chen, C. J., Yang, L. Y., and Cai, H. Z.: A modified empirical mode decomposition method for multiperiod time-series detrending and the application in full-waveform induced polarization data, Geophys. J. Int., 217, 1058–1079, https://doi.org/10.1093/gji/ggz067, 2019. 

Lu, Q. T., Zhang, X. P., Tang, J. T., Jin, S., Liang, L. Z., Wang, X. B., Lin, P. R., Yao, C. L., Gao, W. l., Gu, J. S., Han, L. G., Cai, Y. Z., Zhang, J. C., Liu, B. L., and Zhao, J. H.: Review on advancement in technology and equipment of geophysical exploration for metallic deposits in china, Chinese J. Geophys., 62, 3629–3664, https://doi.org/10.6038/cjg2019N0056, 2019 (in Chinese). 

Meng, Q. X., Pan, H. P., and Luo, M.: A study on the discrete image method for calculation of transient electromagnetic fields in geological media, Appl. Geophys., 12, 493–502, https://doi.org/10.1007/s11770-015-0517-x, 2015. 

Mo, D., Jiang, Q. Y., Li, D. Q., Chen, C. J., Zhang, B. M., and Liu, J. W.: Controlled-source electromagnetic data processing based on gray system theory and robust estimation, Appl. Geophys., 14, 570–580, https://doi.org/10.1007/s11770-017-0646-5, 2017. 

Neelamani, R., Baumstein, A. I., Gillard, D. G., Hadidi, M. T., and Soroka, W. L.: Coherent and random noise attenuation using the curvelet transform, The Leading Edge, 27, 240–248, https://doi.org/10.1190/1.2840373, 2008. 

Pelton, W. H. and Sill, W. R.: Interpretation of complex resistivity and dielectric data, Geophysical Transactions, 29, 297–330, 1983. 

Revil, A., Karaoulis, M., Johnson, T., and Kemna, A.: Review: Some low-frequency electrical methods for subsurface characterization and monitoring in hydrogeology, Hydrogeol. J., 20, 617–658, https://doi.org/10.1007/s10040-011-0819-x 2012. 

Revil, A., Razdan, M., Julien, S., Coperey, A., Abdulsamad, F., Ghorbani, A., and Rossi, M.: Induced polarization response of porous media with metallic particles – Part 9: Influence of permafrost, Geophysics, 84, E337–E355, https://doi.org/10.1190/geo2019-0013.1, 2019. 

Seigel, H. O.: Mathematical formulation and type curves for induced polarization, Geophysics, 24, 547–565, https://doi.org/10.1190/1.1438625, 1959. 

Vinegar, H. J. and Waxman, M. H.: Induced polarization of shaly sands, Geophysics, 49, 1267–1287, https://https://doi.org/10.1190/1.1441755, 1984. 

Wang, Y. B. and He, J. S.: A hybrid coding and its application to the oil and gas fracturing intelligent real time monitoring system based on pseudorandom signal, Geophysical and Geochemical Exploration, 44, 74–80, https://doi.org/10.11720/wtyht.2020.2288, 2020. 

Wait, J. R.: The variable-frequency method, Overvoltage research and geophysical applications, International series of monographs on earth sciences, Pergamon, 29–49, https://doi.org/10.1016/b978-0-08-009272-0.50008-x, 1959. 

Xi, X. L., Yang, H. C., He, L. F., and Chen, R. J.: Chromite mapping using induced polarization method based on spread spectrum technology, Symposium on the Application of Geophysics to Engineering and Environmental Problems (SAGEEP 2013), 17–21 March 2013, Denver, Colorado, USA, EEGS, Expanded Abstracts, 13–19, https://doi.org/10.4133/sageep2013-015.1, 2013. 

Xi, X. L., Yang, H. C., Zhao, X. F., Yao, H. C., Qiu, J. T., Shen, R. J., and Chen, R. J.: Large-scale distributed 2D/3D FDIP system based on ZigBee network and GPS, Symposium on the Application of Geophysics to Engineering and Environmental Problems (SAGEEP 2014), 16–20 March 2014, Boston, Massachusetts, USA, EEGS, Expanded Abstracts, 130–139, https://doi.org/10.1190/SAGEEP.27-055, 2014. 

Zhang, B. Y., He, G., and Wang J.: New High-density Electrical Instrument Measuring System, Instrument Technique and Sensor, 1, 24–26, https://doi.org/10.3969/j.issn.1002-1841.2014.01.009, 2014.  

Zhang, Q. D., Hao, K. X., and Li, M.: Improved correlation identification of subsurface using all phase FFT algorithm, KSII Transactions on Internet & Information Systems, 14, 495–513, https://doi.org/10.3837/tiis.2020.02.002, 2020. 

Zonge, K. L. and Wynn, J. C.: Recent advances and applications in complex resistivity measurements, Geophysics, 40, 851–864, https://doi.org/10.1190/1.1440572, 1975. 

Download
Short summary
We propose an enhanced correlation identification (ECI) algorithm to attenuate background noise. The cross-correlation matching method is helpful for the extraction of useful components of the raw SSIP data and suppression of background noise. Experiments on both synthetic and real SSIP data show that the ECI algorithm is proposed to preserve the valid information of the raw SSIP data to display the actual location and shape of adjacent high-resistivity anomalies.