The forward gated-diode method is used to extract the dielectric oxide thickness and body doping concentration of MOSFETs, especially when both of the variables are unknown previously. First, the dielectric oxide thic...The forward gated-diode method is used to extract the dielectric oxide thickness and body doping concentration of MOSFETs, especially when both of the variables are unknown previously. First, the dielectric oxide thickness and the body doping concentration as a function of forward gated-diode peak recombination-generation （R-G） current are derived from the device physics. Then the peak R-G current characteristics of the MOSFETs with different dielectric oxide thicknesses and body doping concentrations are simulated with ISE-Dessis for parameter extraction. The results from the simulation data demonstrate excellent agreement with those extracted from the forward gated-diode method.展开更多
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.展开更多
In this study,we propose a three-dimensional(3D)forward modeling algorithm of surface-to-borehole transient electromagnetic(SBTEM)fields based on an unstructured vector fi nite-element method to analyze the characteri...In this study,we propose a three-dimensional(3D)forward modeling algorithm of surface-to-borehole transient electromagnetic(SBTEM)fields based on an unstructured vector fi nite-element method to analyze the characteristics of SBTEM responses for complex geoelectrical models.To solve the double-curl diff usion equation for the electric fi eld,we use an unstructured tetrahedral mesh to discretize the model domain and select the unconditionally stable backward Euler scheme to discretize the time derivative.In our numerical experiments,we use a grounded wire as a transmitting source.After validating the algorithm’s eff ectiveness,we first analyze the diffusion characteristics and detectability of the electromagnetic field.After that,we focus our attention on the distribution and the cause of zero bands for Ex and dBy/dt components with the hope of guiding future field surveys.Finally,by simulating diff erent models,we analyze the capability of the SBTEM method in detecting typical mineral veins so that we can provide a reference for mineral resource exploration in the deep earth.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been deve...A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.展开更多
The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of c...The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.展开更多
With the development of distribution automation system, the centralized meter reading system has been adopted more and more extensively, which provides real-time electricity consumption data of end-users, and conseque...With the development of distribution automation system, the centralized meter reading system has been adopted more and more extensively, which provides real-time electricity consumption data of end-users, and consequently lays foundation for operating condition on-line analysis of distribution network. In this paper, a modified back/forward sweep method, which directly uses real-time electricity consumption data acquired from the centralized meter reading system, is proposedto realize voltage analysis based on 24-hour electricity consumption data of a typical transformer district. Furthermore, the calculated line losses are verified through data collected from the energy metering of the distribution transformer, illustrating that the proposed method can be applied in analyzing voltage level and discovering unknown energy losses, which will lay foundation for on-line analysis, calculation and monitoring of power distribution network.展开更多
It is one of the most important part to build an accurate gravity model in geophysical exploration.Traditional gravity modelling is usually based on grid method,such as difference method and finite element method wide...It is one of the most important part to build an accurate gravity model in geophysical exploration.Traditional gravity modelling is usually based on grid method,such as difference method and finite element method widely used.Due to self-adaptability lack of division meshes and the difficulty of high-dimensional calculation.展开更多
This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of com...This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.展开更多
This paper proposes a hybrid method based on the forward-backward method （FBM） and the reciprocity theorem （RT） for evaluating the scattering field from dielectric rough surface with a 2D target above it. Here, th...This paper proposes a hybrid method based on the forward-backward method （FBM） and the reciprocity theorem （RT） for evaluating the scattering field from dielectric rough surface with a 2D target above it. Here, the equivalent electric/magnetic current densities on the rough surface as well as the scattering field from it are numerically calculated by FBM, and the scattered field from the isolated target is obtained utilizing the method of moments （MOM）. Meanwhile, the rescattered coupling interactions between the target and the surface are evaluated employing the combination of FBM and RT. Our hybrid method is first validated by available MOM results. Then, the functional dependences of bistatic and monostatic scattering from the target above rough surface upon the target altitude, incident and scattering angles are numerically simulated and discussed. This study presents a numerical description for the scattering mechanism associated with rescattered coupling interactions between a target and an underlying randomly rough surface.展开更多
The hydrodynamic performance of a high forward-speed ship in obliquely propagating waves is numerically examined to assess both free motions and wave field in comparison with a low forward-speed ship.This numerical mo...The hydrodynamic performance of a high forward-speed ship in obliquely propagating waves is numerically examined to assess both free motions and wave field in comparison with a low forward-speed ship.This numerical model is based on the time-domain potential flow theory and higher-order boundary element method,where an analytical expression is completely expanded to determine the base-unsteady coupling flow imposed on the moving condition of the ship.The ship in the numerical model may possess different advancing speeds,i.e.stationary,low speed,and high speed.The role of the water depth,wave height,wave period,and incident wave angle is analyzed by means of the accurate numerical model.It is found that the resonant motions of the high forward-speed ship are triggered by comparison with the stationary one.More specifically,a higher forward speed generates a V-shaped wave region with a larger elevation,which induces stronger resonant motions corresponding to larger wave periods.The shoaling effect is adverse to the motion of the low-speed ship,but is beneficial to the resonant motion of the high-speed ship.When waves obliquely propagate toward the ship,the V-shaped wave region would be broken due to the coupling effect between roll and pitch motions.It is also demonstrated that the maximum heave motion occurs in beam seas for stationary cases but occurs in head waves for high speeds.However,the variation of the pitch motion with period is hardly affected by wave incident angles.展开更多
A finite difference method is introduced to solve the forward-backward heat equation in two space dimensions. In this procedure, the backward and forward difference scheme in two subdomains and a coarse-mesh second-or...A finite difference method is introduced to solve the forward-backward heat equation in two space dimensions. In this procedure, the backward and forward difference scheme in two subdomains and a coarse-mesh second-order central difference scheme at the middle interface are used. Maximum norm error estimate for the procedure is derived. Then an iterative method based on domain decomposition is presented for the numerical scheme and the convergence of the given method is established. Then numerical experiments are presented to support the theoretical analysis.展开更多
As a fiber sensor, optical time domain L reflectometer becomes more and more popular to measure parameters, such as strain and temperature in structural health monitoring （SHM） simultaneously. Since the accuracy of...As a fiber sensor, optical time domain L reflectometer becomes more and more popular to measure parameters, such as strain and temperature in structural health monitoring （SHM） simultaneously. Since the accuracy of range resolution in optical time domain reflectometer （OTDR） is determined by the pulse width of laser, the range resolution in order of centimeter is achieved by employing of picoseconds lasers which are not commercial. In this paper, to achieve this accuracy with conventional OTDR, Fourier wavelet regularized deconvolution （ForWaRD） method is employed to deconvolve and denoise the detected signal simultaneously. Simulations show that this method improves resolution of conventional OTDR system to the order of several centimeters.展开更多
The airborne electromagnetic （AEM） method has a high sampling rate and survey flexibility. However, traditional numerical modeling approaches must use high-resolution physical grids to guarantee modeling accuracy, e...The airborne electromagnetic （AEM） method has a high sampling rate and survey flexibility. However, traditional numerical modeling approaches must use high-resolution physical grids to guarantee modeling accuracy, especially for complex geological structures such as anisotropic earth. This can lead to huge computational costs. To solve this problem, we propose a spectral-element （SE） method for 3D AEM anisotropic modeling, which combines the advantages of spectral and finite-element methods. Thus, the SE method has accuracy as high as that of the spectral method and the ability to model complex geology inherited from the finite-element method. The SE method can improve the modeling accuracy within discrete grids and reduce the dependence of modeling results on the grids. This helps achieve high-accuracy anisotropic AEM modeling. We first introduced a rotating tensor of anisotropic conductivity to Maxwell＇s equations and described the electrical field via SE basis functions based on GLL interpolation polynomials. We used the Galerkin weighted residual method to establish the linear equation system for the SE method, and we took a vertical magnetic dipole as the transmission source for our AEM modeling. We then applied fourth-order SE calculations with coarse physical grids to check the accuracy of our modeling results against a 1D semi-analytical solution for an anisotropic half-space model and verified the high accuracy of the SE. Moreover, we conducted AEM modeling for different anisotropic 3D abnormal bodies using two physical grid scales and three orders of SE to obtain the convergence conditions for different anisotropic abnormal bodies. Finally, we studied the identification of anisotropy for single anisotropic abnormal bodies, anisotropic surrounding rock, and single anisotropic abnormal body embedded in an anisotropic surrounding rock. This approach will play a key role in the inversion and interpretation of AEM data collected in regions with anisotropic geology.展开更多
基金Project supported by the Key Project of the National Natural Science Foundation of China(No.60936005)the Shenzhen Science & Technology Foundation,China(No.JSA200903160146A)+1 种基金the Industry,Education and Academy Cooperation Program of Guangdong Province,China(No.2009B090300318)the Fundamental Research Project of Shenzhen Science & Technology Foundation,China (No.JC200903160353A)
文摘The forward gated-diode method is used to extract the dielectric oxide thickness and body doping concentration of MOSFETs, especially when both of the variables are unknown previously. First, the dielectric oxide thickness and the body doping concentration as a function of forward gated-diode peak recombination-generation （R-G） current are derived from the device physics. Then the peak R-G current characteristics of the MOSFETs with different dielectric oxide thicknesses and body doping concentrations are simulated with ISE-Dessis for parameter extraction. The results from the simulation data demonstrate excellent agreement with those extracted from the forward gated-diode method.
基金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.
基金This work was supported by the National Natural Science Foundation of China(Grant Nos..42030806,41774125,41804098,41904104)the Key National Research Project of China(Grant No.2018YFC0603300).
文摘In this study,we propose a three-dimensional(3D)forward modeling algorithm of surface-to-borehole transient electromagnetic(SBTEM)fields based on an unstructured vector fi nite-element method to analyze the characteristics of SBTEM responses for complex geoelectrical models.To solve the double-curl diff usion equation for the electric fi eld,we use an unstructured tetrahedral mesh to discretize the model domain and select the unconditionally stable backward Euler scheme to discretize the time derivative.In our numerical experiments,we use a grounded wire as a transmitting source.After validating the algorithm’s eff ectiveness,we first analyze the diffusion characteristics and detectability of the electromagnetic field.After that,we focus our attention on the distribution and the cause of zero bands for Ex and dBy/dt components with the hope of guiding future field surveys.Finally,by simulating diff erent models,we analyze the capability of the SBTEM method in detecting typical mineral veins so that we can provide a reference for mineral resource exploration in the deep earth.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金Projects(41674080,41674079)supported by the National Natural Science Foundation of China
文摘A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.
文摘The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.
文摘With the development of distribution automation system, the centralized meter reading system has been adopted more and more extensively, which provides real-time electricity consumption data of end-users, and consequently lays foundation for operating condition on-line analysis of distribution network. In this paper, a modified back/forward sweep method, which directly uses real-time electricity consumption data acquired from the centralized meter reading system, is proposedto realize voltage analysis based on 24-hour electricity consumption data of a typical transformer district. Furthermore, the calculated line losses are verified through data collected from the energy metering of the distribution transformer, illustrating that the proposed method can be applied in analyzing voltage level and discovering unknown energy losses, which will lay foundation for on-line analysis, calculation and monitoring of power distribution network.
基金provided by China Geological Survey with the project(Nos.DD20190707,DD20190012)the Fundamental Research Funds for China Central public research Institutes with the project(No.JKY202014)
文摘It is one of the most important part to build an accurate gravity model in geophysical exploration.Traditional gravity modelling is usually based on grid method,such as difference method and finite element method widely used.Due to self-adaptability lack of division meshes and the difficulty of high-dimensional calculation.
基金This work was supported by China Postdoctoral Science Foundation(No.2022M723391)the Science and Technology Innovation Project of Higher Education in Shanxi Province(No.2019L0754)+1 种基金the Central Guiding Local Science and Technology Development Fund Project(No.YDZJSX2021B021)Shanxi Province Basic Research Plan General Project(No.202203021221294).
文摘This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.
基金Project supported by the National Natural Science Foundation of China (Grant No 60571058)the National Defense Foundation of China
文摘This paper proposes a hybrid method based on the forward-backward method （FBM） and the reciprocity theorem （RT） for evaluating the scattering field from dielectric rough surface with a 2D target above it. Here, the equivalent electric/magnetic current densities on the rough surface as well as the scattering field from it are numerically calculated by FBM, and the scattered field from the isolated target is obtained utilizing the method of moments （MOM）. Meanwhile, the rescattered coupling interactions between the target and the surface are evaluated employing the combination of FBM and RT. Our hybrid method is first validated by available MOM results. Then, the functional dependences of bistatic and monostatic scattering from the target above rough surface upon the target altitude, incident and scattering angles are numerically simulated and discussed. This study presents a numerical description for the scattering mechanism associated with rescattered coupling interactions between a target and an underlying randomly rough surface.
基金supported by the National Natural Science Foundation of China(Grant Nos.52271278 and 52111530137)the Natural Science Foundation of Jiangsu Province(Grant No.SBK2022020579)the Newton Advanced Fellowships by the Royal Society(Grant No.NAF\R1\180304).
文摘The hydrodynamic performance of a high forward-speed ship in obliquely propagating waves is numerically examined to assess both free motions and wave field in comparison with a low forward-speed ship.This numerical model is based on the time-domain potential flow theory and higher-order boundary element method,where an analytical expression is completely expanded to determine the base-unsteady coupling flow imposed on the moving condition of the ship.The ship in the numerical model may possess different advancing speeds,i.e.stationary,low speed,and high speed.The role of the water depth,wave height,wave period,and incident wave angle is analyzed by means of the accurate numerical model.It is found that the resonant motions of the high forward-speed ship are triggered by comparison with the stationary one.More specifically,a higher forward speed generates a V-shaped wave region with a larger elevation,which induces stronger resonant motions corresponding to larger wave periods.The shoaling effect is adverse to the motion of the low-speed ship,but is beneficial to the resonant motion of the high-speed ship.When waves obliquely propagate toward the ship,the V-shaped wave region would be broken due to the coupling effect between roll and pitch motions.It is also demonstrated that the maximum heave motion occurs in beam seas for stationary cases but occurs in head waves for high speeds.However,the variation of the pitch motion with period is hardly affected by wave incident angles.
基金Supported by National Science Foundation of China(Grant 10871179)the National Basic Research Programme of China(Grant 2008CB717806)the Department of Education of Zhejiang Province(GrantY200803559).
文摘A finite difference method is introduced to solve the forward-backward heat equation in two space dimensions. In this procedure, the backward and forward difference scheme in two subdomains and a coarse-mesh second-order central difference scheme at the middle interface are used. Maximum norm error estimate for the procedure is derived. Then an iterative method based on domain decomposition is presented for the numerical scheme and the convergence of the given method is established. Then numerical experiments are presented to support the theoretical analysis.
文摘As a fiber sensor, optical time domain L reflectometer becomes more and more popular to measure parameters, such as strain and temperature in structural health monitoring （SHM） simultaneously. Since the accuracy of range resolution in optical time domain reflectometer （OTDR） is determined by the pulse width of laser, the range resolution in order of centimeter is achieved by employing of picoseconds lasers which are not commercial. In this paper, to achieve this accuracy with conventional OTDR, Fourier wavelet regularized deconvolution （ForWaRD） method is employed to deconvolve and denoise the detected signal simultaneously. Simulations show that this method improves resolution of conventional OTDR system to the order of several centimeters.
基金financially supported by the Key Program of National Natural Science Foundation of China(No.41530320)China Natural Science Foundation for Young Scientists(No.41404093)+1 种基金Key National Research Project of China(Nos2016YFC0303100 and 2017YFC0601900)China Natural Science Foundation(No.41774125)
文摘The airborne electromagnetic （AEM） method has a high sampling rate and survey flexibility. However, traditional numerical modeling approaches must use high-resolution physical grids to guarantee modeling accuracy, especially for complex geological structures such as anisotropic earth. This can lead to huge computational costs. To solve this problem, we propose a spectral-element （SE） method for 3D AEM anisotropic modeling, which combines the advantages of spectral and finite-element methods. Thus, the SE method has accuracy as high as that of the spectral method and the ability to model complex geology inherited from the finite-element method. The SE method can improve the modeling accuracy within discrete grids and reduce the dependence of modeling results on the grids. This helps achieve high-accuracy anisotropic AEM modeling. We first introduced a rotating tensor of anisotropic conductivity to Maxwell＇s equations and described the electrical field via SE basis functions based on GLL interpolation polynomials. We used the Galerkin weighted residual method to establish the linear equation system for the SE method, and we took a vertical magnetic dipole as the transmission source for our AEM modeling. We then applied fourth-order SE calculations with coarse physical grids to check the accuracy of our modeling results against a 1D semi-analytical solution for an anisotropic half-space model and verified the high accuracy of the SE. Moreover, we conducted AEM modeling for different anisotropic 3D abnormal bodies using two physical grid scales and three orders of SE to obtain the convergence conditions for different anisotropic abnormal bodies. Finally, we studied the identification of anisotropy for single anisotropic abnormal bodies, anisotropic surrounding rock, and single anisotropic abnormal body embedded in an anisotropic surrounding rock. This approach will play a key role in the inversion and interpretation of AEM data collected in regions with anisotropic geology.
