CN107288626A - 一种页岩气地质甜度与工程甜度计算方法 - Google Patents
一种页岩气地质甜度与工程甜度计算方法 Download PDFInfo
- Publication number
- CN107288626A CN107288626A CN201610202327.6A CN201610202327A CN107288626A CN 107288626 A CN107288626 A CN 107288626A CN 201610202327 A CN201610202327 A CN 201610202327A CN 107288626 A CN107288626 A CN 107288626A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- sugariness
- engineering
- shale
- 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.)
- Granted
Links
- 238000000205 computational method Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000003825 pressing Methods 0.000 claims abstract description 11
- 238000010276 construction Methods 0.000 claims abstract description 10
- 238000009826 distribution Methods 0.000 claims abstract description 6
- 235000021185 dessert Nutrition 0.000 claims description 63
- 239000012530 fluid Substances 0.000 claims description 13
- 238000004519 manufacturing process Methods 0.000 claims description 12
- 230000001186 cumulative effect Effects 0.000 claims description 10
- 230000000052 comparative effect Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 239000004576 sand Substances 0.000 claims description 9
- 238000013459 approach Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 238000011161 development Methods 0.000 claims description 5
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 5
- 239000011707 mineral Substances 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 239000003795 chemical substances by application Substances 0.000 claims description 4
- 239000007788 liquid Substances 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 235000009508 confectionery Nutrition 0.000 claims description 3
- 239000005416 organic matter Substances 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 2
- 230000035699 permeability Effects 0.000 claims description 2
- 235000008331 Pinus X rigitaeda Nutrition 0.000 claims 1
- 235000011613 Pinus brutia Nutrition 0.000 claims 1
- 241000018646 Pinus brutia Species 0.000 claims 1
- 241001122767 Theaceae Species 0.000 claims 1
- 206010016256 fatigue Diseases 0.000 abstract 1
- 208000010392 Bone Fractures Diseases 0.000 description 13
- 206010017076 Fracture Diseases 0.000 description 13
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 5
- 229910052799 carbon Inorganic materials 0.000 description 5
- 239000011435 rock Substances 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- CZMRCDWAGMRECN-UGDNZRGBSA-N Sucrose Chemical compound O[C@H]1[C@H](O)[C@@H](CO)O[C@@]1(CO)O[C@@H]1[C@H](O)[C@@H](O)[C@H](O)[C@@H](CO)O1 CZMRCDWAGMRECN-UGDNZRGBSA-N 0.000 description 2
- 229930006000 Sucrose Natural products 0.000 description 2
- 244000269722 Thea sinensis Species 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000005720 sucrose Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 208000006670 Multiple fractures Diseases 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000000418 atomic force spectrum Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
- E21B43/267—Methods for stimulating production by forming crevices or fractures reinforcing fractures by propping
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种页岩气地质甜度与工程甜度的计算方法。该方法选取具有独立性的评价参数,以压后累产量为依据,采用灰色关联度方法计算各参数的权重分配,应用欧式贴进度的计算方法表征页岩“甜度”。本发明成功解决了页岩油气藏中的井进行压裂施工段簇划分针对性不强、有效出气段簇比例不高的难题,为水平井分段压裂段、簇位置优选提供了有力依据,从而最大限度地实现“降本增效”的目标。
Description
技术领域
本发明涉及石油开采技术领域,尤其涉及一种页岩气地质甜度与工程甜度计算方法。
背景技术
目前,关于页岩地质甜点与工程甜点的评价方法有很多种。地质甜点的评价参数主要包括有机碳含量(TOC)、热成熟度(RO)、含气量、孔隙度、天然裂缝特性参数及孔隙压力等;工程甜点的评价参数主要包括岩石矿物组分及岩石力学参数等。但是目前已知的评价方法主要存在以下问题:
1)参与地质甜点及工程甜点计算的评价参数的独立性有待商榷。这是因为把有关联性的参数都放进模型中往往会使计算结果出现偏差。例如,有机碳含量TOC与含气量有一定的相关性,有机碳含量高只能表明页岩储层具有高含气量的物质基础,但不代表含气量一定高,其还取决于热演化程度RO的高低。因此通常选择含气量作为表征页岩含气性的评价参数。
2)甜点的程度没有进一步界定,甜点值的大小与压后产量的相关性不强。为了进一步表征页岩甜点的好坏程度,引入“甜度”的概念来表征甜点的程度大小,并将甜度与压后产量相关联(注:甜度是一个相对值,例如通常以蔗糖作为基准物,一般以10%或15%的蔗糖水溶液在20℃时的甜度为1.0,其他糖的甜度则与之相比较得到)。
3)地质甜点与工程甜点的权重计算方法有待改进。以往大多采用对等的权重分配来计算总的甜点指标,而没有从压后产量的关联性来寻求地质甜点与工程甜点的权重分配。
综上所述,有必要改进现有地质甜点及工程甜点的评价指标及计算方法,以增加水平井分段压裂段、簇位置选择的科学性和针对性,更好地实现页岩气开发的“降本增效”目标。
发明内容
为了解决上述技术问题,本发明提供了一种新的适用于页岩油气藏的地质甜度与工程甜度计算方法,以增加水平井分段压裂段、簇位置选择的科学性和针对性,更好地实现页岩气开发的“降本增效”目标。
一种页岩气地质甜度与工程甜度计算方法,包括以下步骤:
S110,分析页岩地层的各项地质参数之间的相关性,从中选出彼此独立的地质参数作为评价参数;
S120,将各项评价参数的最佳值作为理想地质条件,计算目标井层实际地质参数与理想地质条件的欧式贴近度作为地质甜度指标;
S130,计算近井工程甜点和远井工程甜点,从中确定最大的近井工程甜点和最大的远井工程甜点;
S140,将最大的近井工程甜点和最大的远井工程甜点作为理想井层工程甜点,计算目标井层工程甜点与理想井层工程甜点的欧式贴近度作为工程甜度指标;
S150,以压后累积产量为依据,对地质甜度指标和工程甜度指标分配权重系数,进而获得综合甜度指标;
S160,根据现场页岩气勘探结果修正步骤S150中的权重系数,进而修正综合甜度指标。
根据本发明的实施例,上述步骤S110中,评价参数主要包括总孔隙度、有机质孔隙度、基质渗透率、页岩厚度、天然裂缝发育程度、总含气量、游离气比例、压力系数、热演化程度、脆性矿物含量、杨氏模量、泊松比。
根据本发明的实施例,上述欧式贴近度的计算方法如下:
设A为由n-1个待选的页岩区块A1、A2、A3、……、An-1及理想的页岩区块组成的集合,P是对应于待选的页岩区块A1,A2,A3……An-1及理想的页岩区块的m个特征参数P1、P2、……、Pm组成的集合,由集合A到集合P的一个模糊关系记为R:
其中,rij表示待选的页岩区块或理想的页岩区块Ai具有参数Pj特征的隶属度;
按最大最小法求取集合A到集合P之间的模糊关系R:
式中,x为待选的页岩区块或理想的页岩区块的任一特征参数,
a1为待选的页岩区块或理想的页岩区块的任一特征参数的最小值,
a2为待选的页岩区块或理想的页岩区块的任一特征参数的最大值;
将模糊矩阵R划分为n个次级模糊矩阵R1、R2、……、Rn-1及通过下式计算Rj(j=1、2、……、n-1)与的接近程度:
根据本发明的实施例,上述步骤S130中,基于压裂施工破裂压力曲线,可以通过以下算式计算近井工程甜点:
式中,SE表征近井工程甜点,无量纲;ST0TCBD和ST0TCCD分别为压裂施工破裂压力曲线所包络的脆性覆盖区域的面积和塑性覆盖区域的面积,单位为MPa·min。
根据本发明的实施例,上述步骤S130中,可以通过计算总的加砂量与总的入井压裂液量的比值来表征远井工程甜点。
根据本发明的实施例,上述总的加砂量为不同种类的支撑剂按照平均粒径比值折算成一种主体支撑剂后的总量,和/或所述总的入井压裂液量为不同种类的压裂液按照黏度比折算成一种主体压裂液后的总量。
根据本发明的实施例,上述主体支撑剂粒径可以是40-70目,和/或所述主体压裂液为滑溜水。
根据本发明的实施例,上述步骤S150中,优选以压后累积产量为依据,计算其与地质甜度指标和工程甜度指标的灰色关联度,再对灰色关联度进行归一化处理,获得地质甜度指标和工程甜度指标的权重系数。
根据本发明的实施例,计算压后累积产量与地质甜度指标和工程甜度指标的灰色关联度的方法包括以下步骤:
1)设置j个标准序列{Xo(i)}={Xo(1),Xo(2),Xo(3),…,Xo(n)},比较序列{Xt(i)}={Xt(1),Xt(2),Xt(3),…,Xt(n)},(i=1,2,3,…,n),其中,标准序列为地质甜度指标和工程甜度指标,比较序列为压后若干年的累积产量;
2)采用均值化法对各评价参数进行无量纲化处理:
3)通过下式计算无量纲化处理后的序列{X’o(i)}与{X’t(i)}序列的关联系数Yot(j):
式中,ρ为分辨系数,其取值范围为[0,1];
4)综合关联系数,获得各比较序列与标准序列之间的灰色关联度:
根据本发明的实施例,上述步骤S160中,现场页岩气勘探结果可以是现场页岩气水平井分段压裂后的产气剖面测试结果。
与现有技术相比,本发明的一个或多个实施例可以具有如下优点:
1)该方法选取具有独立性的评价参数,避免因重复考虑而导致计算结果出现较大的偏差。
2)该方法选择评价参数的最佳值组合成理想页岩模型,通过计算页岩实际地质参数与理想地质模型之间的欧式贴进度来表征页岩地质甜度,更加科学合理。
3)该方法综合考虑近井工程甜点和远井工程甜点,获得综合的工程甜度更加科学和具有针对性。
4)该方法以压后累产量为依据,采用灰色关联度方法对地质甜度和工程甜度进行权重分配,从而规避了常规段簇划分中不考虑页岩气井产量的技术风险,为最大限度地实现“降本增效”的目标提供了技术支持。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例共同用于解释本发明,并不构成对本发明的限制。在附图中:
图1是本发明提供的计算方法的流程图;
图2是本发明一示例中计算近井工程甜点的方法的示意图;
图3是本发明一示例中通过本发明提供的计算方法获得的地质甜度,并据此进行段簇划分的结果的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,以下结合附图对本发明作进一步地详细说明。
图1是本发明提供的计算方法的流程图。从图中可以看出,该方法主要包括六大步骤:
S110,确定页岩关键的地质参数。
在该步骤中,主要分析页岩地层的各项地质参数之间的相关性,从中选出彼此独立的地质参数作为评价参数。
例如,与地质甜度关联的独立参数主要包括总孔隙度、有机质孔隙度、基质渗透率、页岩厚度、天然裂缝发育程度(0不发育、0.5一般发育、1极发育)、总含气量、游离气比例、压力系数、热演化程度RO、石英及碳酸盐岩等脆性矿物含量、杨氏模量、泊松比。评价方法包括岩心测试、测井及录井等方法,这些都为常规方法,在此不赘。如背景技术中提及的那样,由于有机碳含量TOC与总含气量不是相互独立,因此在此略去有机碳含量这一参数。
S120,计算页岩地质甜度指标。
将上述12项评价参数的最佳值(注,最佳值不一定是最大值,也可能是最小值或中间值)作为理想地质条件,计算目标井层实际地质参数与理想地质条件的欧式贴近度作为地质甜度指标。
欧式贴近度的计算原理如下:
设A为由n-1个待选的页岩区块A1、A2、A3、……、An-1及理想的页岩区块组成的集合,P是对应于待选的页岩区块A1,A2,A3……An-1及理想的页岩区块的m个特征参数P1、P2、……、Pm组成的集合,由集合A到集合P的一个模糊关系记为R:
其中,rij表示待选的页岩区块或理想的页岩区块Ai具有参数Pj特征的隶属度;
按最大最小法求取集合A到集合P之间的模糊关系R:
式中,x为待选的页岩区块或理想的页岩区块的任一特征参数,
a1为待选的页岩区块或理想的页岩区块的任一特征参数的最小值,
a2为待选的页岩区块或理想的页岩区块的任一特征参数的最大值;
将模糊矩阵R划分为n个次级模糊矩阵R1、R2、……、Rn-1及通过下式计算Rj(j=1、2、……、n-1)与的接近程度:
S130,计算近井工程甜点和远井工程甜点,从中确定最大的近井工程甜点和最大的远井工程甜点。
常规的工程甜点计算方法一般采用脆性矿物含量或岩石力学参数的方法,而本方法的关键的地质参数中已包含了这些内容。在此,本发明的工程参数的计算主要是指近井工程甜点的计算和远井工程甜点的计算。
I、计算近井工程甜点
计算方法主要基于压裂施工破裂压力曲线的形态进行。破裂压力曲线的特征可以精细刻画岩石脆性在宏观上的表现。只要压裂液性质及排量稳定,井口压力在破裂压力峰值后会出现不同程度的下降。压力下降的越快,说明页岩的脆性特征越强,下降得越慢说明页岩的塑性特征越强。因此,基于压裂施工破裂压力曲线能够得出近井工程甜点。
在本实施例中,可以以破裂压力峰值至压力降落平稳时的时间区间为界,破裂压力曲线所包络的面积(破裂压力曲线对时间进行积分)与排量(选择恒定排量时发生的破裂)的乘积,即为页岩破裂期间所消耗的功。显然地,脆性地层破裂后,压力快速下降,消耗的功小,而塑性地层压力下降慢或几乎不降,消耗的功就较多。为简化起见,假设塑性页岩破裂后压力维持不变,则近井脆性指数(也即近井工程甜点)的计算公式为:
式中,SE为工程甜点,无量纲;ST0TCBD和ST0TCCD分别为图2所示的压裂施工破裂压力曲线多边形T0TCBDT0所包络的脆性覆盖区域的面积和多边形T0TCCDT0所包络的塑性覆盖区域的面积,单位为MPa·min。
如果有多次破裂压力显示,说明页岩的脆性更好(工程甜点更好)。按照上述同样的方法计算,只是要把每个排量对应的脆性面积相加,以及每个排量对应的塑性面积相加,然后将结果代入上式进行计算即可。
II、计算远井工程甜点
在该步骤中,通过计算总的加砂量与总的入井压裂液量的比值来表征远井工程甜点。在实际施工中,往往存在多种不同类型的支撑剂和/或多种不同类型的压裂液,因此需要对它们进行折算,以方便计算。例如,将多种支撑剂按照平均粒径比值折算为例如40-70目的主体支撑剂,从而获得总的加砂量,和/或将多种压裂液按照黏度比折算成一种主体压裂液的量,例如滑溜水的体积,从而获得总的入井压裂液量。然后求取折算的总加砂量与折算的总入井压裂液量的比值,也即综合砂液比。显然地,综合砂液比反应了远井地带的页岩可压性大小(工程甜点),其值越大,裂缝导流能力维持时间越长,压后的累产量应当越高。因此,本发明选用它来表征远井工程甜点。
S140,计算页岩工程甜度指标。
根据I和II的计算结果,从中选取最大的近井甜点及最大的远井甜点作为理想井层工程甜点,通过计算目标井层工程甜点与理想井层工程甜点的欧式贴近度来作为页岩的工程甜度指标。
S150,以压后累积产量为依据,对地质甜度指标和工程甜度指标分配权重系数,进而获得综合甜度指标。
具体地,以压后累积产量为依据,计算其与地质甜度指标和工程甜度指标的灰色关联度,再对灰色关联度进行归一化处理,获得地质甜度指标和工程甜度指标的权重系数。具体计算过程如下:
1)设置j个标准序列{Xo(i)}={Xo(1),Xo(2),Xo(3),…,Xo(n)},比较序列{Xt(i)}={Xt(1),Xt(2),Xt(3),…,Xt(n)},(i=1,2,3,…,n),其中,标准序列为地质甜度指标和工程甜度指标,比较序列为压后若干年的累积产量(在本实施例中,以压后三年累产量为依据)。
2)采用均值化法对各评价参数进行无量纲化处理:
3)通过下式计算无量纲化处理后的序列{X’o(i)}与{X’t(i)}序列的关联系数Yot(j):
式中,ρ为分辨系数,其取值范围为[0,1]。
4)综合关联系数,获得各比较序列与标准序列之间的灰色关联度:
5)其中,参数的权重系数可通过该参数的灰色关联度值与各参数灰色关联度的集合之间的比值求得。换言之,通过对灰色关联度进行归一化处理获得地质甜度指标和工程甜度指标的权重系数。
最后进行计算,最终可获得综合地质甜度指标与工程甜度指标的总的甜度指标,也称为综合甜度指标。
S160,检验和修正综合甜度指标。
在该步骤中,根据现场页岩气勘探结果修正步骤S150中的权重系数,进而修正综合甜度指标。
在本实施例中,可以结合现场页岩气水平井分段压裂后的产气剖面测试结果进行验证,由此建立目标区块页岩气的甜度指标规范,以指导后续的水平井分段压裂段、簇位置优选,从而最大限度地实现“降本增效”的目标。
图3所示的实施例是本发明的计算方法在川东南地区某井体积压裂改造施工中得到应用的结果。该井垂深2402.6m,测深4100m,水平段长1530m。通过本发明所提供的方法,计算了沿水平井筒连续分布的地质“甜度”,然后以此为基础进行段簇划分。由最终的实施结果可以得知,该井5个月累产气量达1651×104m3,与同平台另外一口井相比,该井5个月累产气量提高14.2%,取得了显著的经济效益。
以上所述,仅为本发明的具体实施案例,本发明的保护范围并不局限于此,任何熟悉本技术的技术人员在本发明所述的技术规范内,对本发明的修改或替换,都应在本发明的保护范围之内。
Claims (10)
1.一种页岩气地质甜度与工程甜度计算方法,包括以下步骤:
S110,分析页岩地层的各项地质参数之间的相关性,从中选出彼此独立的地质参数作为评价参数;
S120,将各项评价参数的最佳值作为理想地质条件,计算目标井层实际地质参数与理想地质条件的欧式贴近度作为地质甜度指标;
S130,计算近井工程甜点和远井工程甜点,从中确定最大的近井工程甜点和最大的远井工程甜点;
S140,将最大的近井工程甜点和最大的远井工程甜点作为理想井层工程甜点,计算目标井层工程甜点与理想井层工程甜点的欧式贴近度作为工程甜度指标;
S150,以压后累积产量为依据,对地质甜度指标和工程甜度指标分配权重系数,进而获得综合甜度指标;
S160,根据现场页岩气勘探结果修正步骤S150中的权重系数,进而修正综合甜度指标。
2.如权利要求1所述的页岩气地质甜度与工程甜度计算方法,其特征在于:
所述步骤S110中,评价参数包括总孔隙度、有机质孔隙度、基质渗透率、页岩厚度、天然裂缝发育程度、总含气量、游离气比例、压力系数、热演化程度、脆性矿物含量、杨氏模量、泊松比。
3.如权利要求1所述的页岩气地质甜度与工程甜度计算方法,其特征在于,所述欧式贴近度的计算方法如下:
设A为由n-1个待选的页岩区块A1、A2、A3、……、An-1及理想的页岩区块组成的集合,P是对应于待选的页岩区块A1,A2,A3……An-1及理想的页岩区块的m个特征参数P1、P2、……、Pm组成的集合,由集合A到集合P的一个模糊关系记为R:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>R</mi>
<mo>=</mo>
<msub>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>n</mi>
<mo>&times;</mo>
<mi>m</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>...</mo>
<mo>,</mo>
<mi>n</mi>
<mo>;</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>...</mo>
<mo>,</mo>
<mi>m</mi>
</mrow>
其中,rij表示待选的页岩区块或理想的页岩区块Ai具有参数Pj特征的隶属度;
按最大最小法求取集合A到集合P之间的模糊关系R:
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>&mu;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mo>&le;</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>/</mo>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo><</mo>
<mi>x</mi>
<mo><</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mo>&GreaterEqual;</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中,x为待选的页岩区块或理想的页岩区块的任一特征参数,
a1为待选的页岩区块或理想的页岩区块的任一特征参数的最小值,
a2为待选的页岩区块或理想的页岩区块的任一特征参数的最大值;
将模糊矩阵R划分为n个次级模糊矩阵R1、R2、……、Rn-1及通过下式计算Rj(j=1、2、……、n-1)与的接近程度:
<mrow>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
<msubsup>
<mi>R</mi>
<mi>n</mi>
<mo>*</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>m</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</msubsup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>R</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>R</mi>
<mi>n</mi>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>.</mo>
</mrow>
4.如权利要求1所述的页岩气地质甜度与工程甜度计算方法,其特征在于:
所述步骤S130中,基于压裂施工破裂压力曲线,通过以下算式计算近井工程甜点:
<mrow>
<msub>
<mi>S</mi>
<mi>E</mi>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>S</mi>
<mrow>
<mi>T</mi>
<mi>o</mi>
<mi>T</mi>
<mi>C</mi>
<mi>B</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>s</mi>
<mrow>
<mi>T</mi>
<mi>o</mi>
<mi>T</mi>
<mi>C</mi>
<mi>C</mi>
<mi>D</mi>
</mrow>
</msub>
</mfrac>
</mrow>
1
式中,SE表征近井工程甜点,无量纲;ST0TCBD和ST0TCCD分别为压裂施工破裂压力曲线所包络的脆性覆盖区域的面积和塑性覆盖区域的面积,单位为MPa·min。
5.如权利要求1所述的页岩气地质甜度与工程甜度计算方法,其特征在于:
所述步骤S130中,通过计算总的加砂量与总的入井压裂液量的比值来表征远井工程甜点。
6.如权利要求5所述的页岩气地质甜度与工程甜度计算方法,其特征在于:
所述总的加砂量为不同种类的支撑剂按照平均粒径比值折算成一种主体支撑剂后的总量,和/或所述总的入井压裂液量为不同种类的压裂液按照黏度比折算成一种主体压裂液后的总量。
7.如权利要求6所述的页岩气地质甜度与工程甜度计算方法,其特征在于:
所述主体支撑剂粒径为40-70目,和/或所述主体压裂液为滑溜水。
8.如权利要求1所述的页岩气地质甜度与工程甜度计算方法,其特征在于:
所述步骤S150中,以压后累积产量为依据,计算其与地质甜度指标和工程甜度指标的灰色关联度,再对灰色关联度进行归一化处理,获得地质甜度指标和工程甜度指标的权重系数。
9.如权利要求8所述的页岩气地质甜度与工程甜度计算方法,其特征在于,以压后累积产量为依据,计算其与地质甜度指标和工程甜度指标的灰色关联度的方法包括以下步骤:
1)设置j个标准序列{Xo(i)}={Xo(1),Xo(2),Xo(3),…,Xo(n)},比较序列{Xt(i)}={Xt(1),Xt(2),Xt(3),…,Xt(n)},(i=1,2,3,…,n),其中,标准序列为地质甜度指标和工程甜度指标,比较序列为压后若干年的累积产量;
2)采用均值化法对各评价参数进行无量纲化处理:
<mrow>
<msub>
<msup>
<mi>X</mi>
<mo>&prime;</mo>
</msup>
<mi>o</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>X</mi>
<mi>o</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mover>
<mi>X</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<mover>
<mi>X</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<msub>
<mi>X</mi>
<mi>o</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
3)通过下式计算无量纲化处理后的序列{X’o(i)}与{X’t(i)}序列的关联系数Yot(j):
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mi>o</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<munder>
<mi>min</mi>
<mi>o</mi>
</munder>
<munder>
<mi>min</mi>
<mi>t</mi>
</munder>
<mo>|</mo>
<msub>
<msup>
<mi>X</mi>
<mo>&prime;</mo>
</msup>
<mi>o</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<msup>
<mi>X</mi>
<mo>&prime;</mo>
</msup>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>+</mo>
<mi>&rho;</mi>
<munder>
<mi>max</mi>
<mi>o</mi>
</munder>
<munder>
<mi>max</mi>
<mi>t</mi>
</munder>
<mo>|</mo>
<msub>
<msup>
<mi>X</mi>
<mo>&prime;</mo>
</msup>
<mi>o</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<msup>
<mi>X</mi>
<mo>&prime;</mo>
</msup>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<msup>
<mi>X</mi>
<mo>&prime;</mo>
</msup>
<mi>o</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<msup>
<mi>X</mi>
<mo>&prime;</mo>
</msup>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>+</mo>
<mi>&rho;</mi>
<munder>
<mi>max</mi>
<mi>o</mi>
</munder>
<munder>
<mi>max</mi>
<mi>t</mi>
</munder>
<mo>|</mo>
<msub>
<msup>
<mi>X</mi>
<mo>&prime;</mo>
</msup>
<mi>o</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<msup>
<mi>X</mi>
<mo>&prime;</mo>
</msup>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
式中,ρ为分辨系数,其取值范围为[0,1];
4)综合关联系数,获得各比较序列与标准序列之间的灰色关联度:
<mrow>
<msub>
<mi>r</mi>
<mrow>
<mi>o</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mi>Y</mi>
<mi>o</mi>
<mi>t</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
10.根据权利要求1所述的页岩气地质甜度与工程甜度计算方法,其特征在于:
所述步骤S160中,现场页岩气勘探结果为现场页岩气水平井分段压裂后的产气剖面测试结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610202327.6A CN107288626B (zh) | 2016-04-01 | 2016-04-01 | 一种页岩气地质甜度与工程甜度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610202327.6A CN107288626B (zh) | 2016-04-01 | 2016-04-01 | 一种页岩气地质甜度与工程甜度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107288626A true CN107288626A (zh) | 2017-10-24 |
CN107288626B CN107288626B (zh) | 2020-08-21 |
Family
ID=60088249
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610202327.6A Active CN107288626B (zh) | 2016-04-01 | 2016-04-01 | 一种页岩气地质甜度与工程甜度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107288626B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107766662A (zh) * | 2017-10-26 | 2018-03-06 | 中国石油化工股份有限公司 | 一种页岩气水平井测试分段评价方法 |
CN109298449A (zh) * | 2018-09-26 | 2019-02-01 | 中国石油化工股份有限公司 | 一种确定陆相页岩油勘探有利单元的方法 |
CN109444334A (zh) * | 2018-10-23 | 2019-03-08 | 青海省第四地质矿产勘查院((青海省煤炭地质勘查院)) | 一种陆相页岩气评价方法 |
CN110288233A (zh) * | 2019-06-26 | 2019-09-27 | 重庆科技学院 | 一种基于模糊灰色关联法的深层页岩气可压性评价方法 |
CN110806600A (zh) * | 2018-08-06 | 2020-02-18 | 中国石油化工股份有限公司 | 一种定量评价页岩气甜点的方法 |
CN111027818A (zh) * | 2019-11-22 | 2020-04-17 | 中国石油天然气股份有限公司 | 一种页岩油分类评价方法 |
CN111460725A (zh) * | 2019-12-27 | 2020-07-28 | 中国地质调查局成都地质调查中心 | 基于多层次模糊识别的页岩气甜点预测 |
CN112418598A (zh) * | 2020-10-14 | 2021-02-26 | 中国石油化工股份有限公司 | 一种页岩气井压后各段产气贡献率评价方法 |
CN112459777A (zh) * | 2020-12-03 | 2021-03-09 | 西南石油大学 | 一种储层工程甜点系数的计算方法 |
CN112814643A (zh) * | 2021-01-07 | 2021-05-18 | 中国石油天然气股份有限公司 | 页岩油气藏甜点区的布井方法 |
CN113027444A (zh) * | 2021-04-25 | 2021-06-25 | 中国石油大学(北京) | 一种河湖相页岩油的甜点区预测方法及装置 |
CN115048744A (zh) * | 2022-07-01 | 2022-09-13 | 成都创源油气技术开发有限公司 | 一种基于地震数据识别岩石工程甜点方法与装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104502971A (zh) * | 2014-12-22 | 2015-04-08 | 中国石油天然气集团公司 | 页岩储层总有机碳含量和脆性空间展布的检测方法 |
CA2867580A1 (en) * | 2013-10-15 | 2015-04-15 | Korea Institute Of Geoscience And Mineral Resources (Kigam) | Fracking method at fracking intervals of a horizontal drilling zone in a sweet spot range determined based on measurements of resistivity and density in the horizontal drilling zone |
CN104750884A (zh) * | 2013-12-26 | 2015-07-01 | 中国石油化工股份有限公司 | 基于多因素非线性回归的页岩油气富集指数定量评价方法 |
CN105221140A (zh) * | 2014-06-20 | 2016-01-06 | 中国石油化工股份有限公司 | 一种确定页岩地层可压裂性指数的方法 |
CN105426650A (zh) * | 2015-06-08 | 2016-03-23 | 中国石油集团川庆钻探工程有限公司 | 利用测井资料优选页岩气井压裂改造井段的方法 |
-
2016
- 2016-04-01 CN CN201610202327.6A patent/CN107288626B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2867580A1 (en) * | 2013-10-15 | 2015-04-15 | Korea Institute Of Geoscience And Mineral Resources (Kigam) | Fracking method at fracking intervals of a horizontal drilling zone in a sweet spot range determined based on measurements of resistivity and density in the horizontal drilling zone |
CN104750884A (zh) * | 2013-12-26 | 2015-07-01 | 中国石油化工股份有限公司 | 基于多因素非线性回归的页岩油气富集指数定量评价方法 |
CN105221140A (zh) * | 2014-06-20 | 2016-01-06 | 中国石油化工股份有限公司 | 一种确定页岩地层可压裂性指数的方法 |
CN104502971A (zh) * | 2014-12-22 | 2015-04-08 | 中国石油天然气集团公司 | 页岩储层总有机碳含量和脆性空间展布的检测方法 |
CN105426650A (zh) * | 2015-06-08 | 2016-03-23 | 中国石油集团川庆钻探工程有限公司 | 利用测井资料优选页岩气井压裂改造井段的方法 |
Non-Patent Citations (1)
Title |
---|
卞晓冰 等: ""基于施工曲线的页岩气井压后评估新方法"", 《天然气工业》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107766662A (zh) * | 2017-10-26 | 2018-03-06 | 中国石油化工股份有限公司 | 一种页岩气水平井测试分段评价方法 |
CN110806600A (zh) * | 2018-08-06 | 2020-02-18 | 中国石油化工股份有限公司 | 一种定量评价页岩气甜点的方法 |
CN109298449A (zh) * | 2018-09-26 | 2019-02-01 | 中国石油化工股份有限公司 | 一种确定陆相页岩油勘探有利单元的方法 |
CN109444334A (zh) * | 2018-10-23 | 2019-03-08 | 青海省第四地质矿产勘查院((青海省煤炭地质勘查院)) | 一种陆相页岩气评价方法 |
CN110288233A (zh) * | 2019-06-26 | 2019-09-27 | 重庆科技学院 | 一种基于模糊灰色关联法的深层页岩气可压性评价方法 |
CN110288233B (zh) * | 2019-06-26 | 2023-08-08 | 重庆科技学院 | 一种基于模糊灰色关联法的深层页岩气可压性评价方法 |
CN111027818B (zh) * | 2019-11-22 | 2022-05-10 | 中国石油天然气股份有限公司 | 一种页岩油分类评价方法 |
CN111027818A (zh) * | 2019-11-22 | 2020-04-17 | 中国石油天然气股份有限公司 | 一种页岩油分类评价方法 |
CN111460725A (zh) * | 2019-12-27 | 2020-07-28 | 中国地质调查局成都地质调查中心 | 基于多层次模糊识别的页岩气甜点预测 |
CN111460725B (zh) * | 2019-12-27 | 2023-08-18 | 中国地质调查局成都地质调查中心 | 基于多层次模糊识别的页岩气甜点预测 |
CN112418598A (zh) * | 2020-10-14 | 2021-02-26 | 中国石油化工股份有限公司 | 一种页岩气井压后各段产气贡献率评价方法 |
CN112418598B (zh) * | 2020-10-14 | 2024-04-19 | 中国石油化工股份有限公司 | 一种页岩气井压后各段产气贡献率评价方法 |
CN112459777A (zh) * | 2020-12-03 | 2021-03-09 | 西南石油大学 | 一种储层工程甜点系数的计算方法 |
CN112459777B (zh) * | 2020-12-03 | 2022-04-15 | 西南石油大学 | 一种储层工程甜点系数的计算方法 |
CN112814643A (zh) * | 2021-01-07 | 2021-05-18 | 中国石油天然气股份有限公司 | 页岩油气藏甜点区的布井方法 |
CN113027444A (zh) * | 2021-04-25 | 2021-06-25 | 中国石油大学(北京) | 一种河湖相页岩油的甜点区预测方法及装置 |
CN115048744A (zh) * | 2022-07-01 | 2022-09-13 | 成都创源油气技术开发有限公司 | 一种基于地震数据识别岩石工程甜点方法与装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107288626B (zh) | 2020-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107288626A (zh) | 一种页岩气地质甜度与工程甜度计算方法 | |
CN109594968B (zh) | 一种页岩气多段压裂水平井压后裂缝参数评价方法及系统 | |
CN107044277B (zh) | 低渗透非均质油藏水平井重复压裂增产潜力评价方法 | |
Gherabati et al. | The impact of pressure and fluid property variation on well performance of liquid-rich Eagle Ford shale | |
CN104199121B (zh) | 一种页岩气藏建产有利区的综合判别方法 | |
CN105986817B (zh) | 一种用于识别页岩地层工程甜点的方法 | |
CN104389594A (zh) | 页岩气井产能评价预测方法 | |
CN104989357B (zh) | 一种选择致密油气藏重复压裂目标井的方法 | |
CN105221140A (zh) | 一种确定页岩地层可压裂性指数的方法 | |
CN105626025A (zh) | 一种页岩储层压裂的可压性评价方法 | |
CN109577929B (zh) | 一种超低渗致密油藏水平井建立有效驱替的定量评价方法 | |
CN106651158B (zh) | 一种超低渗致密油藏水平井注水开发见效程度定量评价方法 | |
CN103308433A (zh) | 一种基于孔隙演化分析评价致密砂岩储层成岩相的方法 | |
CN103645519A (zh) | 一种基于火山岩风化壳识别及划分标准的油气勘探方法 | |
CN112983368B (zh) | 一种注采优化协同化学驱实现原油均衡驱替的方法 | |
CN114186440B (zh) | 一种地质-工程“双轨制”页岩可压性综合评价方法 | |
CN106680172A (zh) | 评价致密油储层的裂缝的方法 | |
CN108661628B (zh) | 一种基于参数优选的工程甜点定量评价方法 | |
CN104483706A (zh) | 一种基于煤岩力学参数的煤体结构测井定量识别方法 | |
CN109522619A (zh) | 一种低渗油田周期采油方法 | |
CN105719065A (zh) | 一种复杂油藏储量品质分类综合评价方法 | |
CN106845685B (zh) | 一种降低水平井试油压裂作业成本的参数优化方法 | |
CN105625992A (zh) | 针对屋脊断块油藏的合理提液方法 | |
CN108661629B (zh) | 一种针对页岩地层的工程甜点定量评价方法 | |
CN109113690A (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 |