The quasi-monochromatic,continuously energy-tunable,and high-brightness gamma rays that are produced by an inverse Compton scattering(ICS)light source provide an ideal probe for gamma-ray imaging.However,owing to the ...The quasi-monochromatic,continuously energy-tunable,and high-brightness gamma rays that are produced by an inverse Compton scattering(ICS)light source provide an ideal probe for gamma-ray imaging.However,owing to the influence of the intrinsic energy-angle correlation spectrum of this type of light source,monochromatic computed tomography(CT),especially in the gamma-ray energy region,can only be realized in a low-efficiency manner,similar to first-generation CT.A dual-energy scan scheme with a large imaging field of view(FOV)was developed in this study to improve the imaging efficiency.The effectiveness of this scheme was demonstrated based on the beam parameters of a typical ICS light source using Monte Carlo simulations.By leveraging the principle of basis material decomposition,the influence of the energyangle correlation spectrum on CT reconstruction was corrected,and a monochromatic CT image of the imaging object was accurately reconstructed.Furthermore,the electron density and effective atomic number of the imaging object could be obtained simultaneously.展开更多
An inverse analysis is presented to estimate line heat source in two-dimensional steady-state and transient heat transfer problems.A constant heat source is considered in the steady-state heat transfer problem(a param...An inverse analysis is presented to estimate line heat source in two-dimensional steady-state and transient heat transfer problems.A constant heat source is considered in the steady-state heat transfer problem(a parameter estimation problem)and a time-varying heat source is considered in the transient heat transfer problem(a function estimation problem).Since a general irregular 2D heat conducting body is considered,a body-fitted grid generation is used to mesh the domain.Then governing equations and associated boundary and initial conditions are transformed from the physical domain to the computational domain and finite difference method is used to solve the governing equations to obtain the temperature distribution in the body.Using an efficient,accurate,and very easy to implement sensitivity analysis incorporated in a gradient based minimization method(here,steepest descentmethod),the unknown heat source is estimated accurately.In the function estimation part,it is assumed that there is no prior information on the functional form of the heat source and the estimation process can be performed with a reasonable initial guess for the heat source.The main advantage of the proposed inverse analysis is that the sensitivity matrix(and hence,the objective function gradient with respect to the unknown variables)can be computed during the direct heat transfer solution through newyet simple explicit expressions with no need to solve extra equations such as the sensitivity and adjoint problems and impose additional computational costs comparable to the direct problem solution ones.Some test cases are presented to investigate the accuracy,efficiency,and effect of measurement error on the estimated parameter and function for the line heat source.展开更多
The complexity of the seismicity pattern for the subduction zone along the oceanic plate triggered the outer rise events and revealed cyclic tectonic deformation conditions along the plate subduction zones.The outer r...The complexity of the seismicity pattern for the subduction zone along the oceanic plate triggered the outer rise events and revealed cyclic tectonic deformation conditions along the plate subduction zones.The outer rise earthquakes have been observed along the Sunda arc,following the estimated rupture area of the 2005 M_(W)8.6 Nias earthquakes.Here,we used kinematic waveform inversion(KIWI)to obtain the source parameters of the 14 May 2021 M_(W)6.6 event off the west coast of northern Sumatra and to define the fault plane that triggered this outer rise event.The KIWI algorithm allows two types of seismic source to be configured:the moment tensor model to describe the type of shear with six moment tensor components and the Eikonal model for the rupture of pure double-couple sources.This method was chosen for its flexibility to be applied for different sources of seismicity and also for the automated full-moment tensor solution with real-time monitoring.We used full waveform traces from 8 broadband seismic stations within 1000 km epicentral distances sourced from the Incorporated Research Institutions for Seismology(IRIS-IDA)and Geofon GFZ seismic record databases.The initial origin time and hypocenter values are obtained from the IRIS-IDA.The synthetic seismograms used in the inversion process are based on the existing regional green function database model and were accessed from the KIWI Tools Green's Function Database.The obtained scalar seismic moment value is 1.18×10^(19)N·m,equivalent to a moment magnitude M_(W)6.6.The source parameters are 140°,44°,and−99°for the strike,dip,and rake values at a centroid depth of 10.2 km,indicating that this event is a normal fault earthquake that occurred in the outer rise area.The outer rise events with normal faults typically occur at the shallow part of the plate,with nodal-plane dips predominantly in the range of 30°-60°on the weak oceanic lithosphere due to hydrothermal alteration.The stress regime around the plate subduction zone varies both temporally and spatially due to the cyclic influences of megathrust earthquakes.Tensional outer rise earthquakes tend to occur after the megathrust events.The relative timing of these events is not known due to the viscous relaxation of the down going slab and poroelastic response in the trench slope region.The occurrence of the 14 May 2021 earthquake shows the seismicity in the outer rise region in the strongly coupled Sunda arc subduction zone due to elastic bending stress within the duration of the seismic cycle.展开更多
With the continuous improvement of the accuracy of geodetic deformation data,the inversion of seismic source parameters puts forward a higher demand for nonlinear inversion algorithms.In this research,an improved Spar...With the continuous improvement of the accuracy of geodetic deformation data,the inversion of seismic source parameters puts forward a higher demand for nonlinear inversion algorithms.In this research,an improved Sparrow Search Algorithm(SSA)is proposed for the seismic source parameter inversion problem.By replacing the original population generation in the improved algorithm with Latin hypercubic sampling,the Sparrow Search Algorithm reduces the repetition of samples in the population initialization.Subsequently,the algorithm introduces adaptive weights in the discoverer generation phase of the sparrow algorithm and combines the Levy flight strategy to make the algorithm more comprehensive and improve the search accuracy during the whole iteration process.Therefore,the improved Latin hypercube-based sparrow search algorithm(ILHSSA)has better advantages in terms of iterative convergence speed and stability.In order to verify the performance of ILHSSA,the basic genetic algorithm(GA)and sparrow search algorithm(SSA)are examined and compared with ILHSSA by simulated earthquakes of two different earthquake types.The simulation experiments show that the improved algorithm ILHSSA outperforms SSA in accuracy and stability.Compared with the GA algorithm,ILHSSA can achieve the same inversion accuracy as GA,and it even surpasses GA in inversion speed and the inversion results of some parameters,demonstrating better stability.Finally,the improved algorithm is used for the 2017 Bodrum-Cos earthquake and the 2016 Amatrice earthquake in Italy.The inversion results all reflect the practicality and reliability of the improved algorithm.展开更多
This paper proposes a current control scheme for a grid-connected pulse width modulator(PWM) voltage source converter(GC-VSC) under imbalanced and distorted supply voltage conditions.The control scheme is implemented ...This paper proposes a current control scheme for a grid-connected pulse width modulator(PWM) voltage source converter(GC-VSC) under imbalanced and distorted supply voltage conditions.The control scheme is implemented in the positive synchronously rotating reference frame and composed of a single proportional integral(PI) regulator and multi-frequency resonant controllers tuned at the frequencies of 2ω and 6ω,respectively.The experimental results,with the target of eliminating the active power oscillations and current harmonics on a prototype GC-VSC system,validate the feasibility of the proposed current control scheme during supply voltage imbalance and distortion.展开更多
Two cardiac functional models are constructed in this paper. One is a single current model and the other is a current multipole model. Parameters denoting the properties of these two models are calculated by a least-s...Two cardiac functional models are constructed in this paper. One is a single current model and the other is a current multipole model. Parameters denoting the properties of these two models are calculated by a least-square fit to the measurements using a simulated annealing algorithm. The measured signals are detected at 36 observation nodes by a superconducting quantum interference device （SQUID）. By studying the trends of position, orientation and magnitude of the single current dipole model and the current multipole model in the QRS complex during one time span and comparing the reconstructed magnetocardiography （MCG） of these two cardiac models, we find that the current multipole model is a more appropriate model to represent cardiac electrophysiological activity.展开更多
N-soliton solutions of the hierarchy of non-isospectral mKdV equation with self-consistent sources andthe hierarchy of non-isospectral sine-Gordon equation with self-consistent sources are obtained via the inverse sca...N-soliton solutions of the hierarchy of non-isospectral mKdV equation with self-consistent sources andthe hierarchy of non-isospectral sine-Gordon equation with self-consistent sources are obtained via the inverse scatteringtransform.展开更多
In reactor physics tests, it is important to monitor sub-criticality continuously during criticality approach. Reactivity measurements by the inverse kinetics method are widely used during the operation of a nuclear r...In reactor physics tests, it is important to monitor sub-criticality continuously during criticality approach. Reactivity measurements by the inverse kinetics method are widely used during the operation of a nuclear reactor. This technique is successfully applied at sufficiently high power level or to a core without an external neutron source where the neutron source term in point reactor kinetics equations may be neglected. For operation at low power levels or in the sub-critical domain, the increase in the fluctuation of the neutron signal may cause difficulties in the evaluation of reactivity and the effect of direct emission from the external neutron source may not be neglected. Therefore, contribution of the neutron source must be taken into account and this implies knowledge of a quantity proportional to the source strength, which calls the source term and then it should be determined. The research work has been conducted to measure reactivity with source term using a dedicated reactivity measurement system by the Least Square Inverse Kinetics Method (LSIKM). Application to a simulator of HANARO research reactor, Korea Atomic Energy Research Institute (KAERI), with known source strength and reactivity worth has showed consistent and satisfactory agreement.展开更多
The elastic plate vibration model is studied under the external force. The size of the source term by the given mode of the source and some observations from the body of the plate is determined over a time interval, w...The elastic plate vibration model is studied under the external force. The size of the source term by the given mode of the source and some observations from the body of the plate is determined over a time interval, which is referred to be an inverse source problem of a plate equation. The uniqueness theorem for this problem is stated, and the fundamental solution to the plate equation is derived. In the case that the plate is driven by the harmonic load, the fundamental solution method （FSM） and the Tikhonov regularization technique axe used to calculate the source term. Numerical experiments of the Euler-Bernoulli beam and the Kirchhoff-Love plate show that the FSM can work well for practical use, no matter the source term is smooth or piecewise.展开更多
In this article,we consider the diffusion equation with multi-term time-fractional derivatives.We first derive,by a subordination principle for the solution,that the solution is positive when the initial value is non-...In this article,we consider the diffusion equation with multi-term time-fractional derivatives.We first derive,by a subordination principle for the solution,that the solution is positive when the initial value is non-negative.As an application,we prove the uniqueness of solution to an inverse problem of determination of the temporally varying source term by integral type information in a subdomain.Finally,several numerical experiments are presented to show the accuracy and efficiency of the algorithm.展开更多
The inverse estimation of a source location of pollutant released into a turbulent flow is a probability problem instead of a deterministic one, as the turbulent flow is chaotic and irreversible. However, researches c...The inverse estimation of a source location of pollutant released into a turbulent flow is a probability problem instead of a deterministic one, as the turbulent flow is chaotic and irreversible. However, researches can be conducted to provide helpful instructions to the possible source location with corresponding uncertainty. This study aims to propose a method of inverse estimation of a passive-scalar source location. Experimental investigation of the dye plume characteristics released into a fully-developed turbulent flow is performed in a water channel. A planar laser-induced fluorescence (PLIF) technique is used to obtain two-dimensional images of spreading dye plumes at a bulk Reynolds number of 20,000. The distributions of high concentration areas in the PLIF images are chosen as features that characterize the traveling (diffusion) distance or time from the dye source. Graphical analysis is used to extract these high concentration areas. The procedure of graphical analysis has three steps: 1) binarization using a threshold to extract high concentration dye patches;2) labeling individual high-concentration dye patches in the binarized images;and 3) pixel-counting to measure the area and perimeter of each dye patch. We examine the variations of fractal dimension of patches, and the fractal dimension is observed to be almost constant irrespective of the distance from the source. The kurtosis of the probability density function curve of the logarithm dimensionless dye patch areas is found to be related with the downstream diffusion distance, based on which an inverse estimation method to locate a passive-scalar point source is proposed and evaluated.展开更多
The inverse scattering problems are to detect the property of obstacles from the measurements outside the obstacles. One of important research areas in this topic is the recovery of boundary property for impenetrable ...The inverse scattering problems are to detect the property of obstacles from the measurements outside the obstacles. One of important research areas in this topic is the recovery of boundary property for impenetrable obstacles. In this paper, we would like to give a brief review about the recently developed singular source methods. There are three different methods in this category, namely, linear sampling method, pointsource method and probe method. We also present some recent new results about the probe method.展开更多
As a high quality seismic imaging method, full waveform inversion （FWI） can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing i...As a high quality seismic imaging method, full waveform inversion （FWI） can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing is computationally expensive, especially for the three-dimension complex medium inversion. Introducing blended source technology into the frequency-domain FWI can greatly reduce the computational burden and improve the efficiency of the inversion. However, this method has two issues： first, crosstalk noise is caused by interference between the sources involved in the encoding, resulting in an inversion result with some artifacts; second, it is more sensitive to ambient noise compared to conventional FWI, therefore noisy data results in a poor inversion. This paper introduces a frequency-group encoding method to suppress crosstalk noise, and presents a frequency- domain auto-adapting FWI based on source-encoding technology. The conventional FWI method and source-encoding based FWI method are combined using an auto-adapting mechanism. This improvement can both guarantee the quality of the inversion result and maximize the inversion efficiency.展开更多
Multiturn coils is an effective transmitter for transient electromagnetic method(TEM) used in narrow space and complex terrain at presently. However, its high mutual inductance coupling and long turn-off time affect t...Multiturn coils is an effective transmitter for transient electromagnetic method(TEM) used in narrow space and complex terrain at presently. However, its high mutual inductance coupling and long turn-off time affect the quality of later data processing and interpretation. Compared with multiturn coils, the new conical source has low mutual inductance and short turn-off time. Based on the superposition principle, we use Hankel transform and numerical filtering method for forward modelling of the conical source field in the layered-media and explore TEM characteristics excited by this source. We apply improved damped least square inversion to integrated transient electromagnetic(TEM) data. We first invert the induced voltage into similar resistivity and apparent depth, and then use the inverted results as input parameters in the initial model and transform the apparent resistivity data into the frequency domain. Then, damped least square inversion is performed in the frequency domain using the initial model. Subsequently, we use automated model building to search for the extremes and inflection points in the resistivity–depth data that are treated as critical layer parameters. The inversion of theoretical and observed data suggests that the method modifies the resistivity and depth and yields a model of the underground layers.展开更多
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.展开更多
We conducted rapid inversions of rupture process for the 2023 earthquake doublet occurred in SE Türkiye,the first with a magnitude of M_(W)7.8 and the second with a magnitude of M_(W)7.6,using teleseismic and str...We conducted rapid inversions of rupture process for the 2023 earthquake doublet occurred in SE Türkiye,the first with a magnitude of M_(W)7.8 and the second with a magnitude of M_(W)7.6,using teleseismic and strong-motion data.The teleseismic rupture models of the both events were obtained approximately 88 and 55 minutes after their occurrences,respectively.The rupture models indicated that the first event was an asymmetric bilateral event with ruptures mainly propagating to the northeast,while the second one was a unilateral event with ruptures propagating to the west.This information could be useful in locating the meizoseismal areas.Compared with teleseismic models,the strong-motion models showed relatively higher resolution.A noticeable difference was found for the M_(W)7.6 earthquake,for which the strong-motion models shows a bilateral event,rather than a unilateral event,but the dominant rupture direction is still westward.Nevertheless,all strong-motion models are consistent with the teleseismic models in terms of magnitudes,durations,and dominant rupture directions.This suggests that both teleseismic and strong-motion data can be used for fast determination of major source characteristics.In contrast,the strong-motion data would be preferable in future emergency responses since they are recorded earlier and have a better resolution ability on the source ruptures.展开更多
The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault par...The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault parameters inversion has nonlinear characteristics, and the gradient-based optimizer(GBO) has the characteristics of fast convergence speed and falling into local optimum hardly. This paper applies GBO algorithm to simulated earthquakes and real LuShan earthquakes in the nonlinear inversion of the Okada model to obtain the source parameters. The simulated earthquake experiment results show that the algorithm is stable, and the seismic source parameters obtained by GBO are slightly closer to the true value than the multi peak particle swarm optimization(MPSO). In the 2013 LuShan earthquake experiment, the root mean square error between the deformation after forwarding of fault parameters obtained by the introduced GBO algorithm and the surface observation deformation was 3.703 mm, slightly better than 3.708 mm calculated by the MPSO. Moreover, the inversion result of GBO algorithm is better than MPSO algorithm in stability. The above results show that the introduced GBO algorithm has a certain practical application value in seismic fault source parameters inversion.展开更多
Based on Fourier transform,the angular spectrum equation and its discrete model of the inverse source problem of the two-dimensional Helmholtz equation are derived.By transforming the two-dimensional inverse source pr...Based on Fourier transform,the angular spectrum equation and its discrete model of the inverse source problem of the two-dimensional Helmholtz equation are derived.By transforming the two-dimensional inverse source problem in the spatial domain into the inverse source problem in the one-dimensional Fourier domain,the source intensity angular spectrum vector in the Fourier domain is calculated by the one-dimensional method of moments.Then,using the discrete Fourier inverse transform,a multi-frequency angular spectrum method of moments is proposed to reconstruct the inverse source problem of the two-dimensional Helmholtz equation.The reconstruction accuracy and stability of this method are improved by optimizing the excitation wavenumber sequence.Simulation calculations show that the multi-frequency angular spectrum method of moments has high reconstruction accuracy and anti-noise ability,and can also design the reconstruction model according to the effective frequency band of the acoustic sensor to achieve the highest resolution.In addition,since the calculation of the coefficient matrix of the proposed method is independent of the distribution of the source intensity function and does not require regularization processing,when the size and resolution of the required source intensity distribution area are determined,after the coefficient matrix is calculated and saved,the source intensity function corresponding to the region of interest can be reconstructed according to different sound pressure detection data,with high calculation accuracy and efficiency.展开更多
A new source location method using wave-equation based traveltime inversion is proposed to locate microseismic events accurately. With a sourceindependent strategy, microseismic events can be located independently reg...A new source location method using wave-equation based traveltime inversion is proposed to locate microseismic events accurately. With a sourceindependent strategy, microseismic events can be located independently regardless of the accuracy of the source signature and the origin time. The traveltime-residuals-based misfit function has robust performance when the velocity model is inaccurate. The new Fréchet derivatives of the misfit function with respect to source location are derived directly based on the acoustic wave equation, accounting for the influence of geometrical perturbation and spatial velocity variation. Unlike the mostly used traveltime inversion methods, no traveltime picking or ray tracing is needed.Additionally, the improved scattering-integral method is applied to reduce the computational cost. Numerical tests show the validity of the proposed method.展开更多
plane with the strike of 127°, the dip of 79° and the rake of 171°. The rupture process inversion result of MS=7.9 earthquake shows that the total rupture duration is about 37 s, the scalar moment tenso...plane with the strike of 127°, the dip of 79° and the rake of 171°. The rupture process inversion result of MS=7.9 earthquake shows that the total rupture duration is about 37 s, the scalar moment tensor is M0=0.97 × 1020 N·m. Rupture mainly occurred on the shallow area with 110 km long and 30 km wide, the location in which the rupture initiated is not where the main rupture took place, and the area with slip greater than 0.5 m basically lies within 35 km deep middle-crust under the earth surface. The maximum static slip is 3.6 m. There are two distinct areas with slip larger than 2.0 m. We noticed that when the rupture propagated towards northwest and closed to the area around the MS=7.3 hypocenter, the slip decreased rapidly, which may indicate that the rupture process was stopped by barriers. The consistence of spatial distribution of slip on the fault plane with the distribution of aftershocks also supports that the rupture is a heterogeneous process owing to the presence of barriers.展开更多
基金supported by the National Natural Science Foundation of China(Nos.12375157,12027902,and 11905011)。
文摘The quasi-monochromatic,continuously energy-tunable,and high-brightness gamma rays that are produced by an inverse Compton scattering(ICS)light source provide an ideal probe for gamma-ray imaging.However,owing to the influence of the intrinsic energy-angle correlation spectrum of this type of light source,monochromatic computed tomography(CT),especially in the gamma-ray energy region,can only be realized in a low-efficiency manner,similar to first-generation CT.A dual-energy scan scheme with a large imaging field of view(FOV)was developed in this study to improve the imaging efficiency.The effectiveness of this scheme was demonstrated based on the beam parameters of a typical ICS light source using Monte Carlo simulations.By leveraging the principle of basis material decomposition,the influence of the energyangle correlation spectrum on CT reconstruction was corrected,and a monochromatic CT image of the imaging object was accurately reconstructed.Furthermore,the electron density and effective atomic number of the imaging object could be obtained simultaneously.
文摘An inverse analysis is presented to estimate line heat source in two-dimensional steady-state and transient heat transfer problems.A constant heat source is considered in the steady-state heat transfer problem(a parameter estimation problem)and a time-varying heat source is considered in the transient heat transfer problem(a function estimation problem).Since a general irregular 2D heat conducting body is considered,a body-fitted grid generation is used to mesh the domain.Then governing equations and associated boundary and initial conditions are transformed from the physical domain to the computational domain and finite difference method is used to solve the governing equations to obtain the temperature distribution in the body.Using an efficient,accurate,and very easy to implement sensitivity analysis incorporated in a gradient based minimization method(here,steepest descentmethod),the unknown heat source is estimated accurately.In the function estimation part,it is assumed that there is no prior information on the functional form of the heat source and the estimation process can be performed with a reasonable initial guess for the heat source.The main advantage of the proposed inverse analysis is that the sensitivity matrix(and hence,the objective function gradient with respect to the unknown variables)can be computed during the direct heat transfer solution through newyet simple explicit expressions with no need to solve extra equations such as the sensitivity and adjoint problems and impose additional computational costs comparable to the direct problem solution ones.Some test cases are presented to investigate the accuracy,efficiency,and effect of measurement error on the estimated parameter and function for the line heat source.
基金supported by the National Natural Science Foundation of China(Grant No.42130312)。
文摘The complexity of the seismicity pattern for the subduction zone along the oceanic plate triggered the outer rise events and revealed cyclic tectonic deformation conditions along the plate subduction zones.The outer rise earthquakes have been observed along the Sunda arc,following the estimated rupture area of the 2005 M_(W)8.6 Nias earthquakes.Here,we used kinematic waveform inversion(KIWI)to obtain the source parameters of the 14 May 2021 M_(W)6.6 event off the west coast of northern Sumatra and to define the fault plane that triggered this outer rise event.The KIWI algorithm allows two types of seismic source to be configured:the moment tensor model to describe the type of shear with six moment tensor components and the Eikonal model for the rupture of pure double-couple sources.This method was chosen for its flexibility to be applied for different sources of seismicity and also for the automated full-moment tensor solution with real-time monitoring.We used full waveform traces from 8 broadband seismic stations within 1000 km epicentral distances sourced from the Incorporated Research Institutions for Seismology(IRIS-IDA)and Geofon GFZ seismic record databases.The initial origin time and hypocenter values are obtained from the IRIS-IDA.The synthetic seismograms used in the inversion process are based on the existing regional green function database model and were accessed from the KIWI Tools Green's Function Database.The obtained scalar seismic moment value is 1.18×10^(19)N·m,equivalent to a moment magnitude M_(W)6.6.The source parameters are 140°,44°,and−99°for the strike,dip,and rake values at a centroid depth of 10.2 km,indicating that this event is a normal fault earthquake that occurred in the outer rise area.The outer rise events with normal faults typically occur at the shallow part of the plate,with nodal-plane dips predominantly in the range of 30°-60°on the weak oceanic lithosphere due to hydrothermal alteration.The stress regime around the plate subduction zone varies both temporally and spatially due to the cyclic influences of megathrust earthquakes.Tensional outer rise earthquakes tend to occur after the megathrust events.The relative timing of these events is not known due to the viscous relaxation of the down going slab and poroelastic response in the trench slope region.The occurrence of the 14 May 2021 earthquake shows the seismicity in the outer rise region in the strongly coupled Sunda arc subduction zone due to elastic bending stress within the duration of the seismic cycle.
基金funded by the National Natural Science Foundation of China(42174011).
文摘With the continuous improvement of the accuracy of geodetic deformation data,the inversion of seismic source parameters puts forward a higher demand for nonlinear inversion algorithms.In this research,an improved Sparrow Search Algorithm(SSA)is proposed for the seismic source parameter inversion problem.By replacing the original population generation in the improved algorithm with Latin hypercubic sampling,the Sparrow Search Algorithm reduces the repetition of samples in the population initialization.Subsequently,the algorithm introduces adaptive weights in the discoverer generation phase of the sparrow algorithm and combines the Levy flight strategy to make the algorithm more comprehensive and improve the search accuracy during the whole iteration process.Therefore,the improved Latin hypercube-based sparrow search algorithm(ILHSSA)has better advantages in terms of iterative convergence speed and stability.In order to verify the performance of ILHSSA,the basic genetic algorithm(GA)and sparrow search algorithm(SSA)are examined and compared with ILHSSA by simulated earthquakes of two different earthquake types.The simulation experiments show that the improved algorithm ILHSSA outperforms SSA in accuracy and stability.Compared with the GA algorithm,ILHSSA can achieve the same inversion accuracy as GA,and it even surpasses GA in inversion speed and the inversion results of some parameters,demonstrating better stability.Finally,the improved algorithm is used for the 2017 Bodrum-Cos earthquake and the 2016 Amatrice earthquake in Italy.The inversion results all reflect the practicality and reliability of the improved algorithm.
基金supported by the National Natural Science Foundation of China(No.50907057)the National High-Tech Research and Development Program (863) of China(No.2007AA05Z419)
文摘This paper proposes a current control scheme for a grid-connected pulse width modulator(PWM) voltage source converter(GC-VSC) under imbalanced and distorted supply voltage conditions.The control scheme is implemented in the positive synchronously rotating reference frame and composed of a single proportional integral(PI) regulator and multi-frequency resonant controllers tuned at the frequencies of 2ω and 6ω,respectively.The experimental results,with the target of eliminating the active power oscillations and current harmonics on a prototype GC-VSC system,validate the feasibility of the proposed current control scheme during supply voltage imbalance and distortion.
基金Project supported by the State Key Development Program for Basic Research of China (Grant No 2006CB601007)the National Natural Science Foundation of China (Grant No 10674006)the National High Technology Research and Development Program of China (Grant No 2007AA03Z238)
文摘Two cardiac functional models are constructed in this paper. One is a single current model and the other is a current multipole model. Parameters denoting the properties of these two models are calculated by a least-square fit to the measurements using a simulated annealing algorithm. The measured signals are detected at 36 observation nodes by a superconducting quantum interference device （SQUID）. By studying the trends of position, orientation and magnitude of the single current dipole model and the current multipole model in the QRS complex during one time span and comparing the reconstructed magnetocardiography （MCG） of these two cardiac models, we find that the current multipole model is a more appropriate model to represent cardiac electrophysiological activity.
基金Supported by the National Natural Science Foundation of China under Grant Nos.10371070,10671121the Foundation of Shanghai Education Committee for Shanghai Prospective Excellent Young Teachers+1 种基金Shanghai Leading Academic Discipline Project under Grant No.J50101 the President Foundation of East China Institute of Technology under Grant No.DHXK0810
文摘N-soliton solutions of the hierarchy of non-isospectral mKdV equation with self-consistent sources andthe hierarchy of non-isospectral sine-Gordon equation with self-consistent sources are obtained via the inverse scatteringtransform.
文摘In reactor physics tests, it is important to monitor sub-criticality continuously during criticality approach. Reactivity measurements by the inverse kinetics method are widely used during the operation of a nuclear reactor. This technique is successfully applied at sufficiently high power level or to a core without an external neutron source where the neutron source term in point reactor kinetics equations may be neglected. For operation at low power levels or in the sub-critical domain, the increase in the fluctuation of the neutron signal may cause difficulties in the evaluation of reactivity and the effect of direct emission from the external neutron source may not be neglected. Therefore, contribution of the neutron source must be taken into account and this implies knowledge of a quantity proportional to the source strength, which calls the source term and then it should be determined. The research work has been conducted to measure reactivity with source term using a dedicated reactivity measurement system by the Least Square Inverse Kinetics Method (LSIKM). Application to a simulator of HANARO research reactor, Korea Atomic Energy Research Institute (KAERI), with known source strength and reactivity worth has showed consistent and satisfactory agreement.
文摘The elastic plate vibration model is studied under the external force. The size of the source term by the given mode of the source and some observations from the body of the plate is determined over a time interval, which is referred to be an inverse source problem of a plate equation. The uniqueness theorem for this problem is stated, and the fundamental solution to the plate equation is derived. In the case that the plate is driven by the harmonic load, the fundamental solution method （FSM） and the Tikhonov regularization technique axe used to calculate the source term. Numerical experiments of the Euler-Bernoulli beam and the Kirchhoff-Love plate show that the FSM can work well for practical use, no matter the source term is smooth or piecewise.
基金supported by National Natural Science Foundation of China(12271277)the Open Research Fund of Key Laboratory of Nonlinear Analysis&Applications(Central China Normal University),Ministry of Education,China.
文摘In this article,we consider the diffusion equation with multi-term time-fractional derivatives.We first derive,by a subordination principle for the solution,that the solution is positive when the initial value is non-negative.As an application,we prove the uniqueness of solution to an inverse problem of determination of the temporally varying source term by integral type information in a subdomain.Finally,several numerical experiments are presented to show the accuracy and efficiency of the algorithm.
文摘The inverse estimation of a source location of pollutant released into a turbulent flow is a probability problem instead of a deterministic one, as the turbulent flow is chaotic and irreversible. However, researches can be conducted to provide helpful instructions to the possible source location with corresponding uncertainty. This study aims to propose a method of inverse estimation of a passive-scalar source location. Experimental investigation of the dye plume characteristics released into a fully-developed turbulent flow is performed in a water channel. A planar laser-induced fluorescence (PLIF) technique is used to obtain two-dimensional images of spreading dye plumes at a bulk Reynolds number of 20,000. The distributions of high concentration areas in the PLIF images are chosen as features that characterize the traveling (diffusion) distance or time from the dye source. Graphical analysis is used to extract these high concentration areas. The procedure of graphical analysis has three steps: 1) binarization using a threshold to extract high concentration dye patches;2) labeling individual high-concentration dye patches in the binarized images;and 3) pixel-counting to measure the area and perimeter of each dye patch. We examine the variations of fractal dimension of patches, and the fractal dimension is observed to be almost constant irrespective of the distance from the source. The kurtosis of the probability density function curve of the logarithm dimensionless dye patch areas is found to be related with the downstream diffusion distance, based on which an inverse estimation method to locate a passive-scalar point source is proposed and evaluated.
文摘The inverse scattering problems are to detect the property of obstacles from the measurements outside the obstacles. One of important research areas in this topic is the recovery of boundary property for impenetrable obstacles. In this paper, we would like to give a brief review about the recently developed singular source methods. There are three different methods in this category, namely, linear sampling method, pointsource method and probe method. We also present some recent new results about the probe method.
基金financially supported by the National Natural Science Foundation of China(No.41074075/D0409)the National Science and Technology Major Project(No.2011ZX05025-001-04)
文摘As a high quality seismic imaging method, full waveform inversion （FWI） can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing is computationally expensive, especially for the three-dimension complex medium inversion. Introducing blended source technology into the frequency-domain FWI can greatly reduce the computational burden and improve the efficiency of the inversion. However, this method has two issues： first, crosstalk noise is caused by interference between the sources involved in the encoding, resulting in an inversion result with some artifacts; second, it is more sensitive to ambient noise compared to conventional FWI, therefore noisy data results in a poor inversion. This paper introduces a frequency-group encoding method to suppress crosstalk noise, and presents a frequency- domain auto-adapting FWI based on source-encoding technology. The conventional FWI method and source-encoding based FWI method are combined using an auto-adapting mechanism. This improvement can both guarantee the quality of the inversion result and maximize the inversion efficiency.
基金supported by the National Natural Science Foundation of China(Nos.41564001,41674133,41572185,and 41604104)the Distinguished Young Talent Foundation of Jiangxi Province(No.20171BCB23068)
文摘Multiturn coils is an effective transmitter for transient electromagnetic method(TEM) used in narrow space and complex terrain at presently. However, its high mutual inductance coupling and long turn-off time affect the quality of later data processing and interpretation. Compared with multiturn coils, the new conical source has low mutual inductance and short turn-off time. Based on the superposition principle, we use Hankel transform and numerical filtering method for forward modelling of the conical source field in the layered-media and explore TEM characteristics excited by this source. We apply improved damped least square inversion to integrated transient electromagnetic(TEM) data. We first invert the induced voltage into similar resistivity and apparent depth, and then use the inverted results as input parameters in the initial model and transform the apparent resistivity data into the frequency domain. Then, damped least square inversion is performed in the frequency domain using the initial model. Subsequently, we use automated model building to search for the extremes and inflection points in the resistivity–depth data that are treated as critical layer parameters. The inversion of theoretical and observed data suggests that the method modifies the resistivity and depth and yields a model of the underground layers.
基金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 National Key Research and Development Program of China(2022YFF0800603).
文摘We conducted rapid inversions of rupture process for the 2023 earthquake doublet occurred in SE Türkiye,the first with a magnitude of M_(W)7.8 and the second with a magnitude of M_(W)7.6,using teleseismic and strong-motion data.The teleseismic rupture models of the both events were obtained approximately 88 and 55 minutes after their occurrences,respectively.The rupture models indicated that the first event was an asymmetric bilateral event with ruptures mainly propagating to the northeast,while the second one was a unilateral event with ruptures propagating to the west.This information could be useful in locating the meizoseismal areas.Compared with teleseismic models,the strong-motion models showed relatively higher resolution.A noticeable difference was found for the M_(W)7.6 earthquake,for which the strong-motion models shows a bilateral event,rather than a unilateral event,but the dominant rupture direction is still westward.Nevertheless,all strong-motion models are consistent with the teleseismic models in terms of magnitudes,durations,and dominant rupture directions.This suggests that both teleseismic and strong-motion data can be used for fast determination of major source characteristics.In contrast,the strong-motion data would be preferable in future emergency responses since they are recorded earlier and have a better resolution ability on the source ruptures.
基金the National Natural Science Foundation of China(Nos.42174011and 41874001).
文摘The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault parameters inversion has nonlinear characteristics, and the gradient-based optimizer(GBO) has the characteristics of fast convergence speed and falling into local optimum hardly. This paper applies GBO algorithm to simulated earthquakes and real LuShan earthquakes in the nonlinear inversion of the Okada model to obtain the source parameters. The simulated earthquake experiment results show that the algorithm is stable, and the seismic source parameters obtained by GBO are slightly closer to the true value than the multi peak particle swarm optimization(MPSO). In the 2013 LuShan earthquake experiment, the root mean square error between the deformation after forwarding of fault parameters obtained by the introduced GBO algorithm and the surface observation deformation was 3.703 mm, slightly better than 3.708 mm calculated by the MPSO. Moreover, the inversion result of GBO algorithm is better than MPSO algorithm in stability. The above results show that the introduced GBO algorithm has a certain practical application value in seismic fault source parameters inversion.
基金supported by the National Natural Science Foundation of China(51775121,12364055).
文摘Based on Fourier transform,the angular spectrum equation and its discrete model of the inverse source problem of the two-dimensional Helmholtz equation are derived.By transforming the two-dimensional inverse source problem in the spatial domain into the inverse source problem in the one-dimensional Fourier domain,the source intensity angular spectrum vector in the Fourier domain is calculated by the one-dimensional method of moments.Then,using the discrete Fourier inverse transform,a multi-frequency angular spectrum method of moments is proposed to reconstruct the inverse source problem of the two-dimensional Helmholtz equation.The reconstruction accuracy and stability of this method are improved by optimizing the excitation wavenumber sequence.Simulation calculations show that the multi-frequency angular spectrum method of moments has high reconstruction accuracy and anti-noise ability,and can also design the reconstruction model according to the effective frequency band of the acoustic sensor to achieve the highest resolution.In addition,since the calculation of the coefficient matrix of the proposed method is independent of the distribution of the source intensity function and does not require regularization processing,when the size and resolution of the required source intensity distribution area are determined,after the coefficient matrix is calculated and saved,the source intensity function corresponding to the region of interest can be reconstructed according to different sound pressure detection data,with high calculation accuracy and efficiency.
文摘A new source location method using wave-equation based traveltime inversion is proposed to locate microseismic events accurately. With a sourceindependent strategy, microseismic events can be located independently regardless of the accuracy of the source signature and the origin time. The traveltime-residuals-based misfit function has robust performance when the velocity model is inaccurate. The new Fréchet derivatives of the misfit function with respect to source location are derived directly based on the acoustic wave equation, accounting for the influence of geometrical perturbation and spatial velocity variation. Unlike the mostly used traveltime inversion methods, no traveltime picking or ray tracing is needed.Additionally, the improved scattering-integral method is applied to reduce the computational cost. Numerical tests show the validity of the proposed method.
基金The Specialized Funds for National Key Basic Study (G1998040704), the Dual Project of China Earthquake Admini-stration (9691309020301) and National Natural Science Foundation of China (46764010).
文摘plane with the strike of 127°, the dip of 79° and the rake of 171°. The rupture process inversion result of MS=7.9 earthquake shows that the total rupture duration is about 37 s, the scalar moment tensor is M0=0.97 × 1020 N·m. Rupture mainly occurred on the shallow area with 110 km long and 30 km wide, the location in which the rupture initiated is not where the main rupture took place, and the area with slip greater than 0.5 m basically lies within 35 km deep middle-crust under the earth surface. The maximum static slip is 3.6 m. There are two distinct areas with slip larger than 2.0 m. We noticed that when the rupture propagated towards northwest and closed to the area around the MS=7.3 hypocenter, the slip decreased rapidly, which may indicate that the rupture process was stopped by barriers. The consistence of spatial distribution of slip on the fault plane with the distribution of aftershocks also supports that the rupture is a heterogeneous process owing to the presence of barriers.
