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.展开更多
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.展开更多
Simultaneous inversion of scattering and absorption attenuation is of great significance for investigating small-scale inhomogeneities and inelastic properties of the subsurface.However,applying this to complex geophy...Simultaneous inversion of scattering and absorption attenuation is of great significance for investigating small-scale inhomogeneities and inelastic properties of the subsurface.However,applying this to complex geophysical issues is constrained by the costly computational requirements for simulations and inversions using existing methods.The coupling effects between the scattering and absorption coefficients cause parameter crosstalk artifacts in multi-parameter inversion,significantly increasing the nonlinearity of the inverse processes.This paper proposes a robust and effective simultaneous inversion method for scattering and absorption attenuation.The propagation of coda energy is modeled using the finite-element method based on the frequency-domain diffusion equation.We employ the truncated Gauss-Newton technique for the simultaneous estimation of scattering and absorption coefficients to decouple the two attenuation parameters during the inversion procedure.Nevertheless,the inversion may reach a local minimum if the energy frequency is inappropriate due to a phase mismatch between the simulated and measured energies.To mitigate this issue,we provide a frequency selection criterion that considers the substantial spectral shift of the energy density spectrum toward lower frequencies.Numerical examples using synthetic and experimental data indicate that our method can significantly reduce computational complexity and suppress two-parameter crosstalk without requiring a precise initial attenuation model.展开更多
In this paper,obtainment of 3 by simultaneous inversion of configurations of C_(22),C_(23) in 2 via its ditosylate or dimesylate was unsuccessful,leading to elimination products instead.However, monotosylate of 2 can ...In this paper,obtainment of 3 by simultaneous inversion of configurations of C_(22),C_(23) in 2 via its ditosylate or dimesylate was unsuccessful,leading to elimination products instead.However, monotosylate of 2 can be readily obtained with exclusive regioselectivity,which in turn gave a single epoxide 7b Hydrolysis of 7b gave 3 and 2 is a ratio of 2:1,with a combined yield of 80%.The recovered 2 can of course be recycled again.展开更多
Prestack seismic inversion methods adopt approximations of the Zoeppritz equations to describe the relation between reflection coefficients and P-wave velocity, S-wave velocity, and density. However, the error in thes...Prestack seismic inversion methods adopt approximations of the Zoeppritz equations to describe the relation between reflection coefficients and P-wave velocity, S-wave velocity, and density. However, the error in these approximations increases with increasing angle of incidence and variation of the elastic parameters, which increases the number of inversion solutions and minimizes the inversion accuracy. In this study, we explore a method for solving the reflection coefficients by using the Zoeppritz equations. To increase the accuracy of prestack inversion, the simultaneous inversion of P-wave velocity, S-wave velocity, and density by using prestack large-angle seismic data is proposed based on generalized linear inversion theory. Moreover, we reduce the ill posedness and increase the convergence of prestack inversion by using the regularization constraint damping factor and the conjugate gradient algorithm. The proposed prestack inversion method uses prestack large-angle seismic data to obtain accurate seismic elastic parameters that conform to prestack seismic data and are consistent with logging data from wells.展开更多
Simultaneous waveform inversion was used to predict lithofacies and fluid type across the field. Very often, characterizing reservoirs in terms of lithology and fluid type using conventional methods is replete with un...Simultaneous waveform inversion was used to predict lithofacies and fluid type across the field. Very often, characterizing reservoirs in terms of lithology and fluid type using conventional methods is replete with uncertainties, especially in marginal fields. An approach is employed in this study that integrated rock physics and waveform inverse modelling for lithology and fluid-type characterization to appropriately identify potential hydrocarbon saturated zones and their corresponding lithology. Seismic and well-log data were analyzed using Hampson Russel software. The method adopted includes lithofacies and fluid content analysis using rock physics parameters and seismic simultaneous inverse modelling. Rock physics analysis identified 2 broad reservoirs namely: HDZ1 and HDZ2 reservoirs. Results from the inverse modelling showed that low values of acoustic impedance from 19,743 to 20,487 (ft/s)(g/cc) reflect hydrocarbon-bearing reservoirs while medium to high values shows brine and shale respectively, with brine zone ranging from 20,487 to 22,531 (ft/s)(g/cc) and shale above 22,531 (ft/s)(g/cc). Two lithofacies were identified from inversion analysis of Vp/Vs and Mu-Rho, namely: sand and shale with VpVs 1.95 values respectively. Mu-Rho > 12.29 (GPa)(g/cc) and <12.29 (GPa) (g/cc) represent sand and shale respectively. From 3D volume, it was observed that a high accumulation of hydrocarbon was observed to be saturated at the north to the eastern part of the field forming a meandering channel. Sands were mainly distributed around the northeastern to the southwestern part of the field, that tends to be away from Well 029. This was also validated by the volume of rigidity modulus (Mu-Rho) showing high values indicating sands fall within the northeastern part of the field.展开更多
Mesozoic-Palaeozoic marine carbonate rocks are crucial hydrocarbon reservoirs in the Central Uplift area of the South Yellow Sea Basin(SYSB).Due to the scarcity of boreholes and the significant heterogeneity of carbon...Mesozoic-Palaeozoic marine carbonate rocks are crucial hydrocarbon reservoirs in the Central Uplift area of the South Yellow Sea Basin(SYSB).Due to the scarcity of boreholes and the significant heterogeneity of carbonate reservoirs,the distribution of porous carbonate reservoirs and their related key controlling factors remain unclear.In this study,factors affecting the distribution of porous Carboniferous-Early Permian carbonate reservoirs in the SYSB were investigated through seismic inversion and isotope analysis.The log-seismic characteristics of porous carbonate reservoirs,sensitive lithology parameters,and physical property parameters were extracted and analyzed.The pre-stack simultaneous inversion technique was applied to predict the lithology and physical properties of porous carbonate reservoirs.Moreover,the sedimentary of carbonate was analyzed using isotopes of carbon,oxygen,and strontium.The results show that porous carbonate reservoirs are mainly developed in the open platform sediments with porosities of 3%-5%and are mainly distributed in the paleo-highland(Huanglong Formation and Chuanshan Formation)and the slope of paleo-highland(Hezhou Formation).The porous carbonate reservoirs of the Qixia Formation are only locally developed.In addition,the negativeδ13C excursions indicate a warm and humid tropical climate with three sea-level fluctuations in the study area from the Carboniferous to Early Permian.The favorable conditions for developing porous carbonate rocks include the sedimentary environment and diagenetic process.The primary pore tends to form in high-energy environments of the paleo-highland,and the secondary pore is increased by dissolution during the syngenetic or quasi-syngenetic period.According to the hydrocarbon potential analysis,the Late Ordovician Wufeng Formation and Lower Silurian Gaojiabian Formation are the source rocks in the high-maturity-over-maturity stage,the Carboniferous-Lower Permian carbonate is the good reservoirs,and the Late Permian Longtan-Dalong Formation is the stable seal,ensuring a huge hydrocarbon accumulation potential in SYSB.The methods proposed in this study can be applied to other carbonate-dominated strata worldwide.展开更多
Considering Zoeppritz equations, reflections of PP and PS are only the function of ratios of density and velocity. So the inversion results will be the same if the ratios are the same but values of density, velocities...Considering Zoeppritz equations, reflections of PP and PS are only the function of ratios of density and velocity. So the inversion results will be the same if the ratios are the same but values of density, velocities of P- wave and S-wave are different without strict constraint. This paper makes efforts to explore nonlinear simultaneous PP and PS inversion with expectation to reduce the ambiguity of AVO analysis by utilizing the redundancy of multi-component AVO measurements. Accurate estimation of ratio parameters depends on independence of input data. There are only two independent AVO attributes for PP reflectivity （i.e. intercept and gradient） and two for PS reflectivity （i.e. pseudo-intercept and pseudo-gradient or extreme amplitude）, respectively. For individual PP and PS inversion, the values of least-squares objective function do not converge around a large neighborhood of chosen true model parameters. Fortunately for joint PP and PS inversion the values of the least-squares objective function show closed contours with single minima. Finally the power function fitting is used to provide a higher precision AVO attributes than traditional polynomial fitting. By using the four independent fitting attributes （two independent attributes for PP and PS respectively）, the inversion of four ratio parameters （velocities and densities） would be estimated with less errors than that in traditional method.展开更多
As an important indicator parameter of fluid identification,fluid factor has always been a concern for scholars.However,when predicting Russell fluid factor or effective pore-fluid bulk modulus,it is necessary to intr...As an important indicator parameter of fluid identification,fluid factor has always been a concern for scholars.However,when predicting Russell fluid factor or effective pore-fluid bulk modulus,it is necessary to introduce a new rock skeleton parameter which is the dry-rock VP/VS ratio squared(DVRS).In the process of fluid factor calculation or inversion,the existing methods take this parameter as a static constant,which has been estimated in advance,and then apply it to the fluid factor calculation and inversion.The fluid identification analysis based on a portion of the Marmousi 2 model and numerical forward modeling test show that,taking the DVRS as a static constant will limit the identification ability of fluid factor and reduce the inversion accuracy.To solve the above problems,we proposed a new method to regard the DVRS as a dynamic variable varying with depth and lithology for the first time,then apply it to fluid factor calculation and inversion.Firstly,the exact Zoeppritz equations are rewritten into a new form containing the fluid factor and DVRS of upper and lower layers.Next,the new equations are applied to the four parameters simultaneous inversion based on the generalized nonlinear inversion(GNI)method.The testing results on a portion of the Marmousi 2 model and field data show that dynamic DVRS can significantly improve the fluid factor identification ability,effectively suppress illusion.Both synthetic and filed data tests also demonstrate that the GNI method based on Bayesian deterministic inversion(BDI)theory can successfully solve the above four parameter simultaneous inversion problem,and taking the dynamic DVRS as a target inversion parameter can effectively improve the inversion accuracy of fluid factor.All these results completely verified the feasibility and effectiveness of the proposed method.展开更多
Let ξn-1<ξn-2 <ξn-2 <… < ξ1 be the zeros of the the (n -1)-th Legendre polynomial Pn-1(x) and - 1 = xn < xn-1 <… < x1 = 1 the zeros of the polynomial W n(x) =- n(n - 1) Pn-1(t)dt = (1 -x2)P&...Let ξn-1<ξn-2 <ξn-2 <… < ξ1 be the zeros of the the (n -1)-th Legendre polynomial Pn-1(x) and - 1 = xn < xn-1 <… < x1 = 1 the zeros of the polynomial W n(x) =- n(n - 1) Pn-1(t)dt = (1 -x2)P'n-1(x). By the theory of the inverse Pal-Type interpolation, for a function f(x) ∈ C[-1 1], there exists a unique polynomial Rn(x) of degree 2n - 2 (if n is even) satisfying conditions Rn(f,ξk) = f(∈ek)(1≤ k≤ n - 1) ;R'n(f,xk) = f'(xk)(1≤ k≤ n). This paper discusses the simultaneous approximation to a differentiable function f by inverse Pal-Type interpolation polynomial {Rn(f,x)} (n is even) and the main result of this paper is that if f ∈ C'[1,1], r≥2, n≥ + 2> and n is even thenholds uniformly for all x ∈ [- 1,1], where h(x) = 1 +展开更多
In this paper, we consider the simultaneous stabilization of discrete linear timevarying systems in an operator-theoretic framework. We give a criterion of simultaneous stabilization on account of one kind of strong r...In this paper, we consider the simultaneous stabilization of discrete linear timevarying systems in an operator-theoretic framework. We give a criterion of simultaneous stabilization on account of one kind of strong representation and also give a parametrization for all the simultaneous stabilizing controllers.展开更多
The mechanisms of the February 3, 1996 Lijiang main shock, Yunnan Province, are estimated by using the principle to inverse the mechanisms of two point sources simultaneously. The results are that the main shock of Li...The mechanisms of the February 3, 1996 Lijiang main shock, Yunnan Province, are estimated by using the principle to inverse the mechanisms of two point sources simultaneously. The results are that the main shock of Lijiang consists of two large ruptures, the time difference and the distance between the two ruptures are about 12 s (by the inversion) and about 26 km respectively. An extensional normal with strike-slip fault in about the north-south direction was formed by the first rupture, the mechanism of the second rupture is to be further studied. The method to inverse mechanisms of two point sources at the same time and the results obtained by directly analyzing P waveform records of the main shock are introduced, some related problems are also discussed. The Wuding earthquakes of October, 1995 and the Lijiang earthquake are considered to be the manifestation of the same dynamic process at different temporal and spatial points and the occurrence order of the two earthquakes is related to the direction of dynamics transmission.展开更多
A simultaneous inversion of earthquake relocation and three-dimensional crustal structure of P-wave velocity in central-western China (21癗~36癗, 98癊~112癊) were performed in this paper. The crustal P-wave velocity m...A simultaneous inversion of earthquake relocation and three-dimensional crustal structure of P-wave velocity in central-western China (21癗~36癗, 98癊~112癊) were performed in this paper. The crustal P-wave velocity model and earthquake relocation for this region are obtained using Pg and Sg phase readings of 9 988 earthquakes from 1992 to 1999 recorded at 193 seismic stations within central-western China by SPHYPIT90 and SPHREL3D90 programs. A lateral inhomogeneous structure of P-wave velocity in this region was obtained. Ob-vious contrast of P-wave velocities was revealed on both sides of active fault zones. Relocated epicenters of 6 459 events show clear lineation along active faults, which indicated a close correlation between seismicity and the active faults in this region. Focal depths of 82% relocated events ranged from 0 to 20 km, which is in good agreement with that from double-difference earthquake location algorithm.展开更多
基金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 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.
基金supported by the National Key R&D Program of China(Grant No.2023YFF0803203)the National Natural Science Foundation of China(Grant Nos.42288201,42404133)+1 种基金the Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology(Grant No.2022B1212010002)the China Postdoctoral Science Foundation(Grant No.2023T160639)。
文摘Simultaneous inversion of scattering and absorption attenuation is of great significance for investigating small-scale inhomogeneities and inelastic properties of the subsurface.However,applying this to complex geophysical issues is constrained by the costly computational requirements for simulations and inversions using existing methods.The coupling effects between the scattering and absorption coefficients cause parameter crosstalk artifacts in multi-parameter inversion,significantly increasing the nonlinearity of the inverse processes.This paper proposes a robust and effective simultaneous inversion method for scattering and absorption attenuation.The propagation of coda energy is modeled using the finite-element method based on the frequency-domain diffusion equation.We employ the truncated Gauss-Newton technique for the simultaneous estimation of scattering and absorption coefficients to decouple the two attenuation parameters during the inversion procedure.Nevertheless,the inversion may reach a local minimum if the energy frequency is inappropriate due to a phase mismatch between the simulated and measured energies.To mitigate this issue,we provide a frequency selection criterion that considers the substantial spectral shift of the energy density spectrum toward lower frequencies.Numerical examples using synthetic and experimental data indicate that our method can significantly reduce computational complexity and suppress two-parameter crosstalk without requiring a precise initial attenuation model.
文摘In this paper,obtainment of 3 by simultaneous inversion of configurations of C_(22),C_(23) in 2 via its ditosylate or dimesylate was unsuccessful,leading to elimination products instead.However, monotosylate of 2 can be readily obtained with exclusive regioselectivity,which in turn gave a single epoxide 7b Hydrolysis of 7b gave 3 and 2 is a ratio of 2:1,with a combined yield of 80%.The recovered 2 can of course be recycled again.
基金supported by the 973 Program of China(No.2011CB201104 and 2011ZX05009)the National Science and the Technology Major Project(No.2011ZX05006-06)
文摘Prestack seismic inversion methods adopt approximations of the Zoeppritz equations to describe the relation between reflection coefficients and P-wave velocity, S-wave velocity, and density. However, the error in these approximations increases with increasing angle of incidence and variation of the elastic parameters, which increases the number of inversion solutions and minimizes the inversion accuracy. In this study, we explore a method for solving the reflection coefficients by using the Zoeppritz equations. To increase the accuracy of prestack inversion, the simultaneous inversion of P-wave velocity, S-wave velocity, and density by using prestack large-angle seismic data is proposed based on generalized linear inversion theory. Moreover, we reduce the ill posedness and increase the convergence of prestack inversion by using the regularization constraint damping factor and the conjugate gradient algorithm. The proposed prestack inversion method uses prestack large-angle seismic data to obtain accurate seismic elastic parameters that conform to prestack seismic data and are consistent with logging data from wells.
文摘Simultaneous waveform inversion was used to predict lithofacies and fluid type across the field. Very often, characterizing reservoirs in terms of lithology and fluid type using conventional methods is replete with uncertainties, especially in marginal fields. An approach is employed in this study that integrated rock physics and waveform inverse modelling for lithology and fluid-type characterization to appropriately identify potential hydrocarbon saturated zones and their corresponding lithology. Seismic and well-log data were analyzed using Hampson Russel software. The method adopted includes lithofacies and fluid content analysis using rock physics parameters and seismic simultaneous inverse modelling. Rock physics analysis identified 2 broad reservoirs namely: HDZ1 and HDZ2 reservoirs. Results from the inverse modelling showed that low values of acoustic impedance from 19,743 to 20,487 (ft/s)(g/cc) reflect hydrocarbon-bearing reservoirs while medium to high values shows brine and shale respectively, with brine zone ranging from 20,487 to 22,531 (ft/s)(g/cc) and shale above 22,531 (ft/s)(g/cc). Two lithofacies were identified from inversion analysis of Vp/Vs and Mu-Rho, namely: sand and shale with VpVs 1.95 values respectively. Mu-Rho > 12.29 (GPa)(g/cc) and <12.29 (GPa) (g/cc) represent sand and shale respectively. From 3D volume, it was observed that a high accumulation of hydrocarbon was observed to be saturated at the north to the eastern part of the field forming a meandering channel. Sands were mainly distributed around the northeastern to the southwestern part of the field, that tends to be away from Well 029. This was also validated by the volume of rigidity modulus (Mu-Rho) showing high values indicating sands fall within the northeastern part of the field.
基金This study was supported by the project ofthe Science and Technology Innovation Fund of Command Center of Natural Resources Intergrated Survey entitled“Temporal and spatial distribution of paleochannel and origin of organic carbon burial in the Western Bohai Sea since 2.28Ma”(KC20220011)the project entitled“Characterization of Carboniferous-Early Permian heterogeneous porous carbonate reservoirs and hydrocarbon potential analysis in the central uplift of the South Yellow Sea Basin”(KLSG2304)+3 种基金by the Key laboratory of Submarine Science,Ministry of Natural Resources,the project entitled“1∶50000 Marine regional Geological survey in Caofeidian Sea Area,Bohai Sea”(ZD20220602)“1∶250000 Marine regional Geological survey in Weihai Sea Area,North Yellow Sea”(DD20230412)“Geological survey on tectonic and sedimentary conditions of Laoshan uplift”(DD2016015)by the China Geological Survey,and the project entitled“Study on Hydrocarbon Accumulation Failure and Fluid Evolution Reduction of the Permian Reservoir in the Laoshan Uplift,South Yellow Sea”(42076220)organized by the National Natural Science Foundation of China.
文摘Mesozoic-Palaeozoic marine carbonate rocks are crucial hydrocarbon reservoirs in the Central Uplift area of the South Yellow Sea Basin(SYSB).Due to the scarcity of boreholes and the significant heterogeneity of carbonate reservoirs,the distribution of porous carbonate reservoirs and their related key controlling factors remain unclear.In this study,factors affecting the distribution of porous Carboniferous-Early Permian carbonate reservoirs in the SYSB were investigated through seismic inversion and isotope analysis.The log-seismic characteristics of porous carbonate reservoirs,sensitive lithology parameters,and physical property parameters were extracted and analyzed.The pre-stack simultaneous inversion technique was applied to predict the lithology and physical properties of porous carbonate reservoirs.Moreover,the sedimentary of carbonate was analyzed using isotopes of carbon,oxygen,and strontium.The results show that porous carbonate reservoirs are mainly developed in the open platform sediments with porosities of 3%-5%and are mainly distributed in the paleo-highland(Huanglong Formation and Chuanshan Formation)and the slope of paleo-highland(Hezhou Formation).The porous carbonate reservoirs of the Qixia Formation are only locally developed.In addition,the negativeδ13C excursions indicate a warm and humid tropical climate with three sea-level fluctuations in the study area from the Carboniferous to Early Permian.The favorable conditions for developing porous carbonate rocks include the sedimentary environment and diagenetic process.The primary pore tends to form in high-energy environments of the paleo-highland,and the secondary pore is increased by dissolution during the syngenetic or quasi-syngenetic period.According to the hydrocarbon potential analysis,the Late Ordovician Wufeng Formation and Lower Silurian Gaojiabian Formation are the source rocks in the high-maturity-over-maturity stage,the Carboniferous-Lower Permian carbonate is the good reservoirs,and the Late Permian Longtan-Dalong Formation is the stable seal,ensuring a huge hydrocarbon accumulation potential in SYSB.The methods proposed in this study can be applied to other carbonate-dominated strata worldwide.
基金funded by National 973 Basic Research Developments Program of China (No.2005CB422104)863 National High Technique Research Development Project of China (No.2007AA060505)National Natural Science Foundation of China (No.40839901)
文摘Considering Zoeppritz equations, reflections of PP and PS are only the function of ratios of density and velocity. So the inversion results will be the same if the ratios are the same but values of density, velocities of P- wave and S-wave are different without strict constraint. This paper makes efforts to explore nonlinear simultaneous PP and PS inversion with expectation to reduce the ambiguity of AVO analysis by utilizing the redundancy of multi-component AVO measurements. Accurate estimation of ratio parameters depends on independence of input data. There are only two independent AVO attributes for PP reflectivity （i.e. intercept and gradient） and two for PS reflectivity （i.e. pseudo-intercept and pseudo-gradient or extreme amplitude）, respectively. For individual PP and PS inversion, the values of least-squares objective function do not converge around a large neighborhood of chosen true model parameters. Fortunately for joint PP and PS inversion the values of the least-squares objective function show closed contours with single minima. Finally the power function fitting is used to provide a higher precision AVO attributes than traditional polynomial fitting. By using the four independent fitting attributes （two independent attributes for PP and PS respectively）, the inversion of four ratio parameters （velocities and densities） would be estimated with less errors than that in traditional method.
基金the National Natural Science Foundation of China(41904116,41874156,42074167 and 42204135)the Natural Science Foundation of Hunan Province(2020JJ5168)the China Postdoctoral Science Foundation(2021M703629)for their funding of this research.
文摘As an important indicator parameter of fluid identification,fluid factor has always been a concern for scholars.However,when predicting Russell fluid factor or effective pore-fluid bulk modulus,it is necessary to introduce a new rock skeleton parameter which is the dry-rock VP/VS ratio squared(DVRS).In the process of fluid factor calculation or inversion,the existing methods take this parameter as a static constant,which has been estimated in advance,and then apply it to the fluid factor calculation and inversion.The fluid identification analysis based on a portion of the Marmousi 2 model and numerical forward modeling test show that,taking the DVRS as a static constant will limit the identification ability of fluid factor and reduce the inversion accuracy.To solve the above problems,we proposed a new method to regard the DVRS as a dynamic variable varying with depth and lithology for the first time,then apply it to fluid factor calculation and inversion.Firstly,the exact Zoeppritz equations are rewritten into a new form containing the fluid factor and DVRS of upper and lower layers.Next,the new equations are applied to the four parameters simultaneous inversion based on the generalized nonlinear inversion(GNI)method.The testing results on a portion of the Marmousi 2 model and field data show that dynamic DVRS can significantly improve the fluid factor identification ability,effectively suppress illusion.Both synthetic and filed data tests also demonstrate that the GNI method based on Bayesian deterministic inversion(BDI)theory can successfully solve the above four parameter simultaneous inversion problem,and taking the dynamic DVRS as a target inversion parameter can effectively improve the inversion accuracy of fluid factor.All these results completely verified the feasibility and effectiveness of the proposed method.
文摘Let ξn-1<ξn-2 <ξn-2 <… < ξ1 be the zeros of the the (n -1)-th Legendre polynomial Pn-1(x) and - 1 = xn < xn-1 <… < x1 = 1 the zeros of the polynomial W n(x) =- n(n - 1) Pn-1(t)dt = (1 -x2)P'n-1(x). By the theory of the inverse Pal-Type interpolation, for a function f(x) ∈ C[-1 1], there exists a unique polynomial Rn(x) of degree 2n - 2 (if n is even) satisfying conditions Rn(f,ξk) = f(∈ek)(1≤ k≤ n - 1) ;R'n(f,xk) = f'(xk)(1≤ k≤ n). This paper discusses the simultaneous approximation to a differentiable function f by inverse Pal-Type interpolation polynomial {Rn(f,x)} (n is even) and the main result of this paper is that if f ∈ C'[1,1], r≥2, n≥ + 2> and n is even thenholds uniformly for all x ∈ [- 1,1], where h(x) = 1 +
基金Supported by the National Natural Science Foundation of China(Grant Nos.1130104711271059)the Fundamental Research Funds for the Central Universities
文摘In this paper, we consider the simultaneous stabilization of discrete linear timevarying systems in an operator-theoretic framework. We give a criterion of simultaneous stabilization on account of one kind of strong representation and also give a parametrization for all the simultaneous stabilizing controllers.
文摘The mechanisms of the February 3, 1996 Lijiang main shock, Yunnan Province, are estimated by using the principle to inverse the mechanisms of two point sources simultaneously. The results are that the main shock of Lijiang consists of two large ruptures, the time difference and the distance between the two ruptures are about 12 s (by the inversion) and about 26 km respectively. An extensional normal with strike-slip fault in about the north-south direction was formed by the first rupture, the mechanism of the second rupture is to be further studied. The method to inverse mechanisms of two point sources at the same time and the results obtained by directly analyzing P waveform records of the main shock are introduced, some related problems are also discussed. The Wuding earthquakes of October, 1995 and the Lijiang earthquake are considered to be the manifestation of the same dynamic process at different temporal and spatial points and the occurrence order of the two earthquakes is related to the direction of dynamics transmission.
文摘A simultaneous inversion of earthquake relocation and three-dimensional crustal structure of P-wave velocity in central-western China (21癗~36癗, 98癊~112癊) were performed in this paper. The crustal P-wave velocity model and earthquake relocation for this region are obtained using Pg and Sg phase readings of 9 988 earthquakes from 1992 to 1999 recorded at 193 seismic stations within central-western China by SPHYPIT90 and SPHREL3D90 programs. A lateral inhomogeneous structure of P-wave velocity in this region was obtained. Ob-vious contrast of P-wave velocities was revealed on both sides of active fault zones. Relocated epicenters of 6 459 events show clear lineation along active faults, which indicated a close correlation between seismicity and the active faults in this region. Focal depths of 82% relocated events ranged from 0 to 20 km, which is in good agreement with that from double-difference earthquake location algorithm.
