High-resolution bone microstructure imaging based on ultrasonic frequency-domain full-waveform inversion

The main challenge in bone ultrasound imaging is the large acoustic impedance contrast and sound velocity differences between the bone and surrounding soft tissue. It is difficult for conventional pulse-echo modalities to give accurate ultrasound images for irregular bone boundaries and microstructures using uniform sound velocity assumption rather than getting a prior knowledge of sound speed. To overcome these limitations, this paper proposed a frequency-domain fullwaveform inversion(FDFWI) algorithm for bone quantitative imaging utilizing ultrasonic computed tomography(USCT).The forward model was calculated in the frequency domain by solving the full-wave equation. The inverse problem was solved iteratively from low to high discrete frequency components via minimizing a cost function between the modeled and measured data. A quasi-Newton method called the limited-memory Broyden–Fletcher–Goldfarb–Shanno algorithm(L-BFGS) was utilized in the optimization process. Then, bone images were obtained based on the estimation of the velocity and density. The performance of the proposed method was verified by numerical examples, from tubular bone phantom to single distal fibula model, and finally with a distal tibia-fibula pair model. Compared with the high-resolution peripheral quantitative computed tomography(HR-p QCT), the proposed FDFWI can also clearly and accurately presented the wavelength scaled pores and trabeculae in bone images. The results proved that the FDFWI is capable of reconstructing high-resolution ultrasound bone images with sub-millimeter resolution. The parametric bone images may have the potential for the diagnosis of bone disease.

3D variable-grid full-waveform inversion on GPU

Full-waveform inversion(FWI)is a powerful tool to reconstruct subsurface geophysical parameters with high resolution.As3 D surveys become widely implemented,corresponding 3 D processing techniques are required to solve complex geological cases,while a large amount of computation is the most challenging problem.We propose an adaptive variable-grid 3 D FWI on graphics processing unit devices to improve computational efficiency without losing accuracy.The irregular-grid discretization strategy is based on a dispersion relation,and the grid size adapts to depth,velocity,and frequency automatically.According to the transformed grid coordinates,we derive a modified acoustic wave equation and apply it to full wavefield simulation.The 3 D variable-grid modeling is conducted on several 3 D models to validate its feasibility,accuracy and efficiency.Then we apply the proposed modeling method to full-waveform inversion for source and residual wavefield propagation.It is demonstrated that the adaptive variable-grid FWI is capable of decreasing computing time and memory requirements.From the inversion results of the 3 D SEG/EAGE overthrust model,our method retains inversion accuracy when recovering both thrust and channels.

Frequency-Domain Full-Waveform Inversion Based on Tunnel Space Seismic Data

Tunnel seismic detection methods are effective for obtaining the geological structure around the tunnel face,which is critical for safe construction and disaster mitigation in tunnel engineering.However,there is often a lack of accuracy in the acquired geological information and physical properties ahead of the tunnel face in the current tunnel seismic detection methods.Thus,we apply a frequency-domain acoustic full-waveform inversion(FWI)method to obtain high-resolution results for the tunnel structure.We discuss the influence of the frequency group selection strategy and the tunnel observation system settings regarding the inversion results and determine the structural imaging and physical property parameter inversion of abnormal geological bodies ahead of the tunnel face.Based on the conventional strategies of frequency-domain acoustic FWI,we propose a frequency group selection strategy that combines a low-frequency selection covering the vertical wavenumber and a high-frequency selection of antialiasing.This strategy can effectively obtain the spatial structure and physical parameters of the geology ahead of the tunnel face and improve the inversion resolution.In addition,by linearly increasing the side length of the tunnel observation system,we share the influence of the length of the two sides of the observation systems of different tunnels on the inversion results.We found out that the inversion results are the best when the side length is approximately five times the width of the tunnel face,and the influence of increasing the side observation length beyond this range on the inversion results can be ignored.Finally,based on this approach,we invert for the complex multi-stratum model,and an accurate structure and physical property parameters of the complex stratum ahead of the tunnel face are obtained,which verifies the feasibility of the proposed method.

Q full-waveform inversion based on the viscoacoustic equation

Presently, most full-waveform inversion methods are developed for elastic media and ignore the effect of attenuation. The calculation of the quality factor Q is based on velocity parameter inversion under the assumption of a given Q-model that is obtained by tomographic inversion. However, the resolution of the latter is low and cannot reflect the amplitude attenuation and phase distortion during wave propagation in viscoelastic media. Thus, a Q waveform inversion method is proposed. First, we use standard linear body theory to describe attenuation and then we derive the simplified viscoacoustic equation that characterizes amplitude attenuation and phase distortion. In comparison with conventional equations, the simplifi ed equation involves no memory variables and therefore requires less memory during computation. Moreover, the implementations of the attenuation compensation are easier. The adjoint equation and the corresponding gradient equation with respect to either L2-norm or the zero-lag cross-correlation objective function are then derived and the regularization strategy for overcoming the instability during numerical solution of the adjoint equation is proposed. The Q waveform inversion is developed using the limited-memory Broyden–Fletcher– Goldfarb–Shanno (L-BFGS) iteration method for known velocity. To alleviate the dependence of the waveform inversion on the initial model and overcome cycle skipping to some extent, we adopt multiscale analysis. Furthermore, anti-noise property and double-parameter inversion are assessed based on the results of numerical modeling.

ultiscale full-waveform inversion based on shot subsampling

常规全波形反演利用全部炮集参与计算,反演的计算量巨大。针对这一问题,本文分析了不同频率反演对炮数的需求,进而提出一种基于频率多尺度反演方法的加速策略。该方法利用反演所需炮数与频率正相关的特性,在反演低频数据时,每次迭代只抽取一部分炮集参与反演,频率升高时,相应地引入更多的炮集参与运算,两次迭代之间通过组内随机炮采样的方法实现炮集的轮换,避免炮集信息的丢失。该方法通过降低反演炮数从而减少计算量,由于不涉及炮集的串扰,因此不会引入额外的噪声,也不受限于观测系统。模型测试结果表明,该方法在炮集数量较多时可以明显减少计算时间,同时,该方法具有一定的抗噪能力,对含噪声的地震记录也能得到较好的反演结果。

Layer-Valley Hall Effect under Inversion and Time-Reversal Symmetries

Hall effects have been the central paradigms in modern physics,materials science and practical applications,and have led to many exciting breakthroughs,including the discovery of topological Chern invariants and the revolution of metrological resistance standard.To date,the Hall effects have mainly focused on a single degree of freedom(Do F),and most of them require the breaking of spatial-inversion and/or time-reversal symmetries.Here we demonstrate a new type of Hall effect,i.e.,layer-valley Hall effect,based on a combined layer-valley Do F characterized by the product of layer and valley indices.The layer-valley Hall effect has a quantum origin arising from the layer-valley contrasting Berry curvature,and can occur in nonmagnetic centrosymmetric crystals with both spatial-inversion and time-reversal symmetries,transcending the symmetry constraints of single Do F Hall effect based on the constituent layer or valley index.Moreover,the layer-valley Hall effect is highly tunable and shows a W-shaped pattern in response to the out-of-plane electric fields.Additionally,we discuss the potential detection approaches and material-specific design principles of layer-valley Hall effect.Our results demonstrate novel Hall physics and open up exotic paradigms for new research direction of layer-valleytronics that exploits the quantum nature of the coupled layer-valley DoF.

Linearized waveform inversion for vertical transversely isotropic elastic media:Methodology and multi-parameter crosstalk analysis

Seismic migration and inversion are closely related techniques to portray subsurface images and identify hydrocarbon reservoirs.Seismic migration aims at obtaining structural images of subsurface geologic discontinuities.More specifically,seismic migration estimates the reflectivity function(stacked average reflectivity or pre-stack angle-dependent reflectivity)from seismic reflection data.On the other hand,seismic inversion quantitatively estimates the intrinsic rock properties of subsurface formulations.Such seismic inversion methods are applicable to detect hydrocarbon reservoirs that may exhibit lateral variations in the inverted parameters.Although there exist many differences,pre-stack seismic migration is similar with the first iteration of the general linearized seismic inversion.Usually,seismic migration and inversion techniques assume an acoustic or isotropic elastic medium.Unconventional reservoirs such as shale and tight sand formation have notable anisotropic property.We present a linearized waveform inversion(LWI)scheme for weakly anisotropic elastic media with vertical transversely isotropic(VTI)symmetry.It is based on two-way anisotropic elastic wave equation and simultaneously inverts for the localized perturbations(ΔVp_(0)/Vp_(0)/Vs_(0)/Vs_(0)/,Δ∈,Δδ)from the long-wavelength reference model.Our proposed VTI-elastic LWI is an iterative method that requires a forward and an adjoint operator acting on vectors in each iteration.We derive the forward Born approximation operator by perturbation theory and adjoint operator via adjoint-state method.The inversion has improved the quality of the images and reduces the multi-parameter crosstalk comparing with the adjoint-based images.We have observed that the multi-parameter crosstalk problem is more prominent in the inversion images for Thomsen anisotropy parameters.Especially,the Thomsen parameter is the most difficult to resolve.We also analyze the multi-parameter crosstalk using scattering radiation patterns.The linearized waveform inversion for VTI-elastic media presented in this article provides quantitative information of the rock properties that has the potential to help identify hydrocarbon reservoirs.

True-temperature inversion algorithm for a multi-wavelength pyrometer based on fractional-order particle-swarm optimization

Herein,a method of true-temperature inversion for a multi-wavelength pyrometer based on fractional-order particle-swarm optimization is proposed for difficult inversion problems with unknown emissivity.Fractional-order calculus has the inherent advantage of easily jumping out of local extreme values;here,it is introduced into the particle-swarm algorithm to invert the true temperature.An improved adaptive-adjustment mechanism is applied to automatically adjust the current velocity order of the particles and update their velocity and position values,increasing the accuracy of the true temperature values.The results of simulations using the proposed algorithm were compared with three algorithms using typical emissivity models:the internal penalty function algorithm,the optimization function(fmincon)algorithm,and the conventional particle-swarm optimization algorithm.The results show that the proposed algorithm has good accuracy for true-temperature inversion.Actual experimental results from a rocket-motor plume were used to demonstrate that the true-temperature inversion results of this algorithm are in good agreement with the theoretical true-temperature values.

Surface wave inversion with unknown number of soil layers based on a hybrid learning procedure of deep learning and genetic algorithm

Surface wave inversion is a key step in the application of surface waves to soil velocity profiling.Currently,a common practice for the process of inversion is that the number of soil layers is assumed to be known before using heuristic search algorithms to compute the shear wave velocity profile or the number of soil layers is considered as an optimization variable.However,an improper selection of the number of layers may lead to an incorrect shear wave velocity profile.In this study,a deep learning and genetic algorithm hybrid learning procedure is proposed to perform the surface wave inversion without the need to assume the number of soil layers.First,a deep neural network is adapted to learn from a large number of synthetic dispersion curves for inferring the layer number.Then,the shear-wave velocity profile is determined by a genetic algorithm with the known layer number.By applying this procedure to both simulated and real-world cases,the results indicate that the proposed method is reliable and efficient for surface wave inversion.

A novel efficient energy absorber with free inversion of a metal foam-filled circular tube

In this paper, a novel efficient energy absorber with free inversion of a metal foam-filled circular tube(MFFCT) is designed, and the axial compressive behavior of the MFFCT under free inversion is studied analytically and numerically. The theoretical analysis reveals that the energy is mainly dissipated through the radial bending of the metal circular tube, the circumferential expansion of the metal circular tube, and the metal filled-foam compression. The principle of energy conservation is used to derive the theoretical formula for the minimum compressive force of the MFFCT over free inversion under axial loading. Furthermore, the free inversion deformation characteristics of the MFFCT are analyzed numerically. The theoretical steady values are found to be in good agreement with the results of the finite element(FE) analysis. The effects of the average diameter of the metal tube, the wall thickness of the metal tube, and the filled-foam strength on the free inversion deformation of the MFFCT are considered. It is observed that in the steady deformation stage, the load-carrying and energy-absorbing capacities of the MFFCT increase with the increase in the average diameter of the metal tube, the wall thickness of the metal tube, or the filled-foam strength. The specific energy absorption(SEA) of free inversion of the MFFCT is significantly higher than that of the metal tube alone.

Multi-scale data joint inversion of minerals and porosity in altered igneous reservoirs—A case study in the South China Sea

There are abundant igneous gas reservoirs in the South China Sea with significant value of research,and lithology classification,mineral analysis and porosity inversion are important links in reservoir evaluation.However,affected by the diverse lithology,complicated mineral and widespread alteration,conventional logging lithology classification and mineral inversion become considerably difficult.At the same time,owing to the limitation of the wireline log response equation,the quantity and accuracy of minerals can hardly meet the exploration requirements of igneous formations.To overcome those issues,this study takes the South China Sea as an example,and combines multi-scale data such as micro rock slices,petrophysical experiments,wireline log and element cutting log to establish a set of joint inversion methods for minerals and porosity of altered igneous rocks.Specifically,we define the lithology and mineral characteristics through core slices and mineral data,and establish an igneous multi-mineral volumetric model.Then we determine element cutting log correction method based on core element data,and combine wireline log and corrected element cutting log to perform the lithology classification and joint inversion of minerals and porosity.However,it is always difficult to determine the elemental eigenvalues of different minerals in inversion.This paper uses multiple linear regression methods to solve this problem.Finally,an integrated inversion technique for altered igneous formations was developed.The results show that the corrected element cutting log are in good agreement with the core element data,and the mineral and porosity results obtained from the joint inversion based on the wireline log and corrected element cutting log are also in good agreement with the core data from X-ray diffraction.The results demonstrate that the inversion technique is applicable and this study provides a new direction for the mineral inversion research of altered igneous formations.

Three-dimensional constrained gravity inversion of Moho depth and crustal structural characteristics at Mozambique continental margin

Mozambique's continental margin in East Africa was formed during the break-off stage of the east and west Gondwana lands. Studying the geological structure and division of continent-ocean boundary(COB) in Mozambique's continental margin is considered of great significance to rebuild Gondwana land and understand its movement mode. Along these lines, in this work, the initial Moho was fit using the known Moho depth from reflection seismic profiles, and a 3D multi-point constrained gravity inversion was carried out. Thus, highaccuracy Moho depth and crustal thickness in the study area were acquired. According to the crustal structure distribution based on the inversion results, the continental crust at the narrowest position of the Mozambique Channel was detected. According to the analysis of the crustal thickness, the Mozambique ridge is generally oceanic crust and the COB of the whole Mozambique continental margin is divided.

Crustal and uppermost mantle structure of the northeastern Qinghai-Xizang Plateau from joint inversion of surface wave dispersions and receiver functions with P velocity constraints

Lithospheric structure beneath the northeastern Qinghai-Xizang Plateau is of vital significance for studying the geodynamic processes of crustal thickening and expansion of the Qinghai-Xizang Plateau. We conducted a joint inversion of receiver functions and surface wave dispersions with P-wave velocity constraints using data from the Chin Array Ⅱ temporary stations deployed across the Qinghai-Xizang Plateau. Prior to joint inversion, we applied the H-κ-c method(Li JT et al., 2019) to the receiver function data in order to correct for the back-azimuthal variations in the arrival times of Ps phases and crustal multiples caused by crustal anisotropy and dipping interfaces. High-resolution images of vS, crustal thickness, and vP/vSstructures in the Qinghai-Xizang Plateau were simultaneously derived from the joint inversion. The seismic images reveal that crustal thickness decreases outward from the Qinghai-Xizang Plateau. The stable interiors of the Ordos and Alxa blocks exhibited higher velocities and lower crustal vP/vSratios. While, lower velocities and higher vP/vSratios were observed beneath the Qilian Orogen and Songpan-Ganzi terrane(SPGZ), which are geologically active and mechanically weak, especially in the mid-lower crust.Delamination or thermal erosion of the lithosphere triggered by hot asthenospheric flow contributes to the observed uppermost mantle low-velocity zones(LVZs) in the SPGZ. The crustal thickness, vS, and vP/vSratios suggest that whole lithospheric shortening is a plausible mechanism for crustal thickening in the Qinghai-Xizang Plateau, supporting the idea of coupled lithospheric-scale deformation in this region.

Statistical Inversion Based on Nonlinear Weighted Anisotropic Total Variational Model and Its Application in Electrical Impedance Tomography

Electrical impedance tomography (EIT) aims to reconstruct the conductivity distribution using the boundary measured voltage potential. Traditional regularization based method would suffer from error propagation due to the iteration process. The statistical inverse problem method uses statistical inference to estimate unknown parameters. In this article, we develop a nonlinear weighted anisotropic total variation (NWATV) prior density function based on the recently proposed NWATV regularization method. We calculate the corresponding posterior density function, i.e., the solution of the EIT inverse problem in the statistical sense, via a modified Markov chain Monte Carlo (MCMC) sampling. We do numerical experiment to validate the proposed approach.

Quantitative Multi-Layer Electromagnetic Induction Inversion and Full-Waveform Inversion of Crosshole Ground Penetrating Radar Data

Due to the recent system developments for the electromagnetic characterization of the subsurface, fast and easy acquisition is made feasible due to the fast measurement speed, easy coupling with GPS systems, and the availability of multi-channel electromagnetic induction（EMI） and ground penetrating radar（GPR） systems. Moreover, the increasing computer power enables the use of accurate forward modeling programs in advanced inversion algorithms where no approximations are used and the full information content of the measured data can be exploited. Here, recent developments of large-scale quantitative EMI inversion and full-waveform GPR inversion are discussed that yield higher resolution of quantitative medium properties compared to conventional approaches. In both cases a detailed forward model is used in the inversion procedure that is based on Maxwell＇s equations. The multi-channel EMI data that have different sensing depths for the different source-receiver offset are calibrated using a short electrical resistivity tomography（ERT） calibration line which makes it possible to invert for electrical conductivity changes with depth over large areas. The crosshole GPR full-waveform inversion yields significant higher resolution of the permittivity and conductivity images compared to ray-based inversion results.

Source-Independent Amplitude-Semblance Full-Waveform Inversion using a Hybrid Time-and Frequency-Domain Approach

Full-waveform inversion is a promising tool to produce accurate and high-resolution subsurface models.Conventional full-waveform inversion requires an accu-rate estimation of the source wavelet,and its computational cost is high.We develop a novel source-independent full-waveform inversion method using a hybrid time-and frequency-domain scheme to avoid the requirement of source wavelet estimation and to reduce the computational cost.We employ an amplitude-semblance objective function to not only effectively remove the source wavelet effect on full-waveform inver-sion,but also to eliminate the impact of the inconsistency of source wavelets among different shot gathers on full-waveform inversion.To reduce the high computational cost of full-waveform inversion in the time domain,we implement our new algorithm using a hybrid time-and frequency-domain approach.The forward and backward wave propagation operations are conducted in the time domain,while the frequency-domain wavefields are obtained during modeling using the discrete-time Fourier trans-form.The inversion process is conducted in the frequency domain for selected frequen-cies.We verify our method using synthetic seismic data for the Marmousi model.The results demonstrate that our novel source-independent full-waveform inversion pro-duces accurate velocity models even if the source signature is incorrect.In addition,our method can significantly reduce the computational time using the hybrid time-and frequency-domain approach compared to the conventional full-waveform inversion in the time domain.

Rapid source inversions of the 2023 SE Türkiye earthquakes with teleseismic and strong-motion data

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.

GBO algorithm for seismic source parameters inversion

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.

Full-waveform Velocity Inversion Based on the Acoustic Wave Equation

Full-waveform velocity inversion based on the acoustic wave equation in the time domain is investigated in this paper. The inversion is the iterative minimization of the misfit between observed data and synthetic data obtained by a numerical solution of the wave equation. Two inversion algorithms in combination with the CG method and the BFGS method are described respectively. Numerical computations for two models including the benchmark Marmousi model with complex structure are implemented. The inversion results show that the BFGS-based algorithm behaves better in inversion than the CG-based algorithm does. Moreover, the good inversion result for Marmousi model with the BFGS-based algorithm suggests the quasi-Newton methods can provide an important tool for large-scale velocity inversion. More computations demonstrate the correctness and effectives of our inversion algorithms and code.

Moment tensor and stress inversion solutions of acoustic emissions during compression and tensile fracturing in crystalline rocks

We investigate the accuracy and robustness of moment tensor(MT)and stress inversion solutions derived from acoustic emissions(AEs)during the laboratory fracturing of prismatic Barre granite specimens.Pre-cut flaws in the specimens introduce a complex stress field,resulting in a spatial and temporal variation of focal mechanisms.Specifically,we consider two experimental setups:(1)where the rock is loaded in compression to generate primarily shear-type fractures and(2)where the material is loaded in indirect tension to generate predominantly tensile-type fractures.In each test,we first decompose AE moment tensors into double-couple(DC)and non-DC terms and then derive unambiguous normal and slip vectors using k-means clustering and an unstructured damped stress inversion algorithm.We explore temporal and spatial distributions of DC and non-DC events at different loading levels.The majority of the DC and the tensile non-DC events cluster around the pre-cut flaws,where macro-cracks later develop.Results of stress inversion are verified against the stress field from finite element(FE)modeling.A good agreement is found between the experimentally derived and numerically simulated stress orientations.To the best of the authors’knowledge,this work presents the first case where stress inversion methodologies are validated by numerical simulations at laboratory scale and under highly heterogeneous stress distributions.