Controlled-source electromagnetic data processing based on gray system theory and robust estimation

We propose a novel method that combines gray system theory and robust M-estimation method to suppress the interference in controlled-source electromagnetic data. We estimate the standard deviation of the data using a gray model because of the weak dependence of the gray system on data distribution and size. We combine the proposed and threshold method to identify and eliminate outliers. Robust M-estimation is applied to suppress the effect of the outliers and improve the accuracy. We treat the M-estimators of the preserved data as the true data. We use our method to reject the outliers in simulated signals containing noise to verify the feasibility of our proposed method. The processed values are observed to be approximate to the expected values with high accuracy. The maximum relative error is 3.6676%, whereas the minimum is 0.0251%. In processing field data, we observe that the proposed method eliminates outliers, minimizes the root-mean-square error, and improves the reliability of controlled-source electromagnetic data in follow-up processing and interpretation.

Comparison of marine controlled-source electromagnetic data acquisition systems by a reservoir sensitivity index:analyzing the effect of water depths

During the past ten years, a marine controlled source electromagnetic（CSEM） method has been developed rapidly as a technology for hydrocarbon exploration. For shallow water environments, two CSEM data acquisition systems： Seabed Logging（SBL） and towed streamer electromagnetics（TSEM） have been developed in recent years. The purpose is to compare the performance of the SBL and TSEM systems at different water depths. Three different methods for the comparison are presented. The first method is a quick one dimensional sensitivity modelling. As a result, the sensitivity of marine CSEM data increases with water depth for the SBL system. Further, the sensitivity decreases with the increasing water depth for the TSEM system. The two other methods use two dimensional synthetic data from a simple 2-D isotropic model. The second method is a reservoir sensitivity index（RSI） method which has been developed to provide a quick comparison of the two systems. The RSI is calculated as the amplitude of the scattered field dividing by data uncertainty. From the calculations, it is found that with the increasing water depth RSI increases for the SBL system, while it decreases for the TSEM system. The third method uses Occam＇s inversion, and applies an anomaly transverse resistance（ATR） ratio for evaluating the resulting resistivity image. In shallow water environments, the resolution of the CSEM inversion results is good for both the SBL and TSEM systems. In deep water environments, the resolution of the CSEM inversion is better for the SBL system than for the TSEM system. The ATR ratios of the resistivity images show the similar conclusion. The SBL data acquisition system has an advantage in deep water environments. The TSEM system, on the other hand, is preferable for the shallow water environments.

Fast measurement and prediction method for electromagnetic susceptibility of receiver

Aiming at evaluating and predicting rapidly and accurately a high sensitivity receiver’s adaptability in complex electromagnetic environments,a novel testing and prediction method based on dual-channel multi-frequency is proposed to improve the traditional two-tone test.Firstly,two signal generators are used to generate signals at the radio frequency(RF)by frequency scanning,and then a rapid measurement at the intermediate frequency(IF)output port is carried out to obtain a huge amount of sample data for the subsequent analysis.Secondly,the IF output response data are modeled and analyzed to construct the linear and nonlinear response constraint equations in the frequency domain and prediction models in the power domain,which provide the theoretical criteria for interpreting and predicting electromagnetic susceptibility(EMS)of the receiver.An experiment performed on a radar receiver confirms the reliability of the method proposed in this paper.It shows that the interference of each harmonic frequency and each order to the receiver can be identified and predicted with the sensitivity model.Based on this,fast and comprehensive evaluation and prediction of the receiver’s EMS in complex environment can be efficiently realized.

New correction method for controlled-source electromagnetics source effects

A novel method for source effect correction based on integral equation method is proposed.By taking the electrical horizontal field Ex of current source as an example,the correction method is validated using both simulated data and field data.The results show that the correction method is feasible and effective for isotropic media.When the field data are processed,the correction method normalizes the sources with different geological structures,which eliminates the geological difference among sources,and retains the geological difference among receivers.The correction results are in line with the expectation in whole.

Spatio-Temporal Changes of Soil Salinity in Arid Areas of South Xinjiang Using Electromagnetic Induction

The aim of this paper was to research the spatio-temporal changes in total soluble salt content （TS） in a typical arid region of South Xinjiang, China, where the climate is arid and soil salinization happens easily. The total soluble salt content was interpreted by measurements made in the horizontal mode with EM38 and EM31. The electromagnetic induction （EM） surveys were made three times with the apparent soil electrical conductivity （ECa） measurements taken at 3 873 locations in Nov. 2008, 4 807 locations in Apr. 2009 and 6 324 locations in Nov. 2009, respectively. For interpreting the ECa measurements into total soluble salt content, calibtion sites were needed for EM survey of each time, e.g., 66 sites were selected in Nov. 2008 to measure ECa, and soils-core samples were taken by different depth layers of 0-10, 10-20 and 20-40 cm at the same time. On every time duplicate samples were taken at five sites to allevaite the local-scale variability, and soil temperatures in different layers through the profiles were also measured. Factors including TS, pH, water content, bulk density were analyzed by lab experiments. ECa calibration equations were obtained by linear regression analysis, which indicated that soil salinity was one primary concern to ECa with a determination coefficient of 0.792 in 0-10 cm layer, 0.711 in 10-20 cm layer and 0.544 in 20-40 cm layer, respectively. The maps of spatial distribution were predicted by Kriging interpolation, which showed that the high soil salinity was located near the drainage canal, which validated the trend effect caused by the irrigation canal and the drainage canal. And by comparing the soil salinity in different layers, the soluble salt accumulated to the top soil surface only in the area where the soil salinization was serious, and in the other areas, the soil salinity trended to increase from the top soil surface to 40 cm depth. Temporal changes showed that the soil salinity in November was higher than that in April, and the soil salinization trended to aggravate, especially in the top soil layer of 0-10 cm.

Study of electromagnetic and acoustic emission in creep experiments of water-containing rock samples

Based on biaxial shear creep tests conducted on rock samples with different water contents, we present the results of our study on the regularities of electromagnetic and acoustic emission during the process of creep experiments in which we have analyzed the contribution of water to the occurrence of electromagnetic radiation. The result shows that in the creep-fracturing course of rock samples, when the water content increases, the initial frequency and amplitude of electromagnetic and acoustic emission also increases, but at a decreasing growth rate caused by loading stress. This can be used as a criterion for the long-term stability of rock masses under conditions of repeated inundation and discharge of water.

Numerical Method for Solving Electromagnetic Scattering Problem by Many Small Impedance Bodies

In this paper, we study electromagnetic (EM) wave scattering problem by many small impedance bodies. A numerical method for solving this problem is presented. The problem is solved under the physical assumptions ka??1, where a is the characteristic size of the bodies and k is the wave number. This problem is solved asymptotically and numerical experiments are provided to illustrate the idea of the method. Error estimate for the asymptotic solution is also discussed.

Guided and Direct Wave Evaluation of Controlled Source Electromagnetic Survey Using Finite Element Method

Deep target hydrocarbon detection is still challenging and expensive. Direct hydrocarbon indicators (DHIs) in seismic data do not correspond to economical hydrocarbon exploration. Due to unreliability in seismic data for the detection of DHIs, new methods have been investigated. Marine controlled source electromagnet (MCSEM) or Sea bed logging (SBL) is new method for the detection of deep target hydrocarbon reservoir. Sea bed logging has also the potential to reduce the risks of DHIs in deep sea environment. Modelling of real sea environment helps to reduce the further risks before drilling the oil wells. 3D electromagnetic (EM) modelling of seabed logging requires more accurate methods for the detection of hydrocarbon reservoir. Finite element method (FEM) is chosen for the modelling of seabed logging to get more precise EM response from hydrocarbon reservoir below 4000 m from seabed. FEM allows to investigate the total electric and magnetic fields instead of scattered electric and magnetic fields, which shows accurate and precise resistivity contrast below the seabed. From the modelling results, It was investigated that Hz field shows higher magni- tude with 342% than the Ex field. It was observed that 0.125 Hz frequency can be able to show better resistivity contrast of Hz field (31.30%) and Ex field (16.49%) at target depth of 1000 m below seafloor for our proposed model. Hz and Ex field delineation was found to decrease as target depth increased from 1000 m to 4000 m. At the target depth of 4000 m, no field delineation response was seen from the current electromagnetic (EM) antenna used by the industry. New EM antenna has been used to see the EM response for deep target hydrocarbon detection. It was investigated that novel EM antenna shows better delineation at 4000 m target depth for Ex and Hz field up to 10.3% and 15.1% respectively. Novel EM antenna also shows better Hz phase response (128.4%) than the Ex phase response (38.3%) at the target depth of 4000 m below the seafloor.

