基于流体因子叠前地震反演的油水识别方法研究与应用

在双相介质理论指导下构建的流体因子Gassmann流体项f可以实现对油水层的有效判识。本文基于包含流体因子f的反射系数近似公式推导,得到相应的流体弹性阻抗方程,将流体因子纳入地震反演过程并直接参与目标储层的流体识别,实现了从弹性阻抗数据体中直接反演流体因子f的方法技术流程。流体因子f的直接反演有效减少间接计算所带来的累计误差,提高了地震流体识别的精度。在实际资料应用中,结合研究区地质背景与地球物理特征开展的岩石物理分析表明,Gassmann流体因子f对含油和含水特征有较好的区分作用。反演结果表明,流体因子f能为有效地识别目的层的油水特征,为下一步的钻井开发提供可靠的技术支持。      

一种裂缝流体因子的提出及应用

为了解决各向异性下的流体识别问题,将纵波各向异性裂缝预测以及Russell的流体因子融合到直角坐标系中,提出了一种能够同时检测裂缝发育情况以及流体性质的新的裂缝流体因子(Factor of Fluid-filled Fracture,FFF),并通过一组岩性参数检验了裂缝流体因子在裂缝预测及流体识别中的有效性.在理论研究的基础上,选取松辽盆地某地区的火成岩裂缝及流体识别研究为应用实例.通过与测井流体及裂缝信息的对比验证,裂缝流体因子能够较为准确地预测研究区裂缝和流体的分布情况,且裂缝流体因子在单井上的计算结果与单井含气饱和度吻合度较高.此外,根据实际应用效果,指出裂缝流体因子在应用中的局限性:裂缝流体因子在平面成图时受地层厚度影响较大,且无法预测裂缝方向.     

弱各向异性介质中P-SV波弹性阻抗与流体识别因子(英文)

基于Petr Jilek正交各向异性介质中的P-SV波反射系数精确公式,通过近似得到VTI和HTI介质中的P-SV波近似公式,同时结合Kenneth横波弹性阻抗原理,推导建立了VTI和HTI介质中P-SV波的弹性阻抗近似公式;在此基础上构建了弱各向异性介质中的流体识别因子公式,通过理论模型对Castagna和Smith的岩性组合进行了识别;最后详细分析了各向异性参数的取值对P-SV波弹性阻抗及识别因子产生的影响及变化趋势。     

基于叠前弹性阻抗反演的古近系流体识别

对于地震-地质条件十分复杂的古近系地层,Gassmann流体因子相对常规岩石物理参数通常具有更好的油气检测能力.然而,在由叠前反演结果计算流体因子的过程中往往会引入新的中间误差.为此,首先基于饱和流体多孔弹性介质理论推导了 Zoeppritz方程的Russell近似表达式,该方程将纵波反射系数表示为流体因子f、剪切模量μ以及密度ρ的函数.进而结合叠前弹性阻抗反演理论,实现了流体因子f的直接计算.该方法兼顾了弹性阻抗反演算法较强的抗噪能力和流体因子较强的油气识别能力,同时有效避免了传统叠前反演方法在间接计算流体因子过程中产生的累积误差.将该方法应用于渤海油田锦州A构造古近系油气检测,研究结果表明,流体因子f反演结果能够有效识别并清晰刻画出研究区古近系东营组强、中、弱等不同振幅对应的油气层.研究成果指导该区井位部署并得到钻探结果证实,为渤海油田古近系油气勘探开发提供了有力的技术支撑.     

针对深层流体识别的两项弹性阻抗反演与Russell流体因子直接估算方法研究

Russell流体因子是判识深层储层含流体特征的敏感指示参数.针对深层地震资料缺失大角度入射信息以及信噪比普遍较低的特点,开展了基于贝叶斯理论框架的两项弹性阻抗反演与Russell流体因子直接估算方法研究.首先,在饱和孔隙弹性介质理论的指导下,推导了包含Russell流体因子和纵波阻抗参数的两项反射系数近似公式,并且进一步得到了新的两项弹性阻抗方程.其次通过模型的近似精度分析,验证了新弹性阻抗方程满足小角度近似精度要求.最后提出了贝叶斯弹性阻抗反演和Russell流体因子的直接估算方法,通过模型试算表明直接估算方法不仅充分利用了弹性阻抗反演的高抗噪性与实用性,并且有效的避免了间接累计误差对流体因子计算精度的影响,实际应用则进一步验证了直接估算的Russell流体因子可以提高深层储层流体识别的可靠性.     

基于挤喷流效应的Russell流体因子推广及应用

在油气勘探与开发过程中,寻求能够从地震资料中直接识别储层油气的流体指示因子至关重要.由于现存多数流体指示因子都是在Biot理论假设前提下建立起来的,因此在双相孔隙介质条件下不能有效地识别流体.为此,本文基于前人提出的双孔介质统一波动理论,考虑岩石裂缝间挤喷流效应,在经典流体指示因子基础上构建了一种新的流体指示因子.该流体因子能够较好地反映岩石内孔隙流体的变化对波传播产生的影响.在实际资料的应用中,采用多种方法分析、对比该流体识别因子与以Biot理论为基础的传统流体因子的优劣,理论分析与实际资料的验证表明该流体因子对于储层中油气的检测有较高的精度和灵敏度.     

基于AVO反演的频变流体识别方法

研究表明流体引起衰减与频散往往表现为频变AVO现象.一些频散地震属性,例如纵波频散,已经证实为可靠的碳氢指示因子.为了更有效地识别流体,基于f-μ-ρ近似构建了新的流体因子D_f,即频变流体项.该属性的反演首先需要连续小波变换(CWT)谱分解得到不同频带地震数据,通过去相关与先验约束来保证反演结果可靠性.模型试算证实了频变反射系数近似公式的精度可靠性,D_f可以识别出强衰减介质所引起的频散现象.实际数据试算中,D_f可以较好地识别储层孔隙流体,尤其对于气层,具有较好的指示效果.该流体因子将Gassmann流体项的高孔隙流体敏感性与叠前数据丰富的振幅频率信息相结合,反演效果与岩石物理认识相符.此研究有助于利用衰减频散现象借助AVO反演实现流体识别.     

叠前同时反演在砂岩油藏预测中的应用

AVO技术已成为非常重要的烃类检测的工具,而以AVO理论为基础的叠前同时反演则日益成为常用的研究方法.本文首先从已钻井岩石物理分析入手,综合已钻井的AVO正演模拟和流体替代,分析了利用AVO技术对研究区目标层进行岩性和油水识别的可行性.利用叠前同时反演得到了纵波阻抗、横波阻抗和密度数据体,在此基础上计算得到纵横波速度比和泊松阻抗数据体,对研究区目标层进行了岩性和油气预测,取得了良好的应用效果.     

基于叠前地震纵横波模量直接反演的流体检测方法

流体因子是储层流体识别的重要方法,而叠前地震反演是获得流体因子的有效途径之一.本文从流体因子的构建出发,基于多孔弹性介质岩石物理模型,建立了流体因子与纵横波模量之间的直接关系,避免了流体因子计算所需的密度参数无法准确求取的问题,通过推导基于纵横波模量的Zeoppritz近似方程及弹性阻抗方程,探讨了基于弹性阻抗的纵横波模量直接反演方法,模型与实际应用表明,基于弹性阻抗的纵横波模量直接反演方法合理、可靠,减少了常规方法间接计算纵横波模量带来的累积误差,基于纵横波模量的流体因子计算方法有较好的实际应用效果.     

综合利用岩石模量与品质因子信息进行储层流体识别

