3-2孔隙度估算解析

3.2孔隙度估算3.2.1孔隙度地震预测的理论基础3.2.1.1沉积岩速度与孔隙度的关系3.2.1.2沉积岩的总体密度与孔隙度的关系3.2.2储层孔隙度定量解释中的直接公式计算3.2.2.1时间平均方程孔隙度预测3.2.2.2扩展时间平均方程法3.2.2.3密度法及密度速度联合预测法3.2.2.4阻抗法3.2.2.5纵横波速度联合预测法3.2.2.6统计分析法3.2.3地质-地球物理统计法3.2.3.1协克里金（Cokriging）法3.2.3.2多元统计分析3.2孔隙度估算储层的复杂性：孔隙度、渗透率、泥质含量、裂缝密度岩心直接测定法测井评价法地震法：地震波的速度地震波的时差地震波阻抗反演的密度资料优点：在井稀少处提前作出预测得到一个三维分布的数据体3.2.1.1沉积岩速度与孔隙度的关系（1）实验分析3.2.1孔隙度地震预测的理论基础mfpvvv11mfpttt1Wyllie（1956）时间平均方程，适应于10~25%之间的大多数固结的沉积岩；mfmttttpmfmCtttt1考虑疏松砂岩的压实校正3.2.1.1沉积岩速度与孔隙度的关系（1）实验分析3.2.1孔隙度地震预测的理论基础fmpvvv21BAv1mfmshmfmttttMttttRaymer（1980）非线性经验公式适应于更大的孔隙度范围Pickett（1963）纯经验公式，可以是纵横波速度；Angeler（1982）考虑了泥质含量对岩石速度的影响；3.2.1.1沉积岩速度与孔隙度的关系（1）实验分析3.2.1孔隙度地震预测的理论基础MvMvsp89.191.453.318.293.659.5Han（1986）分离了纵横波速度的回归拟合公式3.2.1.1沉积岩速度与孔隙度的关系（1）实验分析3.2.1孔隙度地震预测的理论基础Mvvsp43.056.055.136.116.1spvvpsmpmspsBBtttt速度比对孔隙度计算更为灵敏（Costagna，1985），对于水饱和的碎屑硅酸盐岩石纵横波速度几乎呈线性关系Ikwuakor（1988）纵横波时差与孔隙度成线性关系3.2.1.1沉积岩速度与孔隙度的关系（2）理论分析3.2.1孔隙度地震预测的理论基础Gassman-Biot-Geertsma方程Domenico（1977）结合墨西哥湾的一些具体岩石参数进行了系统研究。孔隙度是决定水饱和岩石纵横波速度的主要因素，但对泊松比的影响不明显；3.2.1.1沉积岩速度与孔隙度的关系（3）纵横波速度比值的意义3.2.1孔隙度地震预测的理论基础水饱和纵横波速度比值2.0：疏松砂岩2.0：固结良好砂岩1.96~1.42：水饱和砂岩1.66~1.30：气饱和砂岩3.2.1.2沉积岩的总体密度与孔隙度的关系3.2.1孔隙度地震预测的理论基础fm1Gregery（1976）提出了沉积岩的总体密度与孔隙度之间良好的线性关系：密度平均方程3.2.2.1时间平均方程孔隙度预测（1）预测流程3.2.2储层孔隙度定量解释中的直接公式计算3.2.2.1时间平均方程孔隙度预测（2）校正泥质校正（层状泥质、结构泥质、分散泥质）流体类型影响校正：校正饱含水、油、气时，时差增大的校正压实校正：对于疏松岩石和欠压实岩层其测量时差偏大的校正3.2.2储层孔隙度定量解释中的直接公式计算3.1~1ppjCC3.2.2.2扩展时间平均方程法3.2.2储层孔隙度定量解释中的直接公式计算岩性时间平均方程物性时间平均方程优点：把岩石放在宏观两相和微观三相的角度下观测，能更为全面地描述孔隙性岩石各组成部分对速度的影响，尤其强调了不能忽视泥质对总体速度的影响；在实际应用可以首先得到泥质含量，然后在求解孔隙度，提高计算的精度；3.2.2.3密度法及密度速度联合预测法3.2.2储层孔隙度定量解释中的直接公式计算shfmMM1mfmshmfmMmshfmshfMMtMtMtt11可同时解出孔隙度和泥质含量两个参数具有泥岩校正的密度平均方程：密度-速度联合求解：3.2.2.4阻抗法（DeBuyl，1986）利用地震反演波阻抗回避速度与密度之间的相互影响及其横向变化的影响3.2.2储层孔隙度定量解释中的直接公式计算mfmfshshmmttvtvMtvM1mshfmshfMMtMtMtv1113.2.2.5纵横波速度联合预测法3.2.2储层孔隙度定量解释中的直接公式计算Mvvsp43.056.055.1psmpmspsBBttttCostagna的实验统计关系Ikwuakor的实验统计关系B是与岩性、有效应力和颗粒接触面积或颗粒结构有关的常数3.2.2.6统计分析法3.2.2储层孔隙度定量解释中的直接公式计算BAv1Pickett（1963）的线性方程tba82.0,42.00078.0t3.0,044.000152.0t8.0,41.016.13.2.2.6统计分析法3.2.2储层孔隙度定量解释中的直接公式计算石灰岩孔隙度估算的精度要比砂岩高得多。Avt~~~~3.2.2.6统计分析法3.2.2储层孔隙度定量解释中的直接公式计算McbawMcbav222111MwMv25.151.315.382.133.637.5砂岩速度降低与孔隙度增加并不是单一的，还受岩性，特别是泥质含量的影响。AmosNur：甘利灯：3.2.3.1协克里金（Cokriging）法Doyen提出的一种地质统计方法利用地震参数和储层参数的空间相关性模拟储层特征的横向变化，用严密的数学方法把精确但稀疏的钻井资料与粗略但密集的地震资料综合在一起，可以改善储层参数平面描述的精度，同时还能提供一个可靠性估计。3.2.3地质-地球物理统计法3.2.3.1协克里金（Cokriging）法3.2.3地质-地球物理统计法tccmttmmjTjncjnniicicxtx11*xyxωxTmncncnccmnnnTxtxtxxyx,,,,,,,,,,,1111yωxmcxx**储层孔隙度在平面上任一位置的估算值可以用该位置周围相关半径以内的n口钻井孔隙度值和m个地震声波时差的线性加权组合来计算。即一种预测滤波的过程。3.2.3.1协克里金（Cokriging）法最小平方法计算加权系数3.2.3地质-地球物理统计法yωxxxxTccc*xx222cEmEωyyωωyxTTTcEEE222mntntnTcTCCCCExxxxxxxxyxs,,,,,11jninttjintjnitjiTCCCCExxxxxxxxyyxRmmrErErCccxxxxmcctttmrEtrErCxxxxmmrtErtErCtcctxxxxttccttmtmrtEtrtErCxxxx估算值与实际值的误差孔隙度的方差最小二乘意义3.2.3.1协克里金（Cokriging）法最小平方法计算加权系数3.2.3地质-地球物理统计法02ωxEsRω21221min2ωsxTE修改与扩展：（1）要有足够的井，平均井距应小于相关半径；（2）对于各向同性和均匀随机场特征的储层参数才适应；（3）储层参数为平稳随机场时其均值为零，对于区域性变化的应用非平稳随机场来模拟；（4）不同的储层参数应选用关系最为紧密的地震参数；孔隙度估算的相对可信度，即最小均方根误差3.2.3地质-地球物理统计法3.2.3地质-地球物理统计法3.2.3地质-地球物理统计法3.2.3地质-地球物理统计法3.2.3地质-地球物理统计法3.2.3.1协克里金（Cokriging）法小结：区域变量理论区域变量：空间从一点到另一点连续变化，但又不能用一个普通的工作函数所表达；空间变异性：空间一个结点上的区域变量样本值与周围结点上的变量有关，这种区域变量的连续变化的性质；半方差图：平均变异性的度量；相关图：在间隔一定距离和角度的两点之间的相似性；对空间随机变量的一种统计估计技术，利用相关性原则对物理测量数据进行内插外推的技术。3.2.3.2基于地震属性参数的多元统计分析算法3.2.3地质-地球物理统计法zyxAzyxAFzyxPm,,,,,,,,1mmAwAwwP110yxfwyxAwyxvwwyxP,,,,3210AWΦΦAAAWT1TN×(M+1)矩阵jjjjjjjjjjjiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiifAvfAfvffAfAvAAvfvAvvfAvN2223210A的选取？评价？它的地质意义和物理意义？