[发明专利]一种用于地震反演的时控储层参数建模方法

摘要

本发明为一种用于地震反演的时控储层参数建模方法；本发明方法利用地震探区内常规测井数据和地震测网信息，在地质年代约束控制下建立标准层位的储层参数模型；与现有技术相比，本发明特别适用于井网密度不高、储层非均质性强的油气田，且具有计算速度快、稳定性好的优点，可直接用于地震资料的反演和解释等工作，为寻找岩性和地层圈闭油气藏和非常规油气藏，如煤层气、页岩气等，提供了重要资料。

主权项

一种用于地震反演的时控储层参数建模方法，其特征在于：所述方法的具体步骤为，步骤1，三维地震数据的预处理步骤：步骤1‑1，获取三维地震数据；步骤1‑2，从所述三维地震数据的卷头及道头中获取所述三维地震数据的范围参数；所述三维地震数据的范围参数包括纵测线的起始线号sl、纵测线的终止线号el、纵测线的测线间隔dl、联络测线的起始线号sc、联络测线的终止线号ec、联络测线的测线间隔dc、每道地震数据的起始时间st、每道地震数据的终止时间et、每道地震数据的时间间隔dt；步骤1‑3，根据所述三维地震数据的范围参数，通过公式(1)获取所述三维地震数据的纵测线数Lnum、联络测线数Cnum以及时间样点数Tnum；Lnum＝(el‑sl)/dlCnum＝(ec‑sc)/dc   (1)；Tnum＝(et‑st)/dt步骤2，获取地质年代层位面：步骤2‑1，获取标准层位线：由每条所述纵测线对应的地震剖面波形图中提取地震反射同相轴，并依次提取所述地震反射同相轴的波峰处相对应的时间值pt，并记录每条所述纵测线的线号pl和联络测线线号pc，由一组所述时间值pt和联络测线线号pc构成一组所述标准层位线；步骤2‑2，获取标准层位面：由所述步骤2‑1中各条所述纵测线上的标准层位线构成标准层位面，所述标准层位面即为地质年代层位面；步骤3，获取所述地质年代层位面上的井点投影坐标：通过钻井数据获取井轨迹数据，根据所述井轨迹数据与所述地质年代层位面的交点坐标，即所述地质年代层位面上的井点投影坐标，并获取所述井点投影坐标的纵测线线号agel、联络测线线号agec以及时间点aget；步骤4，获取所述地质年代层位面上的井点投影坐标的波阻抗值ageAI，其过程是：求取井中波阻抗曲线与所述地质年代层位面的交点，根据交点求得投影坐标的波阻抗值ageAI；步骤5，通过公式(2)获取所述地质年代层位面的相对坐标以及所述井点投影坐标的相对坐标；所述地质年代层位面的相对坐标包括所述地质年代层位面的相对纵测线坐标plr以及所述地质年代层位面的相对联络测线坐标pcr；所述井点投影坐标的相对坐标包括所述井点投影坐标的相对纵测线坐标agele以及所述井点投影坐标的相对联络测线坐标agecr；plr=(pl-sl)/dlpcr=(pc-sc)/dcagelr=(agel-sl)/dlagecr=(agec-sc)/dc---(2);]]>其中，sl,dl,sc,dc为所述步骤1‑2中获取的所述三维地震数据的范围参数；pl,pc为所述步骤2‑1中获取的所述标准层位线中每个样点的纵测线线号和联络测线线号；agel,agec为所述步骤3中获取的所述井点投影坐标的纵测线线号和联络测线线号；步骤6，利用所述步骤4中获取的所述井点投影坐标的波阻抗值ageAI以及所述步骤5中获取的所述井点投影坐标的相对坐标，并根据公式(3)进行插值处理，获取所述地质年代层位面的相对坐标处各样点的波阻抗值ageAIi,j；ageAIi,j=Σk=1n[wi,j·ageAIk]Σi=0LnumΣj=0Cnumwi,j=1---(3);]]>其中，n为所述地质年代层位面上井点投影的个数；ageAIk为所述地质年代层位面上各井点投影处的波阻抗值，k的取值范围是1～n；wi,j为所述地质年代层位面上各样点处的权重，其取值范围是0～1，i的取值范围是0～Lnum、j的取值范围是0～Cnum；步骤7，求取光滑渐变度R，其过程是：根据公式(4)建立储层参数模型，获取所述地质年代层位面的光滑渐变度R；R＝||Rl||2+||Rc||2   (4)；其中，Rl为所述地质年代层位面上纵测线方向的光滑渐变度，Rc为所述地质年代层位面上联络测线方向的光滑渐变度，||Rl||2和||Rc||2表示对Rl和Rc求平方和，然后再开二次方；||Rl||2和||Rc||2的值通过公式(5)获取：||Rl||2=Σj=0CnumΣi=0Lnum[ageAIi+1,j-ageAIi,jageAIi+1,j+ageAIi,j]2||Rc||2=Σi=0LnumΣj=0Cnum[ageAIi,j+1-ageAIi,jageAIi,j+1+ageAIi,j]2---(5);]]>其中，ageAIi,j指第i条纵测线、第j条联络测线处的波阻抗值，ageAIi+1,j指第i+1条纵测线、第j条联络测线处的波阻抗值；ageAIi+1,j‑ageAIi,j指相邻纵测线的波阻抗之差，ageAIi+1,j+ageAIi,j指相邻纵测线的波阻抗之和；ageAIi,j+1指第i条纵测线、第j+1条联络测线处的的波阻抗值，ageAIi,j+1‑ageAIi,j指相邻联络测线的波阻抗之差；ageAIi,j+1+ageAIi,j指相邻联络测线的波阻抗之和；步骤8，建立目标函数：根据公式(6)建立用于迭代反演的目标函数；f(w‾)=||ageAIi,j-ageAIk||2+α·R---(6);]]>其中，为所述地质年代层位面上各样点处的权重，其展开式为wi,j,i＝0～Lnum,j＝0～Cnum；ageAIi,j为所述步骤6中通过所述公式(3)获取的所述地质年代层位面的相对坐标处各点的波阻抗值；ageAIk为所述地质年代层位面上各井点投影处的波阻抗值；α为调整控制点处拟合与插值波阻抗光滑性和渐变性程度的参数，α的取值范围是[0,1]，R为通过所述公式(4)获取的所述地质年代层位面的光滑渐变度；步骤9，获取反演结果，其过程是：采用线性优化算法对所述公式(6)进行迭代反演，获取反演结果，即所述地质年代层位面上各样点处的权重步骤10，获取所述地质年代层位面上的储层参数模型，其过程是：通过所述步骤9获取的反演结果以及公式(3)，获取所述地质年代层位面上各样点处的波阻抗值ageAIi,j，从而获取所述地质年代层位面上的储层参数模型。