以地震资料为主体的储层流体识别技术对现阶段油气勘探具有重要意义。本文利用岩石模量参数与品质因子作为储层流体指示因子,从提高流体指示因子估算可靠性的角度出发,分别研究了基于弹性阻抗反演的弹性模量参数直接估算方法与基于小波域QVO的品质因子估算方法。岩石模量参数直接估算方法在充分利用弹性阻抗反演高抗噪性与实用性的同时,有效地避免了间接累计误差的影响。品质因子的估算方法基于零相位子波假设,通过对小波域叠前品质因子的归零处理消除叠加平均效应。通过模型试算验证两种方法的流体判识可行性;并且结合实际资料进一步表明岩石模量与品质因子信息的综合应用可以提高储层流体的检测质量。     

On the Aerosol factor in the twilight method)

Summary It is proved by means of measuring series for skylight polarization during twilight carried out at different places that the zenith polarization is affected by the tropospheric turbidity still during the middle phase of twilight. This result is important for the inverse problem of the twilight theory. 3 approaches which explain this effect are examined. — Aerosol in the mesopause region is the cause for the well-known dip in the polarization curve during the middle phase of twilight. This interpretation can be supported by the distinct spectral differences.

---

Zusammenfassung An Meßreihen der Himmelslicht-Polarisation in der Dämmerung von verschiedenen Meßorten wird nachgewiesen, daß die Zenitpolarisation von der troposphärischen Trübung noch in der mittleren Dämmerungsphase beeinflußt wird. Dieses Ergebnis ist für die inverse Aufgabe der Dämmerungstheorie wichtig. 3 Möglichkeiten der Erklärung dieses Effekts werden untersucht. —Aerosol im Mesopausenbereich ist die Ursache für die bekannte Stufe in der Polarisationskurve während der mittleren Dämmerungsphase. Diese Deutung kann durch die deutlichen spektralen Unterschiede gestützt werden.

---

---

This work has been presented at the IUGG/WMO Symposium on Radiation including Satellite Techniques (Bergen/Norway, August 22–28, 1968).     

Scale factor for lengths of the geopotential model

Summary Using the data in[1], the scale factor for lengths is derived of the geopotential model R ₀ =GM/W ₀ (W ₀ is the potential on a generalized geoid). The resulting value, R ₀ ==6 363 672.9 m, which is2 m less than the original value[5], is practically the same as that in[6].     

On the period and quality factor of the Chandler wobble

The motion of the Earth’s geophysical pole is modeled by numerical integration of the Liouville equations. The dependence of variable parameters of these equations (the period of free nutation T and the quality factor of the shell Q) is examined as a function of time and initial data sets used for constructing geophysical perturbing functions. In particular, we used oceanic perturbing functions constructed from TOPEX/POSEIDON altimetry data. The annual and Chandler components of the pole motion were also modeled. Apart from the subtraction of quadratic trends, initial data were not filtered or averaged. The performed analysis provided model values of T = 425–440 days and Q = 20–60, and Q was found to be somewhat unstable with time.     

Estimation of the behavior factor of existing RC-MRF buildings

In recent years, several research groups have studied a new generation of analysis methods for seismic response assessment of existing buildings. Nevertheless, many important developments are still needed in order to define more reliable and effective assessment procedures. Moreover, regarding existing buildings, it should be highlighted that due to the low knowledge level, the linear elastic analysis is the only analysis method allowed. The same codes (such as NTC2008, EC8) consider the linear dynamic analysis with behavior factor as the reference method for the evaluation of seismic demand. This type of analysis is based on a linear-elastic structural model subject to a design spectrum, obtained by reducing the elastic spectrum through a behavior factor. The behavior factor (reduction factor or q factor in some codes) is used to reduce the elastic spectrum ordinate or the forces obtained from a linear analysis in order to take into account the non-linear structural capacities. The behavior factors should be defined based on several parameters that influence the seismic nonlinear capacity, such as mechanical materials characteristics, structural system, irregularity and design procedures. In practical applications, there is still an evident lack of detailed rules and accurate behavior factor values adequate for existing buildings. In this work, some investigations of the seismic capacity of the main existing RC-MRF building types have been carried out. In order to make a correct evaluation of the seismic force demand, actual behavior factor values coherent with force based seismic safety assessment procedure have been proposed and compared with the values reported in the Italian seismic code, NTC08.     

