We applied double-difference tomography to relocate seismic events and determine the lithospheric velocity structure beneath the New Britain Island arc and the South Bismarck Sea Basin,based on the local P wave arriva...We applied double-difference tomography to relocate seismic events and determine the lithospheric velocity structure beneath the New Britain Island arc and the South Bismarck Sea Basin,based on the local P wave arrival time dataset collected by the International Seismological Centre.Results of the seismic relocation and velocity inversion show that the subduction of Solomon Sea Plate along the New Britain Trench is spatially different above 150 km,and the subduction angle of the slab on the west side is higher than that on the east side.The relocated earthquakes also show that there are double seismic zones at the depths of about 30–90km beneath the New Britain Island Arc.The velocity structure shows that the dehydration of the subducting slab caused the low-velocity anomalies in mantle wedge above the slab,which are associated with the magmatic activities around the New Guinea-New Britain Island arc.Moreover,it shows that there is another low-velocity anomaly zone beneath the Bismarck mid-oceanic ridge with spatial variation.Beneath the west of the Bismarck mid-oceanic ridge,the low-velocity anomaly is weakly connected to the subducted Solomon Sea slab.Conversely,the low-velocity anomaly beneath the Manus Sea Basin is highly intertwined to the subducting slab and its mantle wedge,indicating that the subduction of the Solomon Sea Plate might be a key deep dynamic factor that drives the spreading of the Manus Sea Basin and the separation of the Bismarck Plate.展开更多
West Java in the western part of the Sunda Arc has a relatively high seismicity due to subduction activity and faults.In this study,double-difference tomography was used to obtain the 3D velocity tomograms of P and S ...West Java in the western part of the Sunda Arc has a relatively high seismicity due to subduction activity and faults.In this study,double-difference tomography was used to obtain the 3D velocity tomograms of P and S waves beneath the western part of Java.To infer the geometry of the structure beneath the study area,precise earthquake hypo・center determination was first performed before tomographic imaging.For this,earthquake waveform data were extracted from the regional Meteorological,Climatological,Geophysical Agency(BMKG)network of Indonesia from South Sumatra to Central Java.The P and S arrival times for about 1,000 events in the period April 2009 to July 2016 were selected,the key features being events of magnitude>3,azimuthal gap<210°and number of phases>8.A nonlinear method using the oct-tree sampling algorithm from the NonLinLoc program was employed to determine the earthquake hypocenters.The hypocenter locations were then relocated using double-difference tomography(tomoDD).A significant reduction of travel-time(root mean square basis)and a better clustering of earthquakes were achieved which correlated well with the geological structure in West Java.Double-difference tomography was found to give a clear velocity structure,especially beneath the volcanic arc area,i.e.,under Mt Anak Krakatau,Mt Salak and the mountains complex in the southern part of West Java.Low velocity anomalies for the P and S waves as well as the vp/vs ratio below the volcanoes indicated possible partial melting of the upper mantle which ascended from the subducted slab beneath the volcanic arc.展开更多
Electrical capacitance tomography(ECT)has been applied to two-phase flow measurement in recent years.Image reconstruction algorithms play an important role in the successful applications of ECT.To solve the ill-posed ...Electrical capacitance tomography(ECT)has been applied to two-phase flow measurement in recent years.Image reconstruction algorithms play an important role in the successful applications of ECT.To solve the ill-posed and nonlinear inverse problem of ECT image reconstruction,a new ECT image reconstruction method based on fast linearized alternating direction method of multipliers(FLADMM)is proposed in this paper.On the basis of theoretical analysis of compressed sensing(CS),the data acquisition of ECT is regarded as a linear measurement process of permittivity distribution signal of pipe section.A new measurement matrix is designed and L1 regularization method is used to convert ECT inverse problem to a convex relaxation problem which contains prior knowledge.A new fast alternating direction method of multipliers which contained linearized idea is employed to minimize the objective function.Simulation data and experimental results indicate that compared with other methods,the quality and speed of reconstructed images are markedly improved.Also,the dynamic experimental results indicate that the proposed algorithm can ful fill the real-time requirement of ECT systems in the application.展开更多
A trust region method is proposed to solve the problem of microwave tomography,which is very difficult to be solved for its ill-posedness and nonlinearity. Compared with the Levenberg-Marquardt method, this method int...A trust region method is proposed to solve the problem of microwave tomography,which is very difficult to be solved for its ill-posedness and nonlinearity. Compared with the Levenberg-Marquardt method, this method introduces more a priori knowledge and might obtain better results, though the two methods are equal in some cases.展开更多
Crosswell seismic tomography can be used to study the lateral variation of reservoirs, reservoir properties and the dynamic movement of fluids. In view of the instability of crosswell seismic tomography, the gradient ...Crosswell seismic tomography can be used to study the lateral variation of reservoirs, reservoir properties and the dynamic movement of fluids. In view of the instability of crosswell seismic tomography, the gradient method was improved by introducing regularization, and a gradient regularization method is presented in this paper. This method was verified by processing numerical simulation data and physical model data.展开更多
Solution to impedance distribution in electrical impedance tomography （EIT） is an ill-posed nonlinear inverse problem. It is especially difficult to reconstruct an EIT image in the center area of a measured object. ...Solution to impedance distribution in electrical impedance tomography （EIT） is an ill-posed nonlinear inverse problem. It is especially difficult to reconstruct an EIT image in the center area of a measured object. Tikhonov regularization with some prior information is a sound regnlarization method for static electrical impedance tomography under the condition that some true impedance distribution information is known a priori. This paper presents a direct search method （DSM） as pretreatment of image reconstruction through which one not only can construct a regularization matrix which may locate in areas of impedance change, but also can obtain an initial impedance distribution more similar to the true impedance distribution, as well as better current modes which can better distinguish the initial distribution and the true distribution. Simulation results indicate that, by using DSM, resolution in the center area of the measured object can be improved significantly.展开更多
This paper presents a gradient-descent travel time tomography method for solving the acoustictype velocity model inversion problem.Similarly to the adjoint-state method,the proposed method is based on the Eikonal equa...This paper presents a gradient-descent travel time tomography method for solving the acoustictype velocity model inversion problem.Similarly to the adjoint-state method,the proposed method is based on the Eikonal equation,enabling simultaneous calculation of contributions from all common-source receivers to the gradient.This overcomes the inefficiency inherent in conventional travel time tomography methods,which rely on a two-point ray tracing process.By directly calculating Fréchet derivatives,our method avoids the complex derivation processes associated with the adjoint-state method.The key to calculating the Fréchet derivatives is to calculate a so-called ray-path term.Consequently,compared to the adjoint-state method,the proposed method can explicitly obtain the ray paths,resulting in a more concise and intuitive derivation process.Furthermore,our method retains the benefits of the adjoint-state method,such as speed,low memory usage,and robustness.This paper focuses on elucidating the principles and algorithms for calculating the raypath term based on the fast sweeping method.The algorithms could be further speeded up by using parallel computational techniques.Synthetic tests demonstrate that our proposed travel time tomographic method accurately calculates ray paths,regardless of the complexity of the model and recording geometry.展开更多
In this paper, a novel reconstruction method is presented for Near Infrared (NIR) 2-D imaging to recover optical absorption coefficients from laboratory phantom data. The main body of this work validates a new generat...In this paper, a novel reconstruction method is presented for Near Infrared (NIR) 2-D imaging to recover optical absorption coefficients from laboratory phantom data. The main body of this work validates a new generation of highly efficient reconstruction algorithms called “Globally Convergent Method” (GCM) based upon actual measurements taken from brain-shape phantoms. It has been demonstrated in earlier studies using computer-simulated data that this type of reconstructions is stable for imaging complex distributions of optical absorption. The results in this paper demonstrate the excellent capability of GCM in working with experimental data measured from optical phantoms mimicking a rat brain with stroke.展开更多
To reduce the shielding effect of hardened layers on electrical resistivity tomography,a ratio method based on the distortion correction principle and the isolation coefficient is proposed.The effects of the resistivi...To reduce the shielding effect of hardened layers on electrical resistivity tomography,a ratio method based on the distortion correction principle and the isolation coefficient is proposed.The effects of the resistivity and thickness of hardened concrete layers on the detection of target objects are explored.Both numerical simulations and indoor tank tests indicate that when the ratio method is employed to correct the original collected data,the maximum allowable error for the isolation coefficient should not exceed 1%.Notably,when the ratio of hardened layer thickness to electrode spacing does not exceed 1,correction through this method significantly enhances the recognition capability of target objects.However,when the hardened layer thickness is greater than the electrode spacing by a factor of 2 or more,the ratio method cannot achieve satisfactory results.The case study of flood control engineering detection in the Zhangxi section of the Huangpen River in Dongzhi County demonstrates that the detection effect after correction by the ratio method is comparable to that for the adjacent unhardened pavement,and the influence of the hardened layer is obviously weakened,resulting in more reliable results.展开更多
This paper proposes a new level-set-based shape recovery approach that can be applied to a wide range of binary tomography reconstructions.In this technique,we derive generic evolution equations for shape reconstructi...This paper proposes a new level-set-based shape recovery approach that can be applied to a wide range of binary tomography reconstructions.In this technique,we derive generic evolution equations for shape reconstruction in terms of the underlying level-set parameters.We show that using the appropriate basis function to parameterize the level-set function results in an optimization problem with a small number of parameters,which overcomes many of the problems associated with the traditional level-set approach.More concretely,in this paper,we use Gaussian functions as a basis function placed at sparse grid points to represent the parametric level-set function and provide more flexibility in the binary representation of the reconstructed image.In addition,we suggest a convex optimization method that can overcome the problem of the local minimum of the cost function by successfully recovering the coefficients of the basis function.Finally,we illustrate the performance of the proposed method using synthetic images and real X-ray CT projection data.We show that the proposed reconstruction method compares favorably to various state-of-the-art reconstruction techniques for limited-data tomography,and it is also relatively stable in the presence of modest amounts of noise.Furthermore,the shape representation using a compact Gaussian radial basis function works well.展开更多
A P-wave velocity model was built in the central southern of the Tanlu Fault based on double-difference tomography.The results suggest the presence of a low-velocity anomaly extending from the surface to a depth of 25...A P-wave velocity model was built in the central southern of the Tanlu Fault based on double-difference tomography.The results suggest the presence of a low-velocity anomaly extending from the surface to a depth of 25 km around the Tanlu and Feixi Faults,representing fault-related fluids caused by partial melting.The relocated earthquakes indicate a significant concentration of seismic activity above 20 km around the Tanlu and Feixi Faults,suggesting that prominent fault systems possibly serve as conduits for the upward migration of deep minerals.The proposed geodynamic model,supported by geological and geophysical data,suggests that the migration of deep mineralized materials extends along the Tanlu Fault.The obtained results serve as a crucial foundation for elucidating the intricate process of mineralization in the central southern segment of the Tanlu Fault,thereby enhancing comprehension regarding the interaction among ore body formation,fault fluids,localized melting,and seismic activity.展开更多
The Longmenshan(LMS)fault zone is located at the junction of the eastern Tibetan Plateau and the Sichuan Basin and is of great significance for studying regional tectonics and earthquake hazards.Although regional velo...The Longmenshan(LMS)fault zone is located at the junction of the eastern Tibetan Plateau and the Sichuan Basin and is of great significance for studying regional tectonics and earthquake hazards.Although regional velocity models are available for the LMS fault zone,high-resolution velocity models are lacking.Therefore,a dense array of 240 short-period seismometers was deployed around the central segment of the LMS fault zone for approximately 30 days to monitor earthquakes and characterize fine structures of the fault zone.Considering the large quantity of observed seismic data,the data processing workflow consisted of deep learning-based automatic earthquake detection,phase arrival picking,and association.Compared with the earthquake catalog released by the China Earthquake Administration,many more earthquakes were detected by the dense array.Double-difference seismic tomography was adopted to determine V_(p),V_(s),and V_(p)/V_(s)models as well as earthquake locations.The checkerboard test showed that the velocity models have spatial resolutions of approximately 5 km in the horizontal directions and 2 km at depth.To the west of the Yingxiu–Beichuan Fault(YBF),the Precambrian Pengguan complex,where most of earthquakes occurred,is characterized by high velocity and low V_(p)/V_(s)values.In comparison,to the east of the YBF,the Upper Paleozoic to Jurassic sediments,where few earthquakes occurred,show low velocity and high V_(p)/V_(s)values.Our results suggest that the earthquake activity in the LMS fault zone is controlled by the strength of the rock compositions.When the high-resolution velocity models were combined with the relocated earthquakes,we were also able to delineate the fault geometry for different faults in the LMS fault zone.展开更多
Accurate prediction of coal and gas outburst(CGO)hazards is paramount in gas disaster prevention and control.This paper endeavors to overcome the constraints posed by traditional prediction indexes when dealing with C...Accurate prediction of coal and gas outburst(CGO)hazards is paramount in gas disaster prevention and control.This paper endeavors to overcome the constraints posed by traditional prediction indexes when dealing with CGO incidents under low gas pressure conditions.In pursuit of this objective,we have studied and established a mechanical model of the working face under abnormal stress and the excitation energy conditions of CGO,and proposed a method for predicting the risk of CGO under abnormal stress.On site application verification shows that when a strong outburst hazard level prediction is issued,there is a high possibility of outburst disasters occurring.In one of the three locations where we predicted strong outburst hazards,a small outburst occurred,and the accuracy of the prediction was higher than the traditional drilling cuttings index S and drilling cuttings gas desorption index q.Finally,we discuss the mechanism of CGO under the action of stress anomalies.Based on the analysis of stress distribution changes and energy accumulation characteristics of coal under abnormal stress,this article believes that the increase in outburst risk caused by high stress abnormal gradient is mainly due to two reasons:(1)The high stress abnormal gradient leads to an increase in the plastic zone of the coal seam.After the working face advances,it indirectly leads to an increase in the gas expansion energy that can be released from the coal seam before reaching a new stress equilibrium.(2)Abnormal stress leads to increased peak stress of coal body in front of working face.When coal body in elastic area transforms to plastic area,its failure speed is accelerated,which induces accelerated gas desorption and aggravates the risk of outburst.展开更多
Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. ...Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm （SAGA） and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A＇nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.展开更多
To solve the problem that metal artifacts severely damage the clarity of the organization structure in computed tomography（CT） images, a sinogram fusion-based metal artifact correction method is proposed. First, the...To solve the problem that metal artifacts severely damage the clarity of the organization structure in computed tomography（CT） images, a sinogram fusion-based metal artifact correction method is proposed. First, the metal image is segmented from the original CT image by the pre-set threshold. The original CT image and metal image are forward projected into the original projection sinogram and metal projection sinogram, respectively. The interpolation-based correction method and mean filter are used to correct the original CT image and preserve the edge of the corrected CT image, respectively. The filtered CT image is forward projected into the filtered image sinogram. According to the position of the metal sinogram in the original sinogram and filtered image sinogram, the corresponding sinograms PM^D （ in the original sinogram） and PM^C （ in the filtered image sinogram）can be acquired from the original sinogram and filtered image sinogram, respectively. Then, PM^D and PM^C are fused into the fused metal sinogram PM^F according to a certain proportion.The final sinogram can be acquired by fusing PM^F , PM^D and the original sinogram P^O. Finally, the final sinogram is reconstructed into the corrected CT image and metal information is compensated into the corrected CT image.Experiments on clinical images demonstrate that the proposed method can effectively reduce metal artifacts. A comparison with classical metal artifacts correction methods shows that the proposed metal artifacts correction method performs better in metal artifacts suppression and tissue feature preservation.展开更多
First arrival travel time tomography has achieved wide application. However, tomographic resolution is insufficient because geometry constraints cause rays to be unevenly distributed in the velocity model. The variabl...First arrival travel time tomography has achieved wide application. However, tomographic resolution is insufficient because geometry constraints cause rays to be unevenly distributed in the velocity model. The variable damping constraint method adopts uneven priori information to match uneven data distribution which can lessen the correlation between velocity correction values and ray coverage density. In this paper, we combine the variable damping constraint with a smoothness constraint which is added into the regularization equations in velocity inversion to avoid instability caused by only using the variable damping constraint method. The alpha-trimmed-mean filter is used to smooth and denoise intermediate results in the velocity inversion process. We use the LSQR algorithm to enhance the convergence rate and suppress error propagation in solving linear equations. In this paper, we apply the proposed tomographic method to perform velocity inversion using VSP data. The application in recovery test of the checkerboard model and velocity inversion of real VSP data show that the variable damping constraint method can improve tomographic quality because it can solve the effects of uneven ray coverage. In addition, the examples show that the tomographic result near geophones is much more reliable than other areas in the velocity model.展开更多
Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique （SIRT） inversion ...Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique （SIRT） inversion is used to deduce the stress redistribution around the longwall mining panel. The mining-induced microseismic events were recorded by mounting an array of receivers on the surface, above the active panel. After processing and filtering the seismic data, the three-dimensional tomography images of the p-wave velocity variations by SIRT passive seismic velocity tomography were provided. To display the velocity changes on coal seam level and subsequently to infer the stress redistribution, these three-dimensional tomograms into the coal seam level were sliced. In addition, the boundary element method （BEM） was used to simulate the stress redistribution. The results show that the inferred stresses from the passive seismic tomograms are conformed to numerical models and theoretical concept of the stress redistribution around the longwall panel. In velocity tomograms, the main zones of the stress redistribution arotmd the panel, including front and side abutment pressures, and gob stress are obvious and also the movement of stress zones along the face advancement is evident. Moreover, the effect of the advance rate of the face on the stress redistribution is demonstrated in tomography images. The research result proves that the SIRT passive seismic velocity tomography has an ultimate potential for monitoring the changes of stress redistribution around the longwall mining panel continuously and subsequently to improve safety of mining operations.展开更多
Ebstein's anomaly(EA) is a rare cardiac congenital malformation with displacement of septal and posterior tricuspid leaflets,resulting in atrialization of the right ventricle.We report a case of EA in which the et...Ebstein's anomaly(EA) is a rare cardiac congenital malformation with displacement of septal and posterior tricuspid leaflets,resulting in atrialization of the right ventricle.We report a case of EA in which the etiology of a malfunctioning prosthetic tricuspid valve is depicted on cardiac computed tomography to be as a result of thrombus lodged in the valve.展开更多
基金supported by the National Natural Science Foundation of China(Nos.41906048 and 91858215).
文摘We applied double-difference tomography to relocate seismic events and determine the lithospheric velocity structure beneath the New Britain Island arc and the South Bismarck Sea Basin,based on the local P wave arrival time dataset collected by the International Seismological Centre.Results of the seismic relocation and velocity inversion show that the subduction of Solomon Sea Plate along the New Britain Trench is spatially different above 150 km,and the subduction angle of the slab on the west side is higher than that on the east side.The relocated earthquakes also show that there are double seismic zones at the depths of about 30–90km beneath the New Britain Island Arc.The velocity structure shows that the dehydration of the subducting slab caused the low-velocity anomalies in mantle wedge above the slab,which are associated with the magmatic activities around the New Guinea-New Britain Island arc.Moreover,it shows that there is another low-velocity anomaly zone beneath the Bismarck mid-oceanic ridge with spatial variation.Beneath the west of the Bismarck mid-oceanic ridge,the low-velocity anomaly is weakly connected to the subducted Solomon Sea slab.Conversely,the low-velocity anomaly beneath the Manus Sea Basin is highly intertwined to the subducting slab and its mantle wedge,indicating that the subduction of the Solomon Sea Plate might be a key deep dynamic factor that drives the spreading of the Manus Sea Basin and the separation of the Bismarck Plate.
基金the Directorate General of Resources for Science Technologythe Higher Education of the Republic of Indonesia for granting a PMDSU scholarship to SR
文摘West Java in the western part of the Sunda Arc has a relatively high seismicity due to subduction activity and faults.In this study,double-difference tomography was used to obtain the 3D velocity tomograms of P and S waves beneath the western part of Java.To infer the geometry of the structure beneath the study area,precise earthquake hypo・center determination was first performed before tomographic imaging.For this,earthquake waveform data were extracted from the regional Meteorological,Climatological,Geophysical Agency(BMKG)network of Indonesia from South Sumatra to Central Java.The P and S arrival times for about 1,000 events in the period April 2009 to July 2016 were selected,the key features being events of magnitude>3,azimuthal gap<210°and number of phases>8.A nonlinear method using the oct-tree sampling algorithm from the NonLinLoc program was employed to determine the earthquake hypocenters.The hypocenter locations were then relocated using double-difference tomography(tomoDD).A significant reduction of travel-time(root mean square basis)and a better clustering of earthquakes were achieved which correlated well with the geological structure in West Java.Double-difference tomography was found to give a clear velocity structure,especially beneath the volcanic arc area,i.e.,under Mt Anak Krakatau,Mt Salak and the mountains complex in the southern part of West Java.Low velocity anomalies for the P and S waves as well as the vp/vs ratio below the volcanoes indicated possible partial melting of the upper mantle which ascended from the subducted slab beneath the volcanic arc.
基金Supported by the National Natural Science Foundation of China(61203021)the Key Science and Technology Program of Liaoning Province(2011216011)+1 种基金the Natural Science Foundation of Liaoning Province(2013020024)the Program for Liaoning Excellent Talents in Universities(LJQ2015061)
文摘Electrical capacitance tomography(ECT)has been applied to two-phase flow measurement in recent years.Image reconstruction algorithms play an important role in the successful applications of ECT.To solve the ill-posed and nonlinear inverse problem of ECT image reconstruction,a new ECT image reconstruction method based on fast linearized alternating direction method of multipliers(FLADMM)is proposed in this paper.On the basis of theoretical analysis of compressed sensing(CS),the data acquisition of ECT is regarded as a linear measurement process of permittivity distribution signal of pipe section.A new measurement matrix is designed and L1 regularization method is used to convert ECT inverse problem to a convex relaxation problem which contains prior knowledge.A new fast alternating direction method of multipliers which contained linearized idea is employed to minimize the objective function.Simulation data and experimental results indicate that compared with other methods,the quality and speed of reconstructed images are markedly improved.Also,the dynamic experimental results indicate that the proposed algorithm can ful fill the real-time requirement of ECT systems in the application.
文摘A trust region method is proposed to solve the problem of microwave tomography,which is very difficult to be solved for its ill-posedness and nonlinearity. Compared with the Levenberg-Marquardt method, this method introduces more a priori knowledge and might obtain better results, though the two methods are equal in some cases.
文摘Crosswell seismic tomography can be used to study the lateral variation of reservoirs, reservoir properties and the dynamic movement of fluids. In view of the instability of crosswell seismic tomography, the gradient method was improved by introducing regularization, and a gradient regularization method is presented in this paper. This method was verified by processing numerical simulation data and physical model data.
文摘Solution to impedance distribution in electrical impedance tomography （EIT） is an ill-posed nonlinear inverse problem. It is especially difficult to reconstruct an EIT image in the center area of a measured object. Tikhonov regularization with some prior information is a sound regnlarization method for static electrical impedance tomography under the condition that some true impedance distribution information is known a priori. This paper presents a direct search method （DSM） as pretreatment of image reconstruction through which one not only can construct a regularization matrix which may locate in areas of impedance change, but also can obtain an initial impedance distribution more similar to the true impedance distribution, as well as better current modes which can better distinguish the initial distribution and the true distribution. Simulation results indicate that, by using DSM, resolution in the center area of the measured object can be improved significantly.
基金supported by 14th Five-Year Plan major science and technology projects(no.KJGG2022-0201)。
文摘This paper presents a gradient-descent travel time tomography method for solving the acoustictype velocity model inversion problem.Similarly to the adjoint-state method,the proposed method is based on the Eikonal equation,enabling simultaneous calculation of contributions from all common-source receivers to the gradient.This overcomes the inefficiency inherent in conventional travel time tomography methods,which rely on a two-point ray tracing process.By directly calculating Fréchet derivatives,our method avoids the complex derivation processes associated with the adjoint-state method.The key to calculating the Fréchet derivatives is to calculate a so-called ray-path term.Consequently,compared to the adjoint-state method,the proposed method can explicitly obtain the ray paths,resulting in a more concise and intuitive derivation process.Furthermore,our method retains the benefits of the adjoint-state method,such as speed,low memory usage,and robustness.This paper focuses on elucidating the principles and algorithms for calculating the raypath term based on the fast sweeping method.The algorithms could be further speeded up by using parallel computational techniques.Synthetic tests demonstrate that our proposed travel time tomographic method accurately calculates ray paths,regardless of the complexity of the model and recording geometry.
文摘In this paper, a novel reconstruction method is presented for Near Infrared (NIR) 2-D imaging to recover optical absorption coefficients from laboratory phantom data. The main body of this work validates a new generation of highly efficient reconstruction algorithms called “Globally Convergent Method” (GCM) based upon actual measurements taken from brain-shape phantoms. It has been demonstrated in earlier studies using computer-simulated data that this type of reconstructions is stable for imaging complex distributions of optical absorption. The results in this paper demonstrate the excellent capability of GCM in working with experimental data measured from optical phantoms mimicking a rat brain with stroke.
基金National Key Research and Development Program of China(No.2021YFC3000103).
文摘To reduce the shielding effect of hardened layers on electrical resistivity tomography,a ratio method based on the distortion correction principle and the isolation coefficient is proposed.The effects of the resistivity and thickness of hardened concrete layers on the detection of target objects are explored.Both numerical simulations and indoor tank tests indicate that when the ratio method is employed to correct the original collected data,the maximum allowable error for the isolation coefficient should not exceed 1%.Notably,when the ratio of hardened layer thickness to electrode spacing does not exceed 1,correction through this method significantly enhances the recognition capability of target objects.However,when the hardened layer thickness is greater than the electrode spacing by a factor of 2 or more,the ratio method cannot achieve satisfactory results.The case study of flood control engineering detection in the Zhangxi section of the Huangpen River in Dongzhi County demonstrates that the detection effect after correction by the ratio method is comparable to that for the adjacent unhardened pavement,and the influence of the hardened layer is obviously weakened,resulting in more reliable results.
基金This work was supported by JST-CREST Grant Number JPMJCR1765,Japan.
文摘This paper proposes a new level-set-based shape recovery approach that can be applied to a wide range of binary tomography reconstructions.In this technique,we derive generic evolution equations for shape reconstruction in terms of the underlying level-set parameters.We show that using the appropriate basis function to parameterize the level-set function results in an optimization problem with a small number of parameters,which overcomes many of the problems associated with the traditional level-set approach.More concretely,in this paper,we use Gaussian functions as a basis function placed at sparse grid points to represent the parametric level-set function and provide more flexibility in the binary representation of the reconstructed image.In addition,we suggest a convex optimization method that can overcome the problem of the local minimum of the cost function by successfully recovering the coefficients of the basis function.Finally,we illustrate the performance of the proposed method using synthetic images and real X-ray CT projection data.We show that the proposed reconstruction method compares favorably to various state-of-the-art reconstruction techniques for limited-data tomography,and it is also relatively stable in the presence of modest amounts of noise.Furthermore,the shape representation using a compact Gaussian radial basis function works well.
基金supported by the National Natural Science Foundation of China(Nos.42574119,42274083,41974049)partly supported by the Urban Geological Survey Project of Linyi,Shandong Province,China(No.SDGP371300202102000468).
文摘A P-wave velocity model was built in the central southern of the Tanlu Fault based on double-difference tomography.The results suggest the presence of a low-velocity anomaly extending from the surface to a depth of 25 km around the Tanlu and Feixi Faults,representing fault-related fluids caused by partial melting.The relocated earthquakes indicate a significant concentration of seismic activity above 20 km around the Tanlu and Feixi Faults,suggesting that prominent fault systems possibly serve as conduits for the upward migration of deep minerals.The proposed geodynamic model,supported by geological and geophysical data,suggests that the migration of deep mineralized materials extends along the Tanlu Fault.The obtained results serve as a crucial foundation for elucidating the intricate process of mineralization in the central southern segment of the Tanlu Fault,thereby enhancing comprehension regarding the interaction among ore body formation,fault fluids,localized melting,and seismic activity.
基金supported by the Scientific Research Foundation for High-level Talents of Anhui University of Science and Technology under Grant 2024yjrc64the National Key R&D Program of China under Grant 2018YFC1504102。
文摘The Longmenshan(LMS)fault zone is located at the junction of the eastern Tibetan Plateau and the Sichuan Basin and is of great significance for studying regional tectonics and earthquake hazards.Although regional velocity models are available for the LMS fault zone,high-resolution velocity models are lacking.Therefore,a dense array of 240 short-period seismometers was deployed around the central segment of the LMS fault zone for approximately 30 days to monitor earthquakes and characterize fine structures of the fault zone.Considering the large quantity of observed seismic data,the data processing workflow consisted of deep learning-based automatic earthquake detection,phase arrival picking,and association.Compared with the earthquake catalog released by the China Earthquake Administration,many more earthquakes were detected by the dense array.Double-difference seismic tomography was adopted to determine V_(p),V_(s),and V_(p)/V_(s)models as well as earthquake locations.The checkerboard test showed that the velocity models have spatial resolutions of approximately 5 km in the horizontal directions and 2 km at depth.To the west of the Yingxiu–Beichuan Fault(YBF),the Precambrian Pengguan complex,where most of earthquakes occurred,is characterized by high velocity and low V_(p)/V_(s)values.In comparison,to the east of the YBF,the Upper Paleozoic to Jurassic sediments,where few earthquakes occurred,show low velocity and high V_(p)/V_(s)values.Our results suggest that the earthquake activity in the LMS fault zone is controlled by the strength of the rock compositions.When the high-resolution velocity models were combined with the relocated earthquakes,we were also able to delineate the fault geometry for different faults in the LMS fault zone.
基金supported by the National Natural Science Foundation of China(52174162)the Fundamental Research Funds for the Central Universities(FRF-TP-20-002A3).
文摘Accurate prediction of coal and gas outburst(CGO)hazards is paramount in gas disaster prevention and control.This paper endeavors to overcome the constraints posed by traditional prediction indexes when dealing with CGO incidents under low gas pressure conditions.In pursuit of this objective,we have studied and established a mechanical model of the working face under abnormal stress and the excitation energy conditions of CGO,and proposed a method for predicting the risk of CGO under abnormal stress.On site application verification shows that when a strong outburst hazard level prediction is issued,there is a high possibility of outburst disasters occurring.In one of the three locations where we predicted strong outburst hazards,a small outburst occurred,and the accuracy of the prediction was higher than the traditional drilling cuttings index S and drilling cuttings gas desorption index q.Finally,we discuss the mechanism of CGO under the action of stress anomalies.Based on the analysis of stress distribution changes and energy accumulation characteristics of coal under abnormal stress,this article believes that the increase in outburst risk caused by high stress abnormal gradient is mainly due to two reasons:(1)The high stress abnormal gradient leads to an increase in the plastic zone of the coal seam.After the working face advances,it indirectly leads to an increase in the gas expansion energy that can be released from the coal seam before reaching a new stress equilibrium.(2)Abnormal stress leads to increased peak stress of coal body in front of working face.When coal body in elastic area transforms to plastic area,its failure speed is accelerated,which induces accelerated gas desorption and aggravates the risk of outburst.
基金supported by the National Natural Science Foundation of China (Grant Nos.40334040 and 40974033)the Promoting Foundation for Advanced Persons of Talent of NCWU
文摘Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm （SAGA） and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A＇nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.
基金Open Research Fund of the Key Laboratory of Computer Netw ork and Information Integration of Ministry of Education of Southeast University(No.K93-9-2014-10C)the Scientific Research Foundation of Education Department of Anhui Province(No.KJ2014A186,SK2015A433)the National Basic Research Program of China(973 Program)(No.2010CB732503)
文摘To solve the problem that metal artifacts severely damage the clarity of the organization structure in computed tomography（CT） images, a sinogram fusion-based metal artifact correction method is proposed. First, the metal image is segmented from the original CT image by the pre-set threshold. The original CT image and metal image are forward projected into the original projection sinogram and metal projection sinogram, respectively. The interpolation-based correction method and mean filter are used to correct the original CT image and preserve the edge of the corrected CT image, respectively. The filtered CT image is forward projected into the filtered image sinogram. According to the position of the metal sinogram in the original sinogram and filtered image sinogram, the corresponding sinograms PM^D （ in the original sinogram） and PM^C （ in the filtered image sinogram）can be acquired from the original sinogram and filtered image sinogram, respectively. Then, PM^D and PM^C are fused into the fused metal sinogram PM^F according to a certain proportion.The final sinogram can be acquired by fusing PM^F , PM^D and the original sinogram P^O. Finally, the final sinogram is reconstructed into the corrected CT image and metal information is compensated into the corrected CT image.Experiments on clinical images demonstrate that the proposed method can effectively reduce metal artifacts. A comparison with classical metal artifacts correction methods shows that the proposed metal artifacts correction method performs better in metal artifacts suppression and tissue feature preservation.
基金supported by the China Important National Science and Technology Specific Projects (No2011ZX05024-001-02)
文摘First arrival travel time tomography has achieved wide application. However, tomographic resolution is insufficient because geometry constraints cause rays to be unevenly distributed in the velocity model. The variable damping constraint method adopts uneven priori information to match uneven data distribution which can lessen the correlation between velocity correction values and ray coverage density. In this paper, we combine the variable damping constraint with a smoothness constraint which is added into the regularization equations in velocity inversion to avoid instability caused by only using the variable damping constraint method. The alpha-trimmed-mean filter is used to smooth and denoise intermediate results in the velocity inversion process. We use the LSQR algorithm to enhance the convergence rate and suppress error propagation in solving linear equations. In this paper, we apply the proposed tomographic method to perform velocity inversion using VSP data. The application in recovery test of the checkerboard model and velocity inversion of real VSP data show that the variable damping constraint method can improve tomographic quality because it can solve the effects of uneven ray coverage. In addition, the examples show that the tomographic result near geophones is much more reliable than other areas in the velocity model.
文摘Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique （SIRT） inversion is used to deduce the stress redistribution around the longwall mining panel. The mining-induced microseismic events were recorded by mounting an array of receivers on the surface, above the active panel. After processing and filtering the seismic data, the three-dimensional tomography images of the p-wave velocity variations by SIRT passive seismic velocity tomography were provided. To display the velocity changes on coal seam level and subsequently to infer the stress redistribution, these three-dimensional tomograms into the coal seam level were sliced. In addition, the boundary element method （BEM） was used to simulate the stress redistribution. The results show that the inferred stresses from the passive seismic tomograms are conformed to numerical models and theoretical concept of the stress redistribution around the longwall panel. In velocity tomograms, the main zones of the stress redistribution arotmd the panel, including front and side abutment pressures, and gob stress are obvious and also the movement of stress zones along the face advancement is evident. Moreover, the effect of the advance rate of the face on the stress redistribution is demonstrated in tomography images. The research result proves that the SIRT passive seismic velocity tomography has an ultimate potential for monitoring the changes of stress redistribution around the longwall mining panel continuously and subsequently to improve safety of mining operations.
文摘Ebstein's anomaly(EA) is a rare cardiac congenital malformation with displacement of septal and posterior tricuspid leaflets,resulting in atrialization of the right ventricle.We report a case of EA in which the etiology of a malfunctioning prosthetic tricuspid valve is depicted on cardiac computed tomography to be as a result of thrombus lodged in the valve.
