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.展开更多
A three-dimensional(3D)electromagnetic(EM)inversion algorithm based on the nonlinear conjugate gradient(NLCG)method and a two-color plane Gauss-Seidel(GS)multigrid(MG)forward solver is developed to improve inversion e...A three-dimensional(3D)electromagnetic(EM)inversion algorithm based on the nonlinear conjugate gradient(NLCG)method and a two-color plane Gauss-Seidel(GS)multigrid(MG)forward solver is developed to improve inversion efficiency.The results indicate that the computational efficiency of each inversion can be improved by approximately a factor of three by using the proposed MG solver.First,the accuracy of the MG solver is validated through a test on a synthetic model.Next,the numerical performance of the inversion algorithm is evaluated using this model.Finally,the inversion algorithm is applied to a field EM data collected at the Beiya gold polymetallic ore district.A 3D resistivity model is obtained,and the formation process of the metal ore is analyzed.展开更多
Co-assembling chiral molecules with achiral compounds via non-covalent interactions like areneperfluoroarene(AP) interactions offers an effective approach for fabricating chiral functional materials.Herein,chiral mole...Co-assembling chiral molecules with achiral compounds via non-covalent interactions like areneperfluoroarene(AP) interactions offers an effective approach for fabricating chiral functional materials.Herein,chiral molecules L/D-PF1 and L/D-PF2 with pyrene groups were synthesized and its chiroptical properties upon co-assembly with achiral compound octafluoronaphthalene(OFN) through AP interaction were systemically studied.The co-assembly of L/D-PF1/OFN and L/D-PF2/OFN exhibited distinct chiroptical properties such as circular dichroism(CD) and circularly polarized luminescence(CPL) signals.Chirality transfer from the chirality center of L/D-PF1 and L/D-PF2 to the achiral OFN and chiral amplification were successfully achieved.Besides,no significant CPL signal was observed in the self-assembly of L/DPF1 or L/D-PF2 while co-assembly with OFN exhibited obvious CPL amplification induced by AP interaction.Notably,a reversal CD signal and CPL signal could be observed in L/D-PF2/OFN when the molar ratio changed from 1:1 to 1:2 while not found in L/D-PF1/OFN,indicating that that minor structural changes of molecules could cause large changes in assembly.In addition,a series of computational calculations were conducted to verify the AP interaction between L-PF1/L-PF2 and OFN.This work demonstrated that arene-perfluoroarene interaction could drive chiral transfer,chiral amplification and chiral inversion and provided a new method for the preparation of chiroptical materials.展开更多
Diabetes is a widespread disease affecting millions of people,making it one of the leading causes of death in the world.It is a leading cause of cardiovascular disease and end-stage renal disease.Despite advancements ...Diabetes is a widespread disease affecting millions of people,making it one of the leading causes of death in the world.It is a leading cause of cardiovascular disease and end-stage renal disease.Despite advancements in treatment,including insulin therapy and glucose monitoring devices,diabetes continues to significantly impact quality of life and current modalities do not reverse the end-organ damage associated with its progression.While traditionally indicated for type 1 diabetes,recent clinical practice refinements have made pancreas transplants available to select type 2 diabetics meeting specific criteria.These transplants are usually a part of a simultaneous kidney-pancreas transplant.However,although less frequently performed,transplants of pancreas alone or pancreas after kidney transplant are still available.For selected diabetic patients,pancreas transplants offer significant survival benefits and the improvement of cardiovascular and metabolic complications;however,they are not without risks.Complications such as bleeding,vascular thrombosis,infection,organ leak,and rejection are possible.Another challenge to pancreas transplantation is the decreasing number of procedures being performed due to decline in the volume of available highquality allografts and resource constraints of transplant centers.Advancements in monitoring and treatment of diabetes are contributing to the decline in pancreas transplants nowadays.展开更多
Owing to intensified globalization and informatization,the structures of the urban scale hierarchy and urban networks between cities have become increasingly intertwined,resulting in different spatial effects.Therefor...Owing to intensified globalization and informatization,the structures of the urban scale hierarchy and urban networks between cities have become increasingly intertwined,resulting in different spatial effects.Therefore,this paper analyzes the spatial interaction between urban scale hierarchy and urban networks in China from 2019 to 2023,drawing on Baidu migration data and employing a spatial simultaneous equation model.The results reveal a significant positive spatial correlation between cities with higher hierarchy and those with greater network centrality.Within a static framework,we identify a positive interaction between urban scale hierarchy and urban network centrality,while their spatial cross-effects manifest as negative neighborhood interactions based on geographical distance and positive cross-scale interactions shaped by network connections.Within a dynamic framework,changes in urban scale hierarchy and urban networks are mutually reinforcing,thereby widening disparities within the urban hierarchy.Furthermore,an increase in a city’s network centrality had a dampening effect on the population growth of neighboring cities and network-connected cities.This study enhances understanding of the spatial organisation of urban systems and offers insights for coordinated regional development.展开更多
The coexistence of emerging containments,such as antibiotic resistant bacteria(ARB),antibiotic-resistant genes(ARGs)and antibiotics,potentially influence elimination efficiencies in UV light-emitting diode(UV-LED)alon...The coexistence of emerging containments,such as antibiotic resistant bacteria(ARB),antibiotic-resistant genes(ARGs)and antibiotics,potentially influence elimination efficiencies in UV light-emitting diode(UV-LED)alone and UV-LED/H_(2)O_(2) system as their complex interactions.Tetracycline(TC)degradation efficiency(kF)correlated closely with its UV molar absorbance(R^(2)=0.831)in UV-LED alone system and with·OH yield(R^(2)=0.999)in UV-LED/H_(2)O_(2) system across studied wavelengths(265,280 and 310 nm).The kF values for intracellular DNA(i-ARGs)also exhibited a high correlation with UV-LED wavelengths in both systems(R^(2)=0.997-0.999).The coexistence of TC and ARB/ARGs resulted in a mutual inhibition of their degradation efficiencies due to competition for photons and·OH,along with the consequent reduction in intracellular ROS within ARB,with their degradation efficiencies exhibiting marked dependence on wavelength in both systems.Notably,the UV-LED/H_(2)O_(2) system at 265 nm effectively achieved the simultaneous removal of TC,ARB and ARGs with minimal energy consumption,and successfully fragmented ARGs.The degradation pathway of TC was analyzed,and the biotoxicity of its degradation intermediates demonstrated the environmental friendliness and safety of UV-LED/H_(2)O_(2) technology.This study elucidated the competitive interactions between antibiotics and ARB/ARGs within UV-LED/H_(2)O_(2) system,providing a promising approach for their simultaneous removal while ensuring energy efficiency.展开更多
To address the challenges of complex fluvial sandbody distribution and difficult remaining oil recovery in mature continental oilfields,this study focuses on key issues in reservoir identification such as ambiguous na...To address the challenges of complex fluvial sandbody distribution and difficult remaining oil recovery in mature continental oilfields,this study focuses on key issues in reservoir identification such as ambiguous narrow-channel boundaries and subdivision of multi-stage superimposed sandbodies.Taking the Upper Cretaceous continental sandstone in the Sazhong Oilfield of the Daqing Placanticline as an example,a technical system integrating OVT high-resolution processing,multi-attribute fusion,and varible-scale inversion was developed to establish a complete workflow from seismic processing to reservoir prediction and remaining oil recovery.The following results are obtained.First,the Offset Vector Tile(OVT)seismic processing technology is extended,for the first time,from fracture imaging to sandbody prediction,in order to address the weak seismic responses from boundaries of narrow and thin sandbodies.A geology-oriented OVT partitioning method is developed to significantly improve the imaging accuracy,enabling identification of channel sandbodies as narrow as 50 m.Second,an amplitude-coherence dual-attribute fusion method is proposed for predicting narrow channel boundaries between wells.Constrained by a sedimentary unit-level sequence chronostratigraphic framework,this method accurately delineates 800-2000 m long subaqueous distributary channels with bifurcation-convergence features.Third,considering the superimposition of multi-stage channels,a three-level variable-scale stratigraphic model(sandstone groups,sublayers,sedimentary units)is constructed to overcome single-scale modeling limitations,successfully characterizing key sedimentary features like meandering river“cut-offs”through 3D seismic inversion.Based on these advances,a direct link between seismic prediction and remaining oil recovery is established.The horizontal wells deployed using narrow-channel predictions encountered oil-bearing sandstones in the horizontal section by 97%,and achieved initial daily production of 12.5 t per well.Precise identification of individual channel boundaries within 17 composite sandbodies guided recovery processes in 135 wells,yielding an average daily increase of 2.8 t per well and a cumulative increase of 13.6×10^(4)t.展开更多
The forward model of optical fiber strain induced by fractures,together with the associated model resolution matrix,is used to demonstrate the interpretability of fracture parameters once the fracture intersects the f...The forward model of optical fiber strain induced by fractures,together with the associated model resolution matrix,is used to demonstrate the interpretability of fracture parameters once the fracture intersects the fiber.A regularized inversion framework for fracture parameters is established to evaluate the influence of measured data quality on the accuracy of iterative regularized inversion.An interpretation approach for both fracture width and height is proposed,and the synthetic forward data with measurement error and field examples are employed to validate the accuracy of the simultaneous inversion of fracture width and height.The results indicate that,after the fracture contacts the fiber,the strain response is strongly sensitive only to the fracture parameters at the intersection location,whereas the interpretability of parameters at other locations remains limited.The iterative regularized inversion method effectively suppresses the impact of measurement error and exhibits high computational efficiency,showing clear advantages for inversion applications.When incorporating the first-order regularization with a Neumann boundary constraint on the tip width,the inverted fracture-width distribution becomes highly sensitive to fracture height;thus,combined with a bisection strategy,simultaneous inversion of fracture width and height can be achieved.Examination using the model resolution matrix,noisy synthetic data,and field data confirms that the iterative regularized inversion model for fracture width and height provides high interpretive accuracy and can be applied to the calculation and analysis of fracture width,fracture height,net pressure and other parameters.展开更多
The high-resolution nonlinear simultaneous inversion of petrophysical parameters is based on Bayesian statistics and combines petrophysics with geostatistical a priori information. We used the fast Fourier transform–...The high-resolution nonlinear simultaneous inversion of petrophysical parameters is based on Bayesian statistics and combines petrophysics with geostatistical a priori information. We used the fast Fourier transform–moving average（FFT–MA） and gradual deformation method（GDM） to obtain a reasonable variogram by using structural analysis and geostatistical a priori information of petrophysical parameters. Subsequently, we constructed the likelihood function according to the statistical petrophysical model. Finally, we used the Metropolis algorithm to sample the posteriori probability density and complete the inversion of the petrophysical parameters. We used the proposed method to process data from an oil fi eld in China and found good match between inversion and real data with high-resolution. In addition, the direct inversion of petrophysical parameters avoids the error accumulation and decreases the uncertainty, and increases the computational effi ciency.展开更多
A new prestack AVA simultaneous inversion using particle swarm optimization algorithm is proposed, which can obtain the elastic parameters such as P-wave and S-wave impedance from P-wave reflection data simultaneously...A new prestack AVA simultaneous inversion using particle swarm optimization algorithm is proposed, which can obtain the elastic parameters such as P-wave and S-wave impedance from P-wave reflection data simultaneously. Compared with the conventional AVA inversion based on generalized linear technique, this method does not depend on the initial model and can reach the global minimum. In order to increase the stability of the inversion, low-frequency trends of P-wave and S-wave impedances are built into the inversion. This method has been successfully applied to synthetic and field data. The estimated P-wave and S-wave impedances can be combined to derive other elastic parameters, which are sensitive for lithology identification and fluid prediction.展开更多
Inversion of Young’s modulus,Poisson’s ratio and density from pre-stack seismic data has been proved to be feasible and effective.However,the existing methods do not take full advantage of the prior information.With...Inversion of Young’s modulus,Poisson’s ratio and density from pre-stack seismic data has been proved to be feasible and effective.However,the existing methods do not take full advantage of the prior information.Without considering the lateral continuity of the inversion results,these methods need to invert the reflectivity first.In this paper,we propose multi-gather simultaneous inversion for pre-stack seismic data.Meanwhile,the total variation（TV）regularization,L1 norm regularization and initial model constraint are used.In order to solve the objective function contains L1norm,TV norm and L2 norm,we develop an algorithm based on split Bregman iteration.The main advantages of our method are as follows：（1）The elastic parameters are calculated directly from objective function rather than from their reflectivity,therefore the stability and accuracy of the inversion process can be ensured.（2）The inversion results are more in accordance with the prior geological information.（3）The lateral continuity of the inversion results are improved.The proposed method is illustrated by theoretical model data and experimented with a 2-D field data.展开更多
Variable-depth streamer seismic data are characterized by low and high frequencies and can be used to obtain high-quality and resolution images of complex subsurface structures. Taking advantage of the frequency range...Variable-depth streamer seismic data are characterized by low and high frequencies and can be used to obtain high-quality and resolution images of complex subsurface structures. Taking advantage of the frequency range in the variable-depth streamer data, we propose the simultaneous inversion of prestack data from variable-offset stack gathers to obtain the P-wave impedance, S-wave impedance, and density. Next, we validate the method by using model and actual variable-depth streamer data from the Huizhou block. The results suggest that the broadband data recorded by variable-depth streamers improve the signal-to-noise ratio and quality of the inversion results and outperform the constant-depth streamer data in delineating the underground stratigraphy.展开更多
In order to interpret the vertical electrical sounding data more reliably and effectively in the case of lacking proper priori information, two inverse schemes are proposed to invert combined re- sistivity and induced...In order to interpret the vertical electrical sounding data more reliably and effectively in the case of lacking proper priori information, two inverse schemes are proposed to invert combined re- sistivity and induced polarization data by using particle swarm optimization technique. Based on the computational formula of induced polarization, the inversion for chargeability/polarizability data can be transformed into inverting equivalent resistivity data. Then, the inversion for combined data can be decomposed into two procedures： inverting resistivity data and inverting equivalent resistivity data. A sequential inversion scheme is presented to run the two procedures sequentially. Contrast to the se- quential scheme, a simultaneous one is proposed to invert resistivity and induced polarization data si- multaneously. Both the sequential and simultaneous schemes are performed via centered-progressive particle swarm optimization algorithm for more exploratory purpose. Numerical experiments show that both the designed inversion algorithms can invert resistivity and induced polarization data suc- cessfully with fast convergence and high accuracy, even performed in a large search space. The inverse results are comparable to the results from generalized linear method. As an approximate importance sampler, the particle swarm optimization based algorithm can provide posterior analysis conveniently. We employ the posterior probability distributions of inverted model parameters to evaluate the per- formance and uncertainty of inversion. The posterior analysis and further field data testing show that the proposed inversion algorithms perform good sampling of the equivalence region and make sure that the global optimum can locate in the high probability areas.展开更多
Using the data of P-wave network and Zhejiang and travel time recorded at the Shanxi-reservoir seismological Fujian local networks, we implemented a simultaneous inversion of earthquake relocation and velocity struct...Using the data of P-wave network and Zhejiang and travel time recorded at the Shanxi-reservoir seismological Fujian local networks, we implemented a simultaneous inversion of earthquake relocation and velocity structure and determined the new locations of earthquakes in the Shanxi-reservoir. The results show that： （1） the overall epicenter distribution is NW directed, and the Shanxi reservoir induced seismicity has a close relationship to the Shuangxi-Jiaoxiyang fault; （2） the focal depth of the Shanxi reservoir induced seismicity is 5.4km in average, less than the average focal depth in the South China earthquake zone; （3） the focal depth is shallower on the reservoir shore and deeper in the reservoir inundation area. At the beginning of the reservoir induced seismicity, the focal depth increased gradually. This may be due to the gradual penetration of water into a larger depth that induced deeper earthquakes; and （4） there is a low P-wave velocity anomaly in the study area, located at the intersection of multiple faults in the reservoir inundation area. The Shanxi reservoir induced seismicity mostly occurred in this lowvelocity anomaly zone. This may be related to water penetration.展开更多
It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regio...It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regional earthquake tomography, but seldom are all three classes of model parameters updated simultaneously. This is mainly due to the trade-off between the different types of model parameters and the lack of different seismic phases to constrain the model parameters. Using a spherical-coordinate ray tracing algorithm for first and later（primary reflected） arrival tracing algorithm in combination with a popular linearized inversion solver, it is possible to simultaneously recover the three classes of model parameters in regional or global tomographic studies. In this paper we incorporate the multistage irregular shortest-path ray tracing algorithm（in a spherical coordinate system） with a subspace inversion solver to formulate a simultaneous inversion algorithm for triple model parameters updating using direct and later arrival time information.Comparison tests for two sets of data（noise free and added noise） indicate that the new triple-class parameter inversion algorithm is capable of obtaining nearly the same results as the double-class parameter inversion scheme. Furthermore,the proposed multi-parameter type inversion method is not sensitive to a modest level of picking error in the traveltime data, and also performs well with a relatively large uncertainty in earthquake hypocentral locations. This shows it to be a feasible and promising approach in regional or global tomographic applications.展开更多
This paper deals with the simultaneous inversion of velocity distribution and interface positions in a laterally heterogeneous medium using seismic travel time data. Due to the application of variable ...This paper deals with the simultaneous inversion of velocity distribution and interface positions in a laterally heterogeneous medium using seismic travel time data. Due to the application of variable step length and linear interpolation techniques (for both theoretical travel times and the partial derivative matrix), the computing speed is increased greatly. As compared with other domestic computing programs of this kind, it is five times or more faster than others. Besides reflected waves, the transmitted waves can also be used to form the partial derivative matrix of travel times with respect to interface positions. So this method can make full use of the information about interface positions contained in various transmitted waves and speed up the convergence rate. Digital tests and processing results of real data indicate the feasibility and effectiveness of the method and its program.展开更多
The structure-coupled joint inversion method of gravity and magnetic data is a powerful tool for?developing improved physical property models with high resolution and compatible features;?however, the conventional pro...The structure-coupled joint inversion method of gravity and magnetic data is a powerful tool for?developing improved physical property models with high resolution and compatible features;?however, the conventional procedure is inefficient due to the truncated singular values decomposition?(SVD) process at each iteration. To improve the algorithm, a technique using damped leastsquares?is adopted to calculate the structural term of model updates, instead of the truncated SVD. This?produces structural coupled density and magnetization images with high efficiency. A so-called?coupling factor is introduced to regulate the tuning of the desired final structural similarity level.?Synthetic examples show that the joint inversion results are internally consistent and achieve?higher?resolution than separated. The acceptable runtime performance of the damped least squares?technique used in joint inversion indicates that it is more suitable for practical use than the truncated SVD method.展开更多
The static correction of a near-surface model may be improved by using travel time tomographic inversion.We discuss unfavorable factors in the inversion of surface seismic waves that have been analyzed by the first br...The static correction of a near-surface model may be improved by using travel time tomographic inversion.We discuss unfavorable factors in the inversion of surface seismic waves that have been analyzed by the first break.These factors show that sources and geophones arranged on the surface,or close to the surface,give a first break that only includes the direct wave and the up going wave from the down going to up going transition.These up going waves have weak directivity when they arrive at a geophone and so the rays passing through the grids have small directional differences and a narrow azimuth.Drawing lessons from the advantages of Vertical Seismic Profiling(VSP) acquisition mode we describe a pseudo well-surface simultaneous travel time tomographic inversion of a near-surface model.The well depth should be increased in the surface seismic study to produce a pure up going wave,to enhance the verticality of the rays and to increase the azimuth and shorten path length of the rays.Simulations of the effect of well depth on a pseudo well-surface simultaneous travel time tomographic inversion model are reported.The results show that the static corrections are improved significantly when the well depth extends below the weathered or sub-weathered layers.The root mean square error of the statics is 1.14 or 0.93 ms for these two situations,respectively.展开更多
A method for simultaneous determination of mixed model parameters,which have different physical dimensions or different responses to data,is presented.Mixed parameter estimation from observed data within a single mode...A method for simultaneous determination of mixed model parameters,which have different physical dimensions or different responses to data,is presented.Mixed parameter estimation from observed data within a single model space shows instabilities and trade-offs of the solutions. We separate the model space into N-subspaces based on their physical properties or computational convenience and solve the N-subspaces systems by damped least-squares and singular-value decomposition. Since the condition number of each subsystem is smaller than that of the single global system,the approach can greatly increase the stability of the inversion. We also introduce different damping factors into the subsystems to reduce the tradeoffs between the different parameters. The damping factors depend on the conditioning of the subsystems and may be adequately chosen in a range from 0.1 % to 10 % of the largest singular value. We illustrate the method with an example of simultaneous determination of source history,source geometry,and hypocentral location from regional seismograms,although it is applicable to any geophysical inversion.展开更多
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.展开更多
基金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.
基金financially supported by the National Science and Technology Major Project,China(No.2024ZD1002100)the National Natural Science Foundation of China(Nos.42330801,42474112,42504062)+1 种基金the China Postdoctoral Science Foundation(No.2024M761704)Shuimu Tsinghua Scholar Program of Tsinghua University,China(No.2024SM114)。
文摘A three-dimensional(3D)electromagnetic(EM)inversion algorithm based on the nonlinear conjugate gradient(NLCG)method and a two-color plane Gauss-Seidel(GS)multigrid(MG)forward solver is developed to improve inversion efficiency.The results indicate that the computational efficiency of each inversion can be improved by approximately a factor of three by using the proposed MG solver.First,the accuracy of the MG solver is validated through a test on a synthetic model.Next,the numerical performance of the inversion algorithm is evaluated using this model.Finally,the inversion algorithm is applied to a field EM data collected at the Beiya gold polymetallic ore district.A 3D resistivity model is obtained,and the formation process of the metal ore is analyzed.
基金financially supported by the National Natural Science Foundation of China (Nos.22171165 and 22371170)Natural Science Foundation of Shandong Province (No.ZR2022MB080)Scientific and Technological Frontiers in Project of Henan Province(No.242102110192)。
文摘Co-assembling chiral molecules with achiral compounds via non-covalent interactions like areneperfluoroarene(AP) interactions offers an effective approach for fabricating chiral functional materials.Herein,chiral molecules L/D-PF1 and L/D-PF2 with pyrene groups were synthesized and its chiroptical properties upon co-assembly with achiral compound octafluoronaphthalene(OFN) through AP interaction were systemically studied.The co-assembly of L/D-PF1/OFN and L/D-PF2/OFN exhibited distinct chiroptical properties such as circular dichroism(CD) and circularly polarized luminescence(CPL) signals.Chirality transfer from the chirality center of L/D-PF1 and L/D-PF2 to the achiral OFN and chiral amplification were successfully achieved.Besides,no significant CPL signal was observed in the self-assembly of L/DPF1 or L/D-PF2 while co-assembly with OFN exhibited obvious CPL amplification induced by AP interaction.Notably,a reversal CD signal and CPL signal could be observed in L/D-PF2/OFN when the molar ratio changed from 1:1 to 1:2 while not found in L/D-PF1/OFN,indicating that that minor structural changes of molecules could cause large changes in assembly.In addition,a series of computational calculations were conducted to verify the AP interaction between L-PF1/L-PF2 and OFN.This work demonstrated that arene-perfluoroarene interaction could drive chiral transfer,chiral amplification and chiral inversion and provided a new method for the preparation of chiroptical materials.
文摘Diabetes is a widespread disease affecting millions of people,making it one of the leading causes of death in the world.It is a leading cause of cardiovascular disease and end-stage renal disease.Despite advancements in treatment,including insulin therapy and glucose monitoring devices,diabetes continues to significantly impact quality of life and current modalities do not reverse the end-organ damage associated with its progression.While traditionally indicated for type 1 diabetes,recent clinical practice refinements have made pancreas transplants available to select type 2 diabetics meeting specific criteria.These transplants are usually a part of a simultaneous kidney-pancreas transplant.However,although less frequently performed,transplants of pancreas alone or pancreas after kidney transplant are still available.For selected diabetic patients,pancreas transplants offer significant survival benefits and the improvement of cardiovascular and metabolic complications;however,they are not without risks.Complications such as bleeding,vascular thrombosis,infection,organ leak,and rejection are possible.Another challenge to pancreas transplantation is the decreasing number of procedures being performed due to decline in the volume of available highquality allografts and resource constraints of transplant centers.Advancements in monitoring and treatment of diabetes are contributing to the decline in pancreas transplants nowadays.
基金Under the auspices of the National Natural Science Foundation of China(No.42371222,41971167)Fundamental Scientific Research Funds of Central China Normal University(No.CCNU24ZZ120)。
文摘Owing to intensified globalization and informatization,the structures of the urban scale hierarchy and urban networks between cities have become increasingly intertwined,resulting in different spatial effects.Therefore,this paper analyzes the spatial interaction between urban scale hierarchy and urban networks in China from 2019 to 2023,drawing on Baidu migration data and employing a spatial simultaneous equation model.The results reveal a significant positive spatial correlation between cities with higher hierarchy and those with greater network centrality.Within a static framework,we identify a positive interaction between urban scale hierarchy and urban network centrality,while their spatial cross-effects manifest as negative neighborhood interactions based on geographical distance and positive cross-scale interactions shaped by network connections.Within a dynamic framework,changes in urban scale hierarchy and urban networks are mutually reinforcing,thereby widening disparities within the urban hierarchy.Furthermore,an increase in a city’s network centrality had a dampening effect on the population growth of neighboring cities and network-connected cities.This study enhances understanding of the spatial organisation of urban systems and offers insights for coordinated regional development.
基金supported by Major Scientific and Technological Innovation Project of Shandong Province(No.2020CXGC011204)Qingdao Natural Science Foundation(No.23-2-1-234-zyyd-jch).
文摘The coexistence of emerging containments,such as antibiotic resistant bacteria(ARB),antibiotic-resistant genes(ARGs)and antibiotics,potentially influence elimination efficiencies in UV light-emitting diode(UV-LED)alone and UV-LED/H_(2)O_(2) system as their complex interactions.Tetracycline(TC)degradation efficiency(kF)correlated closely with its UV molar absorbance(R^(2)=0.831)in UV-LED alone system and with·OH yield(R^(2)=0.999)in UV-LED/H_(2)O_(2) system across studied wavelengths(265,280 and 310 nm).The kF values for intracellular DNA(i-ARGs)also exhibited a high correlation with UV-LED wavelengths in both systems(R^(2)=0.997-0.999).The coexistence of TC and ARB/ARGs resulted in a mutual inhibition of their degradation efficiencies due to competition for photons and·OH,along with the consequent reduction in intracellular ROS within ARB,with their degradation efficiencies exhibiting marked dependence on wavelength in both systems.Notably,the UV-LED/H_(2)O_(2) system at 265 nm effectively achieved the simultaneous removal of TC,ARB and ARGs with minimal energy consumption,and successfully fragmented ARGs.The degradation pathway of TC was analyzed,and the biotoxicity of its degradation intermediates demonstrated the environmental friendliness and safety of UV-LED/H_(2)O_(2) technology.This study elucidated the competitive interactions between antibiotics and ARB/ARGs within UV-LED/H_(2)O_(2) system,providing a promising approach for their simultaneous removal while ensuring energy efficiency.
基金Supported by the China National Science and Technology Major Project(2025ZD1407000)PetroChina Science and Technology Major Project(2023ZZ22)。
文摘To address the challenges of complex fluvial sandbody distribution and difficult remaining oil recovery in mature continental oilfields,this study focuses on key issues in reservoir identification such as ambiguous narrow-channel boundaries and subdivision of multi-stage superimposed sandbodies.Taking the Upper Cretaceous continental sandstone in the Sazhong Oilfield of the Daqing Placanticline as an example,a technical system integrating OVT high-resolution processing,multi-attribute fusion,and varible-scale inversion was developed to establish a complete workflow from seismic processing to reservoir prediction and remaining oil recovery.The following results are obtained.First,the Offset Vector Tile(OVT)seismic processing technology is extended,for the first time,from fracture imaging to sandbody prediction,in order to address the weak seismic responses from boundaries of narrow and thin sandbodies.A geology-oriented OVT partitioning method is developed to significantly improve the imaging accuracy,enabling identification of channel sandbodies as narrow as 50 m.Second,an amplitude-coherence dual-attribute fusion method is proposed for predicting narrow channel boundaries between wells.Constrained by a sedimentary unit-level sequence chronostratigraphic framework,this method accurately delineates 800-2000 m long subaqueous distributary channels with bifurcation-convergence features.Third,considering the superimposition of multi-stage channels,a three-level variable-scale stratigraphic model(sandstone groups,sublayers,sedimentary units)is constructed to overcome single-scale modeling limitations,successfully characterizing key sedimentary features like meandering river“cut-offs”through 3D seismic inversion.Based on these advances,a direct link between seismic prediction and remaining oil recovery is established.The horizontal wells deployed using narrow-channel predictions encountered oil-bearing sandstones in the horizontal section by 97%,and achieved initial daily production of 12.5 t per well.Precise identification of individual channel boundaries within 17 composite sandbodies guided recovery processes in 135 wells,yielding an average daily increase of 2.8 t per well and a cumulative increase of 13.6×10^(4)t.
基金Supported by the Ministry of Education U40 Program(ZYGXONJSKYCXNLZCXM-E19)National Natural Science Foundation of China(52574078)。
文摘The forward model of optical fiber strain induced by fractures,together with the associated model resolution matrix,is used to demonstrate the interpretability of fracture parameters once the fracture intersects the fiber.A regularized inversion framework for fracture parameters is established to evaluate the influence of measured data quality on the accuracy of iterative regularized inversion.An interpretation approach for both fracture width and height is proposed,and the synthetic forward data with measurement error and field examples are employed to validate the accuracy of the simultaneous inversion of fracture width and height.The results indicate that,after the fracture contacts the fiber,the strain response is strongly sensitive only to the fracture parameters at the intersection location,whereas the interpretability of parameters at other locations remains limited.The iterative regularized inversion method effectively suppresses the impact of measurement error and exhibits high computational efficiency,showing clear advantages for inversion applications.When incorporating the first-order regularization with a Neumann boundary constraint on the tip width,the inverted fracture-width distribution becomes highly sensitive to fracture height;thus,combined with a bisection strategy,simultaneous inversion of fracture width and height can be achieved.Examination using the model resolution matrix,noisy synthetic data,and field data confirms that the iterative regularized inversion model for fracture width and height provides high interpretive accuracy and can be applied to the calculation and analysis of fracture width,fracture height,net pressure and other parameters.
基金sponsored by the National Basic Research Program of China(No.2013CB228604)the Major National Science and Technology Projects(No.2011ZX05009)+1 种基金the Natural Science Foundation of Shandong Province(No.ZR2011DQ013)the National Science Foundation of China(No.41204085)
文摘The high-resolution nonlinear simultaneous inversion of petrophysical parameters is based on Bayesian statistics and combines petrophysics with geostatistical a priori information. We used the fast Fourier transform–moving average（FFT–MA） and gradual deformation method（GDM） to obtain a reasonable variogram by using structural analysis and geostatistical a priori information of petrophysical parameters. Subsequently, we constructed the likelihood function according to the statistical petrophysical model. Finally, we used the Metropolis algorithm to sample the posteriori probability density and complete the inversion of the petrophysical parameters. We used the proposed method to process data from an oil fi eld in China and found good match between inversion and real data with high-resolution. In addition, the direct inversion of petrophysical parameters avoids the error accumulation and decreases the uncertainty, and increases the computational effi ciency.
基金supported by the National Natural Science Foundation of China(Nos.41004096 and 41230318)
文摘A new prestack AVA simultaneous inversion using particle swarm optimization algorithm is proposed, which can obtain the elastic parameters such as P-wave and S-wave impedance from P-wave reflection data simultaneously. Compared with the conventional AVA inversion based on generalized linear technique, this method does not depend on the initial model and can reach the global minimum. In order to increase the stability of the inversion, low-frequency trends of P-wave and S-wave impedances are built into the inversion. This method has been successfully applied to synthetic and field data. The estimated P-wave and S-wave impedances can be combined to derive other elastic parameters, which are sensitive for lithology identification and fluid prediction.
基金supported by the National Natural Science Foundation of China (Nos.61775030,61571096,41301460,61362018,and 41274127)the key projects of Hunan Provincial Department of Education (No.16A174)
文摘Inversion of Young’s modulus,Poisson’s ratio and density from pre-stack seismic data has been proved to be feasible and effective.However,the existing methods do not take full advantage of the prior information.Without considering the lateral continuity of the inversion results,these methods need to invert the reflectivity first.In this paper,we propose multi-gather simultaneous inversion for pre-stack seismic data.Meanwhile,the total variation（TV）regularization,L1 norm regularization and initial model constraint are used.In order to solve the objective function contains L1norm,TV norm and L2 norm,we develop an algorithm based on split Bregman iteration.The main advantages of our method are as follows：（1）The elastic parameters are calculated directly from objective function rather than from their reflectivity,therefore the stability and accuracy of the inversion process can be ensured.（2）The inversion results are more in accordance with the prior geological information.（3）The lateral continuity of the inversion results are improved.The proposed method is illustrated by theoretical model data and experimented with a 2-D field data.
基金sponsored by Comprehensive Research of CNOOC(China) Limited(No.YXKY–2013–SZ–02)
文摘Variable-depth streamer seismic data are characterized by low and high frequencies and can be used to obtain high-quality and resolution images of complex subsurface structures. Taking advantage of the frequency range in the variable-depth streamer data, we propose the simultaneous inversion of prestack data from variable-offset stack gathers to obtain the P-wave impedance, S-wave impedance, and density. Next, we validate the method by using model and actual variable-depth streamer data from the Huizhou block. The results suggest that the broadband data recorded by variable-depth streamers improve the signal-to-noise ratio and quality of the inversion results and outperform the constant-depth streamer data in delineating the underground stratigraphy.
基金supported by the National Natural Science Foundation of China(No.41574123)
文摘In order to interpret the vertical electrical sounding data more reliably and effectively in the case of lacking proper priori information, two inverse schemes are proposed to invert combined re- sistivity and induced polarization data by using particle swarm optimization technique. Based on the computational formula of induced polarization, the inversion for chargeability/polarizability data can be transformed into inverting equivalent resistivity data. Then, the inversion for combined data can be decomposed into two procedures： inverting resistivity data and inverting equivalent resistivity data. A sequential inversion scheme is presented to run the two procedures sequentially. Contrast to the se- quential scheme, a simultaneous one is proposed to invert resistivity and induced polarization data si- multaneously. Both the sequential and simultaneous schemes are performed via centered-progressive particle swarm optimization algorithm for more exploratory purpose. Numerical experiments show that both the designed inversion algorithms can invert resistivity and induced polarization data suc- cessfully with fast convergence and high accuracy, even performed in a large search space. The inverse results are comparable to the results from generalized linear method. As an approximate importance sampler, the particle swarm optimization based algorithm can provide posterior analysis conveniently. We employ the posterior probability distributions of inverted model parameters to evaluate the per- formance and uncertainty of inversion. The posterior analysis and further field data testing show that the proposed inversion algorithms perform good sampling of the equivalence region and make sure that the global optimum can locate in the high probability areas.
基金supported by the National Key Technology R&D Program(2008BAC38B03-01-05)the Earthquake Scientific Research Project(200708020),China
文摘Using the data of P-wave network and Zhejiang and travel time recorded at the Shanxi-reservoir seismological Fujian local networks, we implemented a simultaneous inversion of earthquake relocation and velocity structure and determined the new locations of earthquakes in the Shanxi-reservoir. The results show that： （1） the overall epicenter distribution is NW directed, and the Shanxi reservoir induced seismicity has a close relationship to the Shuangxi-Jiaoxiyang fault; （2） the focal depth of the Shanxi reservoir induced seismicity is 5.4km in average, less than the average focal depth in the South China earthquake zone; （3） the focal depth is shallower on the reservoir shore and deeper in the reservoir inundation area. At the beginning of the reservoir induced seismicity, the focal depth increased gradually. This may be due to the gradual penetration of water into a larger depth that induced deeper earthquakes; and （4） there is a low P-wave velocity anomaly in the study area, located at the intersection of multiple faults in the reservoir inundation area. The Shanxi reservoir induced seismicity mostly occurred in this lowvelocity anomaly zone. This may be related to water penetration.
基金partially supported by the Doctoral Programming Research Fund of Higher Education, Chinese Ministry of Education (No. 20110205110010)
文摘It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regional earthquake tomography, but seldom are all three classes of model parameters updated simultaneously. This is mainly due to the trade-off between the different types of model parameters and the lack of different seismic phases to constrain the model parameters. Using a spherical-coordinate ray tracing algorithm for first and later（primary reflected） arrival tracing algorithm in combination with a popular linearized inversion solver, it is possible to simultaneously recover the three classes of model parameters in regional or global tomographic studies. In this paper we incorporate the multistage irregular shortest-path ray tracing algorithm（in a spherical coordinate system） with a subspace inversion solver to formulate a simultaneous inversion algorithm for triple model parameters updating using direct and later arrival time information.Comparison tests for two sets of data（noise free and added noise） indicate that the new triple-class parameter inversion algorithm is capable of obtaining nearly the same results as the double-class parameter inversion scheme. Furthermore,the proposed multi-parameter type inversion method is not sensitive to a modest level of picking error in the traveltime data, and also performs well with a relatively large uncertainty in earthquake hypocentral locations. This shows it to be a feasible and promising approach in regional or global tomographic applications.
文摘This paper deals with the simultaneous inversion of velocity distribution and interface positions in a laterally heterogeneous medium using seismic travel time data. Due to the application of variable step length and linear interpolation techniques (for both theoretical travel times and the partial derivative matrix), the computing speed is increased greatly. As compared with other domestic computing programs of this kind, it is five times or more faster than others. Besides reflected waves, the transmitted waves can also be used to form the partial derivative matrix of travel times with respect to interface positions. So this method can make full use of the information about interface positions contained in various transmitted waves and speed up the convergence rate. Digital tests and processing results of real data indicate the feasibility and effectiveness of the method and its program.
文摘The structure-coupled joint inversion method of gravity and magnetic data is a powerful tool for?developing improved physical property models with high resolution and compatible features;?however, the conventional procedure is inefficient due to the truncated singular values decomposition?(SVD) process at each iteration. To improve the algorithm, a technique using damped leastsquares?is adopted to calculate the structural term of model updates, instead of the truncated SVD. This?produces structural coupled density and magnetization images with high efficiency. A so-called?coupling factor is introduced to regulate the tuning of the desired final structural similarity level.?Synthetic examples show that the joint inversion results are internally consistent and achieve?higher?resolution than separated. The acceptable runtime performance of the damped least squares?technique used in joint inversion indicates that it is more suitable for practical use than the truncated SVD method.
文摘The static correction of a near-surface model may be improved by using travel time tomographic inversion.We discuss unfavorable factors in the inversion of surface seismic waves that have been analyzed by the first break.These factors show that sources and geophones arranged on the surface,or close to the surface,give a first break that only includes the direct wave and the up going wave from the down going to up going transition.These up going waves have weak directivity when they arrive at a geophone and so the rays passing through the grids have small directional differences and a narrow azimuth.Drawing lessons from the advantages of Vertical Seismic Profiling(VSP) acquisition mode we describe a pseudo well-surface simultaneous travel time tomographic inversion of a near-surface model.The well depth should be increased in the surface seismic study to produce a pure up going wave,to enhance the verticality of the rays and to increase the azimuth and shorten path length of the rays.Simulations of the effect of well depth on a pseudo well-surface simultaneous travel time tomographic inversion model are reported.The results show that the static corrections are improved significantly when the well depth extends below the weathered or sub-weathered layers.The root mean square error of the statics is 1.14 or 0.93 ms for these two situations,respectively.
基金supported by Innovation Project of Chinese Academy of Sciences
文摘A method for simultaneous determination of mixed model parameters,which have different physical dimensions or different responses to data,is presented.Mixed parameter estimation from observed data within a single model space shows instabilities and trade-offs of the solutions. We separate the model space into N-subspaces based on their physical properties or computational convenience and solve the N-subspaces systems by damped least-squares and singular-value decomposition. Since the condition number of each subsystem is smaller than that of the single global system,the approach can greatly increase the stability of the inversion. We also introduce different damping factors into the subsystems to reduce the tradeoffs between the different parameters. The damping factors depend on the conditioning of the subsystems and may be adequately chosen in a range from 0.1 % to 10 % of the largest singular value. We illustrate the method with an example of simultaneous determination of source history,source geometry,and hypocentral location from regional seismograms,although it is applicable to any geophysical inversion.
基金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.