Determination of the geopotential scale factor from satellite altimetry

Summary The geopotential scale factor R ₀ =GM/W ₀ has been determined on the basis of satellite altimetry as R _{0=(6 363 672·5±0·3) m} and/or the geopotential value on the geoid W ₀ =(62 636 256·5±3) m ² s ^–2 . It has been stated that R ₀ and/or W ₀ is independent of the tidal distortion of surface W=W ₀ due to the zero frequency tide.

---

¶rt;a nmu amumuu u ama amnmuaa R ₀ =GM/W ₀ =(6 363 672,5±0,3) m u/uu aunmuaa a nmuu¶rt;a W ₀ =(62 636 256,5±3) m² s^–2. m, m R ₀ u/uu W ₀ auum m nm amu a a nuu ¶rt;au nmu W=W ₀ .

     

Longitudinal factor in the spatial distribution of energetic solar flares

Discrete location of large flares near zero Carrington longitude results in artificial breaks within physically related flare clusters. This effect is eliminated by using the data presentation algorithm, which results in the conclusion that drifting zones of intensified sunspot formation (concentrated cores of activity complexes) and energetic flare generation exist in either hemisphere. The indicated L zone is hypothetically related to a non-asymmetric relic solar magnetic field or to the regions where large-scale convective cells, extending to the convective zone bottom, originate.     

可变点约束叠前流体因子直接提取方法

以Gassmann流体因子(Gassmann Fluid Item,GFI)为目标,提出了一种流体因子直接提取的新方法.首先,以贝叶斯反演框架为基础,将似然函数、先验信息以及Gassmann流体因子近似方程相结合,得到初始的目标函数;其次,进一步在初始目标函数中加入可变数量的点约束信息,并得到最终的目标函数;最后,通过求解该目标函数,就直接提取出了Gassmann流体因子.该方法的主要特点是不需要初始模型的参与,而是通过一个约束模型来控制提取结果的稳定性和准确性,并且可以从约束模型中选定不同数量的约束点进行约束,称为可变点约束.给出并讨论了三种常用的不同点约束模式和原则,并用模型说明了它们不同的约束效果.模型验证和实际应用结果皆以表明,该方法即使在叠前数据信噪比很低的情况下也能较好地提取出Gassmann流体因子,流体因子提取结果客观性高、稳定性好,并且能够与已知的流体解释结果很好地匹配,益于进一步推广应用.     

Inversion of the asymmetry factor for desert areas of China

Deserts,which have high surface albedo and wide area,are important components of the earth system.It is very important for the research of surface radiation and energy balance to understand the anisotropic scattering of desert areas.The emergence and development of multi-angle remote sensing made possible the inversion of the anisotropic scattering of desert areas at the regional or global scale.Firstly,this paper explored the accuracy of the inversion of asymmetry factor using the Hapke model and the simul...     

The ductility reduction factor in the seismic design of buildings

This paper presents new trends in the relationship between the ductility reduction factor and the ductility demand in the seismic design of buildings. A total of 4860 inelastic time-history analyses were carried out to study this relationship using 60 single-degree-of-freedom models excited by an ensemble of 81 earthquake accelerogram records from around the world. The asymmetrical distribution of the results highlighted the inaccuracies associated with assuming a normal distribution simply described by the mean and standard deviation to represent the data. A probability of exceedence approach has been used based on counting the number of occurrences the ductility demand exceeds a specified level. The ductility reduction factors developed in this study are consistent with other studies in the long-period range but are different in the short-period range. The ductility reduction factor for very short period buildings of limited ductility has been found to be greater than previously predicted. © 1998 John Wiley & Sons, Ltd.     

Specifying the values of delta factor for the Dolgoe Ledovo permanent gravity station

Continuous long-term highly precise gravity measurements at the Dolgoe Ledovo gravity station operated by the Schmidt Institute of Physics of the Earth, Russian Academy of Sciences, are carried out. The long-term time series of tidal corrections are obtained, and the values of the delta factor are specified for this station, which enables the correct intercomparison of the absolute gravity instruments.