《基于非线性算法的裂缝及裂缝内流体识别方法.pdf》由会员分享，可在线阅读，更多相关《基于非线性算法的裂缝及裂缝内流体识别方法.pdf（21页珍藏版）》请在专利查询网上搜索。
1、10申请公布号CN104018826A43申请公布日20140903CN104018826A21申请号201410153256622申请日20140416E21B49/0020060171申请人孙赞东地址100000北京市昌平区府学路18号中国石油大学（北京）72发明人孙赞东杨沛田军邓新良74专利代理机构北京国林贸知识产权代理有限公司11001代理人李桂玲54发明名称基于非线性算法的裂缝及裂缝内流体识别方法57摘要一种基于非线性算法的裂缝及裂缝内流体识别方法。本发明基于三项RGER反射系数方程，通过不同方位相减的办法推导出包含RGER方程第三项SIN2TAN2项和THOMSEN各向异性参数的非。
2、线性方程。利用非线性算法粒子群算法，从所推导的方程中求解出裂缝密度和THOMSEN各向异性参数等。利用求解出的裂缝密度、各向异性参数推导出了裂缝流体指示因子表达式。本发明在实际应用中，基于所提出的非线性方程得到的裂缝密度远远比常规线性AVOZ反演得到的裂缝密度接近理论值，反演的各向异性参数的值也非常接近理论值；利用非线性方程的解得到的裂缝流体指示因子也很好地识别出了裂缝流体类型，充分说明了该方法的可行性。51INTCL权利要求书2页说明书9页附图9页19中华人民共和国国家知识产权局12发明专利申请权利要求书2页说明书9页附图9页10申请公布号CN104018826ACN104018826A1/。
3、2页21一种基于非线性算法的裂缝及裂缝内流体识别方法，其特征在于，所述方法包括如下步骤步骤1把分方位共反射点道集分别记为方位1、方位2、方位3方位N，接着把N1个方位的共反射点道集分别减去某一个方位的共反射点道集，得到N1个方位的共反射点振幅差异道集，为非线性反演准备数据；步骤2利用N个方位的共反射点道集进行裂缝反演，算出每个反射点上每个时间采样点上的裂缝方位；步骤3把步骤2中得到的裂缝方位输入求初值的线性方程，得到求解非线性方程所需要的初值，所述初值包括裂缝密度、各向异性参数；步骤4利用步骤1所述的振幅差异道集，步骤2得到的裂缝方位、步骤3得到的裂缝密度和各向异性参数作为初始值，基于所述非线。
4、性方程进行反演，得到裂缝密度、各向异性参数；步骤5利用步骤4得到的裂缝密度、各向异性参数求得裂缝流体指示因子的大小，进行流体类型识别。2根据权利要求1所述的基于非线性算法的裂缝及裂缝内流体识别方法，其特征在于，所述的步骤2中所用的裂缝反演方法是基于两项RGER方程的常规线性AVOZ反演。3根据权利要求1所述的基于非线性算法的裂缝及裂缝内流体识别方法，其特征在于，所述步骤2中的每个反射点上每个时间采样点上的裂缝方位取自RGER反射系数方程。4根据权利要求1所述的基于非线性算法的裂缝及裂缝内流体识别方法，其特征在于，所述的步骤3的初值的求取是基于一个线性方程，而不是随机地给定非线性反演所需要的初始。
5、值。5根据权利要求1所述的基于非线性算法的裂缝及裂缝内流体识别方法，其特征在于，步骤3中所述求初值的线性方程是RPPI,BM0SIN2I5EM1SIN2ITAN2IDM2SIN2ITAN2I8式中，6根据权利要求1所述的基于非线性算法的裂缝及裂缝内流体识别方法，其特征在于，所述的步骤4的实施是基于正演信息的指导，是通过数值模型得到从三项RGER方程出发，通过不同方位相减的办法推导出一个非线性方程，通过求解该方程得到裂缝密度和各向异性参数。7根据权利要求1所述的基于非线性算法的裂缝及裂缝内流体识别方法，其特征在于，所述的步骤4中非线性方程的求解利用的是粒子群算法。8根据权利要求1所述的基于非线性。
6、算法的裂缝及裂缝内流体识别方法，其特征在于，步骤4中所述非线性方程是权利要求书CN104018826A2/2页3式中，ACOS2ICOS20；BSIN2ISIN20CCOS2ICOS20；DSIN2ISIN20BV8GR；EV；DV。9根据权利要求1所述的基于非线性算法的裂缝及裂缝内流体识别方法，其特征在于，所述的步骤5中裂缝流体指示因子，是从1995年SCHOENBERG和SAYERS推导的裂缝流体指示因子GN/T出发，结合各向异性参数和法向柔量N、切向柔量T的关系，推导出裂缝流体指示因子一个表达形式，而且数值模型的应用表明利用非线性方程的解得到的裂缝流体指示因子也很好地识别出了裂缝流体类型。
7、。权利要求书CN104018826A1/9页4基于非线性算法的裂缝及裂缝内流体识别方法技术领域0001本发明属于石油勘探领域，涉及一种裂缝储层综合评价中裂缝以及流体类型预测的方法，尤其涉及一种利用三项RGER方程的裂缝参数非线性反演方法。背景技术0002裂缝既是油气的储集空间，也是流体的运移通道，能够将储层中分散的孔洞储集空间连通起来。因此，识别裂缝和裂缝内流体对油气勘探具有重大意义。储层内发育裂缝时对外表现为各向异性特征，对于垂直排列的平行裂缝，可以将其等效为HTI介质HORIZONTALTRANSVERSEISOTROPY，即具有水平对称轴的横向各向同性介质模型来研究，常用THOMSEN参。
8、数表征HTI介质的各向异性特征THOMSEN，1995。描述这类垂直排列的平行裂缝的参数主要有裂缝密度和裂缝走向等。目前工业界大多用常规AVOZ方法进行裂缝探测，即基于两项RGER反射系数公式RGER，1998。实际上，三项RGER反射系数公式更能表征裂缝特征参数，并且，该公式中的第三项包含THOMSEN各向异性参数。下面给出RGER（1998）推导的反射系数公式0003RPP,ABISOBANICOS2SIN2CSIN2TAN210004其中，000500060007式中、分别为纵波速度和横波速度；V、V、是THOMSEN各向异性参数；为采集方位角；为裂缝走向；为入射角；Z为纵波垂直入射时的。
9、波阻抗；G为横波的切向模量。在（1）式中的第二项（即与SIN2I相乘的项）是AVO梯度，梯度项分为两部分各向同性梯度BISO和各向异性梯度BANI（该参数能够指示裂缝引起的各向异性程度即裂缝密度）。0008基于RGER近似公式的AVOZ反演方法一般只用其前两项，不用第三项含SIN2TAN2项。实际上，在用RGER反射系数近似公式描述地下某个成像点的各向异性特征时非常受第三项的影响，尤其是对入射角相对较大的情况。目前随着地震采集技术的提高，可以获得大偏移距数据，相应地，反射波入射角范围也较大。当入射角范围较大大于30时，RGER公式中的第三项SIN2TAN2对纵波反射系数的贡献也很大，不能忽略。。
10、因此，为了获得较精确的地下裂缝特征参数，如何充分利用RGER方程中的第三项是需要解决的一个问题。0009同时，裂缝流体也是目前地球物理工作者研究的热点和难点。根据SCHOENBERG和SAYERS1995的线性滑动理论，裂缝流体指示因子FLUIDINDICATOR可以写成如下形式说明书CN104018826A2/9页500100011BAKULIN等2000推导出各向异性参数和线性滑动理论中法向柔量N、切向柔量T的关系，如下公式00120013同时，根据RGER（1998）明确V和的关系V。0014据SCHOENBERG和DOUMA1988，在裂缝饱含不同流体类型时，法向柔量N和切向柔量T可以。
11、表示为00151对于干裂缝干的或者充填气，K0，则001600172对于湿裂缝充填油或者水，则有00180019其中即横纵波速度比的平方；E是HUDSON理论中的裂缝体密度。0020基于以上理论，虽然有很多人也对如何识别裂缝内流体类型做了研究，但是大多只是定性描述。发明内容0021本发明提出一种基于非线性算法的裂缝及裂缝内流体识别方法，它是基于三项RGER方程，该方法能够用于实际资料的裂缝以及裂缝内流体类型的识别。通过本发明提供的方法，能够识别出裂缝发育区，并能识别出裂缝内充填的流体类型，为更精细的裂缝储层描述提供依据。0022本发明实现上述目的的具体方案如下0023一种基于非线性算法的裂缝及。
12、裂缝内流体识别方法，其特征在于，所述方法包括如下步骤0024步骤1把分方位共反射点道集分别记为方位1、方位2、方位3方位N，接着把N1个方位的共反射点道集分别减去某一个方位的共反射点道集，得到N1个方位的共反射点振幅差异道集，为非线性反演准备数据；0025步骤2利用N个方位的共反射点道集进行裂缝反演，算出每个反射点上每个时间采样点上的裂缝方位；说明书CN104018826A3/9页60026步骤3把步骤2中得到的裂缝方位输入求初值的线性方程，得到求解非线性方程所需要的初值，这些初值包括裂缝密度、各向异性参数；0027步骤4利用步骤1所述的振幅差异道集，步骤2得到的裂缝方位、步骤3得到的裂缝密度。
13、和各向异性参数作为初始值，基于所述非线性方程进行反演，得到裂缝密度、各向异性参数；0028步骤5利用步骤4得到的裂缝密度、各向异性参数求得裂缝流体指示因子的大小，进行流体类型识别。0029所述的步骤2中所用的裂缝反演方法是基于两项RGER方程的常规线性AVOZ反演。0030所述步骤2中的每个反射点上每个时间采样点上的裂缝方位取自世界著名的RGER反射系数方程。0031所述的步骤3的初值的求取是基于一个线性方程，而不是随机地给定非线性反演所需要的初始值。0032所述的步骤4的实施是基于正演信息的指导，是通过数值模型得到从世界著名的三项RGER方程出发，通过不同方位相减的办法推导出一个新的非线性方。
14、程，通过求解该方程得到裂缝密度和各向异性参数。数值模型的应用表明基于非线性方程得到的裂缝密度远远比常规线性AVOZ反演得到的裂缝密度接近理论值，另外，反演的各向异性参数的值也非常接近理论值。0033所述的步骤4中非线性方程的求解利用的是粒子群算法。0034所述的步骤5中裂缝流体指示因子，是从1995年SCHOENBERG和SAYERS推导的裂缝流体指示因子GN/T出发，结合各向异性参数和法向柔量N、切向柔量T的关系，推导出裂缝流体指示因子一个表达形式，而且数值模型的应用表明利用非线性方程的解得到的裂缝流体指示因子也很好地识别出了裂缝流体类型。0035本发明与已有技术相比产生的有益效果是0036。
15、在反射波入射角较大的情况下30，基于两项RGER反射系数方程的线性AVOZ反演不能够精确描述地下裂缝诱导的各向异性特征。然而，目前随着野外采集技术的提高，地震资料的偏移距越来越大，也意味着入射角越来越大。00371、本发明基于三项RGER反射系数方程，通过不同方位相减的办法推导出包含RGER方程第三项SIN2TAN2项和THOMSEN各向异性参数的非线性方程。利用非线性算法粒子群算法，从所推导的方程中求解出裂缝密度和THOMSEN各向异性参数等。利用求解出的裂缝密度、各向异性参数推导出了裂缝流体指示因子的新的表达形式。00382、数值模型的应用表明基于所提出的非线性方程得到的裂缝密度远远比常规。
16、线性AVOZ反演得到的裂缝密度接近理论值，反演的各向异性参数的值也非常接近理论值；00393、利用非线性方程的解得到的裂缝流体指示因子也很好地识别出了裂缝流体类型。00404、将本发明应用到塔里木哈7井区奥陶系碳酸盐岩储层中，反演求出的裂缝参数与FMI结果吻合，且比常规线性AVOZ反演结果更精确，得到的裂缝流体因子也和井生产情况吻合，充分说明了该方法的可行性。说明书CN104018826A4/9页7附图说明0041图1是基于非线性算法的裂缝及裂缝内流体识别方法流程图；0042图2（A）是方位角道集图；0043图2（B）不同方位做差后的道集图；0044图3（A）是无噪音情况下非线性反演和常规AV。
17、OZ反演的裂缝密度BANI对比图；0045图3（B）是无噪音情况下非线性反演和常规AVOZ反演的各向异性参数V对比图；0046图3（C）是无噪音情况下非线性反演和常规AVOZ反演的各向异性参数V对比图；0047图4表示无噪音情况下第19个CDP点的目标函数值随迭代次数的变化情况图；0048图5（A）是含15噪音情况下非线性反演和常规AVOZ反演的裂缝密度BANI对比图；0049图5（B）是含15噪音情况下非线性反演和常规AVOZ反演的各向异性参数V对比图；0050图5（C）是含15噪音情况下非线性反演和常规AVOZ反演的各向异性参数V对比图；0051图6表示含15噪音情况下第19个CDP点的目。
18、标函数值随迭代次数的变化情况图。0052图7（A）是无噪音情况下裂缝流体因子的理论值、初始值、非线性反演结果的对比图；0053图7（B）是含15噪音情况下裂缝流体因子的理论值、初始值、非线性反演结果的对比图；0054图8是不同方位的偏移距道集图；0055图9是不同方位做差后的偏移距道集图；0056图10是测井储层参数综合评价图；0057图11是过井叠前时间偏移剖面图；0058图12是基于非线性反演的裂缝密度剖面图；0059图13是裂缝流体指示因子图。具体实施方式0060以下结合实施例与附图说明对本发明做进一步说明。0061图1是基于非线性算法的裂缝及裂缝内流体识别方法流程图。0062一种基于非。
19、线性算法的裂缝及裂缝内流体识别方法，包括如下步骤0063步骤1把分方位共反射点道集分别记为方位1、方位2、方位3方位6，接着把5个方位的共反射点道集分别减去某一个方位的共反射点道集，得到5个方位的共反射点振幅差异道集，为非线性反演准备数据；0064步骤2利用6个方位的共反射点道集进行裂缝反演，算出每个反射点上每个时间采样点上的裂缝方位（裂缝方位取自世界著名的RGER反射系数方程（1）；0065步骤3把步骤2中得到的裂缝方位输入本发明所推导的求初值的线性方程（8），说明书CN104018826A5/9页8得到求解非线性方程所需要的初值，这些初值包括裂缝密度BANI、各向异性参数V和V；0066步。
20、骤4利用步骤1的振幅差异道集，步骤2得到的裂缝方位、步骤3得到的裂缝密度和各向异性参数作为初始值，基于本发明推导的非线性方程（6）进行反演，得到裂缝密度BANI、各向异性参数V和V；0067步骤5利用步骤4得到的裂缝密度、各向异性参数求得裂缝流体指示因子的大小，进行流体类型识别。0068所述的步骤2中所用的裂缝反演方法是基于两项RGER方程的常规线性AVOZ反演。0069所述步骤2中的每个反射点上每个时间采样点上的裂缝方位取自世界著名的RGER反射系数方程。0070所述的步骤3的初值的求取是基于一个线性方程，而不是随机地给定非线性反演所需要的初始值。0071所述的步骤4的实施是基于正演信息的指。
21、导，是通过数值模型得到从世界著名的三项RGER方程出发，通过不同方位相减的办法推导出一个新的非线性方程，通过求解该方程得到裂缝密度和各向异性参数。数值模型的应用表明基于非线性方程得到的裂缝密度远远比常规线性AVOZ反演得到的裂缝密度接近理论值，另外，反演的各向异性参数的值也非常接近理论值。0072所述的步骤4中非线性方程的求解利用的是粒子群算法。0073所述的步骤5中裂缝流体指示因子，是从1995年SCHOENBERG和SAYERS推导的裂缝流体指示因子GN/T出发，结合各向异性参数和法向柔量N、切向柔量T的关系，推导出裂缝流体指示因子一个表达形式，而且数值模型的应用表明利用非线性方程的解得到。
22、的裂缝流体指示因子也很好地识别出了裂缝流体类型。0074本发明所提供的基于非线性算法的裂缝及裂缝内流体识别方法的基本原理如下0075从公式1出发，对于两个不同方位的反射地震数据，分别记为RPP,1和RPP,2，对两者做差，把这个差值记为RPP,，并且12，通过一系列的化简，可以得到公式6，00760077式中，0078ACOS2ICOS20；BSIN2ISIN200079CCOS2ICOS20；DSIN2ISIN200080BV8GR；EV；DV0081I，0均表示野外采集测线方位，0同时还表示被减去的方位；S为裂缝方位；表示垂直横纵波速度比的平方；V，V，R是THOMSEN各向异性参数。非线。
23、性方程6中的B表示常规线性AVOZ反演BANI，E，D分别是THOMSEN各向异性参数V和V。说明书CN104018826A6/9页9方程6中B，E，D和S是要求解的未知量。0082要求解方程6就要建立目标函数，本发明中将不同方位的叠前共反射点道集数据作为非线性反演的输入，利用公式6求解计算出的R,与实际观测的误差E作为目标函数00830084式中，ROBS,为实际观测的两个方位反射振幅之差；RCAL,是由公式6计算得到的理论值；I为第I个入射角和透射角的平均角度；J指第J个方位数据的平均方位角；T为迭代次数。通过求解非线性方程7可以得到描述裂缝特征的参数裂缝密度BANI各向异性强度、裂缝方位。
24、S，以及各向异性参数V和V。0085本发明使用一种新型的非线性反演方法粒子群优化算法PSO方法来求解方程7，这就需要B，E，D和S的初值。如果已知裂缝方位S本文中使用的S来自于AVOZ反演，把其代入方程6，那么可以得到线性方程8，如下0086RPPI,BM0SIN2I5EM1SIN2ITAN2IDM2SIN2ITAN2I80087式中，00880089式8可以写为，0090AXR90091这样，通过方程8求得的B，E，D作为求解目标函数7的初值。0092接着，结合方程3和方程6中各向异性参数的表示形式，可以把式2中裂缝流体指示因子简化成如下形式00930094公式10中的E，B，D均可以从方程。
25、7中求出。0095本发明中裂缝探测的适用范围如下对于无裂缝和裂缝走向呈现不同方向的井（多组裂缝），各向异性程度弱，建立在RGER方程基础上裂缝参数非线性反演方法无法预测出井附近裂缝发育区和井中的流体类型；只有对于裂缝走向呈现某一优势方向（单组裂缝）的井，其各向异性程度较强时，才能预测其中的裂缝和裂缝流体类型。0096图2用数值模型正演模拟说明不同方位道集做差后能量差异更明显。首先利用RGER公式（包含第三项）进行模型正演。设计一双层介质模型，上层为均匀各向同性介质层，下层为HTI介质层，模型参数如表1，并且HTI介质中充填气。同时，设定裂缝纵横比为0001，裂缝密度范围0005015，间隔是0。
26、005，即共有30个CDP点。利用HUDSON1981，1986等效介质理论和THOMSEN理论1986，求得各个CDP点的THOMSEN各向异性参数。设定入射角为15、25和35，方位角为0、30、60、90、120、150。利用公式说明书CN104018826A7/9页101得到反射系数，再将30个CDP点的反射系数序列分别与40HZ的零相位雷克子波进行褶积，得到的合成地震记录。图1A是第30个CDP点的方位角道集，图1B是第30个CDP点的0、30、60、120、150方位分别减去90方位后的道集，从图中可以看出做差后的不同道集之间的能量差异明显增大，这为下面进行非线性反演得到比较好的结。
27、果奠定了数学基础。0097表1各向同性/HTI双层单界面模型参数00980099图3是不含噪音的情况下非线性反演和常规AVOZ反演的结果对比情况，非线性反演结果中包括裂缝密度BANI、各向异性参数V和V，充分说明了非线性反演比常规AVOZ反演更精确。求解计算误差函数也是目标函数ET最小的具体步骤如下每个CDP点都以90方位作为被减方位，得到5个方位做差后的道集，作为目标函数的观测值；以合成的地震记录为输入，利用AVOZ反演求得裂缝方位S，求解过程详见SUN等2012A；把步骤得到的S输入方程8求出用于非线性反演的初始值E0，B0，D0；把通过步骤得到的S和步骤得到的E0，B0，D0作为目标函数。
28、7的初值进行PSO反演得到E，B，D和S，PSO反演的基本原理详见SUN等2012B。本发明所展示的均是反射界面处的结果。图3A是通过不同方法求得的裂缝密度BANI对比结果，黑色实线表示理论值，带圈实线表示线性反演的结果，即AVOZ反演的结果，绿色虚线表示初始值，蓝色点表示非线性反演的结果。可以看出，在裂缝内充填流体的情况下，用线性AVOZ反演的BANI值偏离理论值很多，尤其是裂缝密度较大的情况，在比较接近实际值的初始值的约束下，通过非线性反演得到的BANI值比较接近于理论值。图3B和图3C分别是通过非线性反演得到的各向异性参数V和V与理论值、以及初始值的结果，实线表示理论值，虚线表示初始值，。
29、蓝色点表示非线性反演的结果。从图3B和图3C中可以看出，非线性反演得到的V和V都比较接近于理论值，非线性反演结果相对于理论值的偏离程度小于初始值相对于理论值的偏离程度。0100图4表示无噪音情况下第19个CDP点的目标函数值随迭代次数的变化情况，这说明随着迭代次数的增加，观测值和计算值之间的误差越来越小，所求的解越来越接近理论值。0101图5是含15随机噪音的情况下非线性反演和常规AVOZ反演的结果对比情况，非线性反演结果中包括裂缝密度BANI、各向异性参数V和V，这不仅说明了非线性反演比常规AVOZ反演更精确，更说明了该方法有一定的抗噪音能力。其具体求解步骤同不含噪音的情况相同。图5A是通过。
30、不同方法求得的裂缝密度BANI对比结果，黑色实线表示理论值，带圈实线表示线性反演的结果，即AVOZ反演的结果，绿色虚线表示初始值，蓝色点表示非线性反演的结果。可以看出，所求得的裂缝密度BANI的整体变化特征和无噪音时候相同，由于有噪音的存在AVOZ反演结果、非线性反演结果和初始值均表现不同程度的抖动，说明书CN104018826A108/9页11相对于初始BANI，通过非线性反演得到的BANI更接近于理论值。图11和图12分别是通过非线性反演得到的各向异性参数V和V与理论值、以及初始值的结果，实线表示理论值，虚线表示初始值，蓝色点表示非线性反演的结果。从这两张图中可以看出，非线性反演得到的V和。
31、V都比较接近于理论值，并且非线性反演结果相对于理论值的偏离程度小于初始值相对于理论值的偏离程度。0102图6是第19个CDP点的目标函数值随迭代次数的变化情况，可见用本发明提出的非线性反演算法，在噪音存在的情况下目标函数仍然能够快速收敛，误差值也相对较小且合理，通过图4和图6对比可以看出由于噪音的存在，目标函数收敛所需要的迭代次数会相应地增大。本发明只显示了第19个CDP点的目标函数随迭代次数的变化情况，实际上不同的CDP点，其收敛情况是不同的，从数值模型的反演结果上看，用PSO方法求解所推导的非线性方程是可行的，结果也比较接近于理论值。0103图7同样的，在其他模型参数不变的条件下按着上面的。
32、步骤求出了裂缝内充填水包括无噪音和有噪音两种情况的裂缝密度BANI方程6中的B，各向异性参数V方程6中的E和V方程6中的D。通过公式10可以知道如果已知B、E、D就可以求出裂缝流体指示因子。通过公式10和11不难发现当裂缝内充填流体时裂缝流体指示因子近似为0，当裂缝内充填气时裂缝流体指示因子不为0，且其值远大于0，根据此可以用流体指示因子判断裂缝内的流体类型。图7A和图7B均是用非线性方程的解B、E、D求得的裂缝流体指示因子，其中图7A显示的是无噪音的结果，图7B是添加15的随机噪音的结果。可以看出当裂缝内饱含水时，流体指示因子较小接近于0，当裂缝充填气时，流体指示因子的值较大，通过裂缝流体指。
33、示因子二者可以明显被区分开来。通过图7B可看出，虽然受噪音的影响，初始值和非线性反演的结果均在一定程度上偏离理论值，但是这不影响利用流体因子区分裂缝流体类型。0104以下以该方法在新疆某区实际三维数据的一条测线为例具体说明该方法的实际应用效果。0105图8是不同方位的偏移距道集。其中每个玫瑰图中红色表示对应的采集方位。0106图9不同方位做差后的道集。它是由图8中方位2方位6分别减去方位1后得到的5个振幅差异道集。可以看出，相减后的振幅差异道集的信噪比明显降低，随机噪音比重大幅增加，但是在目标层3950MS附近，依然表现出一定的特征，经过相减后的道集之间的差异更明显。把不同方位的振幅差异道集按。
34、不同的入射角范围叠加作为非线性反演的输入。0107图10是研究井的测井储层参数综合评价图。本发明选择一条过井测线进行研究，目标层在一间房组和鹰山组一段。该井钻遇串珠群储层，它是由裂缝连通的溶蚀孔而成的。玫瑰图是由FMI资料统计，由玫瑰图显示该井张开缝发育，并且FMI显示张开缝倾角较陡，各向异性程度强，声成像解释结论表明该井段裂缝发育，测井评价结果也显示目标层内裂缝孔隙度比较高。0108图11表示过井的叠前时间偏移剖面。地震表现为串珠响应特征。0109图12表示利用本发明求出的过井裂缝密度。可以看出，利用本发明求出的裂缝密度在井的周围表现高值，即强各向异性特征，与井资料吻合。0110图13表示利。
35、用本发明求出的裂缝流体指示因子。利用非线性反演得到的各向异说明书CN104018826A119/9页12性参数V、V和裂缝密度BANI值，根据公式13得到流体指示因子。根据公式4和5可知，当裂缝内充填油或水时，流体指示因子的值较低，而充填气时，流体指示因子值较高。生产数据表明该井是口油井，井累计产液52197T，产油51392T，流体指示因子在井位置表现低值，与生产情况吻合，这表明利用本发明求的裂缝流体因子可以较好地识别裂缝流体类型。说明书CN104018826A121/9页13图1图2说明书附图CN104018826A132/9页14图3说明书附图CN104018826A143/9页15图4说明书附图CN104018826A154/9页16图5图6说明书附图CN104018826A165/9页17图7说明书附图CN104018826A176/9页18图8图9说明书附图CN104018826A187/9页19图10说明书附图CN104018826A198/9页20图11图12说明书附图CN104018826A209/9页21图13说明书附图CN104018826A21。