CN107167575B - 一种基于岩心的裂缝连续表征方法 - Google Patents

一种基于岩心的裂缝连续表征方法 Download PDF

Info

Publication number
CN107167575B
CN107167575B CN201710319039.3A CN201710319039A CN107167575B CN 107167575 B CN107167575 B CN 107167575B CN 201710319039 A CN201710319039 A CN 201710319039A CN 107167575 B CN107167575 B CN 107167575B
Authority
CN
China
Prior art keywords
crack
fracture
unit
core
cracks
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710319039.3A
Other languages
English (en)
Other versions
CN107167575A (zh
Inventor
赵文韬
荆铁亚
郜时旺
王金意
张健
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huaneng Clean Energy Research Institute
Huaneng Group Technology Innovation Center Co Ltd
Original Assignee
Huaneng Clean Energy Research Institute
Huaneng Group Technology Innovation Center Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huaneng Clean Energy Research Institute, Huaneng Group Technology Innovation Center Co Ltd filed Critical Huaneng Clean Energy Research Institute
Priority to CN201710319039.3A priority Critical patent/CN107167575B/zh
Publication of CN107167575A publication Critical patent/CN107167575A/zh
Application granted granted Critical
Publication of CN107167575B publication Critical patent/CN107167575B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Food Science & Technology (AREA)
  • Analytical Chemistry (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Medicinal Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

一种基于岩心的裂缝连续表征方法,基于实测岩心裂缝中心深度、长度、开度、充填程度及充填物、倾角及岩心直径等关键参数,确定测量精度,设计基于单位岩心长度和相关参数的裂缝发育程度计算方法,进而对裂缝发育程度进行连续定量表征,提高岩心裂缝密度和面缝率的计算精度,为裂缝定量预测提供可靠的数据资料保证;本发明将实测岩心裂缝参数与实际刻画精度相结合,可连续表征裂缝在垂向上的分布情况和集中程度,能更好地满足油气勘探开发、特别是裂缝型和非常规油气勘探开发的实际需求。

Description

一种基于岩心的裂缝连续表征方法
技术领域
本发明属于油气勘探开发技术中的储层裂缝表征技术,特别涉及一种基于岩心的裂缝连续表征方法。
背景技术
裂缝作为石油和天然气的一种重要储存空间和运输通道,无论是在流体储集还是在油气运移方面都起着至关主要的作用。因此,裂缝评价、尤其是在孔隙度和渗透率都极低的非常规油气储层中的裂缝评价,成为优选油气勘探有利区以及制定下一步开发方案的一项重要工作。裂缝的有效、精确表征是裂缝评价的必要前提,只有建立起一套可方便实施、可定量描述、可纵横对比的裂缝表征方案,才能为后期的储层评价与生产决策提供正确指导。
现有的岩心裂缝表征方法中,最常用的是利用裂缝线密度(即单位岩心长度内裂缝发育的总条数)对裂缝发育程度进行描述。但是,线密度描述方法具有明显的劣势:其一,裂缝线密度仅考虑了裂缝的条数而未考虑裂缝倾角、裂缝性质等重要参数对裂缝表征的影响;其二,裂缝线密度的单位统计长度通常为某地层厚度,长度较长,从而导致计算精度较低。因此,在岩心观测中采用线密度对裂缝进行表征,不能满足油气勘探高精度定量描述的实际需要。裂缝面密度(即单位岩心面积内裂缝发育的总长度)描述方法在一定程度上提高了表征精度,但是,裂缝面密度测量通常是基于单个岩心进行裂缝统计,测量较为繁琐,且受岩心长度这一非地质因素影响较大,计算精度偏低,垂向上也不具备连续性,纵、横向可对比性较差。因此,有必要形成一套方便岩心实测、测量精度统一、纵向可连续表征、横向可量化对比的岩心裂缝表征新方法,以满足油气前期勘探的实际需求。
该方法在国内首次针对岩心层段提出了裂缝连续表征方法,简化了测试流程,保证了参数的完整性,能够在岩心实测的基础上垂向连续、定量、统一地反映裂缝集中程度,提高了裂缝指标的可对比性。
发明内容
为了克服现有方法的不足,本发明的目的在于提供一种基于岩心的裂缝连续表征方法,基于实测岩心裂缝中心深度、长度、开度、充填程度及充填物、倾角及岩心直径等关键参数,确定测量精度,设计基于单位岩心长度和相关参数的裂缝发育程度计算方法,进而对裂缝发育程度进行连续定量表征,提高岩心裂缝面密度和面缝率计算精度,为裂缝定量预测提供可靠的数据资料保证。
为了达到上述目的,本发明的技术方案如下:
一种基于岩心的裂缝连续表征方法,包括以下步骤:
步骤一、裂缝类型识别:裂缝的有效识别是裂缝准确统计的基本前提,主要是在岩心观测中对不同成因和性质的裂缝进行区分;裂缝整体划分为人工裂缝和天然裂缝两类,在裂缝观测时要排除人工裂缝对裂缝表征的干扰;天然裂缝又分为构造裂缝和非构造裂缝两类:构造裂缝主要包括张性裂缝、剪性裂缝和压性裂缝,而非构造裂缝则主要包括溶蚀缝、压实缝、风化缝、层间缝和沉积裂缝;
裂缝类型按以下特点进行识别:
人工裂缝指由于人工施钻所形成的相关裂缝,其形态在很大程度上取决于压裂井附近的现代地应力状态及钻孔方向;在未压裂的井段中,人工裂缝通常仅见于岩心表面,很难贯穿至岩心内部或切穿岩心,且形态与钻头形状具有明显相似性,因此其能够与天然裂缝有效区分;
天然裂缝中,构造裂缝往往与区域构造作用密切相关,断裂或褶皱附近的构造裂缝通常较为发育,成组出现、走向稳定、平直、穿一层或多层,且剪性裂缝的延伸性和平整度显著优于张性裂缝和压性裂缝;非构造裂缝成因复杂,形态多样,尺度跨度和走向变化大,与构造作用无关:溶蚀缝指成岩后地层中的部分物质在酸性流体作用下溶蚀形成的裂缝,通常在裂缝面上可见明显矿物溶蚀特征;而由于地层存在一些封闭的含液成岩缝和气孔,在垂向的压实作用下,液体压力上升会导致在侧向压力较小的方向扩展形成裂缝,这类裂缝则称之为压实缝,往往呈头盖骨“缝合线”状;风化缝指由于表层风化、剥蚀和盖层沉积,在风化壳和不整合面上所形成的裂缝,横截面多呈“V”字型,切割岩层较浅,通常不切穿岩层,走向不稳定;层间缝为在相邻地层间,由于层面间不整合或沉积条件不同所形成的层间裂缝群,通常倾角变化大,与地层面平行,走向离散,不切穿地层;而由于流体的定向流动和冲击作用,砾石、岩屑会呈定向排列,形成液体流动通道,进而形成沉积裂缝,并与沉积构造特征密切相关;
步骤二、裂缝参数统计:基于不同类型裂缝的识别,统计各条裂缝的相关参数,相关参数主要包括裂缝性质、裂缝中心深度Fmid,单位m、裂缝长度Lfrac,单位m、裂缝开度w,单位mm、充填程度及充填物、裂缝倾角dip和岩心直径d,单位m;其中,裂缝中心深度Fmid测量不便时,首先测量裂缝的顶深Ftop,单位m和底深Fbtm,单位m,然后利用公式:
Figure BDA0001289283220000031
间接计算得到裂缝中心深度Fmid;裂缝长度Lfrac指裂缝在岩心上截取的直线长度,通常表示为裂缝椭圆截面的长轴长度;裂缝开度w指裂缝的张开程度,亦即裂缝两侧的间隔距离;充填程度反映裂缝内的充填情况,主要包括全充填、半充填和未充填三类;裂缝倾角dip指裂缝面与岩心横截面间的夹角;裂缝长度Lfrac、裂缝顶深Ftop和底深Fbtm、裂缝倾角dip之间满足关系式:
Fbtm-Ftop=Lfrac·sin(dip) (2)
步骤三、表征精度确定:表征精度指进行裂缝计量的最小岩心单元,用每次裂缝表征的单位岩心长度Lunit表示;表征精度为均匀划分或根据实际计量需求按地层、岩性、沉积相等非均匀划分;单位岩心顶深和底深分别为Utop,单位m和Ubtm,单位m,且满足以下关系式:
Lunit=Ubtm-Utop (3)
步骤四、裂缝密度表征:裂缝密度反映裂缝的集中程度,有线密度、面密度和体密度之分;线密度指一条单位长度的测线上通过的裂缝总条数,面密度指单位面积内所有裂缝的总长度,而体密度指在单位体积内裂缝穿切的总面积;本发明按面密度进行统计,即:
Figure BDA0001289283220000041
式中,Ds代表裂缝面密度,S代表测量面面积,li代表该测量面内各条裂缝的长度;将各单位岩心段的纵截面作为测量面,则上式可变换为:
Figure BDA0001289283220000042
针对每条裂缝,由于裂缝的对应深度与单位岩心段深度并非完全一致,二者存在多种交互关系,而实际仅需统计在该单位岩心段内的裂缝长度,因此上式进一步转换为:
Figure BDA0001289283220000043
Ds(unit)为某一单位岩心内的裂缝面密度;Lreal(unit)代表各条裂缝在该单位岩心上截取的实际长度,其进一步表征为:
Figure BDA0001289283220000051
式中,min(Ubtm,Fbtm)代表某单位岩心底深与裂缝底深二者中的较小值,max(Utop,Ftop)代表某单位岩心顶深与裂缝顶深二者中的较大值;max[0,t]代表对t取正值,而利用max[0,min(Ubtm,Fbtm)-max(Utop,Ftop)]则能够计算各条裂缝在该单位岩心内的裂缝垂向长度,再除以sin(dip)即得该单位岩心内的真实裂缝长度Lreal(unit);根据公式易知,当裂缝与单位岩心段不相交时,Lreal(unit)=0;当裂缝与单位岩心段部分相交时,Lreal(unit)≤Lfrac;而当裂缝完全在单位岩心段内时,Lreal(unit)=Lfrac;当某段岩心的裂缝面密度值Ds(unit)较高时,反映该段裂缝较为发育;而当某段岩心的裂缝面密度值Ds(unit)较低时,则反映该段裂缝发育能力不足;
步骤五、裂缝面缝率表征:裂缝面缝率(φ)指测量面内裂缝所占面积与单位测量面积的比值,计算公式为:
Figure BDA0001289283220000052
公式中,S代表岩心测量面,Si代表测量面内每条裂缝所占面积;由于每个裂缝面积均表示为裂缝长度li与裂缝宽度wi之积,因此上式进一步转换为:
Figure BDA0001289283220000053
将岩心测量中各单位岩心段的纵截面作为测量面,则上述公式表述为:
φ(unit)为某单位岩心段内的裂缝面缝率;(li·wi)代表测量面内每条裂缝所占面积;同裂缝面密度计算方法一致,上式进一步演变为:
Figure BDA0001289283220000061
式中,Sreal(unit)代表在该单位岩心段内各裂缝截面积,其详述为:
Figure BDA0001289283220000062
当某段岩心的裂缝面缝率值φ(unit)较高时,反映该段裂缝能提供的孔隙度较高;而当某段岩心的裂缝面缝率值φ(unit)较低时,则指示该段裂缝所提供的岩石孔隙相对不足;
步骤六、构造/非构造裂缝比:用于反映某井裂缝的主要成因;在此采用构造裂缝面密度与非构造裂缝面密度的比值进行表征,以反映裂缝以构造成因为主还是以非构造成因为主;根据步骤四,构造裂缝面密度表示为:
Figure BDA0001289283220000063
而非构造裂缝面密度则表示为:
Figure BDA0001289283220000064
两式相除,即得:
Figure BDA0001289283220000065
当某段岩心的构造/非构造裂缝比Rstr/nonstr高于1时,反映该段岩心裂缝主要以构造成因为主;而当某段岩心的构造/非构造裂缝比Rstr/nonstr低于1时,反映该段岩心裂缝则主要以非构造成因为主。
本发明基于岩心的裂缝连续表征方法,将实测岩心裂缝参数与实际刻画精度相结合,可连续表征裂缝在垂向上的分布情况和集中程度,能更好地满足油气勘探开发、特别是裂缝型和非常规油气勘探开发的实际需求。
附图说明
图1为本发明中的岩心裂缝连续表征方法流程图。
图2为本发明中的岩心裂缝相关参数示意图。
图3为本发明中的表征精度示意图。
图4为本发明中的裂缝面密度表征示意图。
图5为本发明中的裂缝长度与单位岩心关系示意图,其中:图5a为贯穿单位岩心的裂缝长度示意图,图5b为穿插单位岩心的裂缝长度示意图,图5c为单位岩心内的裂缝长度示意图,图5d为单位岩心外的裂缝长度示意图。
具体实施方式
如图1所示,本发明为一种基于岩心的裂缝连续表征方法,包括以下步骤:
步骤一、裂缝类型识别:裂缝的有效识别是裂缝准确统计的基本前提,主要是在岩心观测中对不同成因和性质的裂缝进行区分。裂缝整体划分为人工裂缝和天然裂缝两类,在裂缝观测时要排除人工裂缝对裂缝表征的干扰。天然裂缝又可以分为构造裂缝和非构造裂缝两类:构造裂缝主要包括张性裂缝、剪性裂缝和压性裂缝,而非构造裂缝则主要包括溶蚀缝、压实缝、风化缝、层间缝和沉积裂缝等。
人工裂缝指由于人工施钻所形成的相关裂缝,其形态在很大程度上取决于压裂井附近的现代地应力状态及钻孔方向。在未压裂的井段中,人工裂缝通常仅见于岩心表面,很难贯穿至岩心内部或切穿岩心,且形态与钻头形状具有明显相似性,因此其能够与天然裂缝有效区分。
天然裂缝中,构造裂缝往往与区域构造作用密切相关,断裂或褶皱附近的构造裂缝通常较为发育,一般成组出现、走向稳定、平直、穿一层或多层,且剪性裂缝的延伸性和平整度显著优于张性裂缝和压性裂缝。非构造裂缝成因复杂,形态多样,尺度跨度和走向变化大,与构造作用无关:溶蚀缝指成岩后地层中的部分物质在酸性流体作用下溶蚀形成的裂缝,通常在裂缝面上可见明显矿物溶蚀特征;而由于地层存在一些封闭的含液成岩缝和气孔,在垂向的压实作用下,液体压力上升会导致在侧向压力较小的方向扩展形成裂缝,这类裂缝则称之为压实缝,往往呈头盖骨“缝合线”状;风化缝指由于表层风化、剥蚀和盖层沉积,在风化壳和不整合面上所形成的裂缝,横截面多呈“V”字型,切割岩层较浅,通常不切穿岩层,走向不稳定;层间缝为在相邻地层间,由于层面间不整合或沉积条件不同所形成的层间裂缝群,通常倾角变化大,与地层面平行,走向离散,不切穿地层;而由于流体的定向流动和冲击作用,砾石、岩屑会呈定向排列,形成液体流动通道,进而可形成沉积裂缝,并与沉积构造特征密切相关。
步骤二、裂缝参数统计:基于不同类型裂缝的识别,统计各条裂缝的相关参数,相关参数主要包括裂缝性质、裂缝中心深度Fmid(m)、裂缝长度Lfrac(m)、裂缝开度w(mm)、充填程度及充填物、裂缝倾角dip(°)和岩心直径d(m)(图2)。其中,裂缝中心深度测量不便时,可首先测量裂缝的顶深Ftop(m)和底深Fbtm(m),然后利用公式:
Figure BDA0001289283220000081
间接计算得到;裂缝长度(Lfrac)指裂缝在岩心上截取的直线长度,通常可表示为裂缝椭圆截面的长轴长度;裂缝开度(w)指裂缝的张开程度,亦即裂缝两侧的间隔距离;充填程度可反映裂缝内的充填情况,主要包括全充填、半充填和未充填三类;裂缝倾角指裂缝面与岩心横截面间的夹角。通常而言,裂缝长度(Lfrac)、裂缝顶深(Ftop)和底深(Fbtm)、裂缝倾角(dip)之间满足关系式:
Fbtm-Ftop=Lfrac·sin(dip). (2)
步骤三、表征精度确定:表征精度指进行裂缝计量的最小岩心单元,可用每次裂缝表征的单位岩心长度(Lunit)表示(图3)。表征精度既可均匀划分,也可根据实际计量需求按地层、岩性、沉积相等非均匀划分。单位岩心顶深和底深分别为Utop(m)和Ubtm(m),且满足以下关系式:
Lunit=Ubtm-Utop. (3)
步骤四、裂缝密度表征:裂缝密度可反映裂缝的集中程度,有线密度、面密度和体密度之分。线密度指一条单位长度的测线上通过的裂缝总条数,面密度指单位面积内所有裂缝的总长度,而体密度指在单位体积内裂缝穿切的总面积。本发明按面密度进行统计,即:
Figure BDA0001289283220000091
式中,Ds代表裂缝面密度,S代表测量面面积,li代表该测量面内各条裂缝的长度。将各单位岩心段的纵截面作为测量面,则上式可变换为(图4):
针对每条裂缝,由于裂缝的对应深度与单位岩心段深度并非完全一致,二者存在多种交互关系,而实际仅需统计在该单位岩心段内的裂缝长度,因此上式可进一步转换为(图5):
Figure BDA0001289283220000093
Ds(unit)为某一单位岩心内的裂缝面密度;Lreal(unit)代表各条裂缝在该单位岩心上截取的实际长度,其可进一步表征为:
Figure BDA0001289283220000101
式中,min(Ubtm,Fbtm)代表某单位岩心底深与裂缝底深二者中的较小值,max(Utop,Ftop)代表某单位岩心顶深与裂缝顶深二者中的较大值。max[0,t]代表对t取正值,而利用max[0,min(Ubtm,Fbtm)-max(Utop,Ftop)]则可计算各条裂缝在该单位岩心内的裂缝垂向长度,再乘以sin(dip)即可得该单位岩心内的真实裂缝长度Lreal(unit)(图5)。根据公式易知,当裂缝与单位岩心段不相交时,Lreal(unit)=0,如图5d;当裂缝与单位岩心段部分相交时,Lreal(unit)≤Lfrac,如图5a、5b;而当裂缝完全在单位岩心段内时,Lreal(unit)=Lfrac,如图5c。当某段岩心的裂缝面密度值Ds(unit)较高时,可反映该段裂缝较为发育;而当某段岩心的裂缝面密度值Ds(unit)较低时,则可反映该段裂缝发育能力不足。
步骤五、裂缝面缝率表征:裂缝面缝率(φ)指测量面内裂缝所占面积与单位测量面积的比值,计算公式为:
Figure BDA0001289283220000102
公式中,S代表岩心测量面,Si代表测量面内每条裂缝所占面积。由于每个裂缝面积均可表示为裂缝长度(li)与裂缝宽度(wi)之积,因此上式可进一步转换为:
Figure BDA0001289283220000103
将岩心测量中各单位岩心段的纵截面作为测量面,则上述公式可表述为:
Figure BDA0001289283220000104
φ(unit)为某单位岩心段内的裂缝面缝率;(li·wi)代表测量面内每条裂缝所占面积。同裂缝面密度计算方法一致,上式可进一步演变为:
式中,Sreal(unit)代表在该单位岩心段内各裂缝截面积,其可详述为:
Figure BDA0001289283220000112
当某段岩心的裂缝面缝率值φ(unit)较高时,可反映该段裂缝能提供的孔隙度较高;而当某段岩心的裂缝面缝率值φ(unit)较低时,则指示该段裂缝所提供的岩石孔隙相对不足。
步骤六、构造/非构造裂缝比:用于反映某井裂缝的主要成因;在此采用构造裂缝面密度与非构造裂缝面密度的比值进行表征,以反映裂缝以构造成因为主还是以非构造成因为主。根据步骤四,构造裂缝面密度可表示为:
Figure BDA0001289283220000113
而非构造裂缝面密度则表示为:
Figure BDA0001289283220000114
两式相除,即可得:
Figure BDA0001289283220000115
当某段岩心的构造/非构造裂缝比Rstr/nonstr高于1时,反映该段岩心裂缝主要以构造成因为主;而当某段岩心的构造/非构造裂缝比Rstr/nonstr低于1时,反映该段岩心裂缝则主要以非构造成因为主。
为了使本领域的普通技术人员可以对本发明专利中的岩心裂缝连续表征方法有更深刻的理解,下面将以A井为例,参照流程图(图1),对本发明中的裂缝连续表征方法做进一步详细介绍。
步骤一、裂缝类型识别:A井为某研究区的页岩气钻井。以A井中3030-3036m深度段作为研究对象,在排除人工裂缝后,主要识别出剪性构造裂缝和层间缝两类裂缝,并以层间缝为主。
步骤二、裂缝参数统计:基于所识别的不同类型裂缝,分别统计岩心直径(d)、每条裂缝的中心深度(Fmid)、长度(Lfrac)、开度(w)、充填程度和充填物以及倾角(dip)等关键参数。通过裂缝类型识别与裂缝参数统计,该井取心层段裂缝主要参数统计如下(表1):
表1 A井主要裂缝参数汇总表
Figure BDA0001289283220000121
A井主要统计了上奥陶统观音桥组和五峰组两套地层的裂缝参数。其中,观音桥组岩性以深灰色、灰黑色生屑泥质灰岩为主,而五峰组岩性则以黑色、灰黑色碳质页岩为主。裂缝类型主要包括构造裂缝和层间缝两种类型,充填类型以未充填为主,充填物均为方解石,裂缝长度介于0.10-0.98m之间,开度为0.1-1.5mm,倾角存在低角度和高角度两个峰值,低角度裂缝主要为层间缝,高角度裂缝主要为张性构造裂缝。
步骤三、表征精度确定:由于统计的深度范围为3030-3036m,共计长度约6m;且鉴于本次表征主要用于反映裂缝在垂向分布上的差异性,因此将表征单位长度Lunit确定为1m,即依次表征3030-3031m、3031-3032m、3032-3033m、3033-3034m、3034-3035m和3035-3036m深度范围的裂缝参数。
步骤四、裂缝密度表征:首先,根据公式(1)-(2)将裂缝中心深度转换为裂缝顶深(Ftop)和底深(Fbtm);其次,根据公式(7),计算每条裂缝在某单位岩心段内的实际统计裂缝长度Lreal(unit),如表2;最后,依据公式(6),将各单位岩心段内的所有实际统计裂缝长度汇总求和,再除以单位测量面面积,即可获得A井垂向上的裂缝密度连续分布值(表2)。
表2 A井裂缝密度统计数据表
Figure BDA0001289283220000141
Figure BDA0001289283220000151
步骤五、裂缝面缝率表征:基于表1内的A井裂缝参数统计数据,首先,根据公式(1)-(2)将裂缝中心深度转换为裂缝顶深(Ftop)和底深(Fbtm);其次,根据公式(12),计算每条裂缝在单位岩心段内实际所占裂缝面积Sreal(unit),如表3;最后,依据公式(11),将各单位岩心段内的所有实际裂缝所占面积汇总求和,再除以单位测量面面积,即可获得A井裂缝面缝率的垂向连续分布值(表3)。
表3 A井裂缝面缝率统计数据表
Figure BDA0001289283220000152
Figure BDA0001289283220000161
Figure BDA0001289283220000171
步骤六、构造/非构造裂缝比:基于表1和表2内的A井裂缝基本参数数据和裂缝密度统计数据,利用公式(13)-(14)分类计算构造裂缝和非构造裂缝的裂缝面密度值(表4和表5);在此基础上,利用公式(15),计算各单位岩心段内构造裂缝密度与非构造裂缝密度的比值,即为构造/非构造裂缝比(表6)。由于本组测试样品中的非构造裂缝均为层间缝,因此该比值结果也为构造/层间裂缝比。岩心段上部构造/层间裂缝比较高,整体大于2.0;下部比值极低,均低于0.5;可反映测试岩心段上部裂缝以构造成因为主,下部裂缝则以非构造成因为主。
表4 A井构造裂缝密度统计数据表
Figure BDA0001289283220000172
表5 A井非构造裂缝密度统计数据表
Figure BDA0001289283220000181
Figure BDA0001289283220000191
表6 A井构造/非构造裂缝比统计数据表
Figure BDA0001289283220000192
本领域的技术人员应当理解,岩心裂缝各参数的准确测量是裂缝连续准确表征的重要前提,而实际取心段岩心的完整性、连续性均可影响后期裂缝密度计算的准确性。因此,为了保证该方法的计算结果可与测井、地震手段的分析结果进行有效对比,在岩心观测之前有必要整理、摆齐岩心,裂缝表征结果才具有较高的精度。

Claims (1)

1.一种基于岩心的裂缝连续表征方法,其特征在于,包括以下步骤:
步骤一、裂缝类型识别:裂缝的有效识别是裂缝准确统计的基本前提,主要是在岩心观测中对不同成因和性质的裂缝进行区分;裂缝整体划分为人工裂缝和天然裂缝两类,在裂缝观测时要排除人工裂缝对裂缝表征的干扰;天然裂缝又分为构造裂缝和非构造裂缝两类:构造裂缝主要包括张性裂缝、剪性裂缝和压性裂缝,而非构造裂缝则主要包括溶蚀缝、压实缝、风化缝、层间缝和沉积裂缝;
裂缝类型按以下特点进行识别:
人工裂缝指由于人工施钻所形成的相关裂缝,其形态取决于压裂井附近的现代地应力状态及钻孔方向;在未压裂的井段中,人工裂缝仅见于岩心表面,很难贯穿至岩心内部或切穿岩心,且形态与钻头形状具有明显相似性,因此其能够与天然裂缝有效区分;
天然裂缝中,构造裂缝与区域构造作用密切相关,断裂或褶皱附近的构造裂缝较为发育,成组出现、走向稳定、平直、穿一层或多层,且剪性裂缝的延伸性和平整度优于张性裂缝和压性裂缝;非构造裂缝成因复杂,形态多样,尺度跨度和走向变化大,与构造作用无关:溶蚀缝指成岩后地层中的部分物质在酸性流体作用下溶蚀形成的裂缝,在裂缝面上可见明显矿物溶蚀特征;而由于地层存在一些封闭的含液成岩缝和气孔,在垂向的压实作用下,液体压力上升会导致在侧向压力较小的方向扩展形成裂缝,这类裂缝则称之为压实缝,呈头盖骨“缝合线”状;风化缝指由于表层风化、剥蚀和盖层沉积,在风化壳和不整合面上所形成的裂缝,横截面多呈“V”字型,切割岩层较浅,通常不切穿岩层,走向不稳定;层间缝为在相邻地层间,由于层面间不整合或沉积条件不同所形成的层间裂缝群,倾角变化大,与地层面平行,走向离散,不切穿地层;而由于流体的定向流动和冲击作用,砾石、岩屑会呈定向排列,形成液体流动通道,进而形成沉积裂缝,并与沉积构造特征密切相关;
步骤二、裂缝参数统计:基于不同类型裂缝的识别,统计各条裂缝的相关参数,相关参数主要包括裂缝性质、裂缝中心深度Fmid,单位m、裂缝长度Lfrac,单位m、裂缝开度w,单位mm、充填程度及充填物、裂缝倾角dip和岩心直径d,单位m;其中,裂缝中心深度Fmid测量不便时,首先测量裂缝的顶深Ftop,单位m和底深Fbtm,单位m,然后利用公式:
Figure FDA0002214023630000021
间接计算得到裂缝中心深度Fmid;裂缝长度Lfrac指裂缝在岩心上截取的直线长度,表示为裂缝椭圆截面的长轴长度;裂缝开度w指裂缝的张开程度,亦即裂缝两侧的间隔距离;充填程度反映裂缝内的充填情况,主要包括全充填、半充填和未充填三类;裂缝倾角dip指裂缝面与岩心横截面间的夹角;裂缝长度Lfrac、裂缝顶深Ftop和底深Fbtm、裂缝倾角dip之间满足关系式:
Fbtm-Ftop=Lfrac·sin(dip) (2)
步骤三、表征精度确定:表征精度指进行裂缝计量的最小岩心单元,用每次裂缝表征的单位岩心长度Lunit表示;表征精度为均匀划分或根据实际计量需求按地层、岩性、沉积相非均匀划分;单位岩心顶深和底深分别为Utop,单位m和Ubtm,单位m,且满足以下关系式:
Lunit=Ubtm-Utop (3)
步骤四、裂缝密度表征:裂缝密度反映裂缝的集中程度,有线密度、面密度和体密度之分;线密度指一条单位长度的测线上通过的裂缝总条数,面密度指单位面积内所有裂缝的总长度,而体密度指在单位体积内裂缝穿切的总面积;以下按面密度进行统计,即:
Figure FDA0002214023630000031
式中,Ds代表裂缝面密度,S代表测量面面积,li代表该测量面内各条裂缝的长度;将各单位岩心段的纵截面作为测量面,则上式可变换为:
Figure FDA0002214023630000032
针对每条裂缝,由于裂缝的对应深度与单位岩心段深度并非完全一致,二者存在多种交互关系,而实际仅需统计在该单位岩心段内的裂缝长度,因此上式进一步转换为:
Figure FDA0002214023630000033
Ds(unit)为某一单位岩心内的裂缝面密度;Lreal(unit)代表各条裂缝在该单位岩心上截取的实际长度,其进一步表征为:
Figure FDA0002214023630000034
式中,min(Ubtm,Fbtm)代表某单位岩心底深与裂缝底深二者中的较小值,max(Utop,Ftop)代表某单位岩心顶深与裂缝顶深二者中的较大值;max[0,t]代表对t取正值,而利用max[0,min(Ubtm,Fbtm)-max(Utop,Ftop)]则能够计算各条裂缝在该单位岩心内的裂缝垂向长度,再除以sin(dip)即得该单位岩心内的真实裂缝长度Lreal(unit);根据公式易知,当裂缝与单位岩心段不相交时,Lreal(unit)=0;当裂缝与单位岩心段部分相交时,Lreal(unit)≤Lfrac;而当裂缝完全在单位岩心段内时,Lreal(unit)=Lfrac;当某段岩心的裂缝面密度值Ds(unit)较高时,反映该段裂缝较为发育;而当某段岩心的裂缝面密度值Ds(unit)较低时,则反映该段裂缝发育能力不足;
步骤五、裂缝面缝率表征:裂缝面缝率(φ)指测量面内裂缝所占面积与单位测量面积的比值,计算公式为:
Figure FDA0002214023630000041
公式中,S代表岩心测量面,Si代表测量面内每条裂缝所占面积;由于每个裂缝面积均表示为裂缝长度li与裂缝宽度wi之积,因此上式进一步转换为:
Figure FDA0002214023630000042
将岩心测量中各单位岩心段的纵截面作为测量面,则上述公式表述为:
Figure FDA0002214023630000043
φ(unit)为某单位岩心段内的裂缝面缝率;(li·wi)代表测量面内每条裂缝所占面积;同裂缝面密度计算方法一致,上式进一步演变为:
Figure FDA0002214023630000051
式中,Sreal(unit)代表在该单位岩心段内各裂缝截面积,其详述为:
Figure FDA0002214023630000052
当某段岩心的裂缝面缝率值φ(unit)较高时,反映该段裂缝能提供的孔隙度较高;而当某段岩心的裂缝面缝率值φ(unit)较低时,则指示该段裂缝所提供的岩石孔隙相对不足;
步骤六、构造/非构造裂缝比:用于反映某井裂缝的主要成因;在此采用构造裂缝面密度与非构造裂缝面密度的比值进行表征,以反映裂缝以构造成因为主还是以非构造成因为主;根据步骤四,构造裂缝面密度表示为:
Figure FDA0002214023630000053
而非构造裂缝面密度则表示为:
两式相除,即得:
当某段岩心的构造/非构造裂缝比Rstr/nonstr高于1时,反映该段岩心裂缝主要以构造成因为主;而当某段岩心的构造/非构造裂缝比Rstr/nonstr低于1时,反映该段岩心裂缝则主要以非构造成因为主。
CN201710319039.3A 2017-05-08 2017-05-08 一种基于岩心的裂缝连续表征方法 Active CN107167575B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710319039.3A CN107167575B (zh) 2017-05-08 2017-05-08 一种基于岩心的裂缝连续表征方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710319039.3A CN107167575B (zh) 2017-05-08 2017-05-08 一种基于岩心的裂缝连续表征方法

Publications (2)

Publication Number Publication Date
CN107167575A CN107167575A (zh) 2017-09-15
CN107167575B true CN107167575B (zh) 2020-01-14

Family

ID=59813641

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710319039.3A Active CN107167575B (zh) 2017-05-08 2017-05-08 一种基于岩心的裂缝连续表征方法

Country Status (1)

Country Link
CN (1) CN107167575B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10947841B2 (en) * 2018-01-30 2021-03-16 Baker Hughes, A Ge Company, Llc Method to compute density of fractures from image logs
CN109238799B (zh) * 2018-09-11 2020-10-30 东北石油大学 裂缝结构精细表征的碳酸盐岩板状岩心制备方法
CN109341528B (zh) * 2018-10-22 2020-06-09 大昌建设集团有限公司 一种基于三维激光扫描的建筑物爆破振动损伤评价方法
CN111272630B (zh) * 2020-02-28 2022-05-10 西南石油大学 致密岩心人工裂缝参数的计算方法
CN111507954A (zh) * 2020-04-14 2020-08-07 中国石油天然气股份有限公司大港油田分公司 用于储层裂缝的统计方法及机器可读存储介质
CN114325854A (zh) * 2022-03-04 2022-04-12 中海油研究总院有限责任公司 一种岩心不同级次裂缝发育程度判断方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103064114A (zh) * 2011-10-18 2013-04-24 中国石油化工股份有限公司 一种裂缝储层的表征方法及装置
CN104989360A (zh) * 2015-07-10 2015-10-21 中国石油天然气股份有限公司 一种低渗透油藏动态裂缝的表征方法
CN105842754A (zh) * 2016-03-24 2016-08-10 中国石油天然气股份有限公司 一种定量确定超深层致密砂岩裂缝有效性的方法
CN106370817A (zh) * 2016-10-12 2017-02-01 西南石油大学 一种基于岩心分析和电成像测井的缝洞定量表征方法
CN107064454A (zh) * 2017-02-16 2017-08-18 中国海洋石油总公司 一种储层裂缝发育程度的岩心确定方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103064114A (zh) * 2011-10-18 2013-04-24 中国石油化工股份有限公司 一种裂缝储层的表征方法及装置
CN104989360A (zh) * 2015-07-10 2015-10-21 中国石油天然气股份有限公司 一种低渗透油藏动态裂缝的表征方法
CN105842754A (zh) * 2016-03-24 2016-08-10 中国石油天然气股份有限公司 一种定量确定超深层致密砂岩裂缝有效性的方法
CN106370817A (zh) * 2016-10-12 2017-02-01 西南石油大学 一种基于岩心分析和电成像测井的缝洞定量表征方法
CN107064454A (zh) * 2017-02-16 2017-08-18 中国海洋石油总公司 一种储层裂缝发育程度的岩心确定方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
巢北地区二叠系页岩裂缝发育特征及主控因素分析;袁雷雷;《中国优秀硕士学位论文全文数据库 基础科学辑》;20150215(第2期);A011-109 *

Also Published As

Publication number Publication date
CN107167575A (zh) 2017-09-15

Similar Documents

Publication Publication Date Title
CN107167575B (zh) 一种基于岩心的裂缝连续表征方法
CN110644980B (zh) 一种超低渗透油藏储层综合分类评价方法
CN105931125B (zh) 一种致密油分段多簇体积压裂水平井产量预测方法
CN107346455A (zh) 一种用于识别页岩气产能的方法
CN112835124B (zh) 基于成像测井和阵列声波测井资料的裂缝有效性评价方法
CN105317435A (zh) 一种水平井裂缝识别方法
CN110320569B (zh) 一种致密砂岩储层单井裂缝发育强度定量评价方法
CN112145165B (zh) 一种微裂缝-孔隙型储层动静态渗透率转换方法
CN100455768C (zh) 确定扩径层段地层密度的方法
CN108593531B (zh) 一种岩体风化程度快速量化评价方法及其应用
CN111206921A (zh) 一种适用于火山岩溢流相有利储层的描述方法
CN112560246B (zh) 一种目标井散点地层压力系数的预测方法
CN117192610A (zh) 一种基于纵横波坐标变换的含油气储层预测方法
Waring et al. Devonian dolomitized reef, D3 reservoir, Leduc field, Alberta, Canada
CN110619490B (zh) 一种致密碳酸盐岩断裂破碎带胶结封闭性的评价方法
CN112147713B (zh) 泥页岩总有机碳含量分段预测方法
CN110764161A (zh) 一种通过钻探资料判识碳酸盐岩断裂破碎带的综合方法
Liu et al. Beyond volumetrics: Petrophysical characterization using rock types to predict dynamic flow behavior in tight gas sands
CN110929364A (zh) 一种致密碎屑岩储层微裂缝定量评价方法
CN113720745A (zh) 含碳屑碎屑岩储层地球物理测井计算孔隙度的方法
CN107622450B (zh) 基于测井资料快速判别水平井具备经济产量的方法
CN112859196B (zh) 井筒内断层断点的精确识别方法
Petrokushyn et al. Reservoirs Type Classification Using Hydraulic Flow Units Approach With An Application To Pivnichno-Korobochkynska Area (Dnipro-Donetsk Depression)
Jia et al. Study On Thin Differential Interpretation Method Of 0.2 m High Resolution Logging Series in Sabei Development Zone
CN106526696A (zh) 煤层顶底板岩性预测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant