一种复杂结构井约束三维密度层序反演方法
【专利摘要】本发明是物探中复杂结构井约束三维密度层序反演方法，利用地震或电法勘探得到标志界面的构造图，网格化后得到初始界面网格数据，将测井声波数据转换为测井密度数据，进行相对密度的3D重力反演得到每个计算点的各个深度层的相对密度值，计算点的每一个深度层的绝对密度背景值，确定井约束系数，实现井约束，计算每个深度计算点反演密度绝对值，迭代反演获得反演拟合异常数据、三维密度数据。本发明实现了井约束，提高了重力资料密度反演的可靠性，反演速度快，实测数据处理时间满足处理需要，实现了方法的实用化。
【专利说明】一种复杂结构井约束三维密度层序反演方法
【技术领域】
[0001]本发明涉及地球物理勘探技术，属于地面重力勘探法，具体是一种复杂结构井约束三维密度反演方法。
【背景技术】
[0002]目前使用的地面重力勘探方法，以石油行业重力勘探为例，执行的标准为《地面重力勘探技术规程(SY/T5819-2002)》，其对重力资料的反演处理技术要求均是2D\2.5D (二度\ 二度半)剖面反演、单界面重力界面反演而进行的，尚没有三维密度反演的规定和成熟推广的技术，更没有涉及复杂结构井约束的三维密度反演。
[0003]现在，三维密度反演技术尚处在研究阶段，姚长利公开了一种重磁遗传算法三维反演的方法(姚长利等，重磁遗传算法三维反演中高速计算及有效存储方法技术，地球物理学报，Vol.46, N0.2，2003.3，P252-258。)，美国犹它大学的日丹诺夫等也公开了类似的方法(Zhdanov M.S.Three-dimensional regularized focusing inversion of gravity tensorcomponent data, Geophysics, 2004.69(4), P925-937。),但是,上述方法使用的模型均是针对矿体的简单模型，反演密度也是相对密度，这样的反演方法和反演结果无法满足石油勘探对复杂三维反演的需要。根据查新了解到，对于适用于石油勘探的复杂地质结构、井约束的三维密度反演尚未见到文献报道。中国专利《基于标准格架的快速三维重力、磁力物性反演的方法》(200810119467.2),公开了一种基于标准格架正演、层序异常分离和层密度(层磁化率)快速近似计算和快速迭代反演的一种重磁反演方法，具有快速计算、处理实测数据的实用性。但是该方法不具备井约束反演功能，也不具备反演逆断层、反转密度(磁化率)的功能，因而它不适用于复杂区重力资料的三维反演。
【发明内容】
[0004]本发明目的在于提供一种利用井约束提高反演精度，解决复杂山前带地区石油重力勘探的复杂结构井约束三维密度层序反演方法。
[0005]本发明通过以下步骤实现:
[0006]I)野外实测采集重力数据，处理得到重力异常平面网格数据；
[0007]步骤I)所述的网格数据为重力异常数据。
[0008]2)利用地震或电法勘探得到标志界面的构造图，网格化后得到初始界面网格数据；
[0009]步骤2)所述的网格数据为构造图界面的深度数据。
[0010]步骤2)所述的初始界面通过探区已知地质资料、已知的地震构造图、电法勘探构造图结合已知钻井资料建立，或直接由反射面的地震构造图转换为需要的初始界面网格数据。
[0011]步骤2)所述的网格化后得到初始界面网格数据，当遇有逆掩断层则网格化为不同的初始界面网格数据文件。[0012]所述的不同的初始界面网格数据文件是断层上盘一侧的构造图界面深度数据网格化成一个文件，而断层下盘一侧的构造图界面深度数据网格化成为另一个文件，在逆掩断层带，两个网格文件的同一位置的数据分别为逆断层的上盘和下盘界面深度数据。
[0013]3)将测井声波数据转换为测井密度数据；
[0014]步骤3)所述的转换测井密度数据由测井声波速度(V)数据按Gardner公式转换获得测井的密度(σ )数据，转换得到的密度数据与测井声波速度数据一一对应，再按重力3D反演纵向剖分的间距对测井密度数据进行取均值。
[0015]步骤3 )所述的测井密度数据由地表向下排列，每一口钻井的密度数据个数与重力3D密度反演的纵向剖分数目相等；
[0016]当测井深度与反演深度一样时，反演纵向剖分为m个层，则测井密度对应求取m个层的密度平均值；当测井深度(η个层)小于反演深度(m个层)时，则测井密度对应求取I层至η层的各个层的密度平均值，(η+1)层至m层的密度数据可以取为该段深度地层的岩石密度值
[0017]所述的Gardner 公式为 o=0.31XVa25。
[0018]4)进行相对密度的3D重力反演，获得相对密度的三维数据体，得到每个计算点的各个深度层的相对密度值；
[0019]步骤4)所述的相对密度的3D重力反演采用基于标准格架的快速三维重力、磁力物性反演的方法进行。
[0020]5)计算点的每一个深度层的绝对密度背景值，先由距离计算点最近的2 口井测井密度数据，根据这2 口井 与计算点的距离反比进行插值，获得计算点与测井深度相同深度的绝对密度背景值，再根 据计算点初始界面深度与最近的I 口测井位置的初始界面深度差，进行计算点各个深度层背景密度值所应计存的深度层的浮动；
[0021]步骤5)所述的计算点各个深度层背景密度值所应计存的深度层的浮动方法是按最靠近计算点的测井初始界面深度值与计算点初始界面值之差进行浮动，计算点初始界面深度低则向下浮动，计算点初始界面深度浅则向上浮动；向上浮动时，如上浮X个剖分层，会引起底部缺少X个剖分层的背景密度，则这X个剖分层的背景密度全部取为第m-x层的数据。向下浮动时，如下浮I个剖分层，会引起顶部缺少I个剖分层的背景密度，则这I个剖分层的背景密度全部取为第y+l层的数据。
[0022]6)确定井约束系数k(r)，实现井约束；
[0023]步骤6)井约束系数k(r)计算如下:
【权利要求】
1. 一种复杂结构井约束三维密度层序反演方法，特点是通过以下步骤实现: 1)野外实测采集重力数据，处理得到重力异常平面网格数据； 2)利用地震或电法勘探得到标志界面的构造图，网格化后得到初始界面网格数据； 3)将测井声波数据转换为测井密度数据； 4)进行相对密度的3D重力反演，获得相对密度的三维数据体，得到每个计算点的各个深度层的相对密度值； 5)计算点的每一个深度层的绝对密度背景值，先由距离计算点最近的2口井测井密度数据，根据这2 口井与计算点的距离反比进行插值，获得计算点与测井深度相同深度的绝对密度背景值，再根据计算点初始界面深度与最近的I 口测井位置的初始界面深度差，进行计算点各个深度层背景密度值所应计存的深度层的浮动； 6)确定井约束系数k(r),实现井约束； 7)采用下式计算每个深度计算点反演密度绝对值: 计算点反演密度的绝对值=该点绝对密度的背景密度值+ (相对密度差值X(l-k(r))),k(r)取值为步骤6)计算得到； 相对密度差值为步骤4)中反演得到的计算点各个深度层的相对密度与最近位置测井的各个深度层的相对密度值之差； 8)迭代反演； 9)计算反演结果，获得反演拟合异常数据、三维密度数据。
2.根据权利要求1所述的方法，特点是步骤I)所述的网格数据为重力异常数据。
3.根据权利要求1所述的方法，特点是步骤2)所述的网格数据为构造图界面的深度数据。
4.根据权利要求1所述的方法，特点是步骤2)所述的初始界面通过探区已知地质资料、已知的地震构造图、电法勘探构造图结合已知钻井资料建立，或直接由反射面的地震构造图转换为需要的初始界面网格数据。
5.根据权利要求1所述的方法，特点是步骤2)所述的网格化后得到初始界面网格数据，当遇有逆掩断层则网格化为不同的初始界面网格数据文件。
6.根据权利要求5所述的方法，特点是所述的不同的初始界面网格数据文件是断层上盘一侧的构造图界面深度数据网格化成一个文件，而断层下盘一侧的构造图界面深度数据网格化成为另一个文件，在逆掩断层带，两个网格文件的同一位置的数据分别为逆断层的上盘和下盘界面深度数据。
7.根据权利要求1所述的方法，特点是步骤3)所述的转换测井密度数据由测井声波速度(V)数据按Gardner公式转换获得测井的密度(σ )数据，转换得到的密度数据与测井声波速度数据一一对应，再按重力3D反演纵向剖分的间距对测井密度数据进行取均值。
8.根据权利要求1所述的方法，特点是步骤3)所述的测井密度数据由地表向下排列，每一口钻井的密度数据个数与重力3D密度反演的纵向剖分数目相等； 当测井深度与反演深度一样时，反演纵向剖分为m个层，则测井密度对应求取m个层的密度平均值；当测井深度(η个层)小于反演深度(m个层)时，则测井密度对应求取I层至η层的各个层的密度平均值，(η+1)层至m层的密度数据可以取为该段深度地层的岩石密度值。
9.根据权利要求7所述的方法，特点是所述的Gardner公式为σ=0.31 XV°_25。
10.根据权利要求1所述的方法，特点是步骤4)所述的相对密度的3D重力反演采用基于标准格架的快速三维重力、磁力物性反演的方法进行。
11.根据权利要求1所述的方法，特点是步骤5)所述的计算点各个深度层背景密度值所应计存的深度层的浮动方法是按最靠近计算点的测井初始界面深度值与计算点初始界面值之差进行浮动，计算点初始界面深度低则向下浮动，计算点初始界面深度浅则向上浮动；向上浮动时，如上浮X个剖分层，会引起底部缺少X个剖分层的背景密度，则这X个剖分层的背景密度全部取为第m-x层的数据。向下浮动时，如下浮y个剖分层，会引起顶部缺少y个剖分层的背景密度，则这I个剖分层的背景密度全部取为第y+l层的数据。
12.根据权利要求1所述的方法，特点是步骤6)井约束系数k(r)计算如下:
13.根据权利要求1所述的方法，特点是步骤8)所述的迭代反演是按照迭代反演结束条件对三维反演绝对密度进行三维正演，获得正演重力异常值，再计算该正演重力异常值与实测重力异常值之差，形成拟合差数据，当该拟合差数据的最大差大于设置的反演结束最大差限差时，再根据拟合差数据进行迭代反演，即把该拟合差数据作为重力异常，按上述步骤4)到步骤8)进行循环计算，直到满足设置的迭代反演结束条件为止。
14.根据权利要求1所述的方法，特点是步骤8)所述的迭代反演结束条件是最大拟合差，等于2倍的采集误差。
【文档编号】G01V11/00GK103576212SQ201210250418
【公开日】2014年2月12日 申请日期:2012年7月19日 优先权日:2012年7月19日
【发明者】刘云祥 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司