Calculation of all-time apparent resistivity of large loop transient electromagnetic method with very fast simulated annealing

In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position information of measuring point.Utilizing very fast simulated annealing(VFSA) to fit the theoretical electromagnetic force(EMF) and measured EMF could obtain the all-time apparent resistivity of the measuring points in rectangular transmitting loop.The selective cope of initial model of VFSA could be confirmed by taking the late time apparent resistivity of transient electromagnetic method as the prior information.For verifying the correctness,the all-time apparent resistivities of the geoelectric models were calculated by VFSA and dichotomy,respectively.The results indicate that the relative differences of apparent resistivities calculated by these two methods are within 3%.The change of measuring point position has little influence on the tracing pattern of all-time apparent resistivity.The first branch of the curve of all-time apparent resistivity is close to the resistivity of the first layer medium and the last branch is close to the resistivity of the last layer medium,which proves the correctness of the arithmetics proposed.

The Effects of Negative Values of Apparent Resistivity in TEM Surveys

The calculation of apparent resistivity in electrical methods is relatively simple and consolidated in the literature. Furthermore, in the case of DC methods, the apparent resistivity values are for the most part, intuitive and enable an initial interpretation of the results. On the other hand, in TEM method the apparent resistivity values are not very intuitive and interpretation based only on the apparent resistivity values cannot be done reasonably. In this way, this paper presents a discussion about the reasons and effects of the negatives values of apparent resistivity in TEM soundings. The main objective of this paper is to clarify the meanings of these negatives values and the possible effects in TEM inversion.

Electric potential and apparent resistivity in rocks containing non-uniformly distributed cracks

In this paper, the formula of electric field distribution and ground apparent resistivity of high resistance rock medium containing low resistance crack are deduced and simulated. The result shows that interstitial parameters, such as buried depth, scale, strike, and real resistivity, etc, have influence on observation and computing result of apparent resistivity. This study provided a useful foundation for earthquake prediction using apparent resistivity method.

A new approach for calculating the apparent resistivity tensor

The apparent resistivity tensor pB is an electromagnetic transfer function,which can be used to analyze and explain the underground electrical structure.Conventional method for obtaining the parameter requires controlled sources and can be easy to be disturbed by cultural noises.We present a new method for calculating the apparent resistivity tensor,the current density is first obtained by measuring the curl operator of the magnetic field on the Earth's surface.This approach is independent of the assumption of a plane wave,and may be used irrespective of source types and field areas.We derived the analytical expressions of the apparent resistivity tensor based on synthetic horizontally layered models with a vertical magnetic dipole source.We then calculate the responses of Pb through numerical modeling examples using both natural sources and controlled sources.Compared to traditional apparent resistivity definitions,our apparent resistivity tensor has the same amplitude value but with more sensitive phases in the far zone,and shows few distortions in the transition zone.And in the near-field zone,it is closer to the resistivity distribution under the ground.The simulation results have demonstrated the effectiveness of the proposed method for calculating the apparent resistivity tensor.

Apparent resistivity of azimuthal anisotropy layered media

The electric field, equations of boundary conditions and calculation formula of apparent resistivity are derived for the azimuthal anisotropy layered media with the DC method based on the anisotropic Ohm's law. Taking the Schlumberger symmetric system as an example and using the recurrence formula of nuclear function, the paper theoretically simulates a model of four layers with the same anisotropy coefficient for each layer. The deep sounding curves of resistivity and the pattern of contours are obtained for the model. The results show that the theoretical formula is correct, and the deep sounding curves not only exhibit the difference of resistivity among layers but also indicate the anisotropy characteristics of layers.

Analysis on Apparent Resistivity Variations of Garzê Station before the 2013 Lushan M_S 7. 0 Earthquake

A trend increase in apparent resistivity has been observed in the N30°E monitoring direction at Garze Seismic Station since July 2011. This increase trend in geo-electric resistivity has been observed in the N60°W direction since 2012. During the period of the increase, the national highway No.317 was expanded in the monitoring area, so the potential electrodes in the N30°E direction had to be moved 10m towards the current electrodes. We interpreted the electric sounding data of Garz6 Seismic Station with a horizontally layered model. Analysis based on this model showed that the shift of potential electrodes can cause a 4 l-l.m rise to the measurements in the N30°E direction. Therefore, apparent resistivity of the two directions increased in the same time in 2012 after offsetting the effects from electrodes shift. Sensitivity coefficients of the two observation directions were also obtained using the model. Sensitivity coefficients of both directions were negative for the shallow layers, which can well explain the unexpected annual variations of Garze Seismic Station. In order to quantitatively analyze the effects from the expansion of the national highway on the observation, we constructed a finite element model based on the electrical structure. Analysis results also suggested that the expansion of the national highway could only cause a 0. 15 Ω·m decrease in the N60°W monitoring direction and 0. 1 Ω· m increase in the N30°E direction. Additionally, the valley values of annual variation of 2013 were distinctively higher than that of other years since 2008, meaning that there was an abnormal rise in apparent resistivity in the two observation directions at Garz~ Seismic Station before the Lushan earthquake. However, the rise was contrary to the decline variation before the Wenchuan earthquake. Therefore, it is still unsure whether or not the rise variation is related to the Lushan earthquake.

Contrast between 2D inversion and 3D inversion based on 2D high-density resistivity data

The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.

Pharmacoepidemiologic study of association between apparent treatment resistant hypertension, cardiovascular disease and interaction effect by sex and age

BACKGROUND A limited number of studies have been conducted to test the magnitudes of the association between apparent treatment resistant hypertension(aTRH)and risk of cardiovascular disease(CVD).AIM To investigate the association between aTRH and risk of CVD and examine whether sex and age modify this association.METHODS We applied an observational analysis study design using data from the United States Antihypertensive and Lipid-Lowering Treatment to Prevent Heart Attack Trial(ALLHAT).ALLHAT recruited participants(n=25516)from 625 primary care settings throughout the United States,Canada,Puerto Rico,and United States Virgin Islands,aged 55 and older with hypertension and at least one additional risk factor for heart disease.aTRH was assessed from the year 2 visit.CVD event was defined as one of the following from the year 2 follow-up visit:Fatal or non-fatal myocardial infarction,coronary revascularization,angina,stroke,heart failure,or peripheral artery disease.Cox proportional hazards regression was used to examine the effect of aTRH on CVD risk.Potential modifications of sex and age on this association were examined on the multiplicative scale by interaction term and additive scale by joint effects and relative excess risk for interaction.RESULTS Of the total study participants(n=25516),5030 experienced a CVD event during a mean of 4.7 years follow-up.aTRH was associated with a 30%increase in risk of CVD compared to non-aTRH[hazards ratio(HR)=1.3,95%CI:1.19-1.42].Sex and age modified this relationship on both multiplicative and additive scales independently.Stratified by sex,aTRH was associated with a 64%increase in risk of CVD(HR=1.64,95%CI:1.43–1.88)in women,and a 13%increase in risk of CVD(HR=1.13,95%CI:1.01–1.27)in men.Stratified by age,aTRH had a stronger impact on the risk of CVD in participants aged<65(HR=1.53,95%CI:1.32–1.77)than it did in those aged≥65(HR=1.18,95%CI:1.05–1.32).Significant two-way interactions of sex and aTRH,and age and aTRH on risk of CVD were observed(P<0.05).The observed joint effect of aTRH and ages≥65 years(HR=1.85,95%CI:1.22–2.48)in males was less than what was expected for both additive and multiplicative models(HR=4.10,95%CI:3.63–4.57 and 4.88,95%CI:3.66–6.31),although three-way interaction of sex,age,and aTRH on the risk of CVD and coronary heart disease did not reach a statistical significance(P>0.05).CONCLUSION aTRH was significantly associated with an increased risk of CVD and this association was modified by both sex and age.Further studies are warranted to test these mechanisms.

Resistivity is used as a tool to evaluate the variability of soil water content

Resistivity is used to evaluate soil water content(SWC),which has the advantages of not causing soil disturbance and in low price.It is an effective way to assess the SWC variability.This paper aims to evaluate the variability of loess slope SWC through the change of resistivity.It provides a simple way for long term SWC monitoring to solve the expensive cost of deploying moisture sensors.In this context,geoelectric and environmental factors such as soil temperature and SWC were monitored for three years.The prediction model of apparent resistivity and SWC was calibrated.The post processing of geoelectric data was introduced.In addition,the SWC collected by Time-Domain Reflectometry(TDR)was used to verify the feasibility of electrical resistivity tomography(ERT)data.The SWC variability in the process of rainfall,the evolution of four seasons,and the alternation of drying and wetting were evaluated.The research results show that:i)the SWC monitored by ERT and TDR can reflect the response and hysteretic effect of water content at 0.5-3.0 m depth.ii)The moisture content monitored by ERT reflects that the soil is relatively wet in summer and autumn and dry in winter and spring.iii)From 2017 to 2020,the SWC increased in August,and the soil became dry in January.iv)Two areas with high SWC and three areas with low SWC on loess slope are reflected by resistivity.The outcome can provide the change information of SWC to a great extent without excavating boreholes.

Inverse Spline Interpolation for All-time Resistivity of Central-Loop TEM

A convenient numerical calculation method (inverse spline interpolation) for all-time apparent resistivity intransient electromagnetic method (TEM) is proposed in this paper. Characteristic of early and late normalized inductiveelectromotive force was investigated. According to the turning point, the transient process is divided into the earlyphase, the turning point, and the late phase. Afterwards, apparent resistivity is obtained through inverse spline interpo-lation in the early and the late phases, respectively. Finally, the resistivities of the early-time and the late-time wereconnected together by the turning point. The result shows that the inverse spline method is feasible and the method alsolays a foundation for initial model construction in the TEM automatic inversion.

Electrical Resistivity Sounding for Groundwater Investigation around Enugu Metropolis and the Environs, Southeast Nigeria

This report evaluates the use of electrical method and borehole data to investigate the subsurface to delineate the groundwater potential in Enugu metropolis and the environs, south-eastern Nigeria other than rely only on resistivity method which could lead to interpretation error. Integrating these 2 data sets is key in this study. The study area is located in the Anambra Basin and is underlain by Nkporo/Enugu Shale which is overlain by the Mamu Formation. It is bounded by Latitudes 6°2 0'00"N to 6°30'00"N and Longitudes 7°25'00"E to 7°35'00"E and covers surface area of about 342 m2. Thirty-one vertical electrical soundings (VES) were carried out across the area using the Schlumberger electrode array with current electrode separation from 2 to 500 m to identify the depths and resistivity values of the identified geo-electric layers. Through data analysis using WinResist software, the apparent resistivity, thicknesses and depths and the thicknesses of the aquifers were generated. The resistivity and depths were modelled to generate resistivity map and depth map. The resistivity of the aquiferous zone within the study area varied from 20.55 - 427.8 ohm-m at depths of between 10.7 - 40.05 m. Depth to the water table appears to be shallow at the south western part of the map. The interpreted geo-electric layers show a sequence of lateritic top soil, shale, sand and shale. The frequency distribution of the VES curves generated shows the presence of 3 to 5 layers with HK type as the highest. Also, a 2D model was generated using the correlation of VES to VES data and borehole data to VES data to show the underlying stratigraphy beneath the study area as well as the direction of ground water flow. Result of the VES curve analysis reveals that the sub-surface is underlain by three lithological layers namely: lateritic top soil, shale, sand and shales with NW direction of groundwater flow from the 2D model. Groundwater prospective zones can be seen along NW, SW and central parts of the study area which have low resistivity values.

Cone-shaped source characteristics and inductance effect of transient electromagnetic method

Small multi-turn coil devices are used with the transient electromagnetic method （TEM） in areas with limited space, particularly in underground environments such as coal mines roadways and engineering tunnels, and for detecting shallow geological targets in environmental and engineering fields. However, the equipment involved has strong mutual inductance coupling, which causes a lengthy turn-off time and a deep “blind zone”. This study proposes a new transmitter device with a conical-shape source and derives the radius formula of each coil and the mutual inductance coefficient of the cone. According to primary field characteristics, results of the two fields created, calculation of the conical-shaped source in a uniform medium using theoretical analysis, and a comparison of the inductance of the new device with that of the multi-turn coil, show that inductance of the multi-turn coil is nine times greater than that of the conical source with the same equivalent magnetic moment of 926.1 A·m2. This indicates that the new source leads to a much shallower “blind zone.” Furthermore, increasing the bottom radius and turn of the cone creates a larger mutual inductance but increasing the cone height results in a lower mutual inductance. Using the superposition principle, the primary and secondary magnetic fields for a conical source in a homogeneous medium are calculated; results indicate that the magnetic behavior of the cone is the same as that of the multi-turn coils, but the transient responses of the secondary field and the total field are more stronger than those of the multi-turn coils. To study the transient response characteristics using a cone-shaped source in a layered earth, a numerical filtering algorithm is then developed using the fast Hankel transform and the improved cosine transform, again using the superposition principle. During development, an average apparent resistivity inverted from the induced electromotive force using each coil is defined to represent the comprehensive resistivity of the conical source. To verify the forward calculation method, the transient responses of H type models and KH type models are calculated, and data are inverted using a “smoke ring” inversion. The results of inversion have good agreement with original models and show that the forward calculation method is effective. The results of this study provide an option for solving the problem of a deep “blind zone” and also provide a theoretical indicator for further research.

Two types of changes in apparent resistivity in earthquake prediction

Two types of changes in apparent resistivity (AR) have been linked to earthquake occurrences. This paper studies the changes and their causes, in detail with the ultimate purpose of developing and assessing a method of earthquake (EQ) prediction. The AR changes of the first type (CFT) are considered to be precursors related to earthquakes (EQs); these appear mostly in the medium-term period before EQs and in the short-term period preceding EQs. The changes of the second type (CST) are characterized by a turning anomaly of a long-trend AR variation or the drastically descending/ascending anomaly superimposed on such a variation; these appear synchronously in large areas, such as the Chinese mainland, and northern and northwestern China, ect. Their spatio-temporal clusters correspond well to high seismicities in the areas and distant great EQs around the Chinese mainland. Based on the behaviors of the two types of changes, the AR changes observed prior to the Ms8.0 Wenchuan EQ of 2008 are studied. The results show that in the medium-term period before the EQ, noticeable anomalies appeared synchronously at four stations around the Songpan-Ganzi active block, but only weak upward changes were observed in the short-term period preceding the EQ, which caused the prediction of the imminent EQ to fail.

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.

Research on application of transient electromagnetic technology in coal mine

Studied the principle of transient electromagnetic method in coalmine and solved the computation of the whole time apparent resistivity and the relation between apparent resistivity and exploration depth and so on. Studied the work method of transient electromagnetic method in coalmine and obtained reasonable arrangement way. Studied data processing and explanation method of transient electromagnetic method and obtained high quality electric section. Finally the purpose to detect water-bearing body and water-bearing structure in front of roadway in advance, and detect the water-bearing property of the roof and floor rock layer of coal face were realized by use of transient electromagnetic method.

Study on electric variations of media in epicentral area by geomagnetic transfer functions

Taking the graph of the apparent resistivity Pyx in Lanzhou area, its geomagnetic transfer function variations andthe investigation of moderate-strong local earthquakes as example and from.the viewpoint of seismomagnetic effect,the authors found that subterranean conductivity values might be applied to the division of Lanzhou areafrom the surface downward into four layers, i. e., the anomalous, the transitional,the non-anomalous and theanomalous layer. Only the transfer function modular |A| and |B|,magnetic azimuth ap with its variance ap,and the total variance yi in the anomalous layers showed obvious anomalies with time changes. Thus, it could beconjectured that the anomalous layer was closely associated with the anisotropic variations of conductivity of thesubterranean media in the seismicity area.This phenomenon might provide a way to the further study of itsmechanism. It is needed to study the inherent connotation deeply.

Some results of modeling D-D seismogenic pattern by the fracture model experiment of large-scale rock samples(I)

sing the natural limestone samples taken from the field with dimension of 500 mm×500 mm×1 000 mm, the D-D (dilatancy-diffusion) seismogeny pattern was modeled under the condition of water injection, which observes the time-space evolutionary features about the relative physics fields of the loaded samples from deformation, formation of microcracks to the occurrence of main rupture. The results of observed apparent resistivity show: ① The process of the deformation from microcrack to main rupture on the loaded rock sample could be characterized by the precursory spatial-temporal changes in the observation of apparent resistivity; ② The precursory temporal changes of observation in apparent resistivity could be divided into several stages, and its spatial distribution shows the difference in different parts of the rock sample; ③ Before the main rupture of the rock sample the obvious ″tendency anomaly′ and ′short-term anomaly″ were observed, and some of them could be likely considered as the ″impending earthquake ″anomaly precursor of apparent resistivity. The changes and distribution features of apparent resistivity show that they are intrinsically related to the dilatancy phenomenon of the loaded rock sample. Finally, this paper discusses the mechanism of resistivity change of loaded rock sample theoretically.

Water abundance of mine floor limestone by simulation experiment

The water abundance of mine floor limestone needs to be solved urgently as the average depth of coal mining in China has increased gradually. A method is presented to discuss water abundance with the numerical and physical layered geoelectrical model being established in the half-space, full-space and full-space with tunnel, respectively. The parameters of water abundance are changed in this study, which includes water quantity, water content and volume of abnormity of water-containing abnormity. Results indicate that the different work fields have different macroscopic influences on the apparent resistivity,and the water abundance parameters of water-containing abnormity have quantitative relationship with the apparent resistivity mean in abnormal regions(three-dimensional space region). The quantitative relationships are shown as following: firstly, the amount of water injection has negative linear correlation with the apparent resistivity mean; secondly, when abnormity is unsaturated, there is a negative power function relationship between water content and apparent resistivity mean; thirdly, the volume of abnormity and apparent resistivity mean behave as a decreasing power function law.

False anomaly recognition with horizontal differential field lines method in high-density electrical technique and its application in Xinlei Quarry,Jiuquan

As an important geophysical tool,high density electrical technique infers the underground geological structures by processing and inverting the apparent resistivity data.Currently,the false anomalies have been frequently occurred in the graph of apparent resistivity pseudo-section or inverted geoelectrical section obtained from high-density electrical technique,and are difficult to remove.In this study,the authors explain the mechanism of the false anomalies and put forward the horizontal differential field method to identify the false anomalies.Based on the analysis of modeling results,this method is applied in the surveying data in Xinlei Quarry of Jiuquan,and the results confirm the effectiveness of the horizontal differential field method.

Earthquake Electromagnetic Precursor Anomalies Detected by a New Ground-based Observation Network

It is a debated topic if there are any observable precursor anomalies prior to the earthquake(EQ hereafter)and if the stronger EQ can be successfully predicted.During last few decades quite a lot of observable electromagnetic(EM)precursors were published by using techniques equipped in either satellites or on ground-based stations.But there are only a few cases that the shortterm precursor anomalies of EM field before earthquakes were observed by using alternate EM fields on ground.This paper will present a new EM observation network built in recent years and show a new finding of EM field with the variation of a one-year cycle observed using the network.As an example,the short-term precursor anomalies of apparent resistivity before the Yangbi EQ(Ms 5.1)occurred on March 27,2017 in Yunnan Province will be studied.The observed anomalous phenomena indicate that the anomaly before the EQ can be captured only if reasonable effective methods including sophisticated analytical techniques are used,and it is believed that continuously observed data on the fixed observation network for a long time is an effective means for studying anomalies that appeared before earthquakes.This network can also play an important role in studying the EM environment from space.