Theoretical Study on Stochastic Modeling of Combined Gravity-Magnetic-Electric-Seismic Inversion and Its Application

As gravity field, magnetic field, electric field and seismic wave field are all physical fields, their object function, reverse function and compound function are certainly infinite continuously differentiable functions which can be expanded into Taylor (Fourier) series within domain of definition and be further reduced into solving stochastic distribution function of series and statistic inference of optimal approximation. This is the basis of combined gravity-magnetic-electric-seismic inversion of stochastic modeling. It is an uncertainty modeling technology of combining gravity-magnetic-electric-seismic inversion built on the basis of separation of field and source gravity-magnetic difference-value (D-value) trend surface, taking distribution-independent fault system as its unit, depths of seismic and electric interfaces of interests as its corresponding bivariate compound reverse function of gravity-magnetic anomalies and using high order polynomial (high order trigonometric function) approximating to its series distribution. The difference from current dominant inversion techniques is that, first, it does not respectively create gravity-seismic, magnetic-seismic deterministic inversion model from theoretical model, but combines gravity-magnetic-electric-seismic stochastic inversion model from stochastic model; second, after the concept of equivalent geological body being introduced, using feature of independent variable of gravity-magnetic field functions, taking density and susceptibility related to gravity-magnetic function as default parameters of model, the deterministic model is established owing to better solution to the contradiction of difficulty in identifying strata and less test analytical data for density and susceptibility in newly explored area; third, under assumption of independent parent distribution, a real modeling by strata, the problem of difficult plane closure arising in profile modeling is avoided. This technology has richer and more detailed fault and strata information than sparse pattern seismic data in newly explored area, successfully inverses and plots structural map of Indosinian discontinuity in Hefei basin with combined gravity-magnetic-electric-seismic inversion. With development of high precision gravity-magnetic and overall geophysical technology, it is certain for introducing new methods of stochastic modeling and computational intelligence and promoting the development of combined gravity-magnetic-electric-seismic inversion to open a new substantial path.

Magnetic Charge Theory—The Unification of Gravity with Electricity and Magnetism

The concept of magnetic charge is further developed. The fundamental equations relating magnetic charge and magnetic fields are defined. A better photon model is developed. The ground state vibrations of particles are investigated and indicate that longitudinal waves in the ether are responsible for gravity while vortex waves can explain the Coulomb interaction. Finally, Einstein’s equation for rest mass energy is derived using electromagnetic theory. A laboratory experiment is proposed to validate the theory.

Gravity and magnetic anomalies field characteristics in the South China Sea and its application for interpretation of igneous rocks

Igneous rocks in the South China Sea have broad prospects for oil and gas exploration.Integrated geophysical methods are important approaches to study the distribution of igneous rocks and to determine and identify igneous rock bodies.Aimed at the characteristics of gravity and magnetic fields in the South China Sea,several potential field processing methods are preferentially selected.Reduction to the pole by variable inclinations in the area of low magnetic latitudes is used to perform reduction processing on magnetic anomalies.The preferential continuation method is used to separate gravity and magnetic anomalies and extract the gravity and magnetic anomaly information of igneous rocks in the shallow part of the South China Sea.The 3D spatial equivalent distribution of igneous rocks in South China Sea is illustrated by the 3 D correlation imaging of magnetic anomalies.Since the local anomaly boundaries are highlighted gravity and magnetic gradients,the distribution characters of different igneous rocks are roughly outlined by gravity and magnetic correlation analysis weighted by gradient.The results show the distribution of igneous rocks is controlled and influenced by deep crustal structure and faulting.

Joint inversion of gravity and magnetic data for a two-layer model

Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose a joint gravity and magnetic inversion methodfor two-layer models by concentrating on the relationship between the change of thicknessI and position of the middle layer and anomaly and discuss the effects of the key parameters. Model tests and application to field data show the validity of this method.

SOLIDIFICATION OF Zn－Pb ALLOY UNDER SIMULANT MICROGRAVITY CONDITION OF ORTHOGONAL ELECTRIC AND MAGNETIC FIELDS

ＳＯＬＩＤＩＦＩＣＡＴＩＯＮＯＦＺｎ－ＰｂＡＬＬＯＹＵＮＤＥＲＳＩＭＵＬＡＮＴＭＩＣＲＯＧＲＡＶＩＴＹＣＯＮＤＩＴＩＯＮＯＦＯＲＴＨＯＧＯＮＡＬＥＬＥＣＴＲＩＣＡＮＤＭＡＧＮＥＴＩＣＦＩＥＬＤＳ￥Ｚｈａｏ，Ｊｉｕｚｈｏｕ；Ｇｕｏ，Ｊｉｎｇｊｉｅ；Ｊｉ...

Geological structure of the yellow sea area from regional gravity and magnetic interpretation

Regional gravity and aeromagnetic data covering the area of 32°- 38° N, 118°-127° E at the scale of 1：1,000,000 are coordinated and integrated in a synthetic study of the South China Yellow Sea and adjacent areas. Depth to magnetic crystalline basement and its structure are determined by magnetic anomaly inversion. Depth to and thickness of the Paleozoic rock are also revealed by gravity anomaly inversion with constrains of the basement and known seismic information from several profiles. Structure units, main faults, basin boundaries, and sub-suppressions are outlined on the basis of gravity data interpretation.

Influence of the substitution of Sm, Gd, and Dy for La in La_(0.7)Sr_(0.3)MnO_3 on its magnetic and electric properties and strengthening effect on room-temperature CMR

A series of La0.7-xSmxSr0.3MnO3, La0.7-xGdxSr0.3MnO3, and La0.7-xDyxSr0.3MnO3 （x=0.00, 0.10, 0.20, 0.30） samples were prepared by the solid-state reaction method. The influence of the substitution of Sm, Gd, and Dy for La on the magnetic and electric properties and on the magnetoresistance （MR） was studied through measurements of M-T curves and p-T curves. The results showed that： lattice distortion induced by substitution of Sm, Gd, and Dy for La and extra magnetism of substitution had great influence on the magnetic and electric properties of pcrovskite manganites; substitution of magnetic rare earth element for La was an effective way to change Curie temperature and to strengthen MR in perovskite manganites; and appropriate substitution proportion would generate large MR near room temperature.

Study of Residual Basin and Tectonolayering Based on Airborne Gravity and Magnetic Data

This paper takes the South Yellow Sea as an example to show a new method for comprehensive geological-geophysical research such as residual basin and tectonolayering using airborne gravity and magnetic data in China. Based on airborne gravity and magnetic data, by measuring and analyzing stratigraphic density and susceptibility, the depths to the pre-Sinian magnetic basement top, the pre-Jurassic top and the Cenozoic bottom, are obtained by forward and inverse methods constrained by seismic and drilling data; and furthermore, the residual thicknesses of the Sinian-Triassic, the Jurassic-Cretaceous and the Cenozoic are calculated. Based on airborne gravity and magnetic anomalies, the faults in the pre-Sinian magnetic basement, the Jurassic-Cretaceous and the Cenozoic are respectively interpreted by the qualitative and quantitative methods. On the basis of the above study, and combining regional important tectonic events, four tectonolayers are divided in the vertical succession in South Yellow Sea, namely the pre-Sinian magnetic basement, the Sinian-Triassic, the Jurassic-Cretaceous and the Cenozoic. The result shows that there are thick Cenozoic and Jurassic-Cretaceous strata and thin residual Sinian-Triassic strata in the Suzhong-Huangnan depression area, but there are thin or only sporadic Mesozoic-Cenozoic terrestrial strata and thick Sinian-Triassic marine strata reserved in Subei-Huangzhong uplift area and Sunan-Wunansha uplift area. The four tectonolayers are very different in structures as well as distributions in plane.

Forward Modeling of Gravity,Gravity Gradients,and Magnetic Anomalies due to Complex Bodies

On the basis of the results of improved analytical expression of computation of gravity anomalies due to a homogeneous polyhedral body composed of polygonal facets, and applying the forward theory with the coordinate transformation of vectors and tensors, we deduced both the analytical expressions for gravity gradient tensors and for magnetic anomalies of a polygon, and obtained new analytical expressions for computing vertical gradients of gravity anomalies and vertical component of magnetic anomalies caused by a polyhedral body. And also we developed explicitly the complete unified expressions for the calculation of gravity anomalies, gravity gradient, and magnetic anomalies due to the homogeneous polyhedron. Furthermore, we deduced new analytical expressions for computing vertical gradients of gravity anomalies due to a finite rectangular prism by applying the newly obtained expressions for gravity gradient tensors due to a polyhedral target body. Comparison with forward calculation of models shows the correctness of these new expressions. It will reduce forward calculation time of gravity-magnetic anomalies and improve computational efficiency by applying our unified expressions for joint forward modeling of gravity-magnetic anomalies due to homogeneous polyhedral bodies.

Crustal structure of the western Indian shield: Model based on regional gravity and magnetic data

Regional surface gravity data and global satellite magnetic data have been utilized to generate a preliminary model of the crustal structure along a southwest-northeast profile （Gadra-Fatehpur） through western Rajasthan.The study area represents the western part of the Indian continental landmass which has undergone several major episodes of repeated subduction/collision,plume traces and rifting from Archaean to recent times.The temporal and spatial relationship between the various geotectonic provinces is quite complex,thereby limiting the emergence of a suitable crustal structure model for this region.Exposures of the Malani Igneous Suite （MIS）,a product of bimodal volcanism （～780 Ma）,and considered to be the third largest felsic magmatic province of the world,is evident along the profile and also to the southwest of the study area.The easternmost part of the profile is close to the DAFB （Delhi Aravalli Fold Belt）,a Proterozoic orogenic belt.This study probes the geometry of the different crustal units in terms of density and susceptibility variations in order to decipher the imprints of the major tectonic processes the region has undergone.In order to decipher the crustal geometry of the Gadra-Fatehpur profile,two NW-SE gravity and magnetic profile vertical sections （A-A＇ in the south and B-B＇ in the north） are modelled on the basis of the constraints provided from previous seismic models.The crustal model of the Gadra-Fatehpur profile is composed of alluvium,Tertiary sediments,MIS,Marwar Supergroup,low-density layers （LDLs） and the middle-lower crustal layers,with a distinct change in configuration from the southwest to northeast.The Moho dips from SW to NE,the MIS in the SW gives way to the thick pile of the Marwar Supergroup to the NE.The evolution of MIS has been suggested to have occurred as a consequence of delamination of the upper mantle.LDLs are incorporated in Gadra-Fatehpur model.In the SW,LDL （2550 kg/m3） lies below the MIS in the NE,another LDL （2604 kg/m3） is depicted below the mid-crustal layer.

3D inversion modeling of joint gravity and magnetic data based on a sinusoidal correlation constraint

Joint inversion based on a correlation constraint utilizes a linear correlation function as a structural constraint.The linear correlation function contains a denominator,which may result in a singularity as the objective function is optimized,leading to an unstable inversion calculation.To improve the robustness of this calculation,this paper proposes a new method in which a sinusoidal correlation function is employed as the structural constraint for joint inversion instead of the conventional linear correlation function.This structural constraint does not contain a denominator,thereby preventing a singularity.Compared with the joint inversion method based on a cross-gradient constraint,the joint inversion method based on a sinusoidal correlation constraint exhibits good performance.An application to actual data demonstrates that this method can process real data.

A 3D Geological Model Constrained by Gravity and Magnetic Inversion and its Exploration Implications for the World-class Zhuxi Tungsten Deposit, South China

The Zhuxi tungsten deposit in Jiangxi Province,South China,contains a total W reserve of about 2.86 Mt at an average grade of 0.54 wt%WO3,representing the largest W deposit in the world.Numerous studies on the metallogeny of the deposit have included its timing,the ore-controlling structures and sedimentary host rocks and their implications for mineral exploration.However,the deep nappe structural style of Taqian-Fuchun metallogenic belt that hosts the W deposit,and the spatial shape and scale of deeply concealed intrusions and their sedimentary host rocks are still poorly defined,which seriously restricts the discovery of new deposits at depth and in surrounding areas of the W deposit.Modern 3 D geological modeling is an important tool for the exploration of concealed orebodies,especially in brownfield environments.There are obvious density contrast and weak magnetic contrast in the ore-controlling strata and granite at the periphery of the deposit,which lays a physical foundation for solving the 3 D spatial problems of the ore-controlling geological body in the deep part of the study area through gravity and magnetic modeling.Gravity data(1:50000)and aeromagnetic data(1:50000)from the latest geophysical surveys of 2016-2018 have been used,firstly,to carry out a potential field separation to obtain residual anomalies for gravity and magnetic interactive inversion.Then,on the basis of the analysis of the relationship between physical properties and lithology,under the constraints of surface geology and borehole data,human-computer interactive gravity and magnetic inversion for 18 cross-sections were completed.Finally,the 3 D geological model of the Zhuxi tungsten deposit and its periphery have been established through these 18 sections,and the spatial shape of the intrusions and strata with a depth of 5 km underground were obtained,initially realizing―transparency‖for ore-controlling bodies.According the analysis of the geophysical,geochemical,and geological characteristics of the Zhuxi tungsten deposit,we discern three principles for prospecting and prediction in the research area,and propose five new exploration targets in its periphery.

Integrated gravity and magnetic study on patterns of petroleum basin occurrence in the China seas and adjacent areas

The China seas and adjacent areas contain numerous petroleum basins.One of the main challenges for future oil and gas exploration is to identify the inherent patterns of petroleum basin distribution.The formation and evolution of petroleum basins along with the migration and accumulation of oil and gas are often closely related to the tectonic environment.The gravity and magnetic fields with high lateral resolution and wide coverage provide important data for regional tectonic research.Based on the gravity data in the Global Satellite Gravity Anomaly Database(V31.1)and magnetic data from the Earth Magnetic Anomaly Grid(2-arc-minute resolution)(V2),this study uses integrated gravity and magnetic field technique to obtain integrated gravity and magnetic field result for the China seas and adjacent areas,and then adopts the normalized vertical derivative of the total horizontal derivative technique to conduct partition.Finally,it identifies the relationship between the partition characteristics and tectonics as well as the patterns of petroleum basin occurrence.The research shows that the partition of gravity and magnetic field integrated result has a good correlation with the Neo-Cathaysian tectonic system and tectonic units.The petroleum basins are characterized according to three blocks arranged from north to south and four zones arranged from east to west.The north-south block structure causes the uneven distribution of oil and gas resources in the mainland area and the differences in the hydrocarbon-bearing strata.Petroleum basins are more abundant in the north than in the south.The ages of the main oil-and gas-bearing strata are“Paleozoic–Mesozoic,Paleozoic–Mesozoic–Cenozoic,and Paleozoic–Mesozoic”,in order from north to south.The difference in the overall type of oil and gas resources in all basins is controlled by the east–west zonation.From east to west,the oil and gas resource type exhibits a wave-like pattern of“oil and gas,gas,oil and gas,gas”.The vertical distribution is characterized by an upper oil(Mesozoic–Cenozoic)and lower gas(Mesozoic–Paleozoic)structure.Within the study area,the Paleozoic marine strata should be the main strata of future natural gas exploration.

On the Tectonic Problems of the Southern East China Sea and Adjacent Regions:Evidence from Gravity and Magnetic Data

In this paper, two sets of gravity and magnetic data were used to study the tectonics of the southern East China Sea and Ryukyu trench-arc system： one data set was from the ＇Geological-geophysical map series of China Seas and adjacent areas＇ database and the other was newly collected by R/VKexue Ⅲ in 2011. Magnetic and gravity data were reorganized and processed using the software MMDP, MGDP and RGIS. In addition to the description of the anomaly patterns in different areas, deep and shallow structure studies were performed by using several kinds of calculation, including a spectrum analysis, upward-continuation of the Bouguer anomaly and horizontal derivatives of the total-field magnetic anomaly. The depth of the Moho and magnetic basement were calculated. Based on the above work, several controversial tectonic problems were discussed. Compared to the shelf area and Ryukyu Arc, the Okinawa Trough has an obviously thinned crust, with the thinnest area having thickness less than 14km in the southern part. The Taiwan-Sinzi belt, which terminates to the south by the NW-SE trending Miyako fault belt, contains the relic volcanic arc formed by the splitting of the paleo Ryukyu volcanic arc as a result of the opening of the Okinawa Trough. As an important tectonic boundary, the strike-slip type Miyako fault belt extends northwestward into the shelf area and consists of several discontinuous segments. A forearc terrace composed of an exotic terrane collided with the Ryukyu Arc following the subduction of the Philippine Sea Plate. Mesozoic strata of varying thicknesses exist beneath the Cenozoic strata in the shelf basin and significantly influence the magnetic pattern of this area. The gravity and magnetic data support the existence of a Great East China Sea, which suggests that the entire southern East China Sea shelf area was a basin in the Mesozoic without alternatively arranged uplifts and depressions, and might have extended southwestward and connected with the northern South China Sea shelf basin.

Rock magnetic property of gravity core CSH1 from the northern Okinawa Trough and the effect of early diagenesis

Detailed rock magnetic investigations were undertaken at 2 -4 cm interval for the gravity core CSHI （with a length of 17.36 m） from the northern Okinawa Trough. Time-scale of the core was constructed by two characteristic tephras and foraminferal assemblages, indicating an age of 50 ka for the bottom of the core. Except for three tephras and abrupt decrease in surface, there are little changes in all kinds of rock magnetic parameters that can be correlated to the climate change for the last 50 ka. Different from the common sediments, most S-ratios （S equals the negative ratio of IRM-0.3T to SIRM, which is an indicator of low coercivity content） of the sediments are smaller than 0. 9, which implies a substantial amount of magnetic minerals with high coercivity. The existence of iron sulphide （ greigite or pyrrhotite） is revealed by representative susceptibility - temperature curves showing 200 350 ℃ Curie temperature in addition to 580 ℃ of magnetite, and also by awful smell during heating and dark heating products. Both the occurrence of authigenic iron sulphide and quick decrease of magnetic parameters near the surface clearly show that sediments from Core CSHI have undergone early diagenesis. The featureless magnetic changes of the whole core except for three tephras mean that the post-depositonal alteration is so strong that most original signals have been destroyed. For the same reason, the organic matter in sediment and sulphate in pore water must have been consumed along with dissolution, precipitation of iron and manganese happening sequencially during the redox reaction series. Great caution must be taken when using these altered chemical parameters for the interpretation of climatic changes.

Basin structure and multiresource potential based on high-precision airborne gravity and magnetic data

In this paper,we use high-precision airborne gravity and magnetic data to study the geophysical characteristics of the western slope of the Songliao Basin and its adjacent areas and evaluate the resource potential.We performed an in-depth analysis of three aspects of the basin characteristics:the characteristics of residual strata,the development characteristics of faults,and the distribution characteristics of magmatic rocks.Next,we analyzed the forming background of organic(oil and gas)resources and inorganic(uranium ore and hot dry rock)resources.The results showed that the new Upper Paleozoic strata have significant differences in different regions of the study area(with a thickness of 0–8000 m),mainly distributed in the eastern and northern regions but absent in the middle eastern and western regions.Furthermore,the thickness and depth of the Mesozoic layer varied between the eastern and western regions;it was thicker and deeper in the middle eastern region but thinner and shallower in the western region,and it is absent in most western regions.The main faults in the region are in the north–northeast(NNE)direction.Faults in the NE–NNE and NW directions jointly controlled the morphology of the secondary structural units.Magmatic rocks are relatively developed in the middle and eastern parts of the region.Most magmatic rocks are distributed along the faults and their sides,clearly reflecting the control of the faults on magmatic activities.The western slope of the Songliao basin and its surroundings have a favorable geological setting for the accumulation(mineralization)of oil,gas,shale oil,hot dry rock,and uranium ore.It is conducive to oil and gas exploration of deep new strata and collaborative exploration of multiple resources.

Evaluation and extraction of weak gravity and magnetic anomalies

In this paper, I introduce what are called weak gravity and magnetic anomalies and propose standards for estimating their reliability. I also introduce new techniques for processing this kind of weak anomaly. These techniques consist of interference elimination and weak signal extraction. Practical applications have proved their effectiveness, Weak gravity and magnetic anomalies will get more attention with the development of targeted exploration.

Structural characteristics of the KPR-CBR triple-junction inferred from gravity and magnetic interpretations,Philippine Sea Plate

The intersection of the Kyushu-Palau Ridge(KPR)and the Central Basin Rift(CBR)of the West Philippine Basin(WPB)is a relic of a trench-trench-rift(TTR)type triple-junction,which preserves some pivotal information on the cessation of the seafloor spreading of the WPB,the emplacement and disintegration of the proto-Izu-Bonin-Mariana(IBM)Arc,and the transition from initial rifting to steady-state spreading of the Parece Vela Basin(PVB).However,the structural characteristics of this triple-junction have not been thoroughly understood.In this paper,using the newly acquired multi-beam bathymetric,gravity,and magnetic data obtained by the Qingdao Institute of Marine Geology,China Geological Survey,the authors depict the topographic,gravity,and magnetic characteristics of the triple-junction and adjacent region.Calculations including the upward continuations and total horizontal derivatives of gravity anomaly are also performed to highlight the major structural features and discontinuities.Based on these works,the morphological and structural features and their formation mechanisms are analyzed.The results show that the last episode amagmatic extension along the CBR led to the formation of a deep rift valley,which extends eastward and incised the KPR.The morphological and structural fabrics of the KPR near and to the south of the triple-junction are consistent with those of the western PVB,manifesting as a series of NNE-SSW-and N-S-trending ridges and troughs,which were produced by the extensional faults associated with the initial rifting of the PVB.The superposition of the above two reasons induced the prominent discontinuity of the KPR in deep and shallow crustal structures between 15°N‒15°30′N and 13°30′N‒14°N.Combined with previous authors’results,we propose that the stress produced by the early spreading of the PVB transmitted westward and promoted the final stage amagmatic extension of the CBR.The eastward propagation of the CBR destroyed the KPR,of which the magmatism had decayed or ceased at that time.The destruction mechanism of the KPR associated with the rifting of the PVB varies along strike the KPR.Adjacent to the triple-junction,the KPR was destroyed mainly due to the oblique intersection of the PVB rifting center.Whereas south of the triple-junction,the KPR was destroyed by the E-W-directional extensional faulting on its whole width.

Monazite Recovery by Magnetic and Gravity Separation of Medium Grade Zircon Concentrate from Senegalese Heavy Mineral Sands Deposit

Gravity, magnetic and electrostatic separation methods allowed to obtain different titanium oxide concentrates (ilmenite, leucoxene, rutile) and different varieties of zircon concentrates (premium zircon, standard zircon, medium grade zircon standard) from Senegal’s heavy mineral sands. During mining separation, monazite, which is a paramagnetic mineral, was found in a non-negligible concentration of 0.57 wt% on average in the medium grade zircon standard which also contains 37.96 wt% zircon and 44.46 wt% titanium oxides. Magnetic and gravity separation tests were carried out on the Medium grade zircon standard (MGZS) to produce a monazite concentrate at Eramet Ideas laboratory. Magnetic separation at 1.5 teslas intensity resulted in the recovery of 94.8% of the monazite from the MGZS. Gravity separation also recovered 76.6% of the monazite from the MGZS. The combination of these two treatment methods can thus produce three concentrates from MGZS (a monazite concentrate, a zircon concentrate, and a titanium oxide concentrate).

Simultaneous Modelling of Gravity and Magnetic Data in a Measured Heat Flux Area to Characterize Geothermal Heat Sources: A Case for Eburru Geothermal Complex, Kenya

Forward modelling of gravity and magnetic data was done simultaneously to show the correlation between gravity and magnetic anomalies on a measured heat flux region. The results were used to characterize the heat source structures in Eburru area. Modelling was done using Oasis montaj geosoft software which is an iteration process where the gravity and magnetic anomalies were calculated and compared to the observed residual anomaly until there was a fit. The start model was constructed based on depths from Euler deconvolution and models constrained using stratigraphy data from the existing wells in the study area. Forward modelling of gravity and magnetic data revealed intrusions within the Earth’s subsurface with depth to the top of the sources ranging from 739 m to 5811 m. The density of the sources ranges between 3.0 g/cm3 and 3.2 g/cm3 while their magnetic susceptibility was zero. This implies that intrusions from the mantle with a magnetic susceptibility of zero have temperatures exceeding the curie temperature of rocks. The density of the intrusions modelled was higher than 2.67 g/cm3, the average crustal density, hence it explains the observed positive gravity anomaly. The results also revealed that areas with high heat flux have shallow heat sources and if the heat sources are deep, then there must be a good heat transfer mechanism to the surface.