Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional fin...Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.展开更多
Overt and harmful diabetes mellitus(DM)has detrimental effects on individuals and,by extension,the community.Among the microvascular DM complications is diabetic retinopathy(DR).DR may cause irreversible vision deteri...Overt and harmful diabetes mellitus(DM)has detrimental effects on individuals and,by extension,the community.Among the microvascular DM complications is diabetic retinopathy(DR).DR may cause irreversible vision deterioration in cases of poor blood glucose regulation.Changes in vascular permeability are key trigger points for diabetic macular edema(DME),a condition characterized by the accumulation of fluid in the macula.The development of vascular endothelial growth factor(VEGF)pathway inhibitors has provided a pathogenesis-based treatment approach for DME.Optical coherence tomography(OCT)provides highresolution imaging of the anatomy,including the aging of DME and its structural damage,in distinct morphologic subtypes of macular edema,thereby supporting the assessment of macular edema treatment.The availability of repeated OCT monitoring provides clinical reassurance through the treatment.OCT angiography(OCTA)provides retinal blood flow maps with high spatial resolution.The ability promotes an understanding of disease pathogenesis and facilitates the implementation of new therapeutic methods.This review compares the potential of OCT and OCTA in the diagnosis and treatment of DME,as well as their respective therapeutic applications.展开更多
Photoacoustic-computed tomography is a novel imaging technique that combines high absorption contrast and deep tissue penetration capability,enabling comprehensive three-dimensional imaging of biological targets.Howev...Photoacoustic-computed tomography is a novel imaging technique that combines high absorption contrast and deep tissue penetration capability,enabling comprehensive three-dimensional imaging of biological targets.However,the increasing demand for higher resolution and real-time imaging results in significant data volume,limiting data storage,transmission and processing efficiency of system.Therefore,there is an urgent need for an effective method to compress the raw data without compromising image quality.This paper presents a photoacoustic-computed tomography 3D data compression method and system based on Wavelet-Transformer.This method is based on the cooperative compression framework that integrates wavelet hard coding with deep learning-based soft decoding.It combines the multiscale analysis capability of wavelet transforms with the global feature modeling advantage of Transformers,achieving high-quality data compression and reconstruction.Experimental results using k-wave simulation suggest that the proposed compression system has advantages under extreme compression conditions,achieving a raw data compression ratio of up to 1:40.Furthermore,three-dimensional data compression experiment using in vivo mouse demonstrated that the maximum peak signal-to-noise ratio(PSNR)and structural similarity index(SSIM)values of reconstructed images reached 38.60 and 0.9583,effectively overcoming detail loss and artifacts introduced by raw data compression.All the results suggest that the proposed system can significantly reduce storage requirements and hardware cost,enhancing computational efficiency and image quality.These advantages support the development of photoacoustic-computed tomography toward higher efficiency,real-time performance and intelligent functionality.展开更多
This comprehensive study investigates the formation and evolution of intermetallic compounds during the solidification process of magnesium alloys using advanced micro X-ray computed tomography.By analyzing both commo...This comprehensive study investigates the formation and evolution of intermetallic compounds during the solidification process of magnesium alloys using advanced micro X-ray computed tomography.By analyzing both common industrial Mg-Al-Zn alloys and a novel rare earth-containing Mg-Ni-Gd-Y alloy,we aim to characterize the nucleation,growth,and distribution of Al-Mn and eutectic intermetallics across various stages of solidification.The non destructive imaging technique employed in this research provides high-resolution,three-dimensional insights into the microstructural development,allowing for a detailed examination of the morphology,spatial arrangement,and interconnectivity of intermetallic phases.This approach overcomes limitations of traditional two-dimensional metallographic methods,offering a more comprehensive understanding of the complex three-dimensional structures formed during solidification.展开更多
Gastric cancer,including stomach tumors,poses significant health challenges due to late-stage diagnosis and limited early detection methods.Accurate imaging and precise tumor margin demarcation are critical for effect...Gastric cancer,including stomach tumors,poses significant health challenges due to late-stage diagnosis and limited early detection methods.Accurate imaging and precise tumor margin demarcation are critical for effective treatment planning and improved patient outcomes.Conventional imaging techniques,such as endoscopy and histopathology,provide valuable diagnostic information but cannot offer real-time assessment.We aim to explore the application of optical coherence tomography(OCT),combined with statistical and machine learning methods,for rapid tumor margin demarcation in gastrointestinal(GI)cancer tissues in exvivo.GI tumor specimens from 17 patients were imaged intraoperatively with OCT.Quantitative features were extracted from the images,and statistical and machine learning analyses were applied to distinguish tissue types.Subsequent histopathological evaluation was used as the reference standard for validation.The combination of OCT imaging and data-driven analysis enables clear differentiation between healthy,marginal and tumor tissues in near real-time.The KNN model achieved the highest classification accuracy(0.921±0.040),with SVM(0.906±0.038)and Extra Trees(0.901±0.034)also demonstrating robust performance in discriminating tissue margins.OCT findings demonstrate strong agreement with histopathology.Integrated with statistical and machine learning techniques,OCT enables rapid assessment of GI tumor margins.This approach provides quantitative objectivity and real-time feedback,closely approximating histopathological analysis and supporting improved surgical decision-making.展开更多
Background:Diabetic macular edema is a prevalent retinal condition and a leading cause of visual impairment among diabetic patients’Early detection of affected areas is beneficial for effective diagnosis and treatmen...Background:Diabetic macular edema is a prevalent retinal condition and a leading cause of visual impairment among diabetic patients’Early detection of affected areas is beneficial for effective diagnosis and treatment.Traditionally,diagnosis relies on optical coherence tomography imaging technology interpreted by ophthalmologists.However,this manual image interpretation is often slow and subjective.Therefore,developing automated segmentation for macular edema images is essential to enhance to improve the diagnosis efficiency and accuracy.Methods:In order to improve clinical diagnostic efficiency and accuracy,we proposed a SegNet network structure integrated with a convolutional block attention module(CBAM).This network introduces a multi-scale input module,the CBAM attention mechanism,and jump connection.The multi-scale input module enhances the network’s perceptual capabilities,while the lightweight CBAM effectively fuses relevant features across channels and spatial dimensions,allowing for better learning of varying information levels.Results:Experimental results demonstrate that the proposed network achieves an IoU of 80.127%and an accuracy of 99.162%.Compared to the traditional segmentation network,this model has fewer parameters,faster training and testing speed,and superior performance on semantic segmentation tasks,indicating its highly practical applicability.Conclusion:The C-SegNet proposed in this study enables accurate segmentation of Diabetic macular edema lesion images,which facilitates quicker diagnosis for healthcare professionals.展开更多
The Urumqi foreland thrust tectonic belt exhibits complex geological structures and strong seismicity.Imaging its shallow crustal structure is of great significance for understanding its tectonic mechanism and seismog...The Urumqi foreland thrust tectonic belt exhibits complex geological structures and strong seismicity.Imaging its shallow crustal structure is of great significance for understanding its tectonic mechanism and seismogenic environment.We obtained a high-resolution S-wave velocity model of the shallow crust at depths of 0–8 km using ambient noise tomography applied to data from a dense seismic array.Sediments are generally thinner in the southeast and thicker in the northwest,with a maximum thickness of more than 8 km.Variations in the velocity structure near the Xishan,Wanyaogou,and Yamalike faults indicate that their formation was related to differences in the physical properties on either side of the fault.In addition,the faults exhibit thrusting of the low-velocity sides towards the high-velocity sides.In the study area,earthquakes rarely occur at depths of less than 3 km and are mostly concentrated in the high-velocity zone in the southern part.Below 3 km depth,more earthquakes were observed,mainly distributed near faults or in relatively high-velocity areas in the southern part.This suggests that high-velocity structures are more prone to stress accumulation,resulting in earthquakes.At 6–8 km depth,the densely distributed earthquakes in the northwestern part of the Bogda mountains are well-aligned with the northwest-oriented low-velocity zone observed in this study,suggesting that this weak zone likely controls seismicity in this area.展开更多
BACKGROUND Early screening,preoperative staging,and diagnosis of lymph node metastasis are crucial for improving the prognosis of gastric cancer(GC).AIM To evaluate the diagnostic value of combined multidetector compu...BACKGROUND Early screening,preoperative staging,and diagnosis of lymph node metastasis are crucial for improving the prognosis of gastric cancer(GC).AIM To evaluate the diagnostic value of combined multidetector computed tomography(MDCT)and gastrointestinal endoscopy for GC screening,preoperative staging,and lymph node metastasis detection,thereby providing a reference for clinical diagnosis and treatment.METHODS In this retrospective study clinical and imaging data of 134 patients with suspected GC who were admitted between January 2023 and October 2024 were initially reviewed.According to the inclusion and exclusion criteria,102 patients were finally enrolled in the analysis.All enrolled patients had undergone both MDCT and gastrointestinal endoscopy examinations prior to surgical intervention.Preoperative clinical staging and lymph node metastasis findings were compared with pathological results.RESULTS The combined use of MDCT and gastrointestinal endoscopy demonstrated a sensitivity of 98.53%,specificity of 97.06%,accuracy of 98.04%,positive predictive value of 98.53%,and negative predictive value of 97.06%for diagnosing GC.These factors were all significantly higher than those of MDCT or endoscopy alone(P<0.05).The accuracy rates of the combined approach for detecting clinical T and N stages were 97.06%and 92.65%,respectively,outperforming MDCT alone(86.76% and 79.41%)and endoscopy alone(85.29% and 70.59%)(P<0.05).Among 68 patients with confirmed GC,50(73.53%)were pathologically diagnosed with lymph node metastasis.The accuracy for detecting lymph node metastasis was 66.00%with endoscopy,76.00%with MDCT,and 92.00% with the combined approach,all with statistically significant differences(P<0.05).CONCLUSION The combined application of MDCT and gastrointestinal endoscopy enhanced diagnostic accuracy for GC,provided greater consistency in preoperative staging,and improved the detection of lymph node metastasis,thereby demonstrating significant clinical utility.展开更多
BACKGROUND:Individualized positive end-expiratory pressure(PEEP)titration is a crucial technique in mechanical ventilation therapy for acute respiratory distress syndrome(ARDS)patients with intra-abdominal hypertensio...BACKGROUND:Individualized positive end-expiratory pressure(PEEP)titration is a crucial technique in mechanical ventilation therapy for acute respiratory distress syndrome(ARDS)patients with intra-abdominal hypertension(IAH).This study aimed to evaluate the eff ectiveness of electrical impedance tomography(EIT)-guided PEEP titration in this population.METHODS:This prospective study enrolled 36 ARDS patients,including 22 patients with IAH and 14 without IAH.All the patients underwent EIT-guided PEEP titration at the intersection point between alveolar overdistension and collapse during a decremental PEEP trial.The changes in pulmonary ventilation distribution,respiratory mechanics and hemodynamics during the titration process were observed.RESULTS:After EIT-guided PEEP titration was performed,the PEEP,peak inspiratory pressure and plateau pressure increased significantly(P<0.05).Furthermore,no significant differences were observed in respiratory system compliance,tidal volume,driving pressure,or the 4*DP+RR index between the two groups(P>0.05).The mechanical power increased in the non-IAH(NIAH)group after PEEP titration(P<0.05).Ventilation in gravity-dependent lung regions significantly increased(P<0.05),and the oxygenation index(PaO2/FiO2)improved signifi cantly(P<0.05)in both groups.However,blood pressure,heart rate,respiratory rate,central venous pressure,and lactate levels did not signifi cantly change.In the IAH group,the PaO2/FiO2 ratio improved less than that in the NIAH group did(P<0.05).CONCLUSION:In our study,individualized PEEP titration guided by EIT improved oxygenation in ARDS patients with concomitant IAH without signifi cantly aff ecting hemodynamics.The presence of IAH may limit the improvement of oxygenation during EIT-guided PEEP titration.展开更多
AIM:To investigate the effects of chronic alcohol consumption on retinal microcirculation by comparing different alcohol-consuming groups using optical coherence tomography(OCT)and OCT angiography(OCTA).METHODS:This o...AIM:To investigate the effects of chronic alcohol consumption on retinal microcirculation by comparing different alcohol-consuming groups using optical coherence tomography(OCT)and OCT angiography(OCTA).METHODS:This observational clinical study utilized a cross-sectional and prospective design,focusing on chronic alcohol consumers alongside a non-consuming control group.OCT/OCTA imaging parameters including central retinal subfield thickness(CST),subfoveal choroidal thickness(SCT),foveal avascular zone(FAZ)and vessel density(VD)in the superficial and deep capillary plexuses in both the macular and optic disc(OD)regions were recorded.Data were analyzed using SPSS 15.0;descriptive statistics were reported,group comparisons were performed with Chisquare,Kruskal–Wallis,and Bonferroni-corrected Mann–Whitney U tests,and relationships were assessed using Spearman correlation,with statistical significance set at P<0.05.RESULTS:A total of 160 eyes of 160 participants(110 females and 50 males with mean age 38.7±9.9y)who don’t smoke were divided into five groups:never,occasional,monthly,weekly and daily drinkers.The mean CST was 216.6±14.2μm and the mean SCT was 358.9±84.5μm.There was no statistically significantly difference in CST and SCT among the groups(P=0.890,0.799).Foveal superficial capillary plexuses(SCPs)VD was higher in monthly drinkers compared to occasional drinkers(P=0.015).Foveal VD in deep capillary plexus was also higher in monthly drinkers than in never and occasional drinkers(P=0.004,0.006).Nasal SCPs VD at the OD was higher in monthly drinkers compared to never drinkers(P=0.005).There was no significant difference FAZ area among the groups(P=0.071).CONCLUSION:Both superficial and deep microvascular structures in the inferior quadrants of macula are positively correlated with frequency of alcohol use.Also in our study results is that the monthly drinker group has uniquely higher VDs in both macula and OD.This leads us to consider moderate alcohol consumption may also have protective effects on retinal microcirculation.展开更多
Although previous studies have demonstrated that transcranial focused ultrasound stimulation protects the ischemic brain,clear criteria for the stimulation time window and intensity are lacking.Electrical impedance to...Although previous studies have demonstrated that transcranial focused ultrasound stimulation protects the ischemic brain,clear criteria for the stimulation time window and intensity are lacking.Electrical impedance tomography enables real-time monitoring of changes in cerebral blood perfusion within the ischemic brain,but investigating the feasibility of using this method to assess post-stroke rehabilitation in vivo remains critical.In this study,ischemic stroke was induced in rats through middle cerebral artery occlusion surgery.Transcranial focused ultrasound stimulation was used to treat the rat model of ischemia,and electrical impedance tomography was used to measure impedance during both the acute stage of ischemia and the rehabilitation stage following the stimulation.Electrical impedance tomography results indicated that cerebral impedance increased after the onset of ischemia and decreased following transcranial focused ultrasound stimulation.Furthermore,the stimulation promoted motor function recovery,reduced cerebral infarction volume in the rat model of ischemic stroke,and induced the expression of brain-derived neurotrophic factor in the ischemic brain.Our results also revealed a significant correlation between the impedance of the ischemic brain post-intervention and improvements in behavioral scores and infarct volume.This study shows that daily administration of transcranial focused ultrasound stimulation for 20 minutes to the ischemic hemisphere 24 hours after cerebral ischemia enhanced motor recovery in a rat model of ischemia.Additionally,our findings indicate that electrical impedance tomography can serve as a valuable tool for quantitatively evaluating rehabilitation after ischemic stroke in vivo.These findings suggest the feasibility of using impedance data collected via electrical impedance tomography to clinically assess the effects of rehabilitatory interventions for patients with ischemic stroke.展开更多
Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of t...Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.展开更多
Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a...Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.展开更多
To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of ...To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.展开更多
During seismic wave propagation on a free surface, a strong material contrast boundary develops in response to interference by P- and S- waves to create a surface-wave phenomenon. To accurately determine the effects o...During seismic wave propagation on a free surface, a strong material contrast boundary develops in response to interference by P- and S- waves to create a surface-wave phenomenon. To accurately determine the effects of this interface on surface-wave propagation, the boundary conditions must be accurately modeled. In this paper, we present a numerical approach based on the dynamic poroelasticity for a space–time-domain staggered-grid finite-difference simulation in porous media that contain a free-surface boundary. We propose a generalized stess mirror formulation of the free-surface boundary for solids and fluids in porous media for the grid mesh on which lays the free-surface plane. Its analog is that used for elastic media, which is suitable for precise and stable Rayleigh-type surface-wave modeling. The results of our analysis of first kind of Rayleigh （R1） waves obtained by this model demonstrate that the discretization of the mesh in a similar way to that for elastic media can realize stable numerical solutions with acceptable precision. We present numerical examples demonstrating the efficiency and accuracy of our proposed method.展开更多
Numerical simulation in the frequency-space domain has inherent advantages, such as： it is possible to simulate wave propagation from multiple sources simultaneously; there are no cumulative errors; only the interest...Numerical simulation in the frequency-space domain has inherent advantages, such as： it is possible to simulate wave propagation from multiple sources simultaneously; there are no cumulative errors; only the interesting frequencies can be selected; and it is more suitable for wave propagation in viscoelastic media. The only obstacle to using the method is the requirement of huge computer storage. We extend the compressed format for storing the coefficient matrix. It can reduce the required computer storage dramatically. We get the optimal coefficients by least-squares method to suppress the numerical dispersion and adopt the perfectly matched layer （PML） boundary conditions to eliminate the artificial boundary reflections. Using larger grid intervals decreases computer storage requirements and provides high computational efficiency. Numerical experiments demonstrate that these means are economic and effective, providing a good basis for elastic wave imaging and inversion.展开更多
The staggered-grid finite-difference （SGFD） method has been widely used in seismic forward modeling. The precision of the forward modeling results directly affects the results of the subsequent seismic inversion and...The staggered-grid finite-difference （SGFD） method has been widely used in seismic forward modeling. The precision of the forward modeling results directly affects the results of the subsequent seismic inversion and migration. Numerical dispersion is one of the problems in this method. The window function method can reduce dispersion by replacing the finite-difference operators with window operators, obtained by truncating the spatial convolution series of the pseudospectral method. Although the window operators have high precision in the low-wavenumber domain, their precision decreases rapidly in the high-wavenumber domain. We develop a least squares optimization method to enhance the precision of operators obtained by the window function method. We transform the SGFD problem into a least squares problem and find the best solution iteratively. The window operator is chosen as the initial value and the optimized domain is set by the error threshold. The conjugate gradient method is also adopted to increase the stability of the solution. Approximation error analysis and numerical simulation results suggest that the proposed method increases the precision of the window function operators and decreases the numerical dispersion.展开更多
The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the M...The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the Maxwell＇s equations, and the magnetic components at the center of each edge of the grid cell are evaluated by applying the staggered-grid finite-difference method. The tipper and its divergence are derived to complete the 3D ZTEM forward modeling algorithm. A synthetic model is then used to compare the responses with those of 2D finite-element forward modeling to verify the accuracy of the algorithm. ZTEM offers high horizontal resolution to both simple and complex distributions of conductivity. This work is the theoretical foundation for the interpretation of ZTEM data and the study of 3D ZTEM inversion.展开更多
Numerical simulation in transverse isotropic media with tilted symmetry axis（TTI） using the standard staggered-grid finite-difference scheme（SSG）results in errors caused by averaging or interpolation. In order to ...Numerical simulation in transverse isotropic media with tilted symmetry axis（TTI） using the standard staggered-grid finite-difference scheme（SSG）results in errors caused by averaging or interpolation. In order to eliminate the errors, a method of rotated staggered-grid finite-difference scheme（RSG） is proposed. However, the RSG brings serious numerical dispersion. The compact staggered-grid finite-difference scheme（CSG） is an implicit difference scheme, which use fewer grid points to suppress dispersion more effectively than the SSG. This paper combines the CSG with the RSG to derive a rotated staggered-grid compact finite-difference scheme（RSGC）. The numerical experiments indicate that the RSGC has weaker numerical dispersion and better accuracy than the RSG.展开更多
To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, t...To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.展开更多
基金supported by National Key Research and Development Program of China(No.2022YFF0800602)Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology(No.2022B1212010002)Shenzhen Science and Technology Program(No.KQTD20170810111725321).
文摘Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.
文摘Overt and harmful diabetes mellitus(DM)has detrimental effects on individuals and,by extension,the community.Among the microvascular DM complications is diabetic retinopathy(DR).DR may cause irreversible vision deterioration in cases of poor blood glucose regulation.Changes in vascular permeability are key trigger points for diabetic macular edema(DME),a condition characterized by the accumulation of fluid in the macula.The development of vascular endothelial growth factor(VEGF)pathway inhibitors has provided a pathogenesis-based treatment approach for DME.Optical coherence tomography(OCT)provides highresolution imaging of the anatomy,including the aging of DME and its structural damage,in distinct morphologic subtypes of macular edema,thereby supporting the assessment of macular edema treatment.The availability of repeated OCT monitoring provides clinical reassurance through the treatment.OCT angiography(OCTA)provides retinal blood flow maps with high spatial resolution.The ability promotes an understanding of disease pathogenesis and facilitates the implementation of new therapeutic methods.This review compares the potential of OCT and OCTA in the diagnosis and treatment of DME,as well as their respective therapeutic applications.
基金supported by the National Key R&D Program of China[Grant No.2023YFF0713600]the National Natural Science Foundation of China[Grant No.62275062]+3 种基金Project of Shandong Innovation and Startup Community of High-end Medical Apparatus and Instruments[Grant No.2023-SGTTXM-002 and 2024-SGTTXM-005]the Shandong Province Technology Innovation Guidance Plan(Central Leading Local Science and Technology Development Fund)[Grant No.YDZX2023115]the Taishan Scholar Special Funding Project of Shandong Provincethe Shandong Laboratory of Advanced Biomaterials and Medical Devices in Weihai[Grant No.ZL202402].
文摘Photoacoustic-computed tomography is a novel imaging technique that combines high absorption contrast and deep tissue penetration capability,enabling comprehensive three-dimensional imaging of biological targets.However,the increasing demand for higher resolution and real-time imaging results in significant data volume,limiting data storage,transmission and processing efficiency of system.Therefore,there is an urgent need for an effective method to compress the raw data without compromising image quality.This paper presents a photoacoustic-computed tomography 3D data compression method and system based on Wavelet-Transformer.This method is based on the cooperative compression framework that integrates wavelet hard coding with deep learning-based soft decoding.It combines the multiscale analysis capability of wavelet transforms with the global feature modeling advantage of Transformers,achieving high-quality data compression and reconstruction.Experimental results using k-wave simulation suggest that the proposed compression system has advantages under extreme compression conditions,achieving a raw data compression ratio of up to 1:40.Furthermore,three-dimensional data compression experiment using in vivo mouse demonstrated that the maximum peak signal-to-noise ratio(PSNR)and structural similarity index(SSIM)values of reconstructed images reached 38.60 and 0.9583,effectively overcoming detail loss and artifacts introduced by raw data compression.All the results suggest that the proposed system can significantly reduce storage requirements and hardware cost,enhancing computational efficiency and image quality.These advantages support the development of photoacoustic-computed tomography toward higher efficiency,real-time performance and intelligent functionality.
基金Project(2023YFB4606200)supported by the National Key Research and Development Program of ChinaProject(2023-SSRF-HZ-503114-2)supported by Shanghai Synchrotron Radiation Facility,Instrument BL16U2,China。
文摘This comprehensive study investigates the formation and evolution of intermetallic compounds during the solidification process of magnesium alloys using advanced micro X-ray computed tomography.By analyzing both common industrial Mg-Al-Zn alloys and a novel rare earth-containing Mg-Ni-Gd-Y alloy,we aim to characterize the nucleation,growth,and distribution of Al-Mn and eutectic intermetallics across various stages of solidification.The non destructive imaging technique employed in this research provides high-resolution,three-dimensional insights into the microstructural development,allowing for a detailed examination of the morphology,spatial arrangement,and interconnectivity of intermetallic phases.This approach overcomes limitations of traditional two-dimensional metallographic methods,offering a more comprehensive understanding of the complex three-dimensional structures formed during solidification.
基金supported by the financial support received from the Indian Council of Medical Research-Department of Health Research(ICMRDHR-CoE-5/3/8/5/2019/I-MDMS)Medical Device and Diagnostics Mission Secretariat(MDMS)and Foundation for Centre for Healthcare Entrepreneurship(CfHE).
文摘Gastric cancer,including stomach tumors,poses significant health challenges due to late-stage diagnosis and limited early detection methods.Accurate imaging and precise tumor margin demarcation are critical for effective treatment planning and improved patient outcomes.Conventional imaging techniques,such as endoscopy and histopathology,provide valuable diagnostic information but cannot offer real-time assessment.We aim to explore the application of optical coherence tomography(OCT),combined with statistical and machine learning methods,for rapid tumor margin demarcation in gastrointestinal(GI)cancer tissues in exvivo.GI tumor specimens from 17 patients were imaged intraoperatively with OCT.Quantitative features were extracted from the images,and statistical and machine learning analyses were applied to distinguish tissue types.Subsequent histopathological evaluation was used as the reference standard for validation.The combination of OCT imaging and data-driven analysis enables clear differentiation between healthy,marginal and tumor tissues in near real-time.The KNN model achieved the highest classification accuracy(0.921±0.040),with SVM(0.906±0.038)and Extra Trees(0.901±0.034)also demonstrating robust performance in discriminating tissue margins.OCT findings demonstrate strong agreement with histopathology.Integrated with statistical and machine learning techniques,OCT enables rapid assessment of GI tumor margins.This approach provides quantitative objectivity and real-time feedback,closely approximating histopathological analysis and supporting improved surgical decision-making.
基金supported by the Guangdong Pharmaceutical University 2024 Higher Education Research Projects(GKP202403,GMP202402)the Guangdong Pharmaceutical University College Students’Innovation and Entrepreneurship Training Programs(Grant No.202504302033,202504302034,202504302036,and 202504302244).
文摘Background:Diabetic macular edema is a prevalent retinal condition and a leading cause of visual impairment among diabetic patients’Early detection of affected areas is beneficial for effective diagnosis and treatment.Traditionally,diagnosis relies on optical coherence tomography imaging technology interpreted by ophthalmologists.However,this manual image interpretation is often slow and subjective.Therefore,developing automated segmentation for macular edema images is essential to enhance to improve the diagnosis efficiency and accuracy.Methods:In order to improve clinical diagnostic efficiency and accuracy,we proposed a SegNet network structure integrated with a convolutional block attention module(CBAM).This network introduces a multi-scale input module,the CBAM attention mechanism,and jump connection.The multi-scale input module enhances the network’s perceptual capabilities,while the lightweight CBAM effectively fuses relevant features across channels and spatial dimensions,allowing for better learning of varying information levels.Results:Experimental results demonstrate that the proposed network achieves an IoU of 80.127%and an accuracy of 99.162%.Compared to the traditional segmentation network,this model has fewer parameters,faster training and testing speed,and superior performance on semantic segmentation tasks,indicating its highly practical applicability.Conclusion:The C-SegNet proposed in this study enables accurate segmentation of Diabetic macular edema lesion images,which facilitates quicker diagnosis for healthcare professionals.
基金supported by the Key Research and Development Program of the Xinjiang Uygur Autonomous Region(No.2020B03006-1)the National Natural Science Foundation of China(Nos.42304069,and 42102275).
文摘The Urumqi foreland thrust tectonic belt exhibits complex geological structures and strong seismicity.Imaging its shallow crustal structure is of great significance for understanding its tectonic mechanism and seismogenic environment.We obtained a high-resolution S-wave velocity model of the shallow crust at depths of 0–8 km using ambient noise tomography applied to data from a dense seismic array.Sediments are generally thinner in the southeast and thicker in the northwest,with a maximum thickness of more than 8 km.Variations in the velocity structure near the Xishan,Wanyaogou,and Yamalike faults indicate that their formation was related to differences in the physical properties on either side of the fault.In addition,the faults exhibit thrusting of the low-velocity sides towards the high-velocity sides.In the study area,earthquakes rarely occur at depths of less than 3 km and are mostly concentrated in the high-velocity zone in the southern part.Below 3 km depth,more earthquakes were observed,mainly distributed near faults or in relatively high-velocity areas in the southern part.This suggests that high-velocity structures are more prone to stress accumulation,resulting in earthquakes.At 6–8 km depth,the densely distributed earthquakes in the northwestern part of the Bogda mountains are well-aligned with the northwest-oriented low-velocity zone observed in this study,suggesting that this weak zone likely controls seismicity in this area.
文摘BACKGROUND Early screening,preoperative staging,and diagnosis of lymph node metastasis are crucial for improving the prognosis of gastric cancer(GC).AIM To evaluate the diagnostic value of combined multidetector computed tomography(MDCT)and gastrointestinal endoscopy for GC screening,preoperative staging,and lymph node metastasis detection,thereby providing a reference for clinical diagnosis and treatment.METHODS In this retrospective study clinical and imaging data of 134 patients with suspected GC who were admitted between January 2023 and October 2024 were initially reviewed.According to the inclusion and exclusion criteria,102 patients were finally enrolled in the analysis.All enrolled patients had undergone both MDCT and gastrointestinal endoscopy examinations prior to surgical intervention.Preoperative clinical staging and lymph node metastasis findings were compared with pathological results.RESULTS The combined use of MDCT and gastrointestinal endoscopy demonstrated a sensitivity of 98.53%,specificity of 97.06%,accuracy of 98.04%,positive predictive value of 98.53%,and negative predictive value of 97.06%for diagnosing GC.These factors were all significantly higher than those of MDCT or endoscopy alone(P<0.05).The accuracy rates of the combined approach for detecting clinical T and N stages were 97.06%and 92.65%,respectively,outperforming MDCT alone(86.76% and 79.41%)and endoscopy alone(85.29% and 70.59%)(P<0.05).Among 68 patients with confirmed GC,50(73.53%)were pathologically diagnosed with lymph node metastasis.The accuracy for detecting lymph node metastasis was 66.00%with endoscopy,76.00%with MDCT,and 92.00% with the combined approach,all with statistically significant differences(P<0.05).CONCLUSION The combined application of MDCT and gastrointestinal endoscopy enhanced diagnostic accuracy for GC,provided greater consistency in preoperative staging,and improved the detection of lymph node metastasis,thereby demonstrating significant clinical utility.
基金PEEP titration in ARDS patients using EIT combined with lung ultrasound,Key Laboratory of Emergency Trauma Research,Ministry of Education (KLET-202201)airway clearance protocol in ICU mechanically ventilated patients based on electrical impedance imaging technology,Natural Science Foundation of Hunan Province (2024JJ9148)effects of end expiratory positive pressure on lung re-expansion in patients with ARDS and intra-abdominal hypertension monitored using lung ultrasound,Natural Science Foundation of Hunan Province (2023JJ60308)
文摘BACKGROUND:Individualized positive end-expiratory pressure(PEEP)titration is a crucial technique in mechanical ventilation therapy for acute respiratory distress syndrome(ARDS)patients with intra-abdominal hypertension(IAH).This study aimed to evaluate the eff ectiveness of electrical impedance tomography(EIT)-guided PEEP titration in this population.METHODS:This prospective study enrolled 36 ARDS patients,including 22 patients with IAH and 14 without IAH.All the patients underwent EIT-guided PEEP titration at the intersection point between alveolar overdistension and collapse during a decremental PEEP trial.The changes in pulmonary ventilation distribution,respiratory mechanics and hemodynamics during the titration process were observed.RESULTS:After EIT-guided PEEP titration was performed,the PEEP,peak inspiratory pressure and plateau pressure increased significantly(P<0.05).Furthermore,no significant differences were observed in respiratory system compliance,tidal volume,driving pressure,or the 4*DP+RR index between the two groups(P>0.05).The mechanical power increased in the non-IAH(NIAH)group after PEEP titration(P<0.05).Ventilation in gravity-dependent lung regions significantly increased(P<0.05),and the oxygenation index(PaO2/FiO2)improved signifi cantly(P<0.05)in both groups.However,blood pressure,heart rate,respiratory rate,central venous pressure,and lactate levels did not signifi cantly change.In the IAH group,the PaO2/FiO2 ratio improved less than that in the NIAH group did(P<0.05).CONCLUSION:In our study,individualized PEEP titration guided by EIT improved oxygenation in ARDS patients with concomitant IAH without signifi cantly aff ecting hemodynamics.The presence of IAH may limit the improvement of oxygenation during EIT-guided PEEP titration.
文摘AIM:To investigate the effects of chronic alcohol consumption on retinal microcirculation by comparing different alcohol-consuming groups using optical coherence tomography(OCT)and OCT angiography(OCTA).METHODS:This observational clinical study utilized a cross-sectional and prospective design,focusing on chronic alcohol consumers alongside a non-consuming control group.OCT/OCTA imaging parameters including central retinal subfield thickness(CST),subfoveal choroidal thickness(SCT),foveal avascular zone(FAZ)and vessel density(VD)in the superficial and deep capillary plexuses in both the macular and optic disc(OD)regions were recorded.Data were analyzed using SPSS 15.0;descriptive statistics were reported,group comparisons were performed with Chisquare,Kruskal–Wallis,and Bonferroni-corrected Mann–Whitney U tests,and relationships were assessed using Spearman correlation,with statistical significance set at P<0.05.RESULTS:A total of 160 eyes of 160 participants(110 females and 50 males with mean age 38.7±9.9y)who don’t smoke were divided into five groups:never,occasional,monthly,weekly and daily drinkers.The mean CST was 216.6±14.2μm and the mean SCT was 358.9±84.5μm.There was no statistically significantly difference in CST and SCT among the groups(P=0.890,0.799).Foveal superficial capillary plexuses(SCPs)VD was higher in monthly drinkers compared to occasional drinkers(P=0.015).Foveal VD in deep capillary plexus was also higher in monthly drinkers than in never and occasional drinkers(P=0.004,0.006).Nasal SCPs VD at the OD was higher in monthly drinkers compared to never drinkers(P=0.005).There was no significant difference FAZ area among the groups(P=0.071).CONCLUSION:Both superficial and deep microvascular structures in the inferior quadrants of macula are positively correlated with frequency of alcohol use.Also in our study results is that the monthly drinker group has uniquely higher VDs in both macula and OD.This leads us to consider moderate alcohol consumption may also have protective effects on retinal microcirculation.
基金supported by the Fundamental Research Funds for the Central Universities,Nos.G2021KY05107,G2021KY05101the National Natural Science Foundation of China,Nos.32071316,32211530049+1 种基金the Natural Science Foundation of Shaanxi Province,No.2022-JM482the Education and Teaching Reform Funds for the Central Universities,No.23GZ230102(all to LL and HH).
文摘Although previous studies have demonstrated that transcranial focused ultrasound stimulation protects the ischemic brain,clear criteria for the stimulation time window and intensity are lacking.Electrical impedance tomography enables real-time monitoring of changes in cerebral blood perfusion within the ischemic brain,but investigating the feasibility of using this method to assess post-stroke rehabilitation in vivo remains critical.In this study,ischemic stroke was induced in rats through middle cerebral artery occlusion surgery.Transcranial focused ultrasound stimulation was used to treat the rat model of ischemia,and electrical impedance tomography was used to measure impedance during both the acute stage of ischemia and the rehabilitation stage following the stimulation.Electrical impedance tomography results indicated that cerebral impedance increased after the onset of ischemia and decreased following transcranial focused ultrasound stimulation.Furthermore,the stimulation promoted motor function recovery,reduced cerebral infarction volume in the rat model of ischemic stroke,and induced the expression of brain-derived neurotrophic factor in the ischemic brain.Our results also revealed a significant correlation between the impedance of the ischemic brain post-intervention and improvements in behavioral scores and infarct volume.This study shows that daily administration of transcranial focused ultrasound stimulation for 20 minutes to the ischemic hemisphere 24 hours after cerebral ischemia enhanced motor recovery in a rat model of ischemia.Additionally,our findings indicate that electrical impedance tomography can serve as a valuable tool for quantitatively evaluating rehabilitation after ischemic stroke in vivo.These findings suggest the feasibility of using impedance data collected via electrical impedance tomography to clinically assess the effects of rehabilitatory interventions for patients with ischemic stroke.
基金This research was supported by the National Nature Science Foundation of China （No. 41074100） and the Program for NewCentury Excellent Talents in the University of the Ministry of Education of China （No. NCET- 10-0812）.
文摘Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.
基金supported by the National Natural Science Foundation of China(No.41474110)Shell Ph.D. Scholarship to support excellence in geophysical research
文摘Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.
文摘To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.
基金sponsed by National Natural Science Foundation of China(NSFC,Grant No.41304077)the Natural Basic Research Program of China(the“973 Project,”Grant No.2013CB733303)Postdoctoral Science Foundation of China(Grant No.2014T70740)
文摘During seismic wave propagation on a free surface, a strong material contrast boundary develops in response to interference by P- and S- waves to create a surface-wave phenomenon. To accurately determine the effects of this interface on surface-wave propagation, the boundary conditions must be accurately modeled. In this paper, we present a numerical approach based on the dynamic poroelasticity for a space–time-domain staggered-grid finite-difference simulation in porous media that contain a free-surface boundary. We propose a generalized stess mirror formulation of the free-surface boundary for solids and fluids in porous media for the grid mesh on which lays the free-surface plane. Its analog is that used for elastic media, which is suitable for precise and stable Rayleigh-type surface-wave modeling. The results of our analysis of first kind of Rayleigh （R1） waves obtained by this model demonstrate that the discretization of the mesh in a similar way to that for elastic media can realize stable numerical solutions with acceptable precision. We present numerical examples demonstrating the efficiency and accuracy of our proposed method.
基金supported by the 863 Program (Grant no.2006AA09Z323)the 973 Program (Grant No.2006CB202402)
文摘Numerical simulation in the frequency-space domain has inherent advantages, such as： it is possible to simulate wave propagation from multiple sources simultaneously; there are no cumulative errors; only the interesting frequencies can be selected; and it is more suitable for wave propagation in viscoelastic media. The only obstacle to using the method is the requirement of huge computer storage. We extend the compressed format for storing the coefficient matrix. It can reduce the required computer storage dramatically. We get the optimal coefficients by least-squares method to suppress the numerical dispersion and adopt the perfectly matched layer （PML） boundary conditions to eliminate the artificial boundary reflections. Using larger grid intervals decreases computer storage requirements and provides high computational efficiency. Numerical experiments demonstrate that these means are economic and effective, providing a good basis for elastic wave imaging and inversion.
基金jointly supported by the NSF(No.41720104006)the Strategic Priority Research Program of the Chinese Academy of Sciences(A)(No.XDA14010303)+2 种基金the National Oil and Gas Project(Nos.2016ZX05002-005-007HZ and 2016ZX05014-001-008HZ)the Shandong Innovation Project(No.2017CXGC1602)the Qingdao Innovation Project(Nos.16-5-1-40-jch and 17CX05011)
文摘The staggered-grid finite-difference （SGFD） method has been widely used in seismic forward modeling. The precision of the forward modeling results directly affects the results of the subsequent seismic inversion and migration. Numerical dispersion is one of the problems in this method. The window function method can reduce dispersion by replacing the finite-difference operators with window operators, obtained by truncating the spatial convolution series of the pseudospectral method. Although the window operators have high precision in the low-wavenumber domain, their precision decreases rapidly in the high-wavenumber domain. We develop a least squares optimization method to enhance the precision of operators obtained by the window function method. We transform the SGFD problem into a least squares problem and find the best solution iteratively. The window operator is chosen as the initial value and the optimized domain is set by the error threshold. The conjugate gradient method is also adopted to increase the stability of the solution. Approximation error analysis and numerical simulation results suggest that the proposed method increases the precision of the window function operators and decreases the numerical dispersion.
基金supported by the Natural Science Foundation of China(No.41374078)Geological Survey Projects of Ministry of Land and Resources of China(No.12120113086100 and 12120113101300)
文摘The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the Maxwell＇s equations, and the magnetic components at the center of each edge of the grid cell are evaluated by applying the staggered-grid finite-difference method. The tipper and its divergence are derived to complete the 3D ZTEM forward modeling algorithm. A synthetic model is then used to compare the responses with those of 2D finite-element forward modeling to verify the accuracy of the algorithm. ZTEM offers high horizontal resolution to both simple and complex distributions of conductivity. This work is the theoretical foundation for the interpretation of ZTEM data and the study of 3D ZTEM inversion.
文摘Numerical simulation in transverse isotropic media with tilted symmetry axis（TTI） using the standard staggered-grid finite-difference scheme（SSG）results in errors caused by averaging or interpolation. In order to eliminate the errors, a method of rotated staggered-grid finite-difference scheme（RSG） is proposed. However, the RSG brings serious numerical dispersion. The compact staggered-grid finite-difference scheme（CSG） is an implicit difference scheme, which use fewer grid points to suppress dispersion more effectively than the SSG. This paper combines the CSG with the RSG to derive a rotated staggered-grid compact finite-difference scheme（RSGC）. The numerical experiments indicate that the RSGC has weaker numerical dispersion and better accuracy than the RSG.
基金National Natural Science Foundation of China (49894190-024) and Geophysical Prospecting Key Laboratory Foundation of China National Petroleum Corporation.
文摘To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.
