Research article 18 Mar 2020
Research article  18 Mar 2020
A comprehensive data quality evaluation method for the currents of marine controlledsource electromagnetic transmitters based on the analytic hierarchy process
 ^{1}China University of Geosciences (Beijing), Beijing, 100083, China
 ^{2}Guangzhou Marine Geological Survey, Guangzhou, 510075, China
 ^{1}China University of Geosciences (Beijing), Beijing, 100083, China
 ^{2}Guangzhou Marine Geological Survey, Guangzhou, 510075, China
Correspondence: Meng Wang (wangmengcugb@qq.com)
Hide author detailsCorrespondence: Meng Wang (wangmengcugb@qq.com)
We present a quality control methodology for the currents of marine controlledsource electromagnetic transmitters . The quality level of the transmitting current directly affects the signaltonoise ratio (SNR) of the electromagneticfield data, as received by a multicomponent electromagnetic receiver from the seabed. Although the transmittingcurrent stability is sufficient under normal circumstances, the SNR of the received signal can change owing to factors such as outside noise. In some emergency cases such as instrument failure or a sudden increase in electromagnetic interference that we are not aware of, the frequency and properties of the transmitting current, such as its size and waveform, may change. The traditional current monitoring and data playback tools fail to detect and evaluate the anomalies well and in a timely manner, which introduces considerable errors in the later dataprocessing procedure. Pertaining to these issues, this paper proposes a comprehensive quality evaluation method for the transmitting current. The proposed algorithm, based on the analytic hierarchy process, is first used to analyze five current stability parameters – current frequency, positive amplitudes, negative amplitudes, discrepancy of ideal waveform, and waveform repetition – and then to define the harmonic energy and calculate the quality of transmitting current (QTC) index of the final data to assess the quality of the transmitting current comprehensively. The results of a marine experiment performed in 2016 show that the algorithm can identify abnormal current data and quantitatively evaluate the current conditions. Under normal circumstances, the QTC index is less than 2 %. The key findings are that the QTC index changes to more than 4 % and some curvilinear features are observed if the transmittingcurrent quality is poor. These results will provide a positive, significant guide for the evaluation and monitoring of transmittingcurrent data in marine experiments.
The marine controlledsource electromagnetic (MCSEM) technique is an effective method of exploring naturalgashydrate reservoirs and petroleum reservoirs (Constable and Srnka, 2007; Wang et al., 2013). With the development of the MCSEM method worldwide, it not only has been used to develop a series of algorithms for forward and inverse calculations (Gribenko and Zhdanov, 2007; Jing et al., 2016) but also has shown great potential for use in practical applications (Cox et al., 1986; Constable, 2010). In realistic marine prospecting work, an MCSEM exploration system is usually composed of a highpower, controlledsource electromagnetic transmitter and a submarine mixedfield source electromagnetic receiver (Chen et al., 2017a, b; Di et al., 2018). Each component remarkably influences the quality, precision, and interference of the electromagneticfield signals. Ensuring the high quality and high signaltonoise ratio (SNR) of the transmitting current are the most important tasks of the transmitting system. In 2007, the EMGS (Electromagnetic Geoservices) company studied the waveform of a transmitting current from the perspective of the harmonicenergy ratio (HER) and proposed a transmittingcurrent waveform that improves the signal quality of the electromagnetic field (Rune and Tor, 2008). Edwards (2005), He et al. (2009) and Luan et al. (2018) state some dataprocessing methods of controlledsource electromagnetic, but they do not undertake too much research on current quality. In this paper, we propose a comprehensive quality evaluation algorithm based on the analytic hierarchy process (Saaty, 1987) for the transmitting current. Considering the impacts of factors such as the frequency stability, positive and negative amplitude stability, discrepancy from the ideal waveform, waveform repetition, and HER, as well as the distribution weight value, a quality of transmitting current (QTC) index can be calculated. This research also provides a reliable path for transmittingcurrent improvement and realization of the realtime monitoring systems.
A personal computer (PC) in a deck unit, used for the realtime monitoring of the transmitting current, monitored the changes in the current amplitude based on simple qualitative observations. Figure 1 shows basic parameters such as the current waveform and frequency components of channel 1 (we use multiple channels to record the data from multiple current sensors). The quality and stability of the transmitting current depend on various parameters such as the transmitting frequency, amplitude, difference from the ideal amplitude, degree of waveform repetition, and harmonicenergy distribution. The MCSEM data preprocessing and inversion are influenced by the transmittingcurrent quality. Therefore, the comprehensive evaluation and feedback of the transmittingcurrent data quality are particularly important.
Earlier, PCs on decks performed realtime monitoring of transmittingcurrent parameters that was limited to observing changes in amplitude, a quality that can only be described qualitatively. To evaluate the properties of current quantitatively, more parameters must be calculated, predominantly the actual work factors such as the frequency, amplitude, idealvalue difference, waveform repetition, and harmonicenergy multiple aspects that influence the quality of the received data. In order to facilitate the analysis of fast Fourier transform (FFT) operations, the data are divided into blocks by fixedperiod number (N), and the corresponding length of time (T) can be obtained by dividing N by the sampling rate (f_{s}) as shown in Eq. (1):
where n is a positive integer, and its value depends on the lengths of the original data (ensuring that there are enough data blocks for analysis) and the transmission frequency ($T\gg \mathrm{1}/f$) to be able to analyze the signals of more cycles. And after that, each data block is calculated by the algorithm. Then, onebyone block analysis of the transmitting data is performed to calculate the QTC index, which reflects the current quality, to evaluate the transmitting current quantitatively.
3.1 Frequency stability
The transmittingsource frequency is one of the core parameters of the MCSEM method and considerably impacts the solution process in the frequency domain. Therefore, the evaluation of the transmitting current must introduce frequencyrelated parameters. Considering the actual work parameters affecting the stability of electromagnetic waves, a frequency stability parameter (a_{i}) was defined and calculated as shown in Eq. (2):
where x_{i} is the frequency of the present data block, x_{j} is the frequency of each block used to calculate the frequency base value (average value), and n is the number of data blocks. The frequency axis is discretized (with step size f_{s}∕N) when the fast Fourier transform is performed; hence, when an x_{i} is input, the actual frequency can be determined by searching the local maximum amplitude. The stability is measured as the ratio of the actual frequency of each data block to the average frequency of all data blocks.
3.2 Positive amplitude stability
The MCSEM method is required to output an alternatingcurrent signal of a certain size while maintaining the stability of the actual transmitting current. Considering an actual transmitting current that deviates slightly numerically, reflecting the positive and negative power supplies of different circuits, the average sizes of the forward and reverse currents, defined as positive and negative amplitudes, respectively, and two parameters, positive and negative amplitude stability, are defined. The formulae for the positive amplitude stability are shown in Eq. (3):
where b_{1} is a given initial value within 1 %, b_{i} is the positive amplitude stability of each data block, and ${J}_{i}^{+}$ is the mean value of the positivesequencecurrent data of data block i.
3.3 Negative amplitude stability
Similarly, the negative amplitude stability is given by Eq. (4):
where c_{1} is a given initial value within 1 %, c_{i} is the negative amplitude stability of each data block, ${J}_{i}^{}$ is the average negativesequencecurrent data of data block i. To evaluate a weighted combination later, the negative sequence is calculated as an absolute value.
3.4 Idealwaveform difference
In general, an MCSEM electricdipole transmitting source is a square wave of a singlefrequency signal or a mixedfrequency signal. However, during an actual transmission, the electromagnetic wave is usually not a standard square wave, for various reasons. Based on this difference, the idealwaveform difference parameter can be defined as shown in Eq. (5):
where J_{k} is the kth value of current data block of the transmitting current, J_{ki} is the ideal transmitting current corresponding to J_{k}, and n is the quantity of data of this block.
3.5 Waveform repetition
The differences between the transmittingcurrent waveforms in adjacent periods can be used to map the stability of the current data with time, and this degree of waveform variation is quantified as the waveform repetition degree parameter, which is given by Eq. (6):
where J_{k} is the kth value of the transmittingcurrent data block, J_{k}−b is the (k−b)th value, n is the number of data to be analyzed, and b is the number of samples per cycle of transmitting waveform.
3.6 HER stability
When frequency domain analysis is performed, the FFT tends to produce higher harmonics, which divide the energy to weaken the signal of fundamental frequency. To introduce a parameter that reflects harmonic energy, the stability of HER is defined as shown in Eq. (7):
where ${h}_{{\mathrm{r}}_{i}}$ is the ratio of the actual assigned current to the theoretical current of the ith fundamentalfrequency point, which can be obtained by the conversion of frequency domain amplitude after the FFT transformation, and n is the number of fundamental frequencies used for synthesis. In this sea trial, the two synthetic fundamentalfrequency points used are 0.5 and 1.5 Hz, respectively, for n=2.
3.7 Evaluation algorithm and comprehensive index
The above parameters are the factors that determine the quality of the transmittingcurrent data. To evaluate the transmittingcurrent quality comprehensively and quantitatively, five general data vectors, a_{i}, b_{i}, c_{i}, d_{i}, and e_{i}, are unified using the analytic hierarchy process (AHP) in the proposed method. Subsequently, the HER stability of all of the data (h_{r}) is combined with w_{6} (an experience weight value) to obtain the composite QTC index. The flowchart of the algorithm is shown in Fig. 2a. The algorithm can be applied if the QTC index (1) can show the difference between the stable and unstable transmitting currents and (2) has a certain ability to detect changes in various stability parameters.
The AHP is a structured technique for organizing and analyzing complex decisions, based on mathematics and psychology. It was developed by Saaty in the 1970s and has been extensively studied and refined since then. It is mainly used for subjective decisionmaking problems under the influence of multiple impact factors. It can be simply divided into the following steps: (1) building a structural model; (2) ranking the importance of the impact factors; (3) comparing and establishing a judgment matrix; (4) calculating the eigenvalues and eigenvectors of the judgment matrix; (5) calculating whether the consistency ratio (CR) is less than 0.10 with the maximum eigenvalue – if so, continue, otherwise return to step 2 and reorder the sorting; and (6) normalizing the eigenvectors of the largest eigenvalues to obtain the weight vector.
The hierarchical analysis model of the five aforementioned parameters is shown in Fig. 2b. The measurement layer includes different channels, and each channel represents data recorded by different sensors. We only have data from two different sensors in this sea trial, but more channels including voltage and other parameters can also use this algorithm in the future. The rule layer contains five general parameters, and the index calculated by the AHP is in the target layer. The judgment matrix, presented in Table 1, is obtained by comparing the five parameter vectors using pairwise comparison of the degree of affecting current mass. For example, if b is the parameter we care the most about and d is the parameter we think has the least impact on current, then $b/d=\mathrm{9}$ and $d/b=\mathrm{1}/\mathrm{9}$; if e is more important than d, then $b/e=\mathrm{7}$, $d/e=\mathrm{1}/\mathrm{2}$, and so on.
A consistency test is performed as follows: the consistency index (CI) is calculated and compared with the random consistency index (RI) to obtain the CR. The formula for the CI is
where λ_{max} is the largest eigenvalue of the judgment matrix and n is the number of factors.
The RI is calculated by randomly constructing a sample matrix, presented in Table 2. The data in this table are the reference data that Saaty obtained through numerous random experiments and can be directly used to calculate the CR.
Subsequently, CR can be calculated as shown in Eq. (9).
When CR < 0.10, the consistency of the judgment matrix is considered acceptable; otherwise, the judgment matrix should be modified appropriately.
When the CR condition is satisfied, the eigenvector corresponding to the maximum eigenvalue of the judgment matrix is obtained. The vector obtained after normalization is the ranking weight (w) of the relative importance of the corresponding factors at the same level as a factor in the previous level. The final QTC index of the ith data block is given by Eq. (10):
where $\mathit{w}=({w}_{\mathrm{1}},{w}_{\mathrm{2}},{w}_{\mathrm{3}},{w}_{\mathrm{4}},{w}_{\mathrm{5}})$ is the corresponding weight vector of five indicators: frequency stability (a_{i}), positive amplitude stability (b_{i}), negative amplitude stability (c_{i}), idealwaveform difference (d_{i}), and waveform repetition (e_{i}). h_{r} is the HER stability, w_{6} is an experience weight, and this formula is calculated in the time domain.
Data examined in this study are obtained from the results of South China sea trials conducted in 2016. The analysis results are shown in Fig. 3. As can be seen from Fig. 3a, the quality of the frequency stability data is high; the positive and negative amplitude changes are less volatile, mostly stable at 1 %; and the idealwaveform difference is about 4 % because the actual transmitting current does not reach the ideal value of 300 A (instead having a maximum of about 290 A). The main effect of it on QTC is to lower the average value. The waveform repetition represents the main reaction of the cycle stability of the output current waveform. And an anomaly occurs when the transmitting current is about to be turned off. It can be seen from Fig. 3b that the results are acceptable and the final QTC index is stable at 1 %. The algorithm is developed for a single frequency, but the transmitting current in the MCSEM method consists of multiple frequencies. To study whether the QTC index is more stable at the component frequencies, the QTC index is calculated multiple times according to different frequencies and a spectrum is plotted. The spectrum diagram (Fig. 3c) mainly shows the changes in the QTC index with frequency and time. Blue indicates smaller stability values and a more stable transmitting current, whereas red indicates the opposite. The main transmission frequencies of 0.5, 1.5 Hz, and higherorder harmonics are relatively stable. Hence, the overall results show that the higher frequencies are more stable than the lower frequencies. And the computing time of the three figures are 1.0410, 1.4447 and 46.1602 s, respectively, in MATLAB.
To confirm the validity of the algorithm and the QTC index of each transmittingcurrent stability parameter detection function, different attributes of squarewave signals are considered as abnormal when waiting to inspect the normal current data. By comparing the differences between the original data and the data after adding abnormal noise, the monitoring current is simulated to analyze the data quality in actual operation.
5.1 Simulation of frequency variation
To simulate the current data fluctuations caused by frequency variations, a squarewave signal with the same sampling rate and peak value as the original data is introduced with an analog signal of frequency 2 Hz. The signal is generated using computer code to simulate the frequencyvarying signals and compared with the original signal to study the ability of the QTC index to detect abnormal frequencies, as shown in Fig. 4. It can be seen from the spectrum analysis results that the amplitude of the data introduced in the signal segment increases significantly. Owing to the influence of odd harmonics, the current size of the frequency point of 6 Hz also increases. Figure 5 shows the changes in the original and analog signal QTC index curve (in blue and red, respectively). It can be seen that the QTC index with the variation frequency band has noticeable changes of about 3 %–4 %. Therefore, the current QTC index can distinctly recognize frequency variations.
5.2 Simulation of amplitude variation
Current amplitude variation is simulated to test the ability of the algorithm to detect transmittingcurrent variations. The middle part of the data is selected for poststack noise processing of the original data (Fig. 6). As the amplitude of only the noise segment is expanded, the overall Fourier transform results show only a slight increase in the amplitude of the base frequency point. In Fig. 7, the blue and red curves represent the QTC indices of the original signal and the data after adding noise, respectively. It is apparent that there are considerable abnormal variation points at the beginning and end of the noise section. Thus, the algorithm can recognize current amplitude variations well.
5.3 Simulation of idealwaveform difference variation
The difference between the actual and ideal current data waveforms is an important factor that affects the transmittingcurrent quality. In the introduced noise section, a 0.5 Hz sine wave is generated to simulate the ability to detect differences from the ideal square wave. It can be observed from the results in Figs. 8 and 9 that, when the current waveform changes, the QTC index also changes, and the simulation data (red) decreases by 2 %–3 % compared to the original data (blue). Consequently, the waveform variation certainly affects the actual QTC index.
5.4 Simulation of waveform repetition variation
Squarewave signals with various frequencies are used as noise sources to simulate waveform repetition degree variation. The results are shown in Figs. 10 and 11. The data from three cycle transformations of square waves with frequencies of 0.5, 2, and 5 Hz, such that each waveform is continuously changing relative to the previous waveform, are used to simulate the transverse variations over time in actual transmitting waveforms. Since the QTC index of the simulated data is smaller than that of the original data, the algorithm also has a certain ability to recognize waveform repetition variations.
5.5 Simulation of HER variation
We define the HER as the ratio of the energy of all the harmonics, except the fundamental wave, to the total energy. The smaller this value is, the larger the expected fundamentalwave energy and the base frequency SNR of the obtained data are. A highfrequency signal is simulated with an 8 Hz frequency square wave as noise, and when the fundamentalfrequency harmonicratio (0.5, 1.5 Hz) energy declines, the harmonic energy increases. The original and simulated data are compared as shown in Fig. 12. It can be seen that the HER of simulated data decreases by 15 %–20 %. The QTC index of the simulated data (red curve in Fig. 13) increases by about 1 % on average compared to that of the original data (blue curve in Fig. 13). Hence, it can be concluded that the harmonic energy affects the QTC index.
5.6 Verification of measured data
The simulation results show that the QTC index has a good variation identification effect on multiple parameters. In particular, it can identify abnormal data and reasons of the transmitting current for transmittingcurrent quality evaluation and monitoring.
Figure 14 shows the transmittingcurrent data from a test. In this experiment, a singlefrequency (8 Hz) squarewave signal and mixedfrequency (0.5, 1.5, 2.5 Hz) signals were sent. The current data in the first three sections of the current in this figure comprise the squarewave signal, and the last part is the mixedfrequency signal. The results obtained after calculation using our algorithm are shown in Fig. 15. In normal stable cases, the QTC index is around 1 %. However, it mutates to 4 %–5 % when a large current signal is transmitted and then returns to the normal value. Since a large current signal input is equivalent to a variation signal, when the algorithm is used to identify the data changes, the QTC index immediately reflects the corresponding variations. Consequently, the algorithm can effectively identify transmittingcurrent data anomalies.
In this paper, an algorithm is proposed to address the lack of a quantitative and comprehensive means of MCSEM transmittingcurrent quality evaluation. After performing calculations using our algorithm, the QTC index can be obtained by combining the frequency, amplitude, waveform difference degree, waveform repetition degree, and harmonicenergy parameters of the transmitting current. The algorithm has the following characteristics:

The calculation process of the algorithm is stable and reliable. It can quantitatively calculate and evaluate the quality of the transmittingcurrent data. Under normal circumstances, the QTC index stability is about 1 %. However, when the current frequency or other properties change, the QTC index increases to more than 4 %. Hence, QTC indices of less than 2 % are normal and those of more than 4 % indicate current data exceptions that require correction in a timely manner.

The algorithm can detect the relationship between transmission current stability and current frequency. The QTC index of the fundamental frequency is smaller than those of the harmonic frequencies, and the frequency is negatively correlated with the QTC index. In other words, the fundamentalfrequency current and relatively high frequency current are more stable.

In addition to enabling quantitative transmission current quality evaluation, the algorithm can also detect abnormal situations in time and the reasons for an abnormal transmitting current. It can be seen from the results of the QTC index simulation of six parameters that an abnormal change in each parameter due to variation certainly influences the composite index. The proposed algorithm can identify the reasons that these abnormal factors affect the transmission current quality. It can also provide a reference for realtime current monitoring of multiple attributes in future research.
There are not publicly available data for this study.
RY wrote the code for the research and analyzed the results. MW provided the data required for the experiment and research ideas. GW provided a lot of help with the experiment. All authors read and provided written comments on the manuscript.
The authors declare that they have no conflict of interest.
We would like to thank China University of Geosciences (Beijing) and Guangzhou Marine Geological Survey for providing us the experiments and data. We would also like to thank Ming Deng and all the coauthors for their suggestions and comments that contributed to improving the paper.
This research has been supported by the National Natural Science Foundation of China (grant nos. 41874142, 61531001, and 41504138) and the Key Technologies R&D Program (grant no. 2016YFC0303100).
This paper was edited by Salvatore Grimaldi and reviewed by Axel Djanni and one anonymous referee.
Chen, K., Jing, J. E., Zhao, Q. X., Luo, X. H., Tu, G. H., and Wang, M.: Ocean bottom EM receiver and application for gashydrate detection, Chinese J. Geophys.Ch., 60, 4262–4272, https://doi.org/10.6038/cjg20171114, 2017a.
Chen, K., Deng, M., Luo, X. H., and Wu, Z. L.: A micro oceanbottom Efield receiver, Geophysics, 82, E233–E241, https://doi.org/10.1190/GEO20160242.1, 2017b.
Constable, S. C.: Ten years of marine CSEM for hydrocarbon exploration, Geophysics, 75, A67–A81, 2010.
Constable, S. C. and Srnka, L. J.: An introduction to marine controlledsource electromagnetic methods for hydrocarbon exploration, Geophysics, 72, WA3–WA12, https://doi.org/10.1190/1.2432483, 2007.
Cox, C. S., Constable, S. C., Chave, A. D., and Webb, S. C.: Controlledsource electromagnetic sounding of the oceanic lithosphere, Nature, 320, 52–54, https://doi.org/10.1038/320052a0, 1986.
Di, Q. Y., Xue, G. Q., Lei, D., Wang, Z. X., Zhang, Y. M., Wang, S., and Zhang, Q. M.: Geophysical survey over molybdenum mines using the newly developed MTEM system, J. Appl. Geophys., 158, 65–70, https://doi.org/10.1016/j.jappgeo.2018.07.008, 2018.
Edwards, R. N.: Marine controlled source electromagnetics: principles, methodologies, future commercial applications, Surv. Geophys., 26, 675–700, https://doi.org/10.1007/s1071200518303, 2005.
Gribenko, A. and Zhdanov, M.: Rigorous 3D inversion of marine CSEM data based on the integral equation method, Geophysics, 72, WA73–WA84, https://doi.org/10.1190/1.2435712, 2007.
He, Z. X., Wang, Z. G., Meng, C. X., Sun, X. M., Hu, X. Y., and Xu, J. H.: Data processing of marine CSEM based on 3D modelling, Chinese J. Geophys.Ch., 52, 2165–2173, https://doi.org/10.3969/j.issn.00015733.2009.08.027, 2009 (in Chinese).
Jing, J. E., Wu, Z. L., Deng, M., Zhao, Q. X., Luo, X. H., Tu, G. H., Chen, K., and Wang, M.: Experiment of marine controlledsource electromagnetic detection in a gas hydrate prospective region of the South China Sea, Chinese J. Geophys.Ch., 59, 2564–2572, https://doi.org/10.6038/cjg20160721, 2016 (in Chinese).
Luan, X. D., Di, Q. Y., Cai, H. Z., Jorgensen, M., and Tang, X. J.: CSAMT Static Shift Recognition and Correction Using Radon Transformation, IEEE Geosci. Remote S., 15, 1001–1005, https://doi.org/10.1109/LGRS.2018.2820743, 2018.
Rune, M. and Tor, S. P.: Shaping optimal transmitter waveforms for marine CSEM surveys, Geophysics, 73, F97, https://doi.org/10.1190/1.2898410, 2008.
Saaty, R. W.: The analytic hierarchy process – what it is and how it is used, Math. Modelling, 9, 161–176, 1987.
Wang, M., Zhang, H. Q., Wu, Z. L., Sheng, Y., Luo, X. H., Jing, J. E., and Chen, K.: Marine controlled source electromagnetic launch system for natural gas hydrate resource exploration, Chinese J. Geophys.Ch., 56, 3708–3717, https://doi.org/10.6038/cjg20131112, 2013 (in Chinese).