We propose a novel workflow for fast forward modeling of well logs in axially symmetric 2D models of the nearwellbore environment.The approach integrates the finite element method with deep residual neural networks to...We propose a novel workflow for fast forward modeling of well logs in axially symmetric 2D models of the nearwellbore environment.The approach integrates the finite element method with deep residual neural networks to achieve exceptional computational efficiency and accuracy.The workflow is demonstrated through the modeling of wireline electromagnetic propagation resistivity logs,where the measured responses exhibit a highly nonlinear relationship with formation properties.The motivation for this research is the need for advanced modeling al-gorithms that are fast enough for use in modern quantitative interpretation tools,where thousands of simulations may be required in iterative inversion processes.The proposed algorithm achieves a remarkable enhancement in performance,being up to 3000 times faster than the finite element method alone when utilizing a GPU.While still ensuring high accuracy,this makes it well-suited for practical applications when reliable payzone assessment is needed in complex environmental scenarios.Furthermore,the algorithm’s efficiency positions it as a promising tool for stochastic Bayesian inversion,facilitating reliable uncertainty quantification in subsurface property estimation.展开更多
Vehicle electrification,an important method for reducing carbon emissions from road transport,has been promoted globally.In this study,we analyze how individuals adapt to this transition in transportation and its subs...Vehicle electrification,an important method for reducing carbon emissions from road transport,has been promoted globally.In this study,we analyze how individuals adapt to this transition in transportation and its subsequent impact on urban structure.Considering the varying travel costs associated with electric and fuel vehicles,we analyze the dynamic choices of households concerning house locations and vehicle types in a two-dimensional monocentric city.A spatial equilibrium is developed to model the interactions between urban density,vehicle age and vehicle type.An agent-based microeconomic residential choice model dynamically coupled with a house rent market is developed to analyze household choices of home locations and vehicle energy types,considering vehicle ages and competition for public charging piles.Key findings from our proposed models show that the proportion of electric vehicles(EVs)peaks at over 50%by the end of the first scrappage period,accompanied by more than a 40%increase in commuting distance and time compared to the scenario with only fuel vehicles.Simulation experiments on a theoretical grid indicate that heterogeneity-induced residential segregation can lead to urban sprawl and congestion.Furthermore,households with EVs tend to be located farther from the city center,and an increase in EV ownership contributes to urban expansion.Our study provides insights into how individuals adapt to EV transitions and the resulting impacts on home locations and land use changes.It offers a novel perspective on the dynamic interactions between EV adoption and urban development.展开更多
A fast and high precision spatial domain algorithm is presented for forward modeling of two-dimensional(2D)body gravity anomalies of arbitrary shape and density distribution.The new algorithm takes advantage of the co...A fast and high precision spatial domain algorithm is presented for forward modeling of two-dimensional(2D)body gravity anomalies of arbitrary shape and density distribution.The new algorithm takes advantage of the convolution properties of the expression for 2D gravity anomalies,uses a rectangular cell as a grid subdivision unit,and then 2D bodies with irregular cross-sections are approximated by a combination of 2D bodies with a rectangular cross section.The closed-form expression is used to calculate the gravitational anomalies of the combination of 2D bodies with a rectangular cross section.To improve computing effi ciency,the new algorithm uses a fast algorithm for the implementation of the Toeplitz matrix and vector multiplication.The synthetic 2D models with rectangular and circular cross-sections and constant and variable densities are designed to evaluate the computational accuracy and speed of the new algorithm.The experiment results show that the computation costs less than 6 s for a grid subdivision with 10000×10000 elements.Compared to the traditional forward modeling methods,the proposed method significantly improved computational effi ciency while guaranteeing computational accuracy.展开更多
From June 2008 to August 2013,approximately 67 kt of CO_(2) was injected into a deep saline formation at the Ketzin pilot CO_(2) storage site.During injection,3D seismic surveys have been performed to monitor the migr...From June 2008 to August 2013,approximately 67 kt of CO_(2) was injected into a deep saline formation at the Ketzin pilot CO_(2) storage site.During injection,3D seismic surveys have been performed to monitor the migration of sequestered CO_(2).Seismic monitoring results are limited by the acquisition and signal-to-noise ratio of the acquired data.The multiphysical reservoir simulation provides information regarding the CO_(2) fluid behavior,and the approximated model should be calibrated with the monitoring results.In this work,property models are delivered from the multiphysical model during 3D repeated seismic surveys.The simulated seismic data based on the models are compared with the real data,and the results validate the effectiveness of the multiphysical inversion method.Time-lapse analysis shows the trend of CO_(2) migration during and after injection.展开更多
The Earth's magnetic field,which has been extensively observed from ground to satellite altitudes over several decades,originates from multiple sources,such as the core dynamo,the conductive mantle,the magnetized ...The Earth's magnetic field,which has been extensively observed from ground to satellite altitudes over several decades,originates from multiple sources,such as the core dynamo,the conductive mantle,the magnetized lithosphere,and the space current systems.Modeling of the lithospheric contribution plays an important role in the geophysical studies and industrial applications.In this paper,we propose a new method for global and regional modeling of the lithospheric magnetic field based on the cubed-sphere.An equivalent dipole source method on a quasi-uniform cubed-sphere grid is employed in the forward modeling.The dipole directions are fixed according to a priori magnetization and the relative intensities are estimated by an inversion procedure of least-squares fitting with minimum model regularization.Several numerical tests are performed to validate the accuracy and efficiency of both forward modeling and inversion procedure.The proposed method is applied to the global and regional modeling based on the latest magnetic data from Swarm Alpha satellite and MSS-1 mission.The model results indicate that the proposed method works quite well for realistic satellite data and MSS-1 data is consistent with the Swarm data in terms of lithospheric field modeling.展开更多
Unmanned aerial vehicle transient electromagnetic(UAV-TEM)is a novel airborne exploration method that offers advantages such as low cost,simple operation,high exploration efficiency and suitability for near-surface ex...Unmanned aerial vehicle transient electromagnetic(UAV-TEM)is a novel airborne exploration method that offers advantages such as low cost,simple operation,high exploration efficiency and suitability for near-surface exploration in complex terrain areas.To improve the accuracy of data interpretation in this method,the authors conducted a systematic three-dimensional(3D)forward modeling and inversion of the UAV-TEM.This study utilized the finite element method based on unstructured tetrahedral elements and employed the second-order backward Euler method for time discretization.This allowed for accurate 3D modeling and accounted for the effects of complex terrain.Based on these,the influence characteristics of flight altitudes and the sizes,burial depths,and resistivities of anomalies are compared and analyzed to explore the UAV-TEM systems’exploration capability.Lastly,four typical geoelectrical models of landslides are designed,and the inversion method based on the Gauss-Newton optimization method is used to image the landslide models and analyze the imaging effect of the UAV-TEM method on landslide geohazards.Numerical results showed that UAV-TEM could have better exploration resolution and fine imaging of nearsurface structures,providing important technical support for monitoring,early warning,and preventing landslides and other geological hazards.展开更多
The superparamagnetic effect arises from the superparamagnetism exhibited by a multitude of nano-sized magnetic mineral particles under an external electric field.This phenomenon manifests in transient electromagnetic...The superparamagnetic effect arises from the superparamagnetism exhibited by a multitude of nano-sized magnetic mineral particles under an external electric field.This phenomenon manifests in transient electromagnetic data primarily as a deceleration in the attenuation rate of late-stage signals,a characteristic difficult to discern directly from airborne transient electromagnetic signals,consequently leading to significant misinterpretations of subterranean electrical structures.This study embarks on 3D forward modeling of airborne electromagnetic responses in the frequency domain,accounting for the superparamagnetic effect,utilizing an unstructured finite element method.Superparamagnetic responses in the time domain were obtained through frequency-time conversion.This investigation explores the influence of various parameters-such as magnetic susceptibility,time constants,and flight altitude-on the superparamagnetic effect by examining the response characteristics of typical targets.Findings indicate that in its late stages,the superparamagnetic effect can induce a relative anomaly of up to 300%.There is a positive correlation between magnetic susceptibility and the strength of the superparamagnetic effect.The influence of the time constant's upper and lower limits on the superparamagnetic effect is minimal;however,the range between these limits significantly affects the effect,showing a negative correlation with its intensity.Higher flight altitudes weaken the superparamagnetic signal.The impact is most pronounced when superparamagnetic minerals are shallowly buried,effectively shielding the underlying geology with the characteristics of a good conductivity anomaly,but this effect diminishes with greater depth.The insights from this study provide a theoretical framework for a deeper understanding of the superparamagnetic effect in transient electromagnetic signals and for more accurate interpretations of subterranean geological and electrical structures.展开更多
The mapping method is a forward-modeling method that transforms the irregular surface to horizontal by mapping the rectangular grid as curved; moreover, the wave field calculations move from the physical domain to the...The mapping method is a forward-modeling method that transforms the irregular surface to horizontal by mapping the rectangular grid as curved; moreover, the wave field calculations move from the physical domain to the calculation domain. The mapping method deals with the irregular surface and the low-velocity layer underneath it using a fine grid. For the deeper high-velocity layers, the use of a fine grid causes local oversampling. In addition, when the irregular surface is transformed to horizontal, the flattened interface below the surface is transformed to curved, which produces inaccurate modeling results because of the presence of ladder-like burrs in the simulated seismic wave. Thus, we propose the mapping method based on the dual-variable finite-difference staggered grid. The proposed method uses different size grid spacings in different regions and locally variable time steps to match the size variability of grid spacings. Numerical examples suggest that the proposed method requires less memory storage capacity and improves the computational efficiency compared with forward modeling methods based on the conventional grid.展开更多
We investigated the effect of microscopic distribution modes of hydrates in porous sediments, and the saturation of hydrates and free gas on the elastic properties of saturated sediments. We simulated the propagation ...We investigated the effect of microscopic distribution modes of hydrates in porous sediments, and the saturation of hydrates and free gas on the elastic properties of saturated sediments. We simulated the propagation of seismic waves in gas hydrate-bearing sediments beneath the seafloor, and obtained the common receiver gathers of compressional waves（P-waves） and shear waves（S-waves）. The numerical results suggest that the interface between sediments containing gas hydrates and free gas produces a large-amplitude bottomsimulating reflector. The analysis of multicomponent common receiver data suggests that ocean-bottom seismometers receive the converted waves of upgoing P- and S-waves, which increases the complexity of the wavefield record.展开更多
Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell＇s equations for an anisotropic mo...Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell＇s equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic （EM） fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.展开更多
AVO forward modeling is based on two-phase medium theory and is considered an effective method for describing reservoir rocks and fluids. However, the method depends on the input matrix mineral bulk modulus and the ra...AVO forward modeling is based on two-phase medium theory and is considered an effective method for describing reservoir rocks and fluids. However, the method depends on the input matrix mineral bulk modulus and the rationality of the two-phase medium model. We used the matrix mineral bulk modulus inversion method and multiple constraints to obtain a two-phase medium model with physical meaning. The proposed method guarantees the reliability of the obtained AVO characteristicsin two-phase media. By the comparative analysis of different lithology of the core sample, the advantages and accuracy of the inversion method can be illustrated. Also, the inversion method can be applied in LH area, and the AVO characteristics can be obtained when the porosity, fluid saturation, and other important lithology parameters are changed. In particular, the reflection coefficient amplitude difference between the fast P wave and S wave as a function of porosity at the same incidence angle, and the difference in the incidence angle threshold can be used to decipher porosity.展开更多
A novel behavioral model using three-layer time-delay feed-forward neural networks （TDFFNN）is adopted to model radio frequency （RF）power amplifiers exhibiting memory nonlinearities. In order to extract the paramet...A novel behavioral model using three-layer time-delay feed-forward neural networks （TDFFNN）is adopted to model radio frequency （RF）power amplifiers exhibiting memory nonlinearities. In order to extract the parameters, the back- propagation algorithm is applied to train the proposed neural networks. The proposed model is verified by the typical odd- order-only memory polynomial model in simulation, and the performance is compared with different numbers of taped delay lines（TDLs） and perceptrons of the hidden layer. For validating the TDFFNN model by experiments, a digital test bench is set up to collect input and output data of power amplifiers at a 60 × 10^6 sample/s sampling rate. The 3.75 MHz 16-QAM signal generated in the vector signal generator（VSG） is chosen as the input signal, when measuring the dynamic AM/AM and AM/PM characteristics of power amplifiers. By comparisons and analyses, the presented model provides a good performance in convergence, accuracy and efficiency, which is approved by simulation results and experimental results in the time domain and frequency domain.展开更多
We discuss the feasibility of using controlled-source electromagnetic （CSEM） in the frequency domain for prospecting marine gas hydrates. Based on the Ocean Drilling Program （ODP） Leg 164 log data, we have establi...We discuss the feasibility of using controlled-source electromagnetic （CSEM） in the frequency domain for prospecting marine gas hydrates. Based on the Ocean Drilling Program （ODP） Leg 164 log data, we have established several 1-D resistivity models which have different gas hydrate concentrations. Meanwhile, we analyzed the electromagnetic response of marine gas hydrates in the frequency domain based on these models. We also studied the relationship between electrical field magnitude or phase and parameters such as receiver-transmitter distance and frequency. Our numerical modeling results provide us with a quantitative reference for exploration and resource evaluation of marine gas hydrates.展开更多
Finite-difference(FD) methods are widely used in seismic forward modeling owing to their computational efficiency but are not readily applicable to irregular topographies. Thus, several FD methods based on the transfo...Finite-difference(FD) methods are widely used in seismic forward modeling owing to their computational efficiency but are not readily applicable to irregular topographies. Thus, several FD methods based on the transformation to curvilinear coordinates using body-fitted grids have been proposed, e.g., stand staggered grid(SSG) with interpolation, nonstaggered grid, rotated staggered grid(RSG), and fully staggered. The FD based on the RSG is somewhat superior to others because it satisfies the spatial distribution of the wave equation without additional memory and computational requirements; furthermore, it is simpler to implement. We use the RSG FD method to transform the firstorder stress–velocity equation in the curvilinear coordinates system and introduce the highprecision adaptive, unilateral mimetic finite-difference(UMFD) method to process the freeboundary conditions of an irregular surface. The numerical results suggest that the precision of the solution is higher than that of the vacuum formalism. When the minimum wavelength is low, UMFD avoids the surface wave dispersion. We compare FD methods based on RSG, SEM, and nonstaggered grid and infer that all simulation results are consistent but the computational efficiency of the RSG FD method is higher than the rest.展开更多
Water quality models are important tools to support the optimization of aquatic ecosystem rehabilitation programs and assess their efficiency. Basing on the flow conditions of the Daqinghe River Mouth of the Dianchi L...Water quality models are important tools to support the optimization of aquatic ecosystem rehabilitation programs and assess their efficiency. Basing on the flow conditions of the Daqinghe River Mouth of the Dianchi Lake, China, a two-dimensional water quality model was developed in the research. The hydrodynamics module was numerically solved by the alternating direction iteration （ADI） method. The parameters of the water quality module were obtained through the in situ experiments and the laboratory analyses that were conducted from 2006 to 2007. The model was calibrated and verified by the observation data in 2007. Among the four modelled key variables, i.e., water level, COD （in CODcr）, NH4＋-N and PO43-P the minimum value of the coefficient of determination （COD） was 0.69, indicating the model performed reasonably well. The developed model was then applied to simulate the water quality changes at a downstream cross-section assuming that the designed restoration programs were implemented. According to the simulated results, the restoration programs could cut down the loads of COD and PO43-P about 15%. Such a load reduction, unfortunately, would have very little effect on the NH4^＋-N removal. Moreover, the water quality at the outlet cross-section would be still in class V （3838-02）, indicating more measures should be taken to further reduce the loads. The study demonstrated the capability of water quality models to support aquatic ecosystem restorations.展开更多
The workload of the 3D magnetotelluric forward modeling algorithm is so large that the traditional serial algorithm costs an extremely large compute time. However, the 3D forward modeling algorithm can process the dat...The workload of the 3D magnetotelluric forward modeling algorithm is so large that the traditional serial algorithm costs an extremely large compute time. However, the 3D forward modeling algorithm can process the data in the frequency domain, which is very suitable for parallel computation. With the advantage of MPI and based on an analysis of the flow of the 3D magnetotelluric serial forward algorithm, we suggest the idea of parallel computation and apply it. Three theoretical models are tested and the execution efficiency is compared in different situations. The results indicate that the parallel 3D forward modeling computation is correct and the efficiency is greatly improved. This method is suitable for large size geophysical computations.展开更多
Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-d...Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-dimensional interior ballistics model of the solid propellant ETC gun(2D-IB-SPETCG)is presented to describe the process of the ETC launch.Both calculated pressure and projectile muzzle velocity accord well with the experimental results.The feasibility of the 2D-IB-SPETCG model is proved.Depending on the experimental data and initial parameters,detailed distribution of the ballistics parameters can be simulated.With the distribution of pressure and temperature of the gas phase and the propellant,the influence of plasma during the ignition process can be analyzed.Because of the radial flowing plasma,the propellant in the area of the DRPG is ignited within 0.01 ms,while all propellant in the chamber is ignited within 0.09 ms.The radial ignition delay time is much less than the axial delay time.During the ignition process,the radial pressure difference is less than 5 MPa at the place 0.025 m away from the breech.The radial ignition uniformity is proved.The temperature of the gas increases from several thousand K(conventional ignition)to several ten thousand K(plasma ignition).Compare the distribution of the density and temperature of the gas,we know that low density and high temperature gas appears near the exits of the DRPG,while high density and low temperature gas appears at the wall near the breech.The simulation of the 2D-IB-SPETCG model is an effective way to investigate the interior ballistics process of the ETC launch.The 2D-IB-SPETC model can be used for prediction and improvement of experiments.展开更多
Mathematical geophone （MG） and equal-time stacking （ETS） principles are used to implement seismic prestack forward modeling with irregular surfaces using the oneway acoustic wave-equation. This method receives sei...Mathematical geophone （MG） and equal-time stacking （ETS） principles are used to implement seismic prestack forward modeling with irregular surfaces using the oneway acoustic wave-equation. This method receives seismic primary reflections from the subsurface using a set of virtual MGs. The receivers can be located anywhere on an irregular observing surface. Moreover, the ETS method utilizes the one-way acoustic wave equation to easily and quickly image and extrapolate seismic reflection data. The method is illustrated using high single-noise ratio common shot gathers computed by numerical forward modeling of two simple models, one with a flat surface and one with an irregular surface, and a complex normal fault model. A prestack depth migration method for irregular surface topography was used to reoroduce the normal fault model with high accuracy.展开更多
基金financially supported by the Russian federal research project No.FWZZ-2022-0026“Innovative aspects of electro-dynamics in problems of exploration and oilfield geophysics”.
文摘We propose a novel workflow for fast forward modeling of well logs in axially symmetric 2D models of the nearwellbore environment.The approach integrates the finite element method with deep residual neural networks to achieve exceptional computational efficiency and accuracy.The workflow is demonstrated through the modeling of wireline electromagnetic propagation resistivity logs,where the measured responses exhibit a highly nonlinear relationship with formation properties.The motivation for this research is the need for advanced modeling al-gorithms that are fast enough for use in modern quantitative interpretation tools,where thousands of simulations may be required in iterative inversion processes.The proposed algorithm achieves a remarkable enhancement in performance,being up to 3000 times faster than the finite element method alone when utilizing a GPU.While still ensuring high accuracy,this makes it well-suited for practical applications when reliable payzone assessment is needed in complex environmental scenarios.Furthermore,the algorithm’s efficiency positions it as a promising tool for stochastic Bayesian inversion,facilitating reliable uncertainty quantification in subsurface property estimation.
基金supported by National Natural Science Foundation of China(72288101,72361137002,and 72101018)the Dutch Research Council(NWO Grant 482.22.01).
文摘Vehicle electrification,an important method for reducing carbon emissions from road transport,has been promoted globally.In this study,we analyze how individuals adapt to this transition in transportation and its subsequent impact on urban structure.Considering the varying travel costs associated with electric and fuel vehicles,we analyze the dynamic choices of households concerning house locations and vehicle types in a two-dimensional monocentric city.A spatial equilibrium is developed to model the interactions between urban density,vehicle age and vehicle type.An agent-based microeconomic residential choice model dynamically coupled with a house rent market is developed to analyze household choices of home locations and vehicle energy types,considering vehicle ages and competition for public charging piles.Key findings from our proposed models show that the proportion of electric vehicles(EVs)peaks at over 50%by the end of the first scrappage period,accompanied by more than a 40%increase in commuting distance and time compared to the scenario with only fuel vehicles.Simulation experiments on a theoretical grid indicate that heterogeneity-induced residential segregation can lead to urban sprawl and congestion.Furthermore,households with EVs tend to be located farther from the city center,and an increase in EV ownership contributes to urban expansion.Our study provides insights into how individuals adapt to EV transitions and the resulting impacts on home locations and land use changes.It offers a novel perspective on the dynamic interactions between EV adoption and urban development.
基金This work is jointly sponsored by the National Natural Science Foundation of China(No.41404106)the Scientific Research Startup Fund for Doctoral Program of Guilin University of Technology,Guangxi Natural Science Foundation Program(No.2018GXNSFBA138049)Guangxi Natural Science Foundation Program for Innovation Research Team(No.2016GXNSFGA380004).
文摘A fast and high precision spatial domain algorithm is presented for forward modeling of two-dimensional(2D)body gravity anomalies of arbitrary shape and density distribution.The new algorithm takes advantage of the convolution properties of the expression for 2D gravity anomalies,uses a rectangular cell as a grid subdivision unit,and then 2D bodies with irregular cross-sections are approximated by a combination of 2D bodies with a rectangular cross section.The closed-form expression is used to calculate the gravitational anomalies of the combination of 2D bodies with a rectangular cross section.To improve computing effi ciency,the new algorithm uses a fast algorithm for the implementation of the Toeplitz matrix and vector multiplication.The synthetic 2D models with rectangular and circular cross-sections and constant and variable densities are designed to evaluate the computational accuracy and speed of the new algorithm.The experiment results show that the computation costs less than 6 s for a grid subdivision with 10000×10000 elements.Compared to the traditional forward modeling methods,the proposed method significantly improved computational effi ciency while guaranteeing computational accuracy.
基金supported by the National Natural Science Foundation of China(Grant No.42025403)the Youth Innovation Promotion Association,Chinese Academy of Sciences(Grant No.2023074).
文摘From June 2008 to August 2013,approximately 67 kt of CO_(2) was injected into a deep saline formation at the Ketzin pilot CO_(2) storage site.During injection,3D seismic surveys have been performed to monitor the migration of sequestered CO_(2).Seismic monitoring results are limited by the acquisition and signal-to-noise ratio of the acquired data.The multiphysical reservoir simulation provides information regarding the CO_(2) fluid behavior,and the approximated model should be calibrated with the monitoring results.In this work,property models are delivered from the multiphysical model during 3D repeated seismic surveys.The simulated seismic data based on the models are compared with the real data,and the results validate the effectiveness of the multiphysical inversion method.Time-lapse analysis shows the trend of CO_(2) migration during and after injection.
基金supported by the National Natural Science Foundation of China(42250101,42250102,42250103,12250013)the Macao Foundation。
文摘The Earth's magnetic field,which has been extensively observed from ground to satellite altitudes over several decades,originates from multiple sources,such as the core dynamo,the conductive mantle,the magnetized lithosphere,and the space current systems.Modeling of the lithospheric contribution plays an important role in the geophysical studies and industrial applications.In this paper,we propose a new method for global and regional modeling of the lithospheric magnetic field based on the cubed-sphere.An equivalent dipole source method on a quasi-uniform cubed-sphere grid is employed in the forward modeling.The dipole directions are fixed according to a priori magnetization and the relative intensities are estimated by an inversion procedure of least-squares fitting with minimum model regularization.Several numerical tests are performed to validate the accuracy and efficiency of both forward modeling and inversion procedure.The proposed method is applied to the global and regional modeling based on the latest magnetic data from Swarm Alpha satellite and MSS-1 mission.The model results indicate that the proposed method works quite well for realistic satellite data and MSS-1 data is consistent with the Swarm data in terms of lithospheric field modeling.
基金Supported by Key Research and Development Project of Guangxi Pr ovince(No.AB21196028).
文摘Unmanned aerial vehicle transient electromagnetic(UAV-TEM)is a novel airborne exploration method that offers advantages such as low cost,simple operation,high exploration efficiency and suitability for near-surface exploration in complex terrain areas.To improve the accuracy of data interpretation in this method,the authors conducted a systematic three-dimensional(3D)forward modeling and inversion of the UAV-TEM.This study utilized the finite element method based on unstructured tetrahedral elements and employed the second-order backward Euler method for time discretization.This allowed for accurate 3D modeling and accounted for the effects of complex terrain.Based on these,the influence characteristics of flight altitudes and the sizes,burial depths,and resistivities of anomalies are compared and analyzed to explore the UAV-TEM systems’exploration capability.Lastly,four typical geoelectrical models of landslides are designed,and the inversion method based on the Gauss-Newton optimization method is used to image the landslide models and analyze the imaging effect of the UAV-TEM method on landslide geohazards.Numerical results showed that UAV-TEM could have better exploration resolution and fine imaging of nearsurface structures,providing important technical support for monitoring,early warning,and preventing landslides and other geological hazards.
文摘The superparamagnetic effect arises from the superparamagnetism exhibited by a multitude of nano-sized magnetic mineral particles under an external electric field.This phenomenon manifests in transient electromagnetic data primarily as a deceleration in the attenuation rate of late-stage signals,a characteristic difficult to discern directly from airborne transient electromagnetic signals,consequently leading to significant misinterpretations of subterranean electrical structures.This study embarks on 3D forward modeling of airborne electromagnetic responses in the frequency domain,accounting for the superparamagnetic effect,utilizing an unstructured finite element method.Superparamagnetic responses in the time domain were obtained through frequency-time conversion.This investigation explores the influence of various parameters-such as magnetic susceptibility,time constants,and flight altitude-on the superparamagnetic effect by examining the response characteristics of typical targets.Findings indicate that in its late stages,the superparamagnetic effect can induce a relative anomaly of up to 300%.There is a positive correlation between magnetic susceptibility and the strength of the superparamagnetic effect.The influence of the time constant's upper and lower limits on the superparamagnetic effect is minimal;however,the range between these limits significantly affects the effect,showing a negative correlation with its intensity.Higher flight altitudes weaken the superparamagnetic signal.The impact is most pronounced when superparamagnetic minerals are shallowly buried,effectively shielding the underlying geology with the characteristics of a good conductivity anomaly,but this effect diminishes with greater depth.The insights from this study provide a theoretical framework for a deeper understanding of the superparamagnetic effect in transient electromagnetic signals and for more accurate interpretations of subterranean geological and electrical structures.
基金financially supported by the National Natural Science Foundation of China(Nos.41104069 and 41274124)the National 973 Project(Nos.2014CB239006 and 2011CB202402)+1 种基金the Shandong Natural Science Foundation of China(No.ZR2011DQ016)Fundamental Research Funds for Central Universities(No.R1401005A)
文摘The mapping method is a forward-modeling method that transforms the irregular surface to horizontal by mapping the rectangular grid as curved; moreover, the wave field calculations move from the physical domain to the calculation domain. The mapping method deals with the irregular surface and the low-velocity layer underneath it using a fine grid. For the deeper high-velocity layers, the use of a fine grid causes local oversampling. In addition, when the irregular surface is transformed to horizontal, the flattened interface below the surface is transformed to curved, which produces inaccurate modeling results because of the presence of ladder-like burrs in the simulated seismic wave. Thus, we propose the mapping method based on the dual-variable finite-difference staggered grid. The proposed method uses different size grid spacings in different regions and locally variable time steps to match the size variability of grid spacings. Numerical examples suggest that the proposed method requires less memory storage capacity and improves the computational efficiency compared with forward modeling methods based on the conventional grid.
基金supported by the National Natural Science Foundation of China(No.41174087,41204089)the National Oil and Gas Major Project(No.2011ZX05005-005)
文摘We investigated the effect of microscopic distribution modes of hydrates in porous sediments, and the saturation of hydrates and free gas on the elastic properties of saturated sediments. We simulated the propagation of seismic waves in gas hydrate-bearing sediments beneath the seafloor, and obtained the common receiver gathers of compressional waves（P-waves） and shear waves（S-waves）. The numerical results suggest that the interface between sediments containing gas hydrates and free gas produces a large-amplitude bottomsimulating reflector. The analysis of multicomponent common receiver data suggests that ocean-bottom seismometers receive the converted waves of upgoing P- and S-waves, which increases the complexity of the wavefield record.
基金financially supported by the National Hi-tech Research and Development Program of China(863 Program)(No.2012AA09A20103)
文摘Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell＇s equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic （EM） fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.
基金supported by the National Natural Science Foundation of China(Grant Nos.41404101,41174114,41274130,and 41404102)
文摘AVO forward modeling is based on two-phase medium theory and is considered an effective method for describing reservoir rocks and fluids. However, the method depends on the input matrix mineral bulk modulus and the rationality of the two-phase medium model. We used the matrix mineral bulk modulus inversion method and multiple constraints to obtain a two-phase medium model with physical meaning. The proposed method guarantees the reliability of the obtained AVO characteristicsin two-phase media. By the comparative analysis of different lithology of the core sample, the advantages and accuracy of the inversion method can be illustrated. Also, the inversion method can be applied in LH area, and the AVO characteristics can be obtained when the porosity, fluid saturation, and other important lithology parameters are changed. In particular, the reflection coefficient amplitude difference between the fast P wave and S wave as a function of porosity at the same incidence angle, and the difference in the incidence angle threshold can be used to decipher porosity.
基金The National Natural Science Foundation of China(No.60621002)the National High Technology Research and Development Pro-gram of China(863 Program)(No.2007AA01Z2B4).
文摘A novel behavioral model using three-layer time-delay feed-forward neural networks （TDFFNN）is adopted to model radio frequency （RF）power amplifiers exhibiting memory nonlinearities. In order to extract the parameters, the back- propagation algorithm is applied to train the proposed neural networks. The proposed model is verified by the typical odd- order-only memory polynomial model in simulation, and the performance is compared with different numbers of taped delay lines（TDLs） and perceptrons of the hidden layer. For validating the TDFFNN model by experiments, a digital test bench is set up to collect input and output data of power amplifiers at a 60 × 10^6 sample/s sampling rate. The 3.75 MHz 16-QAM signal generated in the vector signal generator（VSG） is chosen as the input signal, when measuring the dynamic AM/AM and AM/PM characteristics of power amplifiers. By comparisons and analyses, the presented model provides a good performance in convergence, accuracy and efficiency, which is approved by simulation results and experimental results in the time domain and frequency domain.
基金supported by the Program for New Century Excellent Talents in University(No.NCET-04-0370)
文摘We discuss the feasibility of using controlled-source electromagnetic （CSEM） in the frequency domain for prospecting marine gas hydrates. Based on the Ocean Drilling Program （ODP） Leg 164 log data, we have established several 1-D resistivity models which have different gas hydrate concentrations. Meanwhile, we analyzed the electromagnetic response of marine gas hydrates in the frequency domain based on these models. We also studied the relationship between electrical field magnitude or phase and parameters such as receiver-transmitter distance and frequency. Our numerical modeling results provide us with a quantitative reference for exploration and resource evaluation of marine gas hydrates.
基金supported by the National Nature Science Foundation of China(Nos.41504102 and 41604037)National Science and Technology Major Project(No.2016ZX05015-006)Yangtze University Youth Found(No.2015cqn32)
文摘Finite-difference(FD) methods are widely used in seismic forward modeling owing to their computational efficiency but are not readily applicable to irregular topographies. Thus, several FD methods based on the transformation to curvilinear coordinates using body-fitted grids have been proposed, e.g., stand staggered grid(SSG) with interpolation, nonstaggered grid, rotated staggered grid(RSG), and fully staggered. The FD based on the RSG is somewhat superior to others because it satisfies the spatial distribution of the wave equation without additional memory and computational requirements; furthermore, it is simpler to implement. We use the RSG FD method to transform the firstorder stress–velocity equation in the curvilinear coordinates system and introduce the highprecision adaptive, unilateral mimetic finite-difference(UMFD) method to process the freeboundary conditions of an irregular surface. The numerical results suggest that the precision of the solution is higher than that of the vacuum formalism. When the minimum wavelength is low, UMFD avoids the surface wave dispersion. We compare FD methods based on RSG, SEM, and nonstaggered grid and infer that all simulation results are consistent but the computational efficiency of the RSG FD method is higher than the rest.
基金supported by the National Hi-Tech Research and Development Program (863) of China (No.2007AA06A405, 2005AA6010100401)
文摘Water quality models are important tools to support the optimization of aquatic ecosystem rehabilitation programs and assess their efficiency. Basing on the flow conditions of the Daqinghe River Mouth of the Dianchi Lake, China, a two-dimensional water quality model was developed in the research. The hydrodynamics module was numerically solved by the alternating direction iteration （ADI） method. The parameters of the water quality module were obtained through the in situ experiments and the laboratory analyses that were conducted from 2006 to 2007. The model was calibrated and verified by the observation data in 2007. Among the four modelled key variables, i.e., water level, COD （in CODcr）, NH4＋-N and PO43-P the minimum value of the coefficient of determination （COD） was 0.69, indicating the model performed reasonably well. The developed model was then applied to simulate the water quality changes at a downstream cross-section assuming that the designed restoration programs were implemented. According to the simulated results, the restoration programs could cut down the loads of COD and PO43-P about 15%. Such a load reduction, unfortunately, would have very little effect on the NH4^＋-N removal. Moreover, the water quality at the outlet cross-section would be still in class V （3838-02）, indicating more measures should be taken to further reduce the loads. The study demonstrated the capability of water quality models to support aquatic ecosystem restorations.
基金This research is sponsored by the National Natural Science Foundation of China （No. 40374024）.
文摘The workload of the 3D magnetotelluric forward modeling algorithm is so large that the traditional serial algorithm costs an extremely large compute time. However, the 3D forward modeling algorithm can process the data in the frequency domain, which is very suitable for parallel computation. With the advantage of MPI and based on an analysis of the flow of the 3D magnetotelluric serial forward algorithm, we suggest the idea of parallel computation and apply it. Three theoretical models are tested and the execution efficiency is compared in different situations. The results indicate that the parallel 3D forward modeling computation is correct and the efficiency is greatly improved. This method is suitable for large size geophysical computations.
文摘Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-dimensional interior ballistics model of the solid propellant ETC gun(2D-IB-SPETCG)is presented to describe the process of the ETC launch.Both calculated pressure and projectile muzzle velocity accord well with the experimental results.The feasibility of the 2D-IB-SPETCG model is proved.Depending on the experimental data and initial parameters,detailed distribution of the ballistics parameters can be simulated.With the distribution of pressure and temperature of the gas phase and the propellant,the influence of plasma during the ignition process can be analyzed.Because of the radial flowing plasma,the propellant in the area of the DRPG is ignited within 0.01 ms,while all propellant in the chamber is ignited within 0.09 ms.The radial ignition delay time is much less than the axial delay time.During the ignition process,the radial pressure difference is less than 5 MPa at the place 0.025 m away from the breech.The radial ignition uniformity is proved.The temperature of the gas increases from several thousand K(conventional ignition)to several ten thousand K(plasma ignition).Compare the distribution of the density and temperature of the gas,we know that low density and high temperature gas appears near the exits of the DRPG,while high density and low temperature gas appears at the wall near the breech.The simulation of the 2D-IB-SPETCG model is an effective way to investigate the interior ballistics process of the ETC launch.The 2D-IB-SPETC model can be used for prediction and improvement of experiments.
基金This work was funded by National Natural Science Foundation of China （No. 40474044）.
文摘Mathematical geophone （MG） and equal-time stacking （ETS） principles are used to implement seismic prestack forward modeling with irregular surfaces using the oneway acoustic wave-equation. This method receives seismic primary reflections from the subsurface using a set of virtual MGs. The receivers can be located anywhere on an irregular observing surface. Moreover, the ETS method utilizes the one-way acoustic wave equation to easily and quickly image and extrapolate seismic reflection data. The method is illustrated using high single-noise ratio common shot gathers computed by numerical forward modeling of two simple models, one with a flat surface and one with an irregular surface, and a complex normal fault model. A prestack depth migration method for irregular surface topography was used to reoroduce the normal fault model with high accuracy.
