We propose a new 3D inversion scheme to invert the near-and transition-zone data of CSAMT with topography accurately.In this new method,the earth was discretized into unstructured tetrahedra to fit the ragged topograp...We propose a new 3D inversion scheme to invert the near-and transition-zone data of CSAMT with topography accurately.In this new method,the earth was discretized into unstructured tetrahedra to fit the ragged topography and the vector fi nite-element method was adopted to obtain precise responses and good sensitivity.To simulate the attitude and shape of the transmitter,we divided a long-grounded transmitter into dipoles and integrated these dipoles to obtain good responses in the near-and transition-fi eld zones.Next,we designed an L2 norm-based objective functional and applied a standard quasi-Newton method as the optimization method to solve the inverse problem and guarantee steady convergence.We tested our 3D inversion method first on synthetic data and then on a field dataset acquired from select sites near Changbai Mountain,China.In both tests,the new inversion algorithm achieved excellent fitting between the predicted and observed data,even in near-and transition-fi eld zones,and the inversion results agreed well with the true model.These fi ndings reveal that the proposed algorithm is eff ective for 3D inversion of CSAMT data.展开更多
The comprehensive utilization of wood is the main goal of log cutting,but knot defects increase the diffi-culty of rationally optimizing cutting.Due to the lack of real shape data of knot defects in logs,it is diffi c...The comprehensive utilization of wood is the main goal of log cutting,but knot defects increase the diffi-culty of rationally optimizing cutting.Due to the lack of real shape data of knot defects in logs,it is diffi cult for detection methods to establish a correlation between signal and defect morphology.An image-processing method is proposed for knot inversion based on distance regularized level set segmentation(DRLSE)and spatial vertex clustering,and with the inversion of the defects existing relative board position in the log,an inversion model of the knot defect is established.First,the defect edges of the top and bottom images of the boards are extracted by DRLSE and ellipse fi tting,and the major axes of the ellipses made coplanar by angle correction;second,the coordinate points of the top and bottom ellipse edges are extracted to form a spatial straight line;third,to solve the intersection dispersion of spatial straight lines and the major axis plane,K-medoids clustering is used to locate the vertex.Finally,with the vertex and the large ellipse,a 3D cone model is constructed which can be used to invert the shape of knots in the board.The experiment was conducted on ten defective larch boards,and the experimental results showed that this method can accurately invert the shapes of defects in solid wood boards with the advantages of low cost and easy operation.展开更多
Ultrasound computed tomography(USCT)is a noninvasive biomedical imaging modality that offers insights into acoustic properties such as the sound speed(SS)and acoustic attenuation(AA)of the human body,enhancing diagnos...Ultrasound computed tomography(USCT)is a noninvasive biomedical imaging modality that offers insights into acoustic properties such as the sound speed(SS)and acoustic attenuation(AA)of the human body,enhancing diagnostic accuracy and therapy planning.Full waveform inversion(FWI)is a promising USCT image reconstruction method that optimizes the parameter fields of a wave propagation model via gradient-based optimization.However,twodimensional FWI methods are limited by their inability to account for three-dimensional wave propagation in the elevation direction,resulting in image artifacts.To address this problem,we propose a three-dimensional time-domain full waveform inversion algorithm to reconstruct the SS and AA distributions on the basis of a fractional Laplacian wave equation,adjoint field formulation,and gradient descent optimization.Validated by two sets of simulations,the proposed algorithm has potential for generating high-resolution and quantitative SS and AA distributions.This approach holds promise for clinical USCT applications,assisting early disease detection,precise abnormality localization,and optimized treatment planning,thus contributing to better healthcare outcomes.展开更多
Funded by The National Key Research and Development Program of China,China Deep Exploration(Sinoprobe)and The China Geological Suvery Project on 2009–2019,a large scale magnetotelluric sounding(MT)survey grid(Fig.1)h...Funded by The National Key Research and Development Program of China,China Deep Exploration(Sinoprobe)and The China Geological Suvery Project on 2009–2019,a large scale magnetotelluric sounding(MT)survey grid(Fig.1)has covered whole south China.展开更多
Based on the analysis of the conjugate gradient algorithm, we implement a threedimensional （3D） conjugate gradient inversion algorithm with magnetotelluric impedance data. During the inversion process, the 3D conjug...Based on the analysis of the conjugate gradient algorithm, we implement a threedimensional （3D） conjugate gradient inversion algorithm with magnetotelluric impedance data. During the inversion process, the 3D conjugate gradient inversion algorithm doesn＇ t need to compute and store the Jacobian matrix but directly updates the model from the computation of the Jacobian matrix. Requiring only one forward and four pseudo-forward modeling applications per frequency to produce the model update at each iteration, this algorithm efficiently reduces the computation of the inversion. From a trial inversion with synthetic magnetotelluric data, the validity and stability of the 3D conjugate gradient inversion algorithm is verified.展开更多
How to get the rapid and stable inversion results and reconstruct the clear subsurface resistivity structures is a focus problem in current magnetotelluric inversion. A stable solution of an ill-posed inverse problem ...How to get the rapid and stable inversion results and reconstruct the clear subsurface resistivity structures is a focus problem in current magnetotelluric inversion. A stable solution of an ill-posed inverse problem was obtained by the regularization methods in which some desired structures were imposed to stabilize the inverse problem. By the smoothness-constrained model and approximate sensitivity method, the stable subsurface resistivity structures were reconstructed. The synthetic examples show that the smoothness-constrained regularized inversion method is effective and can be reasonable to reconstruct three-dimensional subsurface resistivity structures.展开更多
Accurate assessment of snowpack volumetric liquid water content and bulk density is essential for understanding snow hydrology,avalanche risk management,and monitoring cryosphere changes.This study presents a novel du...Accurate assessment of snowpack volumetric liquid water content and bulk density is essential for understanding snow hydrology,avalanche risk management,and monitoring cryosphere changes.This study presents a novel dual-parameter inversion framework that integrates synthetic electromagnetic modelling,dimensionality reduction,and machine learning algorithms to extract relative permittivity and log-resistivity from ground-penetrating radar(GPR)data.Traditional snowpack measurements are invasive,labor-intensive,and limited to point observations.To overcome these limitations,we developed a non-invasive,scalable,and data-driven framework that uses synthetic GPR datasets representing diverse snowpack conditions with variable moisture and density profiles.Synthetic 1D time series reflections(A-scans)are generated using finite-difference time-domain simulations in the state-of-the-art electromagnetic simulator gprMax.Principal component analysis(PCA)is applied to compress each A-scan while preserving key features,which significantly improved and enhanced the model training efficiency.Four machine learning models,including random forest,neural network,support vector machine,and eXtreme gradient boosting,are trained on PCA-reduced features.Among these,the neural network model achieved the best performance,with R^(2)>0.97 for permittivity and R 2>0.92 for resistivity.Gaussian noise(signal-to-noise ratio of 6 dB)is introduced to the synthetic data,and then targeted domain adaptation is employed to enhance generalization to field data.The framework is validated on two contrasting GPR transects in the Altay Mountains of the Chinese mainland,representing moist(T750)and wet(G125)snowpack conditions.The neural network model predictions are most consistent with the GPR derived estimates,Snowfork measurements,and snow pit data,achieving volumetric liquid water content deviation of≤1.5% and bulk density error within the range of 30-84 kg m^(-3).The results demonstrate that machine learning-based inversion,supported by realistic simulations and data augmentation enables scalable,non-invasive snowpack characterization with significant applications in hydrological forecasting,snow monitoring,and water resource management.展开更多
Co-assembling chiral molecules with achiral compounds via non-covalent interactions like areneperfluoroarene(AP) interactions offers an effective approach for fabricating chiral functional materials.Herein,chiral mole...Co-assembling chiral molecules with achiral compounds via non-covalent interactions like areneperfluoroarene(AP) interactions offers an effective approach for fabricating chiral functional materials.Herein,chiral molecules L/D-PF1 and L/D-PF2 with pyrene groups were synthesized and its chiroptical properties upon co-assembly with achiral compound octafluoronaphthalene(OFN) through AP interaction were systemically studied.The co-assembly of L/D-PF1/OFN and L/D-PF2/OFN exhibited distinct chiroptical properties such as circular dichroism(CD) and circularly polarized luminescence(CPL) signals.Chirality transfer from the chirality center of L/D-PF1 and L/D-PF2 to the achiral OFN and chiral amplification were successfully achieved.Besides,no significant CPL signal was observed in the self-assembly of L/DPF1 or L/D-PF2 while co-assembly with OFN exhibited obvious CPL amplification induced by AP interaction.Notably,a reversal CD signal and CPL signal could be observed in L/D-PF2/OFN when the molar ratio changed from 1:1 to 1:2 while not found in L/D-PF1/OFN,indicating that that minor structural changes of molecules could cause large changes in assembly.In addition,a series of computational calculations were conducted to verify the AP interaction between L-PF1/L-PF2 and OFN.This work demonstrated that arene-perfluoroarene interaction could drive chiral transfer,chiral amplification and chiral inversion and provided a new method for the preparation of chiroptical materials.展开更多
According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing ...According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing differential equations are derived by the Hamilton's variational principle,and then solved by the iterative Separation-of-Variable(i SOV)method,which are applicable to arbitrary combinations of homogeneous Boundary Conditions(BCs).However,only the simply-support,clamped and cantilever panels are considered in this work for the sake of clarity.With the closed-form eigensolutions,the flutter frequency,flutter mode and flutter boundary are presented,and the effect of shear deformation and aerodynamic damping on flutter frequencies is investigated.Besides,the relation between panel energy and the work of aerodynamic load is discussed.The numerical comparisons reveal the following.(A)The flutter eigenvalues obtained by the present method are accurate,validated by the Finite Element Method(FEM)and the Galerkin method.(B)When the span-chord ratio is larger than 3,simplifying a 3D panel to 2D(two-dimensional)panel is reasonable and the relative differences of the flutter points predicted by the two models are less than one percent.(C)The reciprocal relationship between the mechanical energy of the panel and the work done by aerodynamic load is verified by using the present flutter eigenvalues and modes,further indicating the high accuracy of the present solutions.(D)The coupling of shear deformation and aerodynamic damping prevents frequency coalescing.展开更多
The gravity inversion is to restore genetic density distribution of the underground target to be explored for explaining the internal structure and distribution of the Earth.In this paper,we propose a new 3D gravity i...The gravity inversion is to restore genetic density distribution of the underground target to be explored for explaining the internal structure and distribution of the Earth.In this paper,we propose a new 3D gravity inversion method based on 3D U-Net++.Compared with two-dimensional gravity inversion,three-dimensional(3D)gravity inversion can more precisely describe the density distribution of underground space.However,conventional 3D gravity inversion method input is two-dimensional,the input and output of the network proposed in our method are three-dimensional.In the training stage,we design a large number of diversifi ed simulation model-data pairs by using the random walk method to improve the generalization ability of the network.In the test phase,we verify the network performance by using the model-data pairs generated by the simulation.To further illustrate the eff ectiveness of the algorithm,we apply the method to the inversion of the San Nicolas mining area,and the inversion results are basically consistent with the borehole measurement results.Moreover,the results of the 3D U-Net++inversion and the 3D U-Net inversion are compared.The density models of the 3D U-Net++inversion have higher resolution,more concentrated inversion results,and a clearer boundary of the density model.展开更多
This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zer...This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zero norm solution. The inversion approach mainly employs forward modeling; a depth weight function is introduced into the objective function of the zero norms. Sparse inversion results are obtained by the corresponding optimal mathematical method. To achieve the practical geophysical and geological significance of the results, penalty function is applied to constrain the density values. Results obtained by proposed provide clear boundary depth and density contrast distribution information. The method's accuracy, validity, and reliability are verified by comparing its results with those of synthetic models. To further explain its reliability, a practical gravity data is obtained for a region in Texas, USA is applied. Inversion results for this region are compared with those of previous studies, including a research of logging data in the same area. The depth of salt dome obtained by the inversion method is 4.2 km, which is in good agreement with the 4.4 km value from the logging data. From this, the practicality of the inversion method is also validated.展开更多
Based on three-dimensional joint finite element, this paper discusses the theory and methodology about inversionof geodetic data. The FEM and inversion formula is given in detail; also a related code is developed. By ...Based on three-dimensional joint finite element, this paper discusses the theory and methodology about inversionof geodetic data. The FEM and inversion formula is given in detail; also a related code is developed. By use of theGreen's function about 3-D FEM, we invert geodetic measurementS of coseismic deformation of the 1989 Ms=7. 1Loma Prieta earthquake to datermine itS source mechanism. The result indicates that the slip on the fault plane isvery heterogeneous. The maximum slip and shear stress are located about 10 kin to northwest of the eathquakesource, the stress drop is about more than 1 MPa.展开更多
The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the tim...The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.展开更多
We developed a three-dimensional（3D） conjugate gradient inversion algorithm for in-verting magnetotelluric impedance tensor measurements.In order to show the importance of including diagonal components of magnetotel...We developed a three-dimensional（3D） conjugate gradient inversion algorithm for in-verting magnetotelluric impedance tensor measurements.In order to show the importance of including diagonal components of magnetotelluric impedance tensor in 3D inversion,synthetic data were inverted using the 3D conjugate gradient inversion,and the inversion results were compared and analyzed.The results from the 3D inversion of synthetic data indicate that both the off-diagonal and the diagonal components are required in inversions to obtain better inversion results when there are no enough data sites to recover the target resistivity structure.These examples show that lots of information about 3D structure is also contained in the diagonal components;as a result,diagonal components should be in-cluded in 3D inversions.The inversion algorithm was also used to invert the impedance tensor data ac-quired in the Kayabe area in Japan.Inversions with the synthetic and real data demonstrated the va-lidity and practicability of the inversion algorithm.展开更多
A three-dimensional crustal density model beneath North China is determined using P-wave traveltimes and gravity datum by sequential inversion method in this paper. To improve the uniqueness of the solution, we used a...A three-dimensional crustal density model beneath North China is determined using P-wave traveltimes and gravity datum by sequential inversion method in this paper. To improve the uniqueness of the solution, we used a linear relationship between velocity and density to achieve mutual conversions and constraints between velocity difference and density difference. Algebraic reconstruction technique （ART） was used in density inversion, which highly improved the calculation speed comparing with common least squares method. The inversion results indicate that the crustal density beneath North China is extremely inhomogeneous with its distribution coherent with regional tectonics. The lower crust of Taihang mountain tectonic zone shows an obvious low density characteristic. We proposed that it may be an indicator of upwelling of mantle materials or partial melting of lowermost crust.展开更多
Based on the analysis of impedance tensor data, tipper data, and the conjugate gradient algorithm, we develop a three-dimensional (3D) conjugate gradient algorithm for inverting magnetotelluric full information data d...Based on the analysis of impedance tensor data, tipper data, and the conjugate gradient algorithm, we develop a three-dimensional (3D) conjugate gradient algorithm for inverting magnetotelluric full information data determined from five electric and magnetic field components and discuss the method to use the full information data for quantitative interpretation of 3D inversion results. Results from the 3D inversion of synthetic data indicate that the results from inverting full information data which combine the impedance tensor and tipper data are better than results from inverting only the impedance tensor data (or tipper data) in improving resolution and reliability. The synthetic examples also demonstrate the validity and stability of this 3D inversion algorithm.展开更多
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.展开更多
This study explores the quasi-real time inversion principle and precision estimation of three-dimensional coordinates of the epicenter, trigger time and magnitude of earthquakes with the aim to improve traditional met...This study explores the quasi-real time inversion principle and precision estimation of three-dimensional coordinates of the epicenter, trigger time and magnitude of earthquakes with the aim to improve traditional methods, which are flawed due to missing information or distortion in the seismograph records. The epicenter, trigger time and magnitude from the Lushan earthquake are inverted and analyzed based on high-frequency GNSS data. The inversion results achieved a high precision, which are consistent with the data published by the China Earthquake Administration. Moreover, it has been proven that the inversion method has good theoretical value and excellent application prospects.展开更多
Abstract. Due to the incapability of high frequency surface wave radar （HFSWR） to detect undercurrent parameters, a new algorithm is proposed to apply a three-dimensional （3D） nonlinear barotropic shallow sea and ...Abstract. Due to the incapability of high frequency surface wave radar （HFSWR） to detect undercurrent parameters, a new algorithm is proposed to apply a three-dimensional （3D） nonlinear barotropic shallow sea and continental shelf model in coordinate system to the inversion of undercurrent. The calculation domain of this model is the area detected by HFSWR. Considering the benthal topography of the detected area and the ocean dynamic parameters, such as surface current, wind and wave detected by HFSWR, the relation between surface current and undercurrent is established in this model, Accordingly, the undercurrent parameters of corresponding area are obtained. The inversion results agree with the law of ocean dynamics and reproduce the movement of undercurrent.展开更多
Herein,a three-dimensional(3D)inversion method in the frequency domain based on a time–frequency transformation was developed to improve the efficiency of the 3D inversion of transient electromagnetic(TEM)data.The Fo...Herein,a three-dimensional(3D)inversion method in the frequency domain based on a time–frequency transformation was developed to improve the efficiency of the 3D inversion of transient electromagnetic(TEM)data.The Fourier transform related to the electromagnetic response in the frequency and time domains becomes a sine or cosine transform under the excitation of downward-step current.We established a transformation matrix based on the digital fi ltering calculation for the sine transform,and then the frequency domain projection of the TEM data was determined from the linear transformation system using the smoothing constrained least squares inversion method,in which only the imaginary part was used to maintain the TEM data transformation equivalence in the bidirectional projection.Thus,the time-domain TEM inversion problem was indirectly and effectively solved in the frequency domain.In the 3D inversion of the transformed frequency-domain data,the limited-memory Broyden–Fletcher–Goldfarb–Shannoquasi–Newton(L-BFGS)method was used and modifi ed with a restart strategy to adjust the regularization parameter when the algorithm tended to a local minimum.Synthetic data tests showed that our domain transformation method can stably project the TEM data into the frequency domain with very high accuracy;furthe rmore,the 3D inversion of the transformed frequency-domain data is stable,can be used to recover the real resistivity model with an acceptable effi ciency.展开更多
基金the Strategic Priority Research Program of the Chinese Academy of Sciences(No.XDA14020102)the National Natural Science Foundation of China(Nos.41774125,41530320,41904104)+1 种基金the Key National Research Project of China(No.2018YFC0603300)the S&T Program of Beijing(No.Z181100005718001).
文摘We propose a new 3D inversion scheme to invert the near-and transition-zone data of CSAMT with topography accurately.In this new method,the earth was discretized into unstructured tetrahedra to fit the ragged topography and the vector fi nite-element method was adopted to obtain precise responses and good sensitivity.To simulate the attitude and shape of the transmitter,we divided a long-grounded transmitter into dipoles and integrated these dipoles to obtain good responses in the near-and transition-fi eld zones.Next,we designed an L2 norm-based objective functional and applied a standard quasi-Newton method as the optimization method to solve the inverse problem and guarantee steady convergence.We tested our 3D inversion method first on synthetic data and then on a field dataset acquired from select sites near Changbai Mountain,China.In both tests,the new inversion algorithm achieved excellent fitting between the predicted and observed data,even in near-and transition-fi eld zones,and the inversion results agreed well with the true model.These fi ndings reveal that the proposed algorithm is eff ective for 3D inversion of CSAMT data.
基金supported fi nancially by the China State Forestry Administration“948”projects(2015-4-52),and Hei-longjiang Natural Science Foundation(C2017005).
文摘The comprehensive utilization of wood is the main goal of log cutting,but knot defects increase the diffi-culty of rationally optimizing cutting.Due to the lack of real shape data of knot defects in logs,it is diffi cult for detection methods to establish a correlation between signal and defect morphology.An image-processing method is proposed for knot inversion based on distance regularized level set segmentation(DRLSE)and spatial vertex clustering,and with the inversion of the defects existing relative board position in the log,an inversion model of the knot defect is established.First,the defect edges of the top and bottom images of the boards are extracted by DRLSE and ellipse fi tting,and the major axes of the ellipses made coplanar by angle correction;second,the coordinate points of the top and bottom ellipse edges are extracted to form a spatial straight line;third,to solve the intersection dispersion of spatial straight lines and the major axis plane,K-medoids clustering is used to locate the vertex.Finally,with the vertex and the large ellipse,a 3D cone model is constructed which can be used to invert the shape of knots in the board.The experiment was conducted on ten defective larch boards,and the experimental results showed that this method can accurately invert the shapes of defects in solid wood boards with the advantages of low cost and easy operation.
基金supported by the National Key Research and Development Program of China(2022YFA1404400)the National Natural Science Foundation of China(62122072,12174368,61705216,62405306)+4 种基金Anhui Provincial Department of Science and Technology(202203a07020020,18030801138)Anhui Provincial Natural Science Foundation(2308085QA21,2408085QF187)the USTC Research Funds of the Double First-Class Initiative(YD2090002015)the Institute of Artificial Intelligence at Hefei Comprehensive National Science Center(23YGXT005)the Fundamental Research Funds for the Central Universities(WK2090000083).
文摘Ultrasound computed tomography(USCT)is a noninvasive biomedical imaging modality that offers insights into acoustic properties such as the sound speed(SS)and acoustic attenuation(AA)of the human body,enhancing diagnostic accuracy and therapy planning.Full waveform inversion(FWI)is a promising USCT image reconstruction method that optimizes the parameter fields of a wave propagation model via gradient-based optimization.However,twodimensional FWI methods are limited by their inability to account for three-dimensional wave propagation in the elevation direction,resulting in image artifacts.To address this problem,we propose a three-dimensional time-domain full waveform inversion algorithm to reconstruct the SS and AA distributions on the basis of a fractional Laplacian wave equation,adjoint field formulation,and gradient descent optimization.Validated by two sets of simulations,the proposed algorithm has potential for generating high-resolution and quantitative SS and AA distributions.This approach holds promise for clinical USCT applications,assisting early disease detection,precise abnormality localization,and optimized treatment planning,thus contributing to better healthcare outcomes.
基金co-supported by the China Geological Survey Project(DD20190012 and DD20160082)
文摘Funded by The National Key Research and Development Program of China,China Deep Exploration(Sinoprobe)and The China Geological Suvery Project on 2009–2019,a large scale magnetotelluric sounding(MT)survey grid(Fig.1)has covered whole south China.
基金sponsored by National Natural Science Foundation of China (Grant Nos. 40774029, 40674037, and 40374024)the National Hi-tech Research and Development Program of China (863 Program) (No. 2007AA09Z310)the Program for New Century Excellent Talents in University (NCET).
文摘Based on the analysis of the conjugate gradient algorithm, we implement a threedimensional （3D） conjugate gradient inversion algorithm with magnetotelluric impedance data. During the inversion process, the 3D conjugate gradient inversion algorithm doesn＇ t need to compute and store the Jacobian matrix but directly updates the model from the computation of the Jacobian matrix. Requiring only one forward and four pseudo-forward modeling applications per frequency to produce the model update at each iteration, this algorithm efficiently reduces the computation of the inversion. From a trial inversion with synthetic magnetotelluric data, the validity and stability of the 3D conjugate gradient inversion algorithm is verified.
基金Project(20110162120064)supported by Higher School Doctor Subject Special Scientific Research Foundation of ChinaProject(10JJ6059)supported by the Natural Science Foundation of Hunan Province,China
文摘How to get the rapid and stable inversion results and reconstruct the clear subsurface resistivity structures is a focus problem in current magnetotelluric inversion. A stable solution of an ill-posed inverse problem was obtained by the regularization methods in which some desired structures were imposed to stabilize the inverse problem. By the smoothness-constrained model and approximate sensitivity method, the stable subsurface resistivity structures were reconstructed. The synthetic examples show that the smoothness-constrained regularized inversion method is effective and can be reasonable to reconstruct three-dimensional subsurface resistivity structures.
基金supported by the National Key R&D Program of China(Grant Nos.2023YFC3008300&2023YFC3008305)the National Natural Science Foundation of China(Grant No.42172320)+1 种基金the Key Laboratory of Mountain Hazards and Engineering Resilience,Institute of Mountain Hazards and Environment,Chinese Academy of Sciences(Grant Nos.KLMHER-Z06&KLMHER-T07)the Science and Technology Research Program of Institute of Mountain Hazards and Environment,Chinese Academy of Sciences(Grant No.IMHE-CXTD.04).
文摘Accurate assessment of snowpack volumetric liquid water content and bulk density is essential for understanding snow hydrology,avalanche risk management,and monitoring cryosphere changes.This study presents a novel dual-parameter inversion framework that integrates synthetic electromagnetic modelling,dimensionality reduction,and machine learning algorithms to extract relative permittivity and log-resistivity from ground-penetrating radar(GPR)data.Traditional snowpack measurements are invasive,labor-intensive,and limited to point observations.To overcome these limitations,we developed a non-invasive,scalable,and data-driven framework that uses synthetic GPR datasets representing diverse snowpack conditions with variable moisture and density profiles.Synthetic 1D time series reflections(A-scans)are generated using finite-difference time-domain simulations in the state-of-the-art electromagnetic simulator gprMax.Principal component analysis(PCA)is applied to compress each A-scan while preserving key features,which significantly improved and enhanced the model training efficiency.Four machine learning models,including random forest,neural network,support vector machine,and eXtreme gradient boosting,are trained on PCA-reduced features.Among these,the neural network model achieved the best performance,with R^(2)>0.97 for permittivity and R 2>0.92 for resistivity.Gaussian noise(signal-to-noise ratio of 6 dB)is introduced to the synthetic data,and then targeted domain adaptation is employed to enhance generalization to field data.The framework is validated on two contrasting GPR transects in the Altay Mountains of the Chinese mainland,representing moist(T750)and wet(G125)snowpack conditions.The neural network model predictions are most consistent with the GPR derived estimates,Snowfork measurements,and snow pit data,achieving volumetric liquid water content deviation of≤1.5% and bulk density error within the range of 30-84 kg m^(-3).The results demonstrate that machine learning-based inversion,supported by realistic simulations and data augmentation enables scalable,non-invasive snowpack characterization with significant applications in hydrological forecasting,snow monitoring,and water resource management.
基金financially supported by the National Natural Science Foundation of China (Nos.22171165 and 22371170)Natural Science Foundation of Shandong Province (No.ZR2022MB080)Scientific and Technological Frontiers in Project of Henan Province(No.242102110192)。
文摘Co-assembling chiral molecules with achiral compounds via non-covalent interactions like areneperfluoroarene(AP) interactions offers an effective approach for fabricating chiral functional materials.Herein,chiral molecules L/D-PF1 and L/D-PF2 with pyrene groups were synthesized and its chiroptical properties upon co-assembly with achiral compound octafluoronaphthalene(OFN) through AP interaction were systemically studied.The co-assembly of L/D-PF1/OFN and L/D-PF2/OFN exhibited distinct chiroptical properties such as circular dichroism(CD) and circularly polarized luminescence(CPL) signals.Chirality transfer from the chirality center of L/D-PF1 and L/D-PF2 to the achiral OFN and chiral amplification were successfully achieved.Besides,no significant CPL signal was observed in the self-assembly of L/DPF1 or L/D-PF2 while co-assembly with OFN exhibited obvious CPL amplification induced by AP interaction.Notably,a reversal CD signal and CPL signal could be observed in L/D-PF2/OFN when the molar ratio changed from 1:1 to 1:2 while not found in L/D-PF1/OFN,indicating that that minor structural changes of molecules could cause large changes in assembly.In addition,a series of computational calculations were conducted to verify the AP interaction between L-PF1/L-PF2 and OFN.This work demonstrated that arene-perfluoroarene interaction could drive chiral transfer,chiral amplification and chiral inversion and provided a new method for the preparation of chiroptical materials.
基金support of the National Natural Science Foundation of China(No.12172023)。
文摘According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing differential equations are derived by the Hamilton's variational principle,and then solved by the iterative Separation-of-Variable(i SOV)method,which are applicable to arbitrary combinations of homogeneous Boundary Conditions(BCs).However,only the simply-support,clamped and cantilever panels are considered in this work for the sake of clarity.With the closed-form eigensolutions,the flutter frequency,flutter mode and flutter boundary are presented,and the effect of shear deformation and aerodynamic damping on flutter frequencies is investigated.Besides,the relation between panel energy and the work of aerodynamic load is discussed.The numerical comparisons reveal the following.(A)The flutter eigenvalues obtained by the present method are accurate,validated by the Finite Element Method(FEM)and the Galerkin method.(B)When the span-chord ratio is larger than 3,simplifying a 3D panel to 2D(two-dimensional)panel is reasonable and the relative differences of the flutter points predicted by the two models are less than one percent.(C)The reciprocal relationship between the mechanical energy of the panel and the work done by aerodynamic load is verified by using the present flutter eigenvalues and modes,further indicating the high accuracy of the present solutions.(D)The coupling of shear deformation and aerodynamic damping prevents frequency coalescing.
基金supported by the Key Laboratory of Geological Survey and Evaluation of Ministry of Education (China University of Geosciences)(No. GLAB2020ZR13)
文摘The gravity inversion is to restore genetic density distribution of the underground target to be explored for explaining the internal structure and distribution of the Earth.In this paper,we propose a new 3D gravity inversion method based on 3D U-Net++.Compared with two-dimensional gravity inversion,three-dimensional(3D)gravity inversion can more precisely describe the density distribution of underground space.However,conventional 3D gravity inversion method input is two-dimensional,the input and output of the network proposed in our method are three-dimensional.In the training stage,we design a large number of diversifi ed simulation model-data pairs by using the random walk method to improve the generalization ability of the network.In the test phase,we verify the network performance by using the model-data pairs generated by the simulation.To further illustrate the eff ectiveness of the algorithm,we apply the method to the inversion of the San Nicolas mining area,and the inversion results are basically consistent with the borehole measurement results.Moreover,the results of the 3D U-Net++inversion and the 3D U-Net inversion are compared.The density models of the 3D U-Net++inversion have higher resolution,more concentrated inversion results,and a clearer boundary of the density model.
基金supported by the Development of airborne gravity gradiometer(No.2017YFC0601601)open subject of Key Laboratory of Petroleum Resources Research,Institute of Geology and Geophysics,Chinese Academy of Sciences(No.KLOR2018-8)
文摘This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zero norm solution. The inversion approach mainly employs forward modeling; a depth weight function is introduced into the objective function of the zero norms. Sparse inversion results are obtained by the corresponding optimal mathematical method. To achieve the practical geophysical and geological significance of the results, penalty function is applied to constrain the density values. Results obtained by proposed provide clear boundary depth and density contrast distribution information. The method's accuracy, validity, and reliability are verified by comparing its results with those of synthetic models. To further explain its reliability, a practical gravity data is obtained for a region in Texas, USA is applied. Inversion results for this region are compared with those of previous studies, including a research of logging data in the same area. The depth of salt dome obtained by the inversion method is 4.2 km, which is in good agreement with the 4.4 km value from the logging data. From this, the practicality of the inversion method is also validated.
文摘Based on three-dimensional joint finite element, this paper discusses the theory and methodology about inversionof geodetic data. The FEM and inversion formula is given in detail; also a related code is developed. By use of theGreen's function about 3-D FEM, we invert geodetic measurementS of coseismic deformation of the 1989 Ms=7. 1Loma Prieta earthquake to datermine itS source mechanism. The result indicates that the slip on the fault plane isvery heterogeneous. The maximum slip and shear stress are located about 10 kin to northwest of the eathquakesource, the stress drop is about more than 1 MPa.
基金supported by the Joint Fund of Seismological Science(Grant No.U1839206)the National R&D Program on Monitoring,Early Warning and Prevention of Major Natural Disaster(Grant No.2017YFC1500301)+2 种基金supported by IGGCAS Research Start-up Funds(Grant No.E0515402)National Natural Science Foundation of China(Grant No.E1115401)supported by National Natural Science Foundation of China(Grant No.11971258).
文摘The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.
基金supported by the National Natural Science Foundation of China (Nos. 40774029, 41004028)the Special Fund for Basic Scientific Research of Central Colleges (No. 2010ZY53)the Program for New Century Excellent Talents in University (NCET)
文摘We developed a three-dimensional（3D） conjugate gradient inversion algorithm for in-verting magnetotelluric impedance tensor measurements.In order to show the importance of including diagonal components of magnetotelluric impedance tensor in 3D inversion,synthetic data were inverted using the 3D conjugate gradient inversion,and the inversion results were compared and analyzed.The results from the 3D inversion of synthetic data indicate that both the off-diagonal and the diagonal components are required in inversions to obtain better inversion results when there are no enough data sites to recover the target resistivity structure.These examples show that lots of information about 3D structure is also contained in the diagonal components;as a result,diagonal components should be in-cluded in 3D inversions.The inversion algorithm was also used to invert the impedance tensor data ac-quired in the Kayabe area in Japan.Inversions with the synthetic and real data demonstrated the va-lidity and practicability of the inversion algorithm.
基金financially supported by the National Natural Science Foundation of China(No.90814009)
文摘A three-dimensional crustal density model beneath North China is determined using P-wave traveltimes and gravity datum by sequential inversion method in this paper. To improve the uniqueness of the solution, we used a linear relationship between velocity and density to achieve mutual conversions and constraints between velocity difference and density difference. Algebraic reconstruction technique （ART） was used in density inversion, which highly improved the calculation speed comparing with common least squares method. The inversion results indicate that the crustal density beneath North China is extremely inhomogeneous with its distribution coherent with regional tectonics. The lower crust of Taihang mountain tectonic zone shows an obvious low density characteristic. We proposed that it may be an indicator of upwelling of mantle materials or partial melting of lowermost crust.
基金supported by the National Hi-tech Research and Development Program of China(863Program)(No.2007AA09Z310) National Natural Science Foundation of China(Grant No.40774029 40374024)+1 种基金 the Fundamental Research Funds for the Central Universities(Grant No.2010ZY53) the Program for New Century Excellent Talents in University(NCET)
文摘Based on the analysis of impedance tensor data, tipper data, and the conjugate gradient algorithm, we develop a three-dimensional (3D) conjugate gradient algorithm for inverting magnetotelluric full information data determined from five electric and magnetic field components and discuss the method to use the full information data for quantitative interpretation of 3D inversion results. Results from the 3D inversion of synthetic data indicate that the results from inverting full information data which combine the impedance tensor and tipper data are better than results from inverting only the impedance tensor data (or tipper data) in improving resolution and reliability. The synthetic examples also demonstrate the validity and stability of this 3D inversion algorithm.
基金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.
基金National Natural Science Foundation under Grant No.51574201Opening Fund of State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology)under Grant No.SKLGP2016K017+2 种基金Open Research Fund by Sichuan Engineering Research Center for Emergency Mapping&Disaster Reduction under Grant No.K2015B008The State Administration of Work Safety under Grant No.2014_3335Soft Science Research Projects in Sichuan Province under Grant No.2015zr0049
文摘This study explores the quasi-real time inversion principle and precision estimation of three-dimensional coordinates of the epicenter, trigger time and magnitude of earthquakes with the aim to improve traditional methods, which are flawed due to missing information or distortion in the seismograph records. The epicenter, trigger time and magnitude from the Lushan earthquake are inverted and analyzed based on high-frequency GNSS data. The inversion results achieved a high precision, which are consistent with the data published by the China Earthquake Administration. Moreover, it has been proven that the inversion method has good theoretical value and excellent application prospects.
文摘Abstract. Due to the incapability of high frequency surface wave radar （HFSWR） to detect undercurrent parameters, a new algorithm is proposed to apply a three-dimensional （3D） nonlinear barotropic shallow sea and continental shelf model in coordinate system to the inversion of undercurrent. The calculation domain of this model is the area detected by HFSWR. Considering the benthal topography of the detected area and the ocean dynamic parameters, such as surface current, wind and wave detected by HFSWR, the relation between surface current and undercurrent is established in this model, Accordingly, the undercurrent parameters of corresponding area are obtained. The inversion results agree with the law of ocean dynamics and reproduce the movement of undercurrent.
基金the National Key Research and Development Program of China(No.2016YFC060110403).
文摘Herein,a three-dimensional(3D)inversion method in the frequency domain based on a time–frequency transformation was developed to improve the efficiency of the 3D inversion of transient electromagnetic(TEM)data.The Fourier transform related to the electromagnetic response in the frequency and time domains becomes a sine or cosine transform under the excitation of downward-step current.We established a transformation matrix based on the digital fi ltering calculation for the sine transform,and then the frequency domain projection of the TEM data was determined from the linear transformation system using the smoothing constrained least squares inversion method,in which only the imaginary part was used to maintain the TEM data transformation equivalence in the bidirectional projection.Thus,the time-domain TEM inversion problem was indirectly and effectively solved in the frequency domain.In the 3D inversion of the transformed frequency-domain data,the limited-memory Broyden–Fletcher–Goldfarb–Shannoquasi–Newton(L-BFGS)method was used and modifi ed with a restart strategy to adjust the regularization parameter when the algorithm tended to a local minimum.Synthetic data tests showed that our domain transformation method can stably project the TEM data into the frequency domain with very high accuracy;furthe rmore,the 3D inversion of the transformed frequency-domain data is stable,can be used to recover the real resistivity model with an acceptable effi ciency.