Semisolid slurry of 7A04 aluminum alloy prepared by electromagnetic stirring and Sc, Zr additions

Slurry preparation is one of the most critical steps for semisolid casting, and its primary goal is to prepare slurry with uniformly distributed fine globules. In this work, electromagnetic stirring(EMS) and the addition of Sc and Zr elements were used to prepare semisolid slurry of 7A04 aluminum alloy in a large diameter slurry maker. The effects of different treatments on the microstructure, composition and their radial homogeneity were investigated. The results show that, compared to the slurry without any treatment, large volume slurry with finer and more uniform microstructure can be obtained when treated by EMS, Sc, or Zr additions individually. EMS is more competent in the microstructural and chemical homogenization of the slurry while Sc and Zr additions are more excellent in its microstructural refinement. The combined treatment of EMS, Sc and Zr produces premium 7A04 aluminum alloy slurry with uniformly distributed fine α-Al globules and composition. The interaction mechanism between EMS and Sc and Zr additions was also discussed.

Progress of Expert Systems in Electromagnetic Engineering

It is urgent to solve various problems in electromagnetic （EM） engineering under the increasingly complicated environment. Some expert systems （ES） come into being just to keep up with the demand for solving these problems. Combined with the analysis of development orES technology and the development trend of EM engineering software in recent years, the application orES technology in EM engineering is discussed, and especially the progress of complete ES in electromagnetic compatible （EMC） is introduced.

A nonlinear control method for the electromagnetic suspension system of the maglev train

To deal with the inherent nonlinearity and open-loop instability of the electromagnetic suspension（EMS） system,a new nonlinear control method is proposed.The simulation results show that,for a PID controller,the over-shoot of the system response to an airgap step disturbance is about 3 mm,and the transient time is 6 s;however,for the proposed nonlinear controller,there is no overshoot and transient time within 2 s.The proposed method has a faster response and stronger robustness.With a designed bi-DSP suspension controller,this nonlinear control method was implemented on the Shanghai Urban Maglev Test Line（SUMTL） to validate its effectiveness and feasibility.

Precision of meshfree methods and application to forward modeling of two-dimensional electromagnetic sources

Meshfree method offers high accuracy and computational capability and constructs the shape function without relying on predefined elements. We comparatively analyze the global weak form meshfree methods, such as element-free Galerkin method （EFGM）, the point interpolation method （PIM）, and the radial point interpolation method （RPIM）. Taking two dimensional Poisson equation as an example, we discuss the support-domain dimensionless size, the field nodes, and background element settings with respect to their effect on calculation accuracy of the meshfree method. RPIM and EFGM are applied to controlled- source two-dimensional electromagnetic modeling with fixed shape parameters. The accuracy of boundary conditions imposed directly and by a penalty function are discussed in the case of forward modeling of two-dimensional magnetotellurics in a homogeneous medium model. The coupling algorithm of EFG-PIM and EFG-RPIM are generated by integrating the PIM or RPIM and EFGM. The results of the numerical modeling suggest the following. First, the proposed meshfree method and corresponding coupled methods are well-suited for electromagnetic numerical modeling. The accuracy of the algorithm is the highest when the support-domain dimensionless size is 1.0 and the distribution of field nodes is consistent with the nodes of background elements. Second, the accuracy of PIM and RPIM are lower than that of EFGM for the Poisson equation but higher than EFGM for the homogeneous medium MT response. Third, RPIM overcomes the matrix inversion problem of PIM and has a wider selection of support-domain dimensionless sizes as compared to RPIM.

Electromagnetic anomaly before earthquakes measured by electromagnetic experiments

Three experiments are carried out for earthquake monitoring using electromagnetic （EM） methods in recent years. Some earthquakes occurred in chance of the measurement time period for each experiment and the anomalies were recorded before the shocks. The observation at a site 20 km away from the epicenter of Zhangbei Ms6.2 earthquake in 1998 shows that the apparent resistivity decreases in the strike direction before and/or during the earthquake and the resistivity increases in the decline direction. This anomalous variation in apparent resistivity accounts for about 20%. The apparent resistivities at a site in the epicentral area decrease in the strike and decline directions before and/or during the earthquake and increase after shocks. The experiments using artificial electromagnetic signals with super low frequency carried out in 1999 show that the resolution and stability of electric and magnetic spectra are improved. The spectra of electric and magnetic fields and apparent resistivity at the Baodi station began to anomalously change two days before the Qian＇an earthquake with 120 km distant to the station. The anomalous variation of electric and magnetic spectra is about twice as great as normal variation and the apparent resistivity changes by about 20%. The measurements in active seismic area of Yunnan province in the year 2005 indicate that the electric and magnetic spectra anomalously change by one order before the Taoyuan earthquake about 100 km away from the observatories. But the measurements at the sites in Beijing area 2 000 km away from the epicenter do not show any anomaly.

Numerical and Experimental Investigation on the Attenuation of Electromagnetic Waves in Unmagnetized Plasmas Using Inductively Coupled Plasma Actuator

The attenuation of electromagnetic （EM） waves in unmagnetized plasma generated by an inductively coupled plasma （ICP） actuator has been investigated both theoretically and experimentally. A numerical study is conducted to investigate the propagation of EM waves in multilayer plasma structures which cover a square fiat plate. Experimentally, an ICP actuator with dimensions of 20 cm×20 cm×4 cm is designed to produce a steady plasma slab. The attenuation of EM waves in the plasma generated by the ICP actuator is measured by a reflectivity arch test method at incident waves of 2.3 GHz and 10.1 GHz, respectively. A contrastive analysis of calculated and measured results of these incident wave frequencies is presented, which suggests that the experiment accords well with our theory. As expected, the plasma slab generated by the ICP actuator can effectively attenuate the EM waves, which may have great potential application prospects in aircraft stealth.

Numerical Method for Solving Electromagnetic Wave Scattering by One and Many Small Perfectly Conducting Bodies

In this paper, we investigate the problem of electromagnetic (EM) wave scattering by one and many small perfectly conducting bodies and present a numerical method for solving it. For the case of one body, the problem is solved for a body of arbitrary shape, using the corresponding boundary integral equation. For the case of many bodies, the problem is solved asymptotically under the physical assumptions a d a is the characteristic size of the bodies, d is the minimal distance between neighboring bodies, λ = 2π/k is the wave length and k is the wave number. Numerical results for the cases of one and many small bodies are presented. Error analysis for the numerical method is also provided.

Inverse Kinematic Solution for Robot Manipulator Based on Electromagnetism-like and Modified DFP Algorithms

为计算数字解决方案到机器的操纵者的反的 kinematics 问题的一个新方法在这份报纸被开发。与联合限制，像电磁的方法(他们) 利用吸引力排斥机制很快向最佳答案移动样品点。基于他们给的这个近似答案，一个修改 Davidon-Fletcher-Powell (DFP ) 算法被开发在需要的精确解决这个问题。不同于传统的算法，这个修改 DFP (MDFP ) 算法随机选择在 0 和 1 之间的搜索步尺寸。因此，计算复杂性极大地被减少。试验性的结果基于十个将军测试功能和这个新在真实附近的时间混血儿方法能最好生产的美洲狮 560 机器人表演性能。

A three-dimensional transient electromagnetic data inversion method based on a time–frequency transformation

Herein,a three-dimensional(3D)inversion method in the frequency domain based on a time–frequency transformation was developed to improve the efficiency of the 3D inversion of transient electromagnetic(TEM)data.The Fourier transform related to the electromagnetic response in the frequency and time domains becomes a sine or cosine transform under the excitation of downward-step current.We established a transformation matrix based on the digital fi ltering calculation for the sine transform,and then the frequency domain projection of the TEM data was determined from the linear transformation system using the smoothing constrained least squares inversion method,in which only the imaginary part was used to maintain the TEM data transformation equivalence in the bidirectional projection.Thus,the time-domain TEM inversion problem was indirectly and effectively solved in the frequency domain.In the 3D inversion of the transformed frequency-domain data,the limited-memory Broyden–Fletcher–Goldfarb–Shannoquasi–Newton(L-BFGS)method was used and modifi ed with a restart strategy to adjust the regularization parameter when the algorithm tended to a local minimum.Synthetic data tests showed that our domain transformation method can stably project the TEM data into the frequency domain with very high accuracy;furthe rmore,the 3D inversion of the transformed frequency-domain data is stable,can be used to recover the real resistivity model with an acceptable effi ciency.

Integrating electromagnetic surface and antenna array for reflection suppression and excellent radiation

The electromagnetic surface antenna array(EMSAA)has been proposed for obtaining reflection suppression and excellent radiation simultaneously.The antenna with rectangular radiation patch is used to design anisotropic electromagnetic surface.Preternatural reflection characteristics of the element antenna can be tailored depending on the incident polarizations.EMSAA can be constructed by using single structured element antenna with 90° rotation and orthometric arrangement.This orthometric arrangement of EMSAA is helpful to achieve reflection suppression and excellent radiation.The simulated results show that the reflection of EMSAA is suppressed from 5.0 GHz to 8.0 GHz with peak reduction of 12.3 dB.The linear-and circular-polarized radiation properties of EMSAA are obtained and the maximum gain is 14.3 dBi.The measured results are consistent with the simulation results.The results demonstrate that the reflection suppression and excellent radiation are achieved simultaneously.Such design of EMSAA will open the path for integrating antenna fields and electromagnetic surface(EMS)fields.

Electromagnetic and Thermal Simulations of Human Neurons for SAR Applications

The impact of the electromagnetic waves (EM) on human neurons (HN) has been under investigation for decades, in efforts to understand the impact of cell phones (radiation) on human health, or radiation absorption by HN for medical diagnosis and treatment. Research issues including the wave frequency, power intensity, reflections and scattering, and penetration depths are of important considerations to be incorporated into the research study. In this study, computer simulation for the EM exposure to HN was studied for the purpose of determining the upper limits of the electric and magnetic field intensities, power consumption, reflections and transmissions, and the change in temperature resulting from the power absorption by human neurons. Both high frequency structural simulators (HFSS) from ANSYS software, and COMSOL multi-physics were used for the simulation of the EM transmissions and reflections, and the temperature profile within the cells, respectively. For the temperature profile estimation, the study considers an electrical source of 0.5 watt input power, 64 MHz. The EM simulation was looking into the uniformity of the fields within the sample cells. The size of the waveguide was set to be appropriate for a small animal model to be conducted in the future. The incident power was fully transmitted throughout the waveguide, and less than 1% reflections were observed from the simulation. The minimum reflected power near the sample under investigation was found to be with negligible reflected field strengths. The temperature profile resulting from the COMSOL simulation was found to be near 0.25 m°K, indicating no change in temperature on the neuro cells under the EM exposure. The paper details the simulation results for the EM response determined by HFSS, and temperature profile simulated by COMSOL.

Separation of inclusions from aluminum melt using alternating electromagnetic field

Effects of processing variables such as frequency of imposed magnetic field, imposed magnetic flux density, processing time, diameter of inclusions, and value of r 1/δ on the electromagnetic separating(EMS) removal efficiency were analyzed theoretically. The higher the frequency, the wider the range of r 1/δ will be. Removal efficiency reaches the maximum while r 1/δ ranges from 1.5 to 2. And the experimental results on aluminum melt show that higher frequency and magnetic flux density make for higher removal efficiency, matching well with the theoretical results. When f is 15.6 kHz, B e is 0.1 T , and imposed time is 10 s, more than 80% inclusion particles with 6 μm diameter can be removed.