[发明专利]一种用于致密含气砂岩储层预测的测井曲线恢复方法

摘要

本发明为一种用于致密含气砂岩储层预测的测井曲线恢复方法，所述恢复方法是一种获取恢复曲线的方法；通过将趋势场因子引入目标测井曲线AC与待恢复曲线GR，并采用频谱整合技术对所述目标测井曲线AC与所述待恢复曲线GR的高频成分进行匹配整合，得到具有所述目标测井曲线的量纲特征、保留有高频成分的恢复曲线；本发明比较完整的保留了待恢复曲线的特征，能够明显的反映岩相、岩性的纵向空间变化，具有计算速度快、稳定性好、计算精度高以及抗噪能力强的特点，可直接用于地震资料的反演。

主权项

一种用于致密含气砂岩储层预测的测井曲线恢复方法，其特征在于：所述恢复方法是一种获取恢复曲线的方法；通过将趋势场因子引入目标测井曲线AC与待恢复曲线GR，并采用频谱整合技术对所述目标测井曲线AC与所述待恢复曲线GR的高频成分进行匹配整合，得到具有所述目标测井曲线的量纲特征、保留有高频成分的恢复曲线，所述恢复方法的具体步骤如下：步骤1，输入数据：输入所述目标测井曲线AC与所述待恢复曲线GR；步骤2，对输入数据进行预处理：去掉所述目标测井曲线AC以及待恢复曲线GR中没有物理意义的样点值；步骤3，求取所述目标测井曲线AC的趋势场因子LAC：步骤3‑1，建立最小二乘拟合方程，如公式(1)所示：f(a,b)=Σi=1n[ACi-yi]2=Σi=1n[ACi-a-b·di]2---(1);]]>其中，a和b为待定系数，ACi为所述目标测井曲线AC的第i个样点的测井数值，di为所述目标测井曲线AC的第i个样点的井深，yi为井深为di时的函数值，n为样点数；步骤3‑2，通过牛顿法或共轭梯度法求解最小二乘拟合方程，求取待定系数a和b，待定系数a和b是公式(1)取得最小值时的取值，即待定系数a和b是使目标测井曲线AC与趋势场因子LAC之间的拟合差取得最小时的取值；步骤3‑3，根据所述步骤3‑2获取的待定系数a、b以及di求取所述目标测井曲线AC的趋势场因子LAC，如公式(2)所示：LAC＝a+b·di     (2)；步骤4，求取所述待恢复曲线GR的趋势场因子LGR：步骤4‑1，建立最小二乘拟合方程，如公式(3)所示：f(a′,b′)=Σi=1n[GRi-y′i]2=Σi=1n[GRi-a′-b′·d′i]2---(3);]]>其中，a'和b'为待定系数，GRi为所述待恢复曲线GR的第i个样点的测井数值，d'i为所述待恢复曲线GR的第i个样点的井深，y'i为井深为d'i时的函数值，n为样点数；步骤4‑2，通过牛顿法或共轭梯度法求解最小二乘拟合方程，求取待定系数a'和b'，a'和b'是公式(3)取得最小值时的取值，即待定系数a'和b'是使待恢复曲线GR与趋势场因子LGR之间的拟合差取得最小时的取值；步骤4‑3，根据所述步骤4‑2获取的待定系数a'、b'以及d'i求取所述待恢复曲线GR的趋势场因子LGR，如公式(4)所示：LGR＝a'+b'·d'i     (4)；步骤5，求取所述目标测井曲线AC与待恢复曲线GR的高频成分：步骤5‑1，通过公式(5)求取所述目标测井曲线AC的高频成分HAC:HAC＝AC‑LAC    (5)；其中，LAC为所述步骤3‑3中求取的所述目标测井曲线AC的趋势场因子；步骤5‑2，通过公式(6)求取所述待恢复曲线GR的高频成分HGR：HGR＝GR‑LGR    (6)；其中，LGR为所述步骤4‑3中求取的所述待恢复曲线GR的趋势场因子；步骤6，获取所述目标测井曲线频谱FHAC以及所述待恢复曲线频谱FHGR；步骤6‑1：对所述目标测井曲线AC的高频成分HAC进行傅里叶变换，获取所述目标测井曲线频谱FHAC；步骤6‑2，对所述待恢复曲线GR的高频成分HGR进行傅里叶变换，获取所述待恢复曲线频谱FHGR；步骤7，获取所述目标测井曲线频谱FHAC以及待恢复曲线频谱FHGR的低频信号频谱FAC和FGR，并获取所述目标测井曲线AC的低通滤波频谱FAC的最大值λAC以及所述待恢复曲线GR的低通滤波频谱FGR的最大值λGR；步骤8，匹配整合频谱F：对所述目标测井曲线AC的低通滤波频谱FAC与所述待恢复曲线GR的低通滤波频谱FGR进行匹配整合处理，得到匹配整合频谱F，如公式(7)所示：F=FGR+λACλGRFAC---(7);]]>其中，λAC和λGR分别为所述步骤7获取的所述目标测井曲线AC的低通滤波频谱FAC的最大值以及所述待恢复曲线GR的低通滤波频谱FGR的最大值；步骤9，获取恢复曲线LACGR：对所述步骤8获取的匹配整合频谱F进行反傅里叶变换，得到匹配整合频谱的反傅里叶变换CGR，并根据公式(8)获取恢复曲线LACGR：LACGR＝λ·LAC+(1‑λ)·CGR      (8)；其中，λ为调节趋势场因子LAC和CGR在所述待恢复曲线中的比重，其取值范围为[0,1]。