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.展开更多
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.展开更多
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.展开更多
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 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.展开更多
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 progress of modern industry has given rise to great requirements for network transmission latency and reliability in domains such as smart grid and intelligent driving.To address these challenges,the concept of Ti...The progress of modern industry has given rise to great requirements for network transmission latency and reliability in domains such as smart grid and intelligent driving.To address these challenges,the concept of Time-sensitive networking(TSN)is proposed by IEEE 802.1TSN working group.In order to achieve low latency,Cyclic queuing and forwarding(CQF)mechanism is introduced to schedule Timetriggered(TT)flows.In this paper,we construct a TSN model based on CQF and formulate the flow scheduling problem as an optimization problem aimed at maximizing the success rate of flow scheduling.The problem is tackled by a novel algorithm that makes full use of the characteristics and the relationship between the flows.Firstly,by K-means algorithm,the flows are initially partitioned into subsets based on their correlations.Subsequently,the flows within each subset are sorted by a new special criteria extracted from multiple features of flow.Finally,a flow offset selecting method based on load balance is used for resource mapping,so as to complete the process of flow scheduling.Experimental results demonstrate that the proposed algorithm exhibits significant advantages in terms of scheduling success rate and time efficiency.展开更多
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.展开更多
Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a...Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.展开更多
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.展开更多
Currently, the selection of receiving traces in geometry design is mostly based on the horizontal layered medium hypothesis, which is unable to meet survey requirements in a complex area. This paper estimates the opti...Currently, the selection of receiving traces in geometry design is mostly based on the horizontal layered medium hypothesis, which is unable to meet survey requirements in a complex area. This paper estimates the optimal number of receiving traces in field geometry using a numerical simulation based on a field test conducted in previous research （Zhu et al., 2011）. A mathematical model is established for total energy and average efficiency energy using fixed trace spacing and optimal receiving traces are estimated. Seismic data acquired in a complex work area are used to verify the correctness of the proposed method. Results of model data calculations and actual data processing show that results are in agreement. This indicates that the proposed method is reasonable, correct, sufficiently scientific, and can be regarded as a novel method for use in seismic geometry design in complex geological regions.展开更多
Existing studies indicate that gas hydrate-bearing formations exhibit notable seismic velocity dispersion and attenuation. The Shenhu area of the South China Sea hold significant gas hydrate resource potential;however...Existing studies indicate that gas hydrate-bearing formations exhibit notable seismic velocity dispersion and attenuation. The Shenhu area of the South China Sea hold significant gas hydrate resource potential;however, the relationship between seismic velocity dispersion, attenuation properties, and gas-hydrate saturation remains insufficiently understood. Furthermore, a significant mismatch exists between the real seismic angle gather near a well and the synthetic angle gather generated using the convolution method, and this discrepancy may arise from the seismic velocity dispersion and attenuation characteristics of the gas hydrate-bearing formations. In this paper, we develop a rock physics model that integrates White's and Dvorkin's models, accounting for varied types of gas-hydrate occurrence states,specifically tailored to the gas hydrate-bearing formations in the Shenhu area. This model is calibrated with well log data and employed to investigate how gas-hydrate saturation influences seismic velocity dispersion and attenuation. Numerical analysis reveals the coexistence of two types of gas-hydrate occurrence states in the region: high gas-hydrate saturation formations are dominated by loadbearing-type gas hydrate, and formations containing both gas hydrate and free gas may exhibit either load-bearing or pore-filling types. The seismic velocity dispersion and attenuation properties vary significantly depending on the gas-hydrate occurrence state. We further apply the proposed model to generate seismic velocity and attenuation logs at various frequencies. These logs are used in seismic forward modeling employing both the convolution method and the propagator matrix method. Well tie analysis indicates that the synthetic angle gather incorporating attenuation via the propagator matrix method aligns more closely with the real seismic angle gather than the convolution method. This study provides valuable insights into frequency-dependent amplitude versus offset(AVO) analysis and the seismic interpretation of gas hydrate-bearing formations in the South China Sea.展开更多
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.展开更多
In elastic wave forward modeling, absorbing boundary conditions （ABC） are used to mitigate undesired reflections from the model truncation boundaries. The perfectly matched layer （PML） has proved to be the best av...In elastic wave forward modeling, absorbing boundary conditions （ABC） are used to mitigate undesired reflections from the model truncation boundaries. The perfectly matched layer （PML） has proved to be the best available ABC. However, the traditional splitting PML （SPML） ABC has some serious disadvantages： for example, global SPML ABCs require much more computing memory, although the implementation is easy. The implementation of local SPML ABCs also has some difficulties, since edges and corners must be considered. The traditional non-splitting perfectly matched layer （NPML） ABC has complex computation because of the convolution. In this paper, based on non-splitting perfectly matched layer （NPML） ABCs combined with the complex frequency-shifted stretching function （CFS）, we introduce a novel numerical implementation method for PML absorbing boundary conditions with simple calculation equations, small memory requirement, and easy programming.展开更多
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.展开更多
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.展开更多
Borehole-to-surface electrical imaging （BSEI） uses a line source and a point source to generate a stable electric field in the ground. In order to study the surface potential of anomalies, three-dimensional forward ...Borehole-to-surface electrical imaging （BSEI） uses a line source and a point source to generate a stable electric field in the ground. In order to study the surface potential of anomalies, three-dimensional forward modeling of point and line sources was conducted by using the finite-difference method and the incomplete Cholesky conjugate gradient （ICCG） method. Then, the damping least square method was used in the 3D inversion of the formation resistivity data. Several geological models were considered in the forward modeling and inversion. The forward modeling results suggest that the potentials generated by the two sources have different surface signatures. The inversion data suggest that the low- resistivity anomaly is outlined better than the high-resistivity anomaly. Moreover, when the point source is under the anomaly, the resistivity anomaly boundaries are better outlined than when using a line source.展开更多
The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the M...The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the Maxwell＇s equations, and the magnetic components at the center of each edge of the grid cell are evaluated by applying the staggered-grid finite-difference method. The tipper and its divergence are derived to complete the 3D ZTEM forward modeling algorithm. A synthetic model is then used to compare the responses with those of 2D finite-element forward modeling to verify the accuracy of the algorithm. ZTEM offers high horizontal resolution to both simple and complex distributions of conductivity. This work is the theoretical foundation for the interpretation of ZTEM data and the study of 3D ZTEM inversion.展开更多
In this paper, a mathematical model consisting of forward and backward models is built on parallel genetic algorithms (PGAs) for fault diagnosis in a transmission power system. A new method to reduce the scale of faul...In this paper, a mathematical model consisting of forward and backward models is built on parallel genetic algorithms (PGAs) for fault diagnosis in a transmission power system. A new method to reduce the scale of fault sections is developed in the forward model and the message passing interface (MPI) approach is chosen to parallel the genetic algorithms by global sin-gle-population master-slave method (GPGAs). The proposed approach is applied to a sample system consisting of 28 sections, 84 protective relays and 40 circuit breakers. Simulation results show that the new model based on GPGAs can achieve very fast computation in online applications of large-scale power systems.展开更多
基金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.
基金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.
基金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 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.
文摘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.
基金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.
基金supported by Science and Technology Project of State Grid Corporation Headquarters under Grant 5108-202218280A-2-170-XG(Development and Application of Power Time-Sensitive Network Switching Chip。
文摘The progress of modern industry has given rise to great requirements for network transmission latency and reliability in domains such as smart grid and intelligent driving.To address these challenges,the concept of Time-sensitive networking(TSN)is proposed by IEEE 802.1TSN working group.In order to achieve low latency,Cyclic queuing and forwarding(CQF)mechanism is introduced to schedule Timetriggered(TT)flows.In this paper,we construct a TSN model based on CQF and formulate the flow scheduling problem as an optimization problem aimed at maximizing the success rate of flow scheduling.The problem is tackled by a novel algorithm that makes full use of the characteristics and the relationship between the flows.Firstly,by K-means algorithm,the flows are initially partitioned into subsets based on their correlations.Subsequently,the flows within each subset are sorted by a new special criteria extracted from multiple features of flow.Finally,a flow offset selecting method based on load balance is used for resource mapping,so as to complete the process of flow scheduling.Experimental results demonstrate that the proposed algorithm exhibits significant advantages in terms of scheduling success rate and time efficiency.
基金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.41474110)Shell Ph.D. Scholarship to support excellence in geophysical research
文摘Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.
基金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(No.41304115)National Key S&T Special Projects(No.2016ZX050 24001-003)+2 种基金Open Fund for Sichuan Province Key Laboratory of Natural Gas Geology(No.2015trqdz02)the Research Project,CNPC(No.2016A-33)"Young and Middle-aged Key Teachers"Training Program in Southwest Petroleum University
文摘Currently, the selection of receiving traces in geometry design is mostly based on the horizontal layered medium hypothesis, which is unable to meet survey requirements in a complex area. This paper estimates the optimal number of receiving traces in field geometry using a numerical simulation based on a field test conducted in previous research （Zhu et al., 2011）. A mathematical model is established for total energy and average efficiency energy using fixed trace spacing and optimal receiving traces are estimated. Seismic data acquired in a complex work area are used to verify the correctness of the proposed method. Results of model data calculations and actual data processing show that results are in agreement. This indicates that the proposed method is reasonable, correct, sufficiently scientific, and can be regarded as a novel method for use in seismic geometry design in complex geological regions.
基金supported by National Natural Science Foundation of China(W2431028,42122029)SINOPEC Fundamental Research Program(P24258)CNPC Investigations on fundamental experiments and advanced theoretical methods in geophysical prospecting applications(2022DQ0604-02).
文摘Existing studies indicate that gas hydrate-bearing formations exhibit notable seismic velocity dispersion and attenuation. The Shenhu area of the South China Sea hold significant gas hydrate resource potential;however, the relationship between seismic velocity dispersion, attenuation properties, and gas-hydrate saturation remains insufficiently understood. Furthermore, a significant mismatch exists between the real seismic angle gather near a well and the synthetic angle gather generated using the convolution method, and this discrepancy may arise from the seismic velocity dispersion and attenuation characteristics of the gas hydrate-bearing formations. In this paper, we develop a rock physics model that integrates White's and Dvorkin's models, accounting for varied types of gas-hydrate occurrence states,specifically tailored to the gas hydrate-bearing formations in the Shenhu area. This model is calibrated with well log data and employed to investigate how gas-hydrate saturation influences seismic velocity dispersion and attenuation. Numerical analysis reveals the coexistence of two types of gas-hydrate occurrence states in the region: high gas-hydrate saturation formations are dominated by loadbearing-type gas hydrate, and formations containing both gas hydrate and free gas may exhibit either load-bearing or pore-filling types. The seismic velocity dispersion and attenuation properties vary significantly depending on the gas-hydrate occurrence state. We further apply the proposed model to generate seismic velocity and attenuation logs at various frequencies. These logs are used in seismic forward modeling employing both the convolution method and the propagator matrix method. Well tie analysis indicates that the synthetic angle gather incorporating attenuation via the propagator matrix method aligns more closely with the real seismic angle gather than the convolution method. This study provides valuable insights into frequency-dependent amplitude versus offset(AVO) analysis and the seismic interpretation of gas hydrate-bearing formations in the South China Sea.
基金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.
基金sponsored by the Chinese National Development and Reform Commission(No.[2005]2372)the Innovative Technological Research Foundation of PetroChina Company Limited(No.060511-1-3)
文摘In elastic wave forward modeling, absorbing boundary conditions （ABC） are used to mitigate undesired reflections from the model truncation boundaries. The perfectly matched layer （PML） has proved to be the best available ABC. However, the traditional splitting PML （SPML） ABC has some serious disadvantages： for example, global SPML ABCs require much more computing memory, although the implementation is easy. The implementation of local SPML ABCs also has some difficulties, since edges and corners must be considered. The traditional non-splitting perfectly matched layer （NPML） ABC has complex computation because of the convolution. In this paper, based on non-splitting perfectly matched layer （NPML） ABCs combined with the complex frequency-shifted stretching function （CFS）, we introduce a novel numerical implementation method for PML absorbing boundary conditions with simple calculation equations, small memory requirement, and easy programming.
基金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.
基金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.
基金sponsored by the National Major Project(No.2016ZX05014-001)the National Natural Science Foundation of China(No.41172130 and U1403191)the Fundamental Research Funds for the Central Universities(No.2-9-2015-209)
文摘Borehole-to-surface electrical imaging （BSEI） uses a line source and a point source to generate a stable electric field in the ground. In order to study the surface potential of anomalies, three-dimensional forward modeling of point and line sources was conducted by using the finite-difference method and the incomplete Cholesky conjugate gradient （ICCG） method. Then, the damping least square method was used in the 3D inversion of the formation resistivity data. Several geological models were considered in the forward modeling and inversion. The forward modeling results suggest that the potentials generated by the two sources have different surface signatures. The inversion data suggest that the low- resistivity anomaly is outlined better than the high-resistivity anomaly. Moreover, when the point source is under the anomaly, the resistivity anomaly boundaries are better outlined than when using a line source.
基金supported by the Natural Science Foundation of China(No.41374078)Geological Survey Projects of Ministry of Land and Resources of China(No.12120113086100 and 12120113101300)
文摘The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the Maxwell＇s equations, and the magnetic components at the center of each edge of the grid cell are evaluated by applying the staggered-grid finite-difference method. The tipper and its divergence are derived to complete the 3D ZTEM forward modeling algorithm. A synthetic model is then used to compare the responses with those of 2D finite-element forward modeling to verify the accuracy of the algorithm. ZTEM offers high horizontal resolution to both simple and complex distributions of conductivity. This work is the theoretical foundation for the interpretation of ZTEM data and the study of 3D ZTEM inversion.
基金the National Natural Science Foundation of China (No. 50677062)the New Century Excellent Talents in Uni-versity of China (No. NCET-07-0745)the Natural Science Foundation of Zhejiang Province, China (No. R107062)
文摘In this paper, a mathematical model consisting of forward and backward models is built on parallel genetic algorithms (PGAs) for fault diagnosis in a transmission power system. A new method to reduce the scale of fault sections is developed in the forward model and the message passing interface (MPI) approach is chosen to parallel the genetic algorithms by global sin-gle-population master-slave method (GPGAs). The proposed approach is applied to a sample system consisting of 28 sections, 84 protective relays and 40 circuit breakers. Simulation results show that the new model based on GPGAs can achieve very fast computation in online applications of large-scale power systems.
