Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppr...Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppritz equations to estimate Young's modulus,which can introduce cumulative errors and reduce the accuracy of inversion results.To address these issues,this paper introduces the analytical solution of the Zoeppritz equation into the inversion process.The equation is re-derived and expressed in terms of Young's modulus,Poisson's ratio,and density.Within the Bayesian framework,we construct an objective function for the joint inversion of PP and PS waves.Traditional gradient-based algorithms often suffer from low precision and the computational complexity.In this study,we address limitations of conventional approaches related to low precision and complicated code by using Circle chaotic mapping,Levy flights,and Gaussian mutation to optimize the quantum particle swarm optimization(QPSO),named improved quantum particle swarm optimization(IQPSO).The IQPSO demonstrates superior global optimization capabilities.We test the proposed inversion method with both synthetic and field data.The test results demonstrate the proposed method's feasibility and effectiveness,indicating an improvement in inversion accuracy over traditional methods.展开更多
In Rayleigh wave exploration,the inversion of dispersion curves is a crucial step for obtaining subsurface stratigraphic information,characterized by its multi-parameter and multi-extremum nature.Local optimization al...In Rayleigh wave exploration,the inversion of dispersion curves is a crucial step for obtaining subsurface stratigraphic information,characterized by its multi-parameter and multi-extremum nature.Local optimization algorithms used in dispersion curve inversion are highly dependent on the initial model and are prone to being trapped in local optima,while classical global optimization algorithms often suffer from slow convergence and low solution accuracy.To address these issues,this study introduces the Osprey Optimization Algorithm(OOA),known for its strong global search and local exploitation capabilities,into the inversion of dispersion curves to enhance inversion performance.In noiseless theoretical models,the OOA demonstrates excellent inversion accuracy and stability,accurately recovering model parameters.Even in noisy models,OOA maintains robust performance,achieving high inversion precision under high-noise conditions.In multimode dispersion curve tests,OOA effectively handles higher modes due to its efficient global and local search capabilities,and the inversion results show high consistency with theoretical values.Field data from the Wyoming region in the United States and a landfill site in Italy further verify the practical applicability of the OOA.Comprehensive test results indicate that the OOA outperforms the Particle Swarm Optimization(PSO)algorithm,providing a highly accurate and reliable inversion strategy for dispersion curve inversion.展开更多
Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this cha...Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.展开更多
Deblending is a data processing procedure used to separate the source interferences of blended seismic data,which are obtained by simultaneous sources with random time delays to reduce the cost of seismic acquisition....Deblending is a data processing procedure used to separate the source interferences of blended seismic data,which are obtained by simultaneous sources with random time delays to reduce the cost of seismic acquisition.There are three types of deblending algorithms,i.e.,filtering-type noise suppression algorithm,inversion-based algorithm and deep-learning based algorithm.We review the merits of these techniques,and propose to use a sparse inversion method for seismic data deblending.Filtering-based deblending approach is applicable to blended data with a low blending fold and simple geometry.Otherwise,it can suffer from signal distortion and noise leakage.At present,the deep learning based deblending methods are still under development and field data applications are limited due to the lack of high-quality training labels.In contrast,the inversion-based deblending approaches have gained industrial acceptance.Our used inversion approach transforms the pseudo-deblended data into the frequency-wavenumber-wavenumher(FKK)domain,and a sparse constraint is imposed for the coherent signal estimation.The estimated signal is used to predict the interference noise for subtraction from the original pseudo-deblended data.Via minimizing the data misfit,the signal can be iteratively updated with a shrinking threshold until the signal and interference are fully separated.The used FKK sparse inversion algorithm is very accurate and efficient compared with other sparse inversion methods,and it is widely applied in field cases.Synthetic example shows that the deblending error is less than 1%in average amplitudes and less than-40 dB in amplitude spectra.We present three field data examples of land,marine OBN(Ocean Bottom Nodes)and streamer acquisitions to demonstrate its successful applications in separating the source interferences efficiently and accurately.展开更多
Predictions of fluid distribution,stress field,and natural fracture are essential for exploiting unconventional shale gas reservoirs.Given the high likelihood of tilted fractures in subsurface formations,this study fo...Predictions of fluid distribution,stress field,and natural fracture are essential for exploiting unconventional shale gas reservoirs.Given the high likelihood of tilted fractures in subsurface formations,this study focuses on simultaneous seismic inversion to estimate fluid bulk modulus,effective stress parameter,and fracture density in the tilted transversely isotropic(TTI)medium.In this article,a novel PP-wave reflection coefficient approximation equation is first derived based on the constructed TTI stiffness matrix incorporating fracture density,effective stress parameter,and fluid bulk modulus.The high accuracy of the proposed equation has been demonstrated using an anisotropic two-layer model.Furthermore,a stepwise seismic inversion strategy with the L_(P) quasi-norm sparsity constraint is implemented to obtain the anisotropic and isotropic parameters.Three synthetic model tests with varying signal-to-noise ratios(SNRs)confirm the method's feasibility and noise robustness.Ultimately,the proposed method is applied to a 3D fractured shale gas reservoir in the Sichuan Basin,China.The results have effectively characterized shale gas distribution,stress fields,and tilted natural fractures,with validation from geological structures,well logs,and microseismic events.These findings can provide valuable guidance for hydraulic fracturing development,enabling more reliable predictions of reservoir heterogeneity and completion quality.展开更多
In the 3D inversion modeling of gravity and magnetic potential field data,the model weighting function is often applied to overcome the skin eff ect of inversion results.However,divergence occurs at the the deep area,...In the 3D inversion modeling of gravity and magnetic potential field data,the model weighting function is often applied to overcome the skin eff ect of inversion results.However,divergence occurs at the the deep area,and artificial weak negative anomalies form around the positive anomalies in the horizontal direction,resulting in a reduction in the overall resolution.To fully utilize the model weighting function,this study constructs a combined model weighting function.First,a new depth weighting function is constructed by adding a regulator into the conventional depth weighting function to overcome the skin eff ect and inhibit the divergence at the deep area of the inversion results.A horizontal weighting function is then constructed by extracting information from the observation data;this function can suppress the formation of artificial weak anomalies and improve the horizontal resolution of the inversion results.Finally,these two functions are coupled to obtain the combined model weighting function,which can replace the conventional depth weighting function in 3D inversion.It improves the vertical and horizontal resolution of the inversion results without increasing the algorithm complexity and calculation amount,is easy to operate,and adapts to any 3D inversion method.Two model experiments are designed to verify the effectiveness,practicability,and anti-noise of the combined model weighting function.Then the function is applied to the 3D inversion of the measured aeromagnetic data in the Jinchuan area in China.The obtained inversion results are in good agreement with the known geological data.展开更多
The complex plate collision process led the South Yellow Sea Basin(SYSB)to go through an intensity tectonic inversion during the Early Cenozoic,leading to a regional unconformity surface development.As a petroliferous...The complex plate collision process led the South Yellow Sea Basin(SYSB)to go through an intensity tectonic inversion during the Early Cenozoic,leading to a regional unconformity surface development.As a petroliferous basin,SYSB saw intense denudation and deposition processes,making it hard to characterize their source-to-sink system(S2S),and this study provided a new way to reveal them quantitatively.According to the seismic interpretation,it was found that two types of tectonic inversion led to the strata shortening process,which was classified according to their difference in planar movements:dip-slip faults and strike-slip ones.As for dip-slip faults,the inversion structure was primarily formed by the dip-slip movement,and many fault-related folds developed,which developed in the North Depression Zone of the SYSB.The strike-slip ones,accompanied by some negative flower structures,dominate the South Depression Zone of the SYSB.To reveal its source-to-sink(S2S)system in the tectonic inversion basin,we rebuild the provenance area with detrital zircon U-Pb data and heavy mineral assemblage.The results show,during the Eocene(tectonic inversion stage),the proximal slump or fan delta from the Central Uplift Zone was prominently developed in the North Depression Zone,and the South Depression Zone is filled by sediments from the proximal area(Central Uplift Zone in SYSB and Wunansha Uplift)and the prograding delta long-axis parallel to the boundary faults.Then,calculations were conducted on the coarse sediment content,fault displacements,catchment relief,sediment migration distance,and discussions about the impact factors of the S2S system developed in various strata shortening patterns with a statistical method.It was found that,within the dip-slip faults-dominated zone,the volume of the sediment routing system and the ratio of coarse-grained sediments merely have a relationship with the amount of sediment supply and average faults break displacement.Compared with the strike-slip faults-dominated zone,the source-to-sink system shows a lower level of sandy sediment influx,and its coarse-grained content is mainly determined by the average faults broken displacement.展开更多
Machine learning(ML)techniques have emerged as powerful tools for improving the predictive capabilities of Reynolds-averaged Navier-Stokes(RANS)turbulence models in separated flows.This improvement is achieved by leve...Machine learning(ML)techniques have emerged as powerful tools for improving the predictive capabilities of Reynolds-averaged Navier-Stokes(RANS)turbulence models in separated flows.This improvement is achieved by leveraging complex ML models,such as those developed using field inversion and machine learning(FIML),to dynamically adjust the constants within the baseline RANS model.However,the ML models often overlook the fundamental calibrations of the RANS turbulence model.Consequently,the basic calibration of the baseline RANS model is disrupted,leading to a degradation in the accuracy,particularly in basic wall-attached flows outside of the training set.To address this issue,a modified version of the Spalart-Allmaras(SA)turbulence model,known as Rubber-band SA(RBSA),has been proposed recently.This modification involves identifying and embedding constraints related to basic wall-attached flows directly into the model.It is shown that no matter how the parameters of the RBSA model are adjusted as constants throughout the flow field,its accuracy in wall-attached flows remains unaffected.In this paper,we propose a new constraint for the RBSA model,which better safeguards the law of wall in extreme conditions where the model parameter is adjusted dramatically.The resultant model is called the RBSA-poly model.We then show that when combined with FIML augmentation,the RBSA-poly model effectively preserves the accuracy of simple wall-attached flows,even when the adjusted parameters become functions of local flow variables rather than constants.A comparative analysis with the FIML-augmented original SA model reveals that the augmented RBSA-poly model reduces error in basic wall-attached flows by 50%while maintaining comparable accuracy in trained separated flows.These findings confirm the effectiveness of utilizing FIML in conjunction with the RBSA model,offering superior accuracy retention in cardinal flows.展开更多
In this study,we employed Bayesian inversion coupled with the summation-by-parts and simultaneousapproximation-term(SBP-SAT)forward simulation method to elucidate the mechanisms behind mininginduced seismic events cau...In this study,we employed Bayesian inversion coupled with the summation-by-parts and simultaneousapproximation-term(SBP-SAT)forward simulation method to elucidate the mechanisms behind mininginduced seismic events caused by fault slip and their potential effects on rockbursts.Through Bayesian inversion,it is determined that the sources near fault FQ14 have a significant shear component.Additionally,we analyzed the stress and displacement fields of high-energy events,along with the hypocenter distribution of aftershocks,which aided in identifying the slip direction of the critically stressed fault FQ14.We also performed forward modeling to capture the complex dynamics of fault slip under varying friction laws and shear fracture modes.The selection of specific friction laws for fault slip models was based on their ability to accurately replicate observed slip behavior under various external loading conditions,thereby enhancing the applicability of our findings.Our results suggest that the slip behavior of fault FQ14 can be effectively understood by comparing different scenarios.展开更多
Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods rema...Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods remains challenging,as it requires the estimation of more than eight parameters.Assuming the reservoir to be a weakly anisotropic ORT medium with small contrasts in the background elastic parameters,a new azimuthal elastic impedance equation was first derived using parameter combinations and mathematical approximations.This equation exhibited almost the same accuracy as the original equation and contained only six model parameters:the compression modulus,anisotropic shear modulus,anisotropic compression modulus,density,normal fracture weakness,and tangential fracture weakness.Subsequently,a stepwise inversion method using second-order derivatives of the elastic impedance was developed to estimate these parameters.Moreover,the Thomsen anisotropy parameter,epsilon,was estimated from the inversion results using the ratio of the anisotropic compression modulus to the compression modulus.Synthetic examples with moderate noise and field data examples confirm the feasibility and effectiveness of the inversion method.The proposed method exhibited accuracy similar to that of previous inversion strategies and could predict richer vertical fracture information.Ultimately,the method was applied to a three-dimensional work area,and the predictions were consistent with logging and geological a priori information,confirming the effectiveness of this method.Summarily,the proposed stepwise inversion method can alleviate the uncertainty of multi-parameter inversion in ORT medium,thereby improving the reliability of fracture detection.展开更多
An innovative gradient inversion approach employing the natural element method within the framework of least square regularization was proposed to enhance the quantitative interpretation of self-potential(SP)data orig...An innovative gradient inversion approach employing the natural element method within the framework of least square regularization was proposed to enhance the quantitative interpretation of self-potential(SP)data originating from mineral polarization.The results indicated that the natural element method effectively addressed the challenge of subdividing complex resistivity models and aided in the accurate forward calculation of SP.By applying this approach to synthetic SP data and lab-measured SP data associated with redox electrochemical half-cell reactions of iron−copper metal blocks within the geobattery model,the 3D fine structure of buried orebody models was successfully reconstructed and the spatial distribution of SP current sources was mapped.This study significantly contributes to understanding the quantitative relationship between the polarization process of metal deposits and their corresponding SP responses and provides a valuable reference for delineating metal deposits in both terrestrial and marine environments through SP surveys.展开更多
Achieving high-resolution intracranial imaging in a safe and portable manner is critical for the diagnosis of intracranial diseases,preoperative planning of craniotomies and intraoperative management during craniotomy...Achieving high-resolution intracranial imaging in a safe and portable manner is critical for the diagnosis of intracranial diseases,preoperative planning of craniotomies and intraoperative management during craniotomy procedures.Adaptive waveform inversion(AWI),a variant of full waveform inversion(FWI),has shown potential in intracranial ultrasound imaging.However,the robustness of AWI is affected by the parameterization of the Gaussian penalty matrix and the challenges posed by transcranial scenarios.Conventional AWI struggles to produce accurate images in these cases,limiting its application in critical medical settings.To address these issues,we propose a stabilized adaptive waveform inversion(SAWI)method,which introduces a user-defined zero-lag position for theWiener filter.Numerical experiments demonstrate that SAWI can achieve accurate imaging under Gaussian penalty matrix parameter settings where AWI fails,perform successful transcranial imaging in configurations where AWI cannot,and maintain the same imaging accuracy as AWI.The advantage of this method is that it achieves these advancements without modifying the AWI framework or increasing computational costs,which helps to promote the application of AWI in medical fields,particularly in transcranial scenarios.展开更多
BACKGROUND Back pain and sciatica are common complaints that often require imaging for accurate diagnosis and management.Conventional lumbar magnetic resonance imaging(MRI)protocols typically include sagittal and axia...BACKGROUND Back pain and sciatica are common complaints that often require imaging for accurate diagnosis and management.Conventional lumbar magnetic resonance imaging(MRI)protocols typically include sagittal and axial T1 and T2 sequences;however,these may miss certain pathologies.The addition of coronal short tau inversion recovery(STIR)sequences offers the potential to enhance the detection of both spinal and extra-spinal abnormalities,thereby improving clinical decisionmaking and patient outcomes.AIM To evaluate the impact of adding coronal STIR sequences to routine lumbar MRI in diagnosing back pain and sciatica.METHODS We prospectively analyzed data from patients aged 6 and older presenting with back pain or sciatica who underwent lumbar spine MRI at our institution.The standardized MRI protocol utilized included sagittal and axial T1 and T2 sequences,complemented by a coronal STIR sequence.Data on structural abnormalities were collected,reviewed,and analyzed using counts,percentages,and Fisher's exact test for categorical variables.RESULTS Our cohort comprised 274 patients(115 males,159 females;mean age 44.91 years).Notably,39 patients exhibited abnormalities across all sequences,while 72.63%showed normal findings on the coronal STIR sequence.Importantly,30.29%of cases were diagnosed as normal without the coronal STIR,and 36 patients with normal T1 and T2 sequences presented abnormalities on the coronal STIR.The coronal STIR sequence successfully identified 26 spinal and 10 non-spinal pathologies,including 17 cases of sacroiliitis,with a significant association(P<0.0001)between sacroiliitis diagnosis and abnormalities visible solely on this sequence.CONCLUSION Integrating coronal STIR into routine lumbar MRI enhances detection of hidden spinal and extra-spinal pathologies,improves patient management,and offers a cost-effective,practical upgrade with significant diagnostic and clinical value.展开更多
Previous studies have often focused on monitoring grassland growth as the primary target of remote sensing investigations on grassland ecological restoration in the northern Tibetan Plateau,overlooking the crucial rol...Previous studies have often focused on monitoring grassland growth as the primary target of remote sensing investigations on grassland ecological restoration in the northern Tibetan Plateau,overlooking the crucial role played by gravel in the ecological restoration of these grasslands.This study utilizes supervised classification and segmentation techniques based on machine learning to extract gravel morphology profiles from field-sampled plot images and calculate their characteristic parameters.Employing a multivariate linear approach combined with Principal Component Analysis(PCA),a model for inferring gravel characteristic parameters is constructed.Statistical features,particle size characteristics,and spatial distribution patterns of gravel are analyzed.Results reveal that gravel predominantly exhibit sub-rounded shapes,with 80%classified as fine gravel.The coefficients of determination(R2)between gravel particle size and coverage,perimeter,and area are 0.444,0.724,and 0.557,respectively,indicating linear relationships.The cumulative contribution rate of the top five remote sensing factors is 95.44%,with the first geological factor contributing 77.64%,collectively reflecting the primary information of the 20 factors used.Modeling shows that areas with larger gravel particle sizes correspond to increased perimeter and coverage.Gravels in the Nagqu Prefecture of northern Xizang have a particle size range of 4-8 mm,primarily comprising fine gravel which accounts for 94.61%.These findings provide a scientific basis for extracting gravel characteristic parameters and understanding their spatial distribution variations in the northern Tibetan Plateau.展开更多
The brittleness index(BI)is crucial for predicting engineering sweet spots and designing fracturing operations in shale oil reservoir exploration and development.Seismic amplitude variation with offset(AVO)inversion i...The brittleness index(BI)is crucial for predicting engineering sweet spots and designing fracturing operations in shale oil reservoir exploration and development.Seismic amplitude variation with offset(AVO)inversion is commonly used to obtain the BI.Traditionally,velocity,density,and other parameters are firstly inverted,and the BI is then calculated,which often leads to accumulated errors.Moreover,due to the limited of well-log data in field work areas,AVO inversion typically faces the challenge of limited information,resulting in not high accuracy of BI derived by existing AVO inversion methods.To address these issues,we first derive an AVO forward approximation equation that directly characterizes the BI in P-wave reflection coefficients.Based on this,an intelligent AVO inversion method,which combines the advantages of traditional and intelligent approaches,for directly obtaining the BI is proposed.A TransUnet model is constructed to establish the strong nonlinear mapping relationship between seismic data and the BI.By incorporating a combined objective function that is constrained by both low-frequency parameters and training samples,the challenge of limited samples is effectively addressed,and the direct inversion of the BI is stably achieved.Tests on model data and applications on field data demonstrate the feasibility,advancement,and practicality of the proposed method.展开更多
Full waveform inversion is a precise method for parameter inversion,harnessing the complete wavefield information of seismic waves.It holds the potential to intricately characterize the detailed features of the model ...Full waveform inversion is a precise method for parameter inversion,harnessing the complete wavefield information of seismic waves.It holds the potential to intricately characterize the detailed features of the model with high accuracy.However,due to inaccurate initial models,the absence of low-frequency data,and incomplete observational data,full waveform inversion(FWI)exhibits pronounced nonlinear characteristics.When the strata are buried deep,the inversion capability of this method is constrained.To enhance the accuracy and precision of FWI,this paper introduces a novel approach to address the aforementioned challenges—namely,a fractional-order anisotropic total p-variation regularization for full waveform inversion(FATpV-FWI).This method incorporates fractional-order total variation(TV)regularization to construct the inversion objective function,building upon TV regularization,and subsequently employs the alternating direction multiplier method for solving.This approach mitigates the step effect stemming from total variation in seismic inversion,thereby facilitating the reconstruction of sharp interfaces of geophysical parameters while smoothing background variations.Simultaneously,replacing integer-order differences with fractional-order differences bolsters the correlation among seismic data and diminishes the scattering effect caused by integer-order differences in seismic inversion.The outcomes of model tests validate the efficacy of this method,highlighting its ability to enhance the overall accuracy of the inversion process.展开更多
Based on waveform fitting,full waveform inversion(FWI)is an important inversion method with the ability to reconstruct multi-parameter models in high precision.However,the strong nonlinear equation used in FWI present...Based on waveform fitting,full waveform inversion(FWI)is an important inversion method with the ability to reconstruct multi-parameter models in high precision.However,the strong nonlinear equation used in FWI presents the following challenges,such as low convergence efficiency,high dependence on the initial model,and the energy imbalance in deep region of the inverted model.To solve these inherent problems,we develop a timedomain elastic FWI method based on gradient preconditioning with the following details:(1)the limited memory Broyden Fletcher Goldfarb Shanno method with faster convergence is adopted to im-prove the inversion stability;(2)a multi-scaled inversion strategy is used to alleviate the nonlinear inversion instead of falling into the local minimum;(3)in addition,the pseudo-Hessian preconditioned illumination operator is involved for preconditioning the parameter gradients to improve the illumination equilibrium degree of deep structures.Based on the programming implementation of the new method,a deep depression model with five diffractors is used for testing.Compared with the conventional elastic FWI method,the technique proposed by this study has better effectiveness and accuracy on the inversion effect and con-vergence,respectively.展开更多
Poly(m-phenylene isophthalamide)(PMIA),a key aromatic polyamide,is widely used for its outstanding mechanical strength,high thermal stability,and excellent insulation properties.However,different applications demand v...Poly(m-phenylene isophthalamide)(PMIA),a key aromatic polyamide,is widely used for its outstanding mechanical strength,high thermal stability,and excellent insulation properties.However,different applications demand varying dielectric properties,so tailoring its dielectric per-formance is essential.PMIA was first synthesized in this study,followed by introducing pores and developing porous PMIA films and PMIA-based composites with reduced dielectric constants.Porous PMIA films were fabricated using the wet phase inversion process with N,N-dimethylac-etamide(DMAC)solvent and water as the non-solvent.The impact of casting solution composi-tion and coagulation bath temperature on pore structures was analyzed.A film produced with 18%PMIA and 5%LiCl in a 35℃coagulation bath achieved the lowest dielectric constant of 1.76 at 1 Hz,48%lower than the standard PMIA film,which had a tensile strength of 18.5 MPa and an initial degradation temperature of 320℃.展开更多
Traditional deep learning methods pursue complex and single network architectures without considering the petrophysical relationship between different elastic parameters.The mathematical and statistical significance o...Traditional deep learning methods pursue complex and single network architectures without considering the petrophysical relationship between different elastic parameters.The mathematical and statistical significance of the inversion results may lead to model overfitting,especially when there are a limited number of well logs in a working area.Multitask learning provides an eff ective approach to addressing this issue.Simultaneously,learning multiple related tasks can improve a model’s generalization ability to a certain extent,thereby enhancing the performance of related tasks with an equal amount of labeled data.In this study,we propose an end-to-end multitask deep learning model that integrates a fully convolutional network and bidirectional gated recurrent unit for intelligent prestack inversion of“seismic data to elastic parameters.”The use of a Bayesian homoscedastic uncertainty-based loss function enables adaptive learning of the weight coeffi cients for diff erent elastic parameter inversion tasks,thereby reducing uncertainty during the inversion process.The proposed method combines the local feature perception of convolutional neural networks with the long-term memory of bidirectional gated recurrent networks.It maintains the rock physics constraint relationships among diff erent elastic parameters during the inversion process,demonstrating a high level of prediction accuracy.Numerical simulations and processing results of real seismic data validate the eff ectiveness and practicality of the proposed method.展开更多
Full waveform inversion methods evaluate the properties of subsurface media by minimizing the misfit between synthetic and observed data.However,these methods omit measurement errors and physical assumptions in modeli...Full waveform inversion methods evaluate the properties of subsurface media by minimizing the misfit between synthetic and observed data.However,these methods omit measurement errors and physical assumptions in modeling,resulting in several problems in practical applications.In particular,full waveform inversion methods are very sensitive to erroneous observations(outliers)that violate the Gauss–Markov theorem.Herein,we propose a method for addressing spurious observations or outliers.Specifically,we remove outliers by inverting the synthetic data using the local convexity of the Gaussian distribution.To achieve this,we apply a waveform-like noise model based on a specific covariance matrix definition.Finally,we build an inversion problem based on the updated data,which is consistent with the wavefield reconstruction inversion method.Overall,we report an alternative optimization inversion problem for data containing outliers.The proposed method is robust because it uses uncertainties.This method enables accurate inversion,even when based on noisy models or a wrong wavelet.展开更多
基金supported by Fundamental Research Funds for the Central Universities,CHD300102264715National Key Research and Development Program of China under Grant 2021YFA0716902Natural Science Basic Research Program of Shaanxi 2024JCYBMS-199。
文摘Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppritz equations to estimate Young's modulus,which can introduce cumulative errors and reduce the accuracy of inversion results.To address these issues,this paper introduces the analytical solution of the Zoeppritz equation into the inversion process.The equation is re-derived and expressed in terms of Young's modulus,Poisson's ratio,and density.Within the Bayesian framework,we construct an objective function for the joint inversion of PP and PS waves.Traditional gradient-based algorithms often suffer from low precision and the computational complexity.In this study,we address limitations of conventional approaches related to low precision and complicated code by using Circle chaotic mapping,Levy flights,and Gaussian mutation to optimize the quantum particle swarm optimization(QPSO),named improved quantum particle swarm optimization(IQPSO).The IQPSO demonstrates superior global optimization capabilities.We test the proposed inversion method with both synthetic and field data.The test results demonstrate the proposed method's feasibility and effectiveness,indicating an improvement in inversion accuracy over traditional methods.
基金sponsored by China Geological Survey Project(DD20243193 and DD20230206508).
文摘In Rayleigh wave exploration,the inversion of dispersion curves is a crucial step for obtaining subsurface stratigraphic information,characterized by its multi-parameter and multi-extremum nature.Local optimization algorithms used in dispersion curve inversion are highly dependent on the initial model and are prone to being trapped in local optima,while classical global optimization algorithms often suffer from slow convergence and low solution accuracy.To address these issues,this study introduces the Osprey Optimization Algorithm(OOA),known for its strong global search and local exploitation capabilities,into the inversion of dispersion curves to enhance inversion performance.In noiseless theoretical models,the OOA demonstrates excellent inversion accuracy and stability,accurately recovering model parameters.Even in noisy models,OOA maintains robust performance,achieving high inversion precision under high-noise conditions.In multimode dispersion curve tests,OOA effectively handles higher modes due to its efficient global and local search capabilities,and the inversion results show high consistency with theoretical values.Field data from the Wyoming region in the United States and a landfill site in Italy further verify the practical applicability of the OOA.Comprehensive test results indicate that the OOA outperforms the Particle Swarm Optimization(PSO)algorithm,providing a highly accurate and reliable inversion strategy for dispersion curve inversion.
基金funded by the National Key R&D Program of China(Grant No.2022YFC2903904)the National Natural Science Foundation of China(Grant Nos.51904057 and U1906208).
文摘Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.
基金supported by National Science and Technology Major Project(Grant No.2017ZX05018-001)。
文摘Deblending is a data processing procedure used to separate the source interferences of blended seismic data,which are obtained by simultaneous sources with random time delays to reduce the cost of seismic acquisition.There are three types of deblending algorithms,i.e.,filtering-type noise suppression algorithm,inversion-based algorithm and deep-learning based algorithm.We review the merits of these techniques,and propose to use a sparse inversion method for seismic data deblending.Filtering-based deblending approach is applicable to blended data with a low blending fold and simple geometry.Otherwise,it can suffer from signal distortion and noise leakage.At present,the deep learning based deblending methods are still under development and field data applications are limited due to the lack of high-quality training labels.In contrast,the inversion-based deblending approaches have gained industrial acceptance.Our used inversion approach transforms the pseudo-deblended data into the frequency-wavenumber-wavenumher(FKK)domain,and a sparse constraint is imposed for the coherent signal estimation.The estimated signal is used to predict the interference noise for subtraction from the original pseudo-deblended data.Via minimizing the data misfit,the signal can be iteratively updated with a shrinking threshold until the signal and interference are fully separated.The used FKK sparse inversion algorithm is very accurate and efficient compared with other sparse inversion methods,and it is widely applied in field cases.Synthetic example shows that the deblending error is less than 1%in average amplitudes and less than-40 dB in amplitude spectra.We present three field data examples of land,marine OBN(Ocean Bottom Nodes)and streamer acquisitions to demonstrate its successful applications in separating the source interferences efficiently and accurately.
基金financially supported by the Natural Science Foundation of Sichuan Province(Grant Nos.2023NSFSC0767 and2024NSFSC0809)the China Postdoctoral Science Foundation(Grant No.2024MF750281)the Postdoctoral Fellowship Program of CPSF(Grant No.GZC20230326)。
文摘Predictions of fluid distribution,stress field,and natural fracture are essential for exploiting unconventional shale gas reservoirs.Given the high likelihood of tilted fractures in subsurface formations,this study focuses on simultaneous seismic inversion to estimate fluid bulk modulus,effective stress parameter,and fracture density in the tilted transversely isotropic(TTI)medium.In this article,a novel PP-wave reflection coefficient approximation equation is first derived based on the constructed TTI stiffness matrix incorporating fracture density,effective stress parameter,and fluid bulk modulus.The high accuracy of the proposed equation has been demonstrated using an anisotropic two-layer model.Furthermore,a stepwise seismic inversion strategy with the L_(P) quasi-norm sparsity constraint is implemented to obtain the anisotropic and isotropic parameters.Three synthetic model tests with varying signal-to-noise ratios(SNRs)confirm the method's feasibility and noise robustness.Ultimately,the proposed method is applied to a 3D fractured shale gas reservoir in the Sichuan Basin,China.The results have effectively characterized shale gas distribution,stress fields,and tilted natural fractures,with validation from geological structures,well logs,and microseismic events.These findings can provide valuable guidance for hydraulic fracturing development,enabling more reliable predictions of reservoir heterogeneity and completion quality.
基金jointly funded by the National Natural Science Foundation of China(No.U2244220,No.42004125)the China Geological Survey Projects(No.DD20240119,No.DD20243245,No.DD20230114,No.DD20243244)the China Postdoctoral Science Foundation(No.2020M670601)。
文摘In the 3D inversion modeling of gravity and magnetic potential field data,the model weighting function is often applied to overcome the skin eff ect of inversion results.However,divergence occurs at the the deep area,and artificial weak negative anomalies form around the positive anomalies in the horizontal direction,resulting in a reduction in the overall resolution.To fully utilize the model weighting function,this study constructs a combined model weighting function.First,a new depth weighting function is constructed by adding a regulator into the conventional depth weighting function to overcome the skin eff ect and inhibit the divergence at the deep area of the inversion results.A horizontal weighting function is then constructed by extracting information from the observation data;this function can suppress the formation of artificial weak anomalies and improve the horizontal resolution of the inversion results.Finally,these two functions are coupled to obtain the combined model weighting function,which can replace the conventional depth weighting function in 3D inversion.It improves the vertical and horizontal resolution of the inversion results without increasing the algorithm complexity and calculation amount,is easy to operate,and adapts to any 3D inversion method.Two model experiments are designed to verify the effectiveness,practicability,and anti-noise of the combined model weighting function.Then the function is applied to the 3D inversion of the measured aeromagnetic data in the Jinchuan area in China.The obtained inversion results are in good agreement with the known geological data.
基金sponsored by the National Natural Science Foundation of China-Youth Science Fund(No.42402150)the Major State Science and Technology Research Program(No.2016ZX05024002-002)the Chinese Scholarship Council(CSC)。
文摘The complex plate collision process led the South Yellow Sea Basin(SYSB)to go through an intensity tectonic inversion during the Early Cenozoic,leading to a regional unconformity surface development.As a petroliferous basin,SYSB saw intense denudation and deposition processes,making it hard to characterize their source-to-sink system(S2S),and this study provided a new way to reveal them quantitatively.According to the seismic interpretation,it was found that two types of tectonic inversion led to the strata shortening process,which was classified according to their difference in planar movements:dip-slip faults and strike-slip ones.As for dip-slip faults,the inversion structure was primarily formed by the dip-slip movement,and many fault-related folds developed,which developed in the North Depression Zone of the SYSB.The strike-slip ones,accompanied by some negative flower structures,dominate the South Depression Zone of the SYSB.To reveal its source-to-sink(S2S)system in the tectonic inversion basin,we rebuild the provenance area with detrital zircon U-Pb data and heavy mineral assemblage.The results show,during the Eocene(tectonic inversion stage),the proximal slump or fan delta from the Central Uplift Zone was prominently developed in the North Depression Zone,and the South Depression Zone is filled by sediments from the proximal area(Central Uplift Zone in SYSB and Wunansha Uplift)and the prograding delta long-axis parallel to the boundary faults.Then,calculations were conducted on the coarse sediment content,fault displacements,catchment relief,sediment migration distance,and discussions about the impact factors of the S2S system developed in various strata shortening patterns with a statistical method.It was found that,within the dip-slip faults-dominated zone,the volume of the sediment routing system and the ratio of coarse-grained sediments merely have a relationship with the amount of sediment supply and average faults break displacement.Compared with the strike-slip faults-dominated zone,the source-to-sink system shows a lower level of sandy sediment influx,and its coarse-grained content is mainly determined by the average faults broken displacement.
基金supported by the National Natural Science Foundation of China(Grant Nos.12388101,12372288,U23A2069,and 92152301).
文摘Machine learning(ML)techniques have emerged as powerful tools for improving the predictive capabilities of Reynolds-averaged Navier-Stokes(RANS)turbulence models in separated flows.This improvement is achieved by leveraging complex ML models,such as those developed using field inversion and machine learning(FIML),to dynamically adjust the constants within the baseline RANS model.However,the ML models often overlook the fundamental calibrations of the RANS turbulence model.Consequently,the basic calibration of the baseline RANS model is disrupted,leading to a degradation in the accuracy,particularly in basic wall-attached flows outside of the training set.To address this issue,a modified version of the Spalart-Allmaras(SA)turbulence model,known as Rubber-band SA(RBSA),has been proposed recently.This modification involves identifying and embedding constraints related to basic wall-attached flows directly into the model.It is shown that no matter how the parameters of the RBSA model are adjusted as constants throughout the flow field,its accuracy in wall-attached flows remains unaffected.In this paper,we propose a new constraint for the RBSA model,which better safeguards the law of wall in extreme conditions where the model parameter is adjusted dramatically.The resultant model is called the RBSA-poly model.We then show that when combined with FIML augmentation,the RBSA-poly model effectively preserves the accuracy of simple wall-attached flows,even when the adjusted parameters become functions of local flow variables rather than constants.A comparative analysis with the FIML-augmented original SA model reveals that the augmented RBSA-poly model reduces error in basic wall-attached flows by 50%while maintaining comparable accuracy in trained separated flows.These findings confirm the effectiveness of utilizing FIML in conjunction with the RBSA model,offering superior accuracy retention in cardinal flows.
基金the Graduate Innovation Program of China University of Mining and Technology,the Fundamental Research Funds for the Central Universities(Grant No.2023WLKXJ017)the Postgraduate Research&Practice Innovation Program of Jiangsu Province(Grant No.KYCX23_2776)the Shandong Energy Group(Grant No.SNKJ2022BJ03-R28)。
文摘In this study,we employed Bayesian inversion coupled with the summation-by-parts and simultaneousapproximation-term(SBP-SAT)forward simulation method to elucidate the mechanisms behind mininginduced seismic events caused by fault slip and their potential effects on rockbursts.Through Bayesian inversion,it is determined that the sources near fault FQ14 have a significant shear component.Additionally,we analyzed the stress and displacement fields of high-energy events,along with the hypocenter distribution of aftershocks,which aided in identifying the slip direction of the critically stressed fault FQ14.We also performed forward modeling to capture the complex dynamics of fault slip under varying friction laws and shear fracture modes.The selection of specific friction laws for fault slip models was based on their ability to accurately replicate observed slip behavior under various external loading conditions,thereby enhancing the applicability of our findings.Our results suggest that the slip behavior of fault FQ14 can be effectively understood by comparing different scenarios.
基金sponsorship of the National Natural Science Foundation of China(42430809,42274157,42030103,42404132)the Fund of State Key Laboratory of Deep Oil and Gas,China University of Petroleum(East China)(SKLDOG2024-ZYTS-02)+5 种基金the Postdoctoral Fellowship Program of CPSF(GZB20240850)the Postdoctoral Project of Qingdao(QDBSH20240102082)the Fundamental Research Funds for the Central Universities(24CX07004A,24CX06036A)the CNPC Innovation Fund(2024DQ02-0505,2024DQ02-0136)the Innovation fund project for graduate student of China University of Petroleum(East China)the Fundamental Research Funds for the Central Universities(24CX04002A).
文摘Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods remains challenging,as it requires the estimation of more than eight parameters.Assuming the reservoir to be a weakly anisotropic ORT medium with small contrasts in the background elastic parameters,a new azimuthal elastic impedance equation was first derived using parameter combinations and mathematical approximations.This equation exhibited almost the same accuracy as the original equation and contained only six model parameters:the compression modulus,anisotropic shear modulus,anisotropic compression modulus,density,normal fracture weakness,and tangential fracture weakness.Subsequently,a stepwise inversion method using second-order derivatives of the elastic impedance was developed to estimate these parameters.Moreover,the Thomsen anisotropy parameter,epsilon,was estimated from the inversion results using the ratio of the anisotropic compression modulus to the compression modulus.Synthetic examples with moderate noise and field data examples confirm the feasibility and effectiveness of the inversion method.The proposed method exhibited accuracy similar to that of previous inversion strategies and could predict richer vertical fracture information.Ultimately,the method was applied to a three-dimensional work area,and the predictions were consistent with logging and geological a priori information,confirming the effectiveness of this method.Summarily,the proposed stepwise inversion method can alleviate the uncertainty of multi-parameter inversion in ORT medium,thereby improving the reliability of fracture detection.
基金supported by the National Natural Science Foundation of China(No.42174170)。
文摘An innovative gradient inversion approach employing the natural element method within the framework of least square regularization was proposed to enhance the quantitative interpretation of self-potential(SP)data originating from mineral polarization.The results indicated that the natural element method effectively addressed the challenge of subdividing complex resistivity models and aided in the accurate forward calculation of SP.By applying this approach to synthetic SP data and lab-measured SP data associated with redox electrochemical half-cell reactions of iron−copper metal blocks within the geobattery model,the 3D fine structure of buried orebody models was successfully reconstructed and the spatial distribution of SP current sources was mapped.This study significantly contributes to understanding the quantitative relationship between the polarization process of metal deposits and their corresponding SP responses and provides a valuable reference for delineating metal deposits in both terrestrial and marine environments through SP surveys.
基金supported by the National Natural Science Foundation of China(Grant No.82151302)the National High Level Hospital Clinical Research Funding(Grant No.2022-PUMCH-B-113)+1 种基金the National High Level Hospital Clinical Research Funding(Grant No.2022-PUMCH-A-019)the CAMS Innovation Fund for Medical Sciences(Grant No.2021-12M-1-014).
文摘Achieving high-resolution intracranial imaging in a safe and portable manner is critical for the diagnosis of intracranial diseases,preoperative planning of craniotomies and intraoperative management during craniotomy procedures.Adaptive waveform inversion(AWI),a variant of full waveform inversion(FWI),has shown potential in intracranial ultrasound imaging.However,the robustness of AWI is affected by the parameterization of the Gaussian penalty matrix and the challenges posed by transcranial scenarios.Conventional AWI struggles to produce accurate images in these cases,limiting its application in critical medical settings.To address these issues,we propose a stabilized adaptive waveform inversion(SAWI)method,which introduces a user-defined zero-lag position for theWiener filter.Numerical experiments demonstrate that SAWI can achieve accurate imaging under Gaussian penalty matrix parameter settings where AWI fails,perform successful transcranial imaging in configurations where AWI cannot,and maintain the same imaging accuracy as AWI.The advantage of this method is that it achieves these advancements without modifying the AWI framework or increasing computational costs,which helps to promote the application of AWI in medical fields,particularly in transcranial scenarios.
文摘BACKGROUND Back pain and sciatica are common complaints that often require imaging for accurate diagnosis and management.Conventional lumbar magnetic resonance imaging(MRI)protocols typically include sagittal and axial T1 and T2 sequences;however,these may miss certain pathologies.The addition of coronal short tau inversion recovery(STIR)sequences offers the potential to enhance the detection of both spinal and extra-spinal abnormalities,thereby improving clinical decisionmaking and patient outcomes.AIM To evaluate the impact of adding coronal STIR sequences to routine lumbar MRI in diagnosing back pain and sciatica.METHODS We prospectively analyzed data from patients aged 6 and older presenting with back pain or sciatica who underwent lumbar spine MRI at our institution.The standardized MRI protocol utilized included sagittal and axial T1 and T2 sequences,complemented by a coronal STIR sequence.Data on structural abnormalities were collected,reviewed,and analyzed using counts,percentages,and Fisher's exact test for categorical variables.RESULTS Our cohort comprised 274 patients(115 males,159 females;mean age 44.91 years).Notably,39 patients exhibited abnormalities across all sequences,while 72.63%showed normal findings on the coronal STIR sequence.Importantly,30.29%of cases were diagnosed as normal without the coronal STIR,and 36 patients with normal T1 and T2 sequences presented abnormalities on the coronal STIR.The coronal STIR sequence successfully identified 26 spinal and 10 non-spinal pathologies,including 17 cases of sacroiliitis,with a significant association(P<0.0001)between sacroiliitis diagnosis and abnormalities visible solely on this sequence.CONCLUSION Integrating coronal STIR into routine lumbar MRI enhances detection of hidden spinal and extra-spinal pathologies,improves patient management,and offers a cost-effective,practical upgrade with significant diagnostic and clinical value.
基金funded by the Major R&D and Achievement Transformation Projects of Xizang(CGZH2024000416)Science and Technology Program of Xizang(XZ202402ZD0001)Major R&D and Achievement Transformation Projects of Qinghai(2022-QY-224)。
文摘Previous studies have often focused on monitoring grassland growth as the primary target of remote sensing investigations on grassland ecological restoration in the northern Tibetan Plateau,overlooking the crucial role played by gravel in the ecological restoration of these grasslands.This study utilizes supervised classification and segmentation techniques based on machine learning to extract gravel morphology profiles from field-sampled plot images and calculate their characteristic parameters.Employing a multivariate linear approach combined with Principal Component Analysis(PCA),a model for inferring gravel characteristic parameters is constructed.Statistical features,particle size characteristics,and spatial distribution patterns of gravel are analyzed.Results reveal that gravel predominantly exhibit sub-rounded shapes,with 80%classified as fine gravel.The coefficients of determination(R2)between gravel particle size and coverage,perimeter,and area are 0.444,0.724,and 0.557,respectively,indicating linear relationships.The cumulative contribution rate of the top five remote sensing factors is 95.44%,with the first geological factor contributing 77.64%,collectively reflecting the primary information of the 20 factors used.Modeling shows that areas with larger gravel particle sizes correspond to increased perimeter and coverage.Gravels in the Nagqu Prefecture of northern Xizang have a particle size range of 4-8 mm,primarily comprising fine gravel which accounts for 94.61%.These findings provide a scientific basis for extracting gravel characteristic parameters and understanding their spatial distribution variations in the northern Tibetan Plateau.
基金supposed by the National Nature Science Foundation of China(Grant No.42304131)the Natural Science Foundation of Heilongjiang Province(Grant No.LH2023D012)+1 种基金the Heilongjiang Postdoctoral Fund(Grant No.LBH-Z22092)the Basic Research Fund for Universities in Xinjiang Uygur Autonomous Region(Grant No.XJEDU2023P166)。
文摘The brittleness index(BI)is crucial for predicting engineering sweet spots and designing fracturing operations in shale oil reservoir exploration and development.Seismic amplitude variation with offset(AVO)inversion is commonly used to obtain the BI.Traditionally,velocity,density,and other parameters are firstly inverted,and the BI is then calculated,which often leads to accumulated errors.Moreover,due to the limited of well-log data in field work areas,AVO inversion typically faces the challenge of limited information,resulting in not high accuracy of BI derived by existing AVO inversion methods.To address these issues,we first derive an AVO forward approximation equation that directly characterizes the BI in P-wave reflection coefficients.Based on this,an intelligent AVO inversion method,which combines the advantages of traditional and intelligent approaches,for directly obtaining the BI is proposed.A TransUnet model is constructed to establish the strong nonlinear mapping relationship between seismic data and the BI.By incorporating a combined objective function that is constrained by both low-frequency parameters and training samples,the challenge of limited samples is effectively addressed,and the direct inversion of the BI is stably achieved.Tests on model data and applications on field data demonstrate the feasibility,advancement,and practicality of the proposed method.
基金supported by the China Postdoctoral Science Foundation(Grant No.2024MF750281)the Postdoctoral Fellowship Program of CPSF(Grant No.GZC20230326)+1 种基金the Natural Science Foundation Project of Sichuan Province(Grant No.2025ZNSFSC1170)Sichuan Science and Technology Program(Grant No.2023ZYD0158).
文摘Full waveform inversion is a precise method for parameter inversion,harnessing the complete wavefield information of seismic waves.It holds the potential to intricately characterize the detailed features of the model with high accuracy.However,due to inaccurate initial models,the absence of low-frequency data,and incomplete observational data,full waveform inversion(FWI)exhibits pronounced nonlinear characteristics.When the strata are buried deep,the inversion capability of this method is constrained.To enhance the accuracy and precision of FWI,this paper introduces a novel approach to address the aforementioned challenges—namely,a fractional-order anisotropic total p-variation regularization for full waveform inversion(FATpV-FWI).This method incorporates fractional-order total variation(TV)regularization to construct the inversion objective function,building upon TV regularization,and subsequently employs the alternating direction multiplier method for solving.This approach mitigates the step effect stemming from total variation in seismic inversion,thereby facilitating the reconstruction of sharp interfaces of geophysical parameters while smoothing background variations.Simultaneously,replacing integer-order differences with fractional-order differences bolsters the correlation among seismic data and diminishes the scattering effect caused by integer-order differences in seismic inversion.The outcomes of model tests validate the efficacy of this method,highlighting its ability to enhance the overall accuracy of the inversion process.
基金supported by the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(Qingdao)(Grant No.2021QNLM020001)the National Key R&D Program of China(Grant No.2019YFC0605503C)+2 种基金the Major Scientific and Technological Projects of China National Petroleum Corporation(CNPC)(Grant No.ZD2019-183-003)the National Outstanding Youth Science Foundation(Grant No.41922028)the National Innovation Group Project(Grant No.41821002).
文摘Based on waveform fitting,full waveform inversion(FWI)is an important inversion method with the ability to reconstruct multi-parameter models in high precision.However,the strong nonlinear equation used in FWI presents the following challenges,such as low convergence efficiency,high dependence on the initial model,and the energy imbalance in deep region of the inverted model.To solve these inherent problems,we develop a timedomain elastic FWI method based on gradient preconditioning with the following details:(1)the limited memory Broyden Fletcher Goldfarb Shanno method with faster convergence is adopted to im-prove the inversion stability;(2)a multi-scaled inversion strategy is used to alleviate the nonlinear inversion instead of falling into the local minimum;(3)in addition,the pseudo-Hessian preconditioned illumination operator is involved for preconditioning the parameter gradients to improve the illumination equilibrium degree of deep structures.Based on the programming implementation of the new method,a deep depression model with five diffractors is used for testing.Compared with the conventional elastic FWI method,the technique proposed by this study has better effectiveness and accuracy on the inversion effect and con-vergence,respectively.
文摘Poly(m-phenylene isophthalamide)(PMIA),a key aromatic polyamide,is widely used for its outstanding mechanical strength,high thermal stability,and excellent insulation properties.However,different applications demand varying dielectric properties,so tailoring its dielectric per-formance is essential.PMIA was first synthesized in this study,followed by introducing pores and developing porous PMIA films and PMIA-based composites with reduced dielectric constants.Porous PMIA films were fabricated using the wet phase inversion process with N,N-dimethylac-etamide(DMAC)solvent and water as the non-solvent.The impact of casting solution composi-tion and coagulation bath temperature on pore structures was analyzed.A film produced with 18%PMIA and 5%LiCl in a 35℃coagulation bath achieved the lowest dielectric constant of 1.76 at 1 Hz,48%lower than the standard PMIA film,which had a tensile strength of 18.5 MPa and an initial degradation temperature of 320℃.
基金supported by National Key R&D Program of China(2018YFA0702501)National Natural Science Foundation of China (41974140)+1 种基金Science and Technology Management Department,China National Petroleum Corporation(2022DQ0604-01)China National Petroleum Corporation-China University of Petroleum (Beijing) Strategy。
文摘Traditional deep learning methods pursue complex and single network architectures without considering the petrophysical relationship between different elastic parameters.The mathematical and statistical significance of the inversion results may lead to model overfitting,especially when there are a limited number of well logs in a working area.Multitask learning provides an eff ective approach to addressing this issue.Simultaneously,learning multiple related tasks can improve a model’s generalization ability to a certain extent,thereby enhancing the performance of related tasks with an equal amount of labeled data.In this study,we propose an end-to-end multitask deep learning model that integrates a fully convolutional network and bidirectional gated recurrent unit for intelligent prestack inversion of“seismic data to elastic parameters.”The use of a Bayesian homoscedastic uncertainty-based loss function enables adaptive learning of the weight coeffi cients for diff erent elastic parameter inversion tasks,thereby reducing uncertainty during the inversion process.The proposed method combines the local feature perception of convolutional neural networks with the long-term memory of bidirectional gated recurrent networks.It maintains the rock physics constraint relationships among diff erent elastic parameters during the inversion process,demonstrating a high level of prediction accuracy.Numerical simulations and processing results of real seismic data validate the eff ectiveness and practicality of the proposed method.
基金National Natural Science Foundation of China under Grant 42276055National Key Research and Development Program under Grant 2022YFC2803503Fundamental Research Funds for the Central Universities under Grant 202262008.
文摘Full waveform inversion methods evaluate the properties of subsurface media by minimizing the misfit between synthetic and observed data.However,these methods omit measurement errors and physical assumptions in modeling,resulting in several problems in practical applications.In particular,full waveform inversion methods are very sensitive to erroneous observations(outliers)that violate the Gauss–Markov theorem.Herein,we propose a method for addressing spurious observations or outliers.Specifically,we remove outliers by inverting the synthetic data using the local convexity of the Gaussian distribution.To achieve this,we apply a waveform-like noise model based on a specific covariance matrix definition.Finally,we build an inversion problem based on the updated data,which is consistent with the wavefield reconstruction inversion method.Overall,we report an alternative optimization inversion problem for data containing outliers.The proposed method is robust because it uses uncertainties.This method enables accurate inversion,even when based on noisy models or a wrong wavelet.
