CN105718745A - 一种酸压裂缝初始导流能力计算方法 - Google Patents
一种酸压裂缝初始导流能力计算方法 Download PDFInfo
- Publication number
- CN105718745A CN105718745A CN201610051684.7A CN201610051684A CN105718745A CN 105718745 A CN105718745 A CN 105718745A CN 201610051684 A CN201610051684 A CN 201610051684A CN 105718745 A CN105718745 A CN 105718745A
- Authority
- CN
- China
- Prior art keywords
- overbar
- crack
- epsiv
- rock beam
- flow conductivity
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Medicinal Chemistry (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Food Science & Technology (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种酸压裂缝初始导流能力计算方法,依次包括以下步骤:(A)开展酸刻蚀物理模拟实验,利用三维激光扫描仪对酸刻蚀后的每个岩板粗糙表面进行扫描,获取岩板表面的三维数据;(B)利用步骤(A)中得到的三维数据计算每个岩板表面的横向曲折比和纵向曲折比两个参数;(C)计算裂缝(两个岩板组成)的平均横向曲折比和平均纵向曲折比(D)计算裂缝初始导流能力(kfw)0。本发明操作简便,原理可靠,安全性好,能够为酸压方案设计提供有效指导。
Description
技术领域
本发明涉及石油领域,尤其是酸化压裂领域中的一种酸压裂缝初始导流能力计算方法
背景技术
酸化压裂(简称酸压)是油气增产改造的主要措施,其主要是在高于油气储集层破裂压力或天然裂缝闭合压力条件下,将酸液(或前置液)注入储层,在储集层中形成新的人工缝或使天然裂缝张开,同时注入的酸液对裂缝表面造成非均匀刻蚀。酸液返排后,由于裂缝表面的粗糙不平,表面的凸起部分使裂缝在闭合压力作用下保持一定的开启程度,以此形成一条为流体提供流动空间、具有一定导流能力的裂缝,从而达到改善油气渗流条件和油气增产的目的。
裂缝导流能力是评价流体在裂缝中流动性的指标,裂缝导流能力越高,说明同一种流体在裂缝中的流动性就越好。裂缝导流能力包括初始导流能力和导流能力变化率,初始导流能力是指裂缝在不受闭合压力作用,两个裂缝表面刚好接触但不相互挤压状态下的导流能力,导流能力变化率是指在闭合压力作用下,导流能力随着闭合压力的增大而减少的速度。
酸刻蚀物理模拟实验是指在实验室利用岩石样品(岩板)模拟现场酸压施工条件下酸液与裂缝表面岩石的反应,实验岩板表面即为裂缝表面。实验的目的是获取用作酸蚀裂缝导流能力评价的测试样品(岩板),对酸压效果进行预测,以优化酸压施工设计。
目前,室内预测酸压裂缝导流能力的方法是:利用目的储集层的岩板开展物理模拟实验,通过导流能力测试实验测得岩板的初始导流能力,再利用以下N-K方程(Nierode,D.E.andKruck,K.F.AnEvaluationofacidfluidlossadditives,retardedacids,andacidizedfractureconductivity[C].FallMeetingoftheSocietyofPetroleumEngineersofAIME,LasVegas,Nevada,30September-3October1973.SPE-4549-MS.)计算不同闭合压力作用下的导流能力:
kfw=(kfw)0e-0.006895Cσ
式中kfw—不同闭合压力下的裂缝导流能力,10-3μm2·cm;
(kfw)0—裂缝初始导流能力,cm;
σ—地层闭合压力,MPa;
Sf—岩石嵌入强度,MPa。
由此可见,只要获得室内实验样品的初始导流能力,即可大致预测地层中酸压裂缝的导流能力,对酸压施工设计提供指导。然而,室内初始导流能力测试实验流程复杂,实验设备沉重,操作过程不仅繁索,且存在一定的安全隐患。因此,设计更加简单、易于操作的初始导流能力计算方法将会对酸压设计起到重要推进作用。
发明内容
本发明的目的在于提供一种酸压裂缝初始导流能力计算方法,该方法操作简便,原理可靠,安全性好,能够为酸压方案设计提供有效指导。
为达到以上技术目的,本发明提供以下技术方案。
酸压裂缝表面的凹凸不平形成了流体的流动通道,裂缝表面的几何形态是影响导流能力的重要因素。酸液刻蚀裂缝表面后形成的刻蚀形态一般分为均匀刻蚀(表面溶蚀均匀,无明显凹凸起伏)、点状刻蚀(各点溶蚀有差异,存在无规律的凹凸起伏形态)和沟槽刻蚀(表面存在明显的沟槽形态)三种。由于初始导流能力是闭合压力为零时的导流能力,而闭合压力为零时,导流能力只受到裂缝初始表面几何形态的影响,故本发明首先通过分析裂缝表面初始几何形态,定义了具有明确物理意义的裂缝表面粗糙程度的表征参数,再利用表征参数获取了裂缝初始导流能力计算模型。
一种酸压裂缝初始导流能力计算方法,依次包括以下步骤:
(A)开展酸刻蚀物理模拟实验,利用三维激光扫描仪对酸刻蚀后的每个岩板粗糙表面进行扫描,获取岩板表面的三维数据。
(B)利用步骤(A)中得到的三维数据计算每个岩板表面的横向曲折比和纵向曲折比两个参数。
(C)计算裂缝(两个岩板组成)的平均横向曲折比和平均纵向曲折比
(D)利用步骤(C)中得到的和两个参数代入以下公式计算裂缝初始导流能力:
(kfw)0—裂缝初始导流能力,μm2·cm。
所述步骤(A)中获取岩板表面的三维数据有X、Y、Z三列,其中列X和列Y为岩板表面上各点的平面坐标,列Z为岩板表面各点的高度值。
所述的步骤(B)中参数和的计算采用以下公式:
式中:ΔLx—三维激光扫描仪在X方向上的扫描步长;
ΔLy—三维激光扫描仪在Y方向上的扫描步长;
Nx—三维激光扫描仪在X方向上的扫描步数;
Ny—三维激光扫描仪在Y方向上的扫描步数;
z(i,j)—岩板表面上x=i,y=j所在的点的高度值。
参数的定义和计算公式的推导过程如下:
1)利用岩板表面三维数据在三维坐标系中可生成一个曲面(即岩板的粗糙表面)。其中,y=yj(0≤j≤Ny)对应的是一条沿X方向(横向)曲线。由于相邻两个扫描点之间的距离非常小,可利用勾股定理算得点(xi,yj)和点(xi+1,yj)之间的曲线段长度为:
2)通过累加,可算得y=yj对应的横向曲线长度为:
3)利用X方向的扫描总长度将上式进行无因次化为:
4)上述计算过程只处理了y=yj一条横向曲线,计算整个曲面上的所有横向曲线的平均值,即得到的计算公式为:
的推导过程与的推导过程一致。
之所以定义了和两个参数,是因为经过对表面形态的分析,发现它们具有以下物理意义:
情况一:当和的值都很小表明岩样表面在横向和纵向上起伏变化较小,呈现均匀刻蚀形态。这种情况下,两个裂缝面刚好接触时,缝中空间小,流体流动通道少。
情况二:当较大,较小代表岩样表面在Y轴方向(纵向)上的流动沟槽比较发育,而在X轴方向(横向)上的流动通道受阻。由于流体在裂缝表面是沿横流动,这种情况会阻碍流体流动,减小裂缝导流能力。
情况三:当较小 较大,与情况二恰恰相反:岩样表面在横向上的流动沟槽比较发育,而在纵向上的流动通道受阻。此时,由于流体沿横向沟槽流动,故导流能力较大。
情况四:当和的值都较大,表明岩样表面在横向和纵向上起伏变化都较大,其表面发育微小的“凸起”和“低谷”,呈现严重程度的点状刻蚀形态。
所述的步骤(C)中,设两个岩板表面的横向曲折比分别是和则该裂缝的平均横向曲折比计算为:
同理,设两个岩板表面的横向曲折比分别是和则该裂缝的平均纵向曲折比计算为:
所述步骤(D)中,裂缝初始导流能力计算公式的拟合过程包括:
首先,根据如前所述的四种情况分析可知,随着值的增大,裂缝表面在Y轴方向(纵向)上的流动沟槽越来越发育。由于流体在裂缝表面是沿横流动,这种情况会阻碍流体流动,减小裂缝导流能力。而随着值的增大,岩样表面在横向上的流动沟槽越来越发育,而在纵向上的流动通道受阻。此时,由于流体沿横向沟槽流动,故这种情况会增大裂缝导流能力。
所以,根据初始导流能力随着和值的变化情趋式,拟定了下表中6种拟合模型,利用这6种模型对50对岩板的初始导流能力实验数据进行曲面拟合。结果表明,模型6对实验数据的曲面拟合度高达0.92(最大值为1),故选定模型6作为初始导流能力计算的经验公式,求出待解系数后即为步骤(D)中的裂缝初始导流能力计算公式。
6种模型对50对岩板初始导流能力实验数据的曲面拟合度
具体实施方式
下面对本发明的步骤进行详细说明。
随机选取10对岩板,分别编号为:1#、2#、3#、4#、5#、6#、7#、8#、9#、10#,按本发明的步骤依次计算每对岩板的初始导流能力。
一种酸压裂缝初始导流能力计算方法,依次包括如下步骤:
(A)开展酸刻蚀物理模拟实验,利用三维激光扫描仪对酸刻蚀后的每个岩板粗糙表面进行扫描,获取岩板表面的三维数据。
三维数据包括X、Y、Z三列,其中列X和列Y为岩板表面上各点的平面坐标,列Z为岩板表面各点的高度值。
(B)利用步骤(A)中得到的三维数据计算每个岩板表面的横向曲折比和纵向曲折比两个参数。
(C)计算裂缝(两个岩板组成)的平均横向曲折比和平均纵向曲折比
(D)利用步骤(C)中得到的和两个参数代入以下公式计算裂缝初始导流能力:
(kfw)0—裂缝初始导流能力,μm2·cm
为了验证该计算方法的精度,本实例又利用这10对岩板开展了初始导流能力测试实验,对利用实验数据和以上计算所得的理论数据进行了对比,采用相对误差评价了该方法的精确程度。如下表,(kfw)0i表示实验测得的初始导流能力数据。
上表可看出,本发明的计算值与实验真实值的相对误差不超过9%,最低仅有0.4%,平均为4.7%。说明了本发明具有较高的精度。
Claims (2)
1.一种酸压裂缝初始导流能力计算方法,依次包括以下步骤:
(A)利用三维激光扫描仪对酸刻蚀后的每个岩板粗糙表面进行扫描,获取岩板表面的三维数据,岩板表面的三维数据有X、Y、Z三列,其中列X和列Y为岩板表面上各点的平面坐标,列Z为岩板表面各点的高度值;
(B)利用步骤(A)中得到的三维数据计算每个岩板表面的横向曲折比和纵向曲折比
式中:ΔLx—三维激光扫描仪在X方向上的扫描步长,
ΔLy—三维激光扫描仪在Y方向上的扫描步长,
Nx—三维激光扫描仪在X方向上的扫描步数,
Ny—三维激光扫描仪在Y方向上的扫描步数,
z(i,j)—岩板表面上x=i,y=j所在的点的高度值;
(C)计算两个岩板组成的裂缝的平均横向曲折比和平均纵向曲折比
(D)将和两个参数代入以下公式计算裂缝初始导流能力:
式中:(kfw)0—裂缝初始导流能力,μm2·cm。
2.如权利要求1所述的一种酸压裂缝初始导流能力计算方法,其特征在于,所述的步骤(C)中,裂缝的平均横向曲折比和裂缝的平均纵向曲折比计算如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610051684.7A CN105718745A (zh) | 2016-01-26 | 2016-01-26 | 一种酸压裂缝初始导流能力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610051684.7A CN105718745A (zh) | 2016-01-26 | 2016-01-26 | 一种酸压裂缝初始导流能力计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105718745A true CN105718745A (zh) | 2016-06-29 |
Family
ID=56154921
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610051684.7A Pending CN105718745A (zh) | 2016-01-26 | 2016-01-26 | 一种酸压裂缝初始导流能力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105718745A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106203699A (zh) * | 2016-07-11 | 2016-12-07 | 西南石油大学 | 一种粗糙裂缝初始导流能力的预测方法 |
CN108152183A (zh) * | 2016-12-06 | 2018-06-12 | 中国石油天然气股份有限公司 | 酸蚀裂缝刻蚀形态及导流能力测试方法 |
CN108489809A (zh) * | 2018-03-06 | 2018-09-04 | 西南石油大学 | 利用实验手段计算应力作用下粗糙错位裂缝变形量的方法 |
CN108871241A (zh) * | 2018-05-10 | 2018-11-23 | 西南石油大学 | 一种定量评价粗糙裂缝酸刻蚀后粗糙程度变化的方法 |
CN110909477A (zh) * | 2019-11-28 | 2020-03-24 | 西南石油大学 | 一种酸蚀裂缝粗糙程度定量表征方法 |
CN112502685A (zh) * | 2020-12-03 | 2021-03-16 | 西南石油大学 | 一种考虑热效应的碳酸盐岩储层交替酸压级数优选方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002183305A (ja) * | 2000-12-11 | 2002-06-28 | Eidai Yakuhin Kk | 医薬品の在庫管理方法 |
CN104989364A (zh) * | 2015-06-11 | 2015-10-21 | 西南石油大学 | 酸刻蚀物理模拟实验中岩板表面刻蚀形态的评价方法 |
-
2016
- 2016-01-26 CN CN201610051684.7A patent/CN105718745A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002183305A (ja) * | 2000-12-11 | 2002-06-28 | Eidai Yakuhin Kk | 医薬品の在庫管理方法 |
CN104989364A (zh) * | 2015-06-11 | 2015-10-21 | 西南石油大学 | 酸刻蚀物理模拟实验中岩板表面刻蚀形态的评价方法 |
Non-Patent Citations (1)
Title |
---|
白翔: "基于刻蚀形态数字化表征的酸蚀裂缝导流能力研究", 《中国优秀硕士学位论文全文数据库工程科技I辑》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106203699A (zh) * | 2016-07-11 | 2016-12-07 | 西南石油大学 | 一种粗糙裂缝初始导流能力的预测方法 |
CN106203699B (zh) * | 2016-07-11 | 2019-07-23 | 西南石油大学 | 一种粗糙裂缝初始导流能力的预测方法 |
CN108152183A (zh) * | 2016-12-06 | 2018-06-12 | 中国石油天然气股份有限公司 | 酸蚀裂缝刻蚀形态及导流能力测试方法 |
CN108489809A (zh) * | 2018-03-06 | 2018-09-04 | 西南石油大学 | 利用实验手段计算应力作用下粗糙错位裂缝变形量的方法 |
CN108871241A (zh) * | 2018-05-10 | 2018-11-23 | 西南石油大学 | 一种定量评价粗糙裂缝酸刻蚀后粗糙程度变化的方法 |
CN108871241B (zh) * | 2018-05-10 | 2019-04-26 | 西南石油大学 | 一种定量评价粗糙裂缝酸刻蚀后粗糙程度变化的方法 |
CN110909477A (zh) * | 2019-11-28 | 2020-03-24 | 西南石油大学 | 一种酸蚀裂缝粗糙程度定量表征方法 |
CN110909477B (zh) * | 2019-11-28 | 2020-08-11 | 西南石油大学 | 一种酸蚀裂缝粗糙程度定量表征方法 |
CN112502685A (zh) * | 2020-12-03 | 2021-03-16 | 西南石油大学 | 一种考虑热效应的碳酸盐岩储层交替酸压级数优选方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105718745A (zh) | 一种酸压裂缝初始导流能力计算方法 | |
CN104989364B (zh) | 酸刻蚀物理模拟实验中岩板表面刻蚀形态的评价方法 | |
Williams et al. | Design of acid fracturing treatments | |
Kuhnle et al. | Measured and simulated flow near a submerged spur dike | |
US8798970B2 (en) | Method for oil prediction in fractured reservoirs and recorded medium therefor | |
Zhang et al. | Three-dimensional simulation of meandering river based on 3-D RNG κ-ɛ turbulence model | |
CN110147561B (zh) | 一种含天然裂缝致密油气储层体积压裂缝网预测方法 | |
EP2831804B1 (en) | System and method for automatic local grid refinement in reservoir simulation systems | |
CN105066912B (zh) | 酸刻蚀物理模拟实验中岩板表面扫描数据的步长标定方法 | |
CN103437381B (zh) | 群桩拖带叠加沉降量评估处理方法 | |
Stachowiak et al. | Analysis of the tribocorrosion mechanisms in a pin-on-plate combination on the example of AISI304 steel | |
CN106372320B (zh) | 一种采用亚滤波尺度模型对公路隧道湍流进行大涡模拟的方法 | |
CN108711145A (zh) | 一种基于三维图像分析的腐蚀形貌描述解释方法 | |
Saleh et al. | Minimizing the hydraulic side effects of weirs construction by using labyrinth weirs | |
CN114722578B (zh) | 一种隧道地表沉降计算方法 | |
CN102141385B (zh) | 沥青路面表面形貌曲面的测试方法 | |
Wu et al. | A model for the conductivity and compliance of unpropped and natural fractures | |
CN106203699B (zh) | 一种粗糙裂缝初始导流能力的预测方法 | |
Sadeghnejad et al. | Field scale characterization of geological formations using percolation theory | |
CN105929117B (zh) | 一种粗糙裂缝流道复杂程度的评价方法 | |
Baghlani | Application of a high-resolution scheme in simulation of flow in curved channel using boundary-fitted curvilinear coordinates | |
CN102472690B (zh) | 用于预测轮胎和道路表面之间的物理相互作用效应的方法 | |
CN114969886A (zh) | 一种钢轨滚动接触疲劳预测方法 | |
CN113435075A (zh) | 一种模拟海底能源管线遭受滑坡灾害冲击时引起荷载响应的计算机方法 | |
Bruzzone et al. | An experimental evaluation of an etching simulation model for photochemical machining |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160629 |
|
WD01 | Invention patent application deemed withdrawn after publication |