CN106484958A - 一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型 - Google Patents
一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型 Download PDFInfo
- Publication number
- CN106484958A CN106484958A CN201610832079.3A CN201610832079A CN106484958A CN 106484958 A CN106484958 A CN 106484958A CN 201610832079 A CN201610832079 A CN 201610832079A CN 106484958 A CN106484958 A CN 106484958A
- Authority
- CN
- China
- Prior art keywords
- theta
- sin
- permeability
- crack
- cos
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Remote Sensing (AREA)
- Computer Hardware Design (AREA)
- Computer Graphics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明涉及油气田勘探开发领域,尤其是一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型。本发明基于岩心统计、薄片观测以及测井解释,通过现今地应力分析,确定裂缝地下开度,利用古地磁定向、岩层产状法或地层倾角法确定裂缝的三维展布,进而建立真三维裂缝渗透率计算模型,推导可编程算法,编制相应的程序,准确的计算三维等渗透率面、最大差渗透率面、不同方向的渗透率剖面等。本发明对于真三维空间裂缝渗透率张量及井网设计、确定注水井与采油井的空间位置关系、寻找最优钻井轨迹等油田开发方案的部署具有较高的实用价值,并且预测成本低廉、可操作性强,预测结果可靠。
Description
技术领域
本发明涉及油气田勘探开发领域,尤其是一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型。
背景技术
在裂缝性储层勘探开发过程中,天然裂缝的渗流方向与裂缝的产状、充填性以及现今地应力大小、方向密切相关,裂缝的主渗流方向影响井网部署、注水开发等措施实施。在油藏中,天然裂缝分布极为复杂;在一个单元内,不同期次、不同展布以及不同开度的裂缝在油气导流中所起的作用不一,如何考虑三维的地质因素,输出三维的裂缝渗透率张量是裂缝描述的一大难点。刘月田等(刘月田丁祖鹏屈亚光赵辰军.油藏裂缝方向表征及渗透率各向异性参数计算[J].石油学报,2011,32(5):842-846.)所提的裂缝渗透率计算方法仅考虑了特殊情况下(剪切裂缝,对称,等开度)裂缝渗透张量的计算,并且该方法也无法定量输出主渗透率方向;刘敬寿等(刘敬寿,戴俊生,邹娟,杨海盟,汪必峰,周巨标.裂缝性储层渗透率张量定量计算方法[J].石油与天然气地质,2015,36(6):1022-1029.)所提的裂缝渗透率张量计算方法考虑的三维的因素,只能输出二维的结果,无法准确的计算三维等渗透率面、最大差渗透率面、不同方向的渗透率剖面等。本发明专利基于岩心统计、薄片观测以及测井解释,在应力场模拟的基础上,确定裂缝地下开度,利用古地磁定向、岩层产状法或地层倾角法确定裂缝的三维展布,进而建立真三维裂缝渗透率计算模型,在推导相应算法的基础上,编制相应的程序,实现准确的计算三维等渗透率面、最大差渗透率面、不同方向的渗透率剖面等。
本发明旨在解决上述问题,提供了一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型,它解决了无法准确输出井筒裂缝真三维空间等渗透率面、最大差渗透率面以及不同方向的渗透率剖面的问题。
本发明的技术方案为:一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型,具体步骤如下:
第一步 井筒裂缝发育程度的观测
利用岩心观测、测井资料解释对裂缝的开度测量,确定裂缝的组系特征;依据成像资料、测井或者岩心观察,求取裂缝的线密度,为后面计算裂缝平行渗透率准备。
第二步 裂缝定向
在岩心裂缝观测中,对非定向岩心采用古地磁定向、岩层产状法或地层倾角法确定每条裂缝的产状;在成像资料、声波测井解释裂缝时,利用地层倾角测井或者成像测井资料对裂缝定向。
第三步 应力场模拟确定裂缝地下的开度
利用水力压裂法、波速法、声发射法、应力解除法、光弹性应力法等可以对工区现今地应力方位进行了判断,并通过测井、压裂资料以及物理实验计算了关键井现今地应力的数值,利用地震反演、测井解释、岩石力学实验确定岩石的力学属性。在此基础上,利用软件建立有限元模型,采用多井约束的方法对研究区现今地应力场进行了数值模拟。利用模拟结果,确定井筒裂缝发育处的现今地应力大小、方向;国外专家Willis-Richards等和Jing等同时考虑了正应力及剪应力对裂缝开度的影响,得出在现今地应力场改造后的裂缝开度计算公式:
公式(1)中,b为裂缝的现今地下实际开度,单位:m;b0为裂缝的原始开度,单位:m;σ′n为有效正应力,单位:MPa;bres代表裂缝面承受最大有效正应力时的裂缝开度,单位:m;σnref为使裂缝开度减小90%的有效正应力,单位:MPa,是一个与岩性相关的系数,已有学者[Willis-Richards J,1996;Jing Z,1998;Durham WB,1994;Chen Z,2000;秦积舜;2002;]给出了不同条件下的测试数值。
第四步 建立真三维裂缝渗透率张量计算几何模型
如图2所示,设裂缝面Φ的倾角为δ,倾向为ω,以裂缝为参照物建立O-ABC坐标系(O-ABC),O-ABC坐标系中的三个坐标轴分别对应于裂缝面的法线方向(OA轴)、裂缝走向方向(OB轴)、裂缝面内垂直于裂缝走向线的方向(OC轴)。设任意渗流面ψ的倾角为α,倾向为β,其单位法向矢量为m;以渗流面ψ为参照物建立坐标系(O-XYZ),其中,OX轴为渗流面ψ的法线方向,OY轴、OZ轴位于渗流面ψ内,OP为渗流面ψ的走向线,HQ为渗流面ψ的倾斜线,HR为渗流面ψ倾斜线的水平投影;线OS、ON、OE、OW、HR、OH、OP、OB在水平面内,定义θ为渗流面ψ内OY轴与OP轴的夹角,通过调整θ的大小,可以求取裂缝在θ中不同方向的渗透率。
第五步 真三维裂缝渗透率张量数学模型
在已知裂缝线密度、地下开度的基础上,单组(单条)裂缝的平行渗透率K可以表示为:
公式(2)中,b为裂缝的现今地下实际开度,单位:m;Dlf为单组裂缝的线密度,条/m。
设第i条裂缝的渗透率张量为K,则其在坐标系O-XYZ中的表达式可表示为:
同理,渗透率张量K在坐标系O-ABC中的表达式可表示为:
OA轴在大地坐标系中的三个分量为:
OX轴在大地坐标系中的三个分量为:
设裂缝的渗透率张量为K,其在O-ABC坐标系中的表达式可表示为:
同理,渗透率张量K则其在O-XYZ坐标系中的表达式可表示为:
在O-XYZ坐标系中,每个单元体内单组裂缝的渗透率张量KABC可表示为:
由公式(7)-(9)可以得到:
KXYZi=T·KABC·TT (10)
其中,O-ABC坐标系转换为O-XYZ坐标系的旋转矩阵T可表示为:
公式(11)的旋转矩阵T中各参数分别表示为:
在O-XYZ坐标系中,单元体内每组(每条)裂缝的渗透率张量KXYZi可表示为:
在复杂的地质条件下,储层裂缝经历多期构造应力场演化,每组裂缝的产状、线密度、开度往往不同,依据公式(10)、公式(11)及公式(13),单元体内发育多组裂缝时,渗透率张量KXYZ可表示为:
公式(14)中,k为单元体内裂缝的组数;bi为第i组裂缝的开度,单位:m;Ti为第i组裂缝在O-ABC坐标轴分量转换为O-XYZ坐标轴分量的旋转矩阵,单位:m;Dlfi为第i组裂缝的线密度,单位:条/m。
第六步 编制程序,输出渗透率剖面
利用公式(1)-(14)推导的算法编制相应的程序,通过循环迭代变量α(渗流面ψ的倾角)、β(渗流面ψ的倾向)、θ(OY轴与OP轴的夹角),其中的α范围为[0-90°],β的范围为[0-360°),θ的范围为[0-90°),根据计算精度需求,设置循环迭代步长(⊿sp,例如:5°,2°,1°,0.5°,0.2°,0.1°等)将α、β、θ逐次累加⊿sp,利用公式(12)-(14),计算坐标系(O-XYZ)不同的α、β、θ对应的KXXαβθ(X方向渗透率主值)、KYYαβθ(Y方向渗透率主值)、KZZαβθ(Z方向渗透率主值)。
利用公式(1)-(14),设置α=0°,β为[0-360°)的任意数值;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ(Y方向渗透率主值)、KZZαβθ(Z方向渗透率主值),得到水平面的渗透率变化。
利用公式(1)-(14),设置α=90°,β=0°;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ(Y方向渗透率主值)、KZZαβθ(Z方向渗透率主值),得到东西向剖面的渗透率变化。
利用公式(1)-(14),设置α=90°,β=90°;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ(Y方向渗透率主值)、KZZαβθ(Z方向渗透率主值),得到南北向剖面的渗透率变化。
利用公式(1)-(14),设置α范围为[0-90°],β的范围为[0-360°),步长,⊿sp;计算不同的α、β(渗流面ψ确定)的平面内——即不同的渗流面ψ内,KYYαβθ(Y方向渗透率主值)、KZZαβθ(Z方向渗透率主值),每个渗流面ψ内计算得到的KYYαβθ、KZZαβθ数目为2n,2n个KYYαβθ、KZZαβθ的平均值为Kaver;计算不同的α、β对应的方差函数Wαβ:
公式(15)中,KYYαβθi、KZZαβθi分别为特定的α、β对应的第i个Y方向渗透率主值、Z方向渗透率主值,单位:10-3μm2;n为每个渗流面ψ计算的Y或Z方向渗透率主值的数目,单位:个;Kaver为2n个KYYαβθ、KZZαβθ的平均值,单位:10-3μm2。
通过公式(15),计算不同的α、β对应的方差函数Wαβ,比较不同的α、β的Wαβ。其中,Wαβ最大值对应的渗流面ψ为最大差渗透率面,记录对应的倾角α,倾向β;Wαβ最小值对应的渗流面ψ为等渗透率面,记录对应的倾角α,倾向β。
本发明的有益效果是:本发明专利基于岩心统计、薄片观测以及测井解释,通过现今地应力分析,确定裂缝地下开度,利用古地磁定向、岩层产状法或地层倾角法确定裂缝的三维展布,进而建立真三维裂缝渗透率计算模型,推导可编程算法,编制相应的程序,实现准确的计算三维等渗透率面、最大差渗透率面、不同方向的渗透率剖面等。本发明由严格的数学算法推导组成,对相应的地质信息数字化后,可以利用计算机编程语言开发相应的计算程序,实现井筒裂缝真三维空间渗透率张量定量计算。本发明对于真三维空间裂缝渗透率张量及井网设计、确定注水井与采油井的空间位置关系、寻找最优钻井轨迹等油田开发方案的部署具有较高的实用价值,并且预测成本低廉、可操作性强,预测结果可靠。
附图说明
图1为一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型的流程图。
图2为一种基于井筒裂缝观测的三维裂缝渗透率张量数学模型示意图。
图3为苏北盆地金湖凹陷铜城断裂带天33断块构造位置图。
图4为天33断块阜二段岩心裂缝开度统计分布图。
图5为天33断块阜二段岩心裂缝走向玫瑰花图。
图6为天33断块水平剖面裂缝渗透率变化图。
图7为天33断块东西剖面裂缝渗透率变化图。
图8为天33断块南北剖面裂缝渗透率变化图。
图9为天33断块最大差渗透率面(252°∠46°)变化图。
图10为天33断块剖面等渗透率面(350°∠81°)变化图。
具体实施方式
下面结合附图说明本发明的具体实施方式:
以苏北盆地金湖凹陷铜城断裂带东翼阜宁组二段(简称阜二段)岩心裂缝三维渗透率张量定量计算为例来说明本发明的具体技术方案:
天33断块地处江苏金湖县境内,研究区主要位于苏北盆地-东台坳陷-金湖凹陷-铜城断裂带的东翼。铜城断层是一条典型的走滑断层,南部交汇于金湖凹陷的边界-杨村断层,向北消失于铜城地区(图3)。研究区主力产油层为阜二段,储层物性差,岩性以低渗透粉砂岩为主。裂缝对油气的渗流起到至关重要的作用。结合区域地质资料,对天33断块阜二段岩心裂缝渗透率张量定量计算的步骤如下:
第1步 对天33断块井筒裂缝发育程度的观测,利用岩心观测、测井资料解释对裂缝的开度测量(图4),确定裂缝的组系特征;依据成像资料、测井或者岩心观察,求取裂缝的线密度,为后面计算裂缝平行渗透率准备。
第2步 采用斜井岩心裂缝产状校正方法对工区的裂缝进行定向(图5)。
第3步 通过确定岩石力学参数、井点地应力状态,建立地质力学模型,利用公式(1)计算裂缝在地下的开度。
第4步 如图2所示,建立真三维裂缝渗透率张量计算几何模型。
第5步 建立真三维裂缝渗透率张量数学模型,如公式(2)-公式(14)所示。
第6步 编制程序,输出渗透率剖面,利用公式(1)-(14)推导的算法编制相应的程序,通过循环迭代变量α(渗流面ψ的倾角)、β(渗流面ψ的倾向)、θ(OY轴与OP轴的夹角),其中的α范围为[0-90°],β的范围为[0-360°),θ的范围为[0-90°),根据计算精度需求,设置循环迭代步长⊿sp=1°,将α、β、θ逐次累加⊿sp,利用公式(12)-(14),计算坐标系(O-XYZ)不同的α、β、θ对应的KXXαβθ(X方向渗透率主值)、KYYαβθ(Y方向渗透率主值)、KZZαβθ(Z方向渗透率主值)。
利用公式(1)-(14),设置α=0°,β=0°;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ(Y方向渗透率主值)、KZZαβθ(Z方向渗透率主值),得到水平面的渗透率变化(图6)。
利用公式(1)-(14),设置α=90°,β=0°;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ(Y方向渗透率主值)、KZZαβθ(Z方向渗透率主值),得到东西向剖面的渗透率变化(图7)。
利用公式(1)-(14),设置α=90°,β=90°;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ(Y方向渗透率主值)、KZZαβθ(Z方向渗透率主值),得到南北向剖面的渗透率变化(图8)。
利用公式(1)-(14),设置α范围为[0-90°],β的范围为[0-360°),计算不同的渗流面ψ内计算得到的KYYαβθ、KZZαβθ数目为180,180个KYYαβθ、KZZαβθ的平均值为Kaver;计算不同的α、β对应的方差函数Wαβ;通过公式(15),计算不同的α、β对应的方差函数Wαβ,比较不同的α、β的Wαβ。其中,Wαβ最大值对应的渗流面ψ为最大差渗透率面(图9);Wαβ最小值对应的渗流面ψ为等渗透率面(图10)。
上面以举例方式对本发明进行了说明,但本发明不限于上述具体实施例,凡基于本发明所做的任何改动或变型均属于本发明要求保护的范围。
Claims (2)
1.一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型,所述的步骤如下:
1)利用岩心观测、测井资料解释对裂缝的开度测量,确定裂缝的组系特征;依据成像资料、测井或者岩心观察,求取裂缝的线密度;
2)在岩心裂缝观测中,对非定向岩心采用古地磁定向、岩层产状法或地层倾角法确定每条裂缝的产状;在成像资料、声波测井解释裂缝时,利用地层倾角测井或者成像测井资料对裂缝定向;
3)利用水力压裂法、波速法、声发射法、应力解除法、光弹性应力法等方法对工区现今地应力方位进行了判断,并通过测井、压裂资料以及物理实验计算了关键井现今地应力的数值,利用地震反演、测井解释、岩石力学实验确定岩石的力学属性;在此基础上,利用软件建立有限元模型,采用多井多层段实测地应力的大小、方向为约束条件,对研究区现今地应力场进行了数值模拟;利用模拟结果,确定井筒裂缝发育处的现今地应力大小、方向;利用公式(1)计算裂缝在地下的开度:
公式(1)中,b为裂缝的现今地下实际开度,单位:m;b0为裂缝的原始开度,单位:m;σ′n为有效正应力,单位:MPa;bres代表裂缝面承受最大有效正应力时的裂缝开度,单位:m;σnref为使裂缝开度减小90%的有效正应力,单位:MPa;
4)如图2所示,建立真三维裂缝渗透率张量计算几何模型,设裂缝面Φ的倾角为δ,倾向为ω,以裂缝为参照物建立O-ABC坐标系(O-ABC),O-ABC坐标系中的三个坐标轴分别对应于裂缝面的法线方向(OA轴)、裂缝走向方向(OB轴)、裂缝面内垂直于裂缝走向线的方向(OC轴);设任意渗流面ψ的倾角为α,倾向为β,其单位法向矢量为m;以渗流面ψ为参照物建立坐标系(O-XYZ),其中,OX轴为渗流面ψ的法线方向,OY轴、OZ轴位于渗流面ψ内,OP为渗流面ψ的走向线,HQ为渗流面ψ的倾斜线,HR为渗流面ψ倾斜线的水平投影;线OS、ON、OE、OW、HR、OH、OP、OB在水平面内,定义θ为渗流面ψ内OY轴与OP轴的夹角,通过调整θ的大小,可以求取裂缝在θ中不同方向的渗透率;
5)建立真三维裂缝渗透率张量数学模型:在已知裂缝线密度、地下开度的基础上,单组(单条)裂缝的平行渗透率K可以表示为:
公式(2)中,b为裂缝的现今地下实际开度,单位:m;Dlf为单组裂缝的线密度,条/m;
设第i条裂缝的渗透率张量为K,则其在坐标系O-XYZ中的表达式可表示为:
同理,渗透率张量K在坐标系O-ABC中的表达式可表示为:
OA轴在大地坐标系中的三个分量为:
OX轴在大地坐标系中的三个分量为:
设裂缝的渗透率张量为K,其在O-ABC坐标系中的表达式可表示为:
同理,渗透率张量K则其在O-XYZ坐标系中的表达式可表示为:
在O-XYZ坐标系中,每个单元体内单组裂缝的渗透率张量KABC可表示为:
由公式(7)-(9)可以得到:
KXYZi=T·KABC·TT (10)
其中,O-ABC坐标系转换为O-XYZ坐标系的旋转矩阵T可表示为:
公式(11)的旋转矩阵T中各参数分别表示为:
在O-XYZ坐标系中,单元体内每组(每条)裂缝的渗透率张量KXYZi可表示为:
在复杂的地质条件下,储层裂缝经历多期构造应力场演化,每组裂缝的产状、线密度、开度往往不同,依据公式(10)、公式(11)及公式(13),单元体内发育多组裂缝时,渗透率张量KXYZ可表示为:
公式(14)中,k为单元体内裂缝的组数;bi为第i组裂缝的开度,单位:m;Ti为第i组裂缝在O-ABC坐标轴分量转换为O-XYZ坐标轴分量的旋转矩阵,单位:m;Dlfi为第i组裂缝的线密度,单位:条/m;
6)编制程序,输出渗透率剖面。
2.根据权利要求1所述的一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型,其特征在于:
所述的编制程序,输出渗透率剖面,其思路如下:
利用公式(1)-(14)推导的算法编制相应的程序,通过循环迭代变量α(渗流面ψ的倾角)、β(渗流面ψ的倾向)、θ(OY轴与OP轴的夹角),其中的α范围为[0-90°],β的范围为[0-360°),θ的范围为[0-90°),根据计算精度需求,设置循环迭代步长⊿sp,将α、β、θ逐次累加⊿sp,利用公式(12)-(14),计算坐标系(O-XYZ)不同的α、β、θ对应的KXXαβθ(X方向渗透率主值)、KYYαβθ(Y方向渗透率主值)、KZZαβθ(Z方向渗透率主值);
利用公式(1)-(14),设置α=0°,β为[0-360°)的任意数值;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ(Y方向渗透率主值)、KZZαβθ(Z方向渗透率主值),得到水平面的渗透率变化;
利用公式(1)-(14),设置α=90°,β=0°;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ(Y方向渗透率主值)、KZZαβθ(Z方向渗透率主值),得到东西向剖面的渗透率变化;
利用公式(1)-(14),设置α=90°,β=90°;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ(Y方向渗透率主值)、KZZαβθ(Z方向渗透率主值),得到南北向剖面的渗透率变化;
利用公式(1)-(14),设置α范围为[0-90°],β的范围为[0-360°),步长,⊿sp;计算不同的α、β(渗流面ψ确定)的平面内——即不同的渗流面ψ内,KYYαβθ(Y方向渗透率主值)、KZZαβθ(Z方向渗透率主值),每个渗流面ψ内计算得到的KYYαβθ、KZZαβθ数目为2n,2n个KYYαβθ、KZZαβθ的平均值为Kaver;计算不同的α、β对应的方差函数Wαβ:
公式(15)中,KYYαβθi、KZZαβθi分别为特定的α、β对应的第i个Y方向渗透率主值、Z方向渗透率主值,单位:10-3μm2;n为每个渗流面ψ计算的Y或Z方向渗透率主值的数目,单位:个;Kaver为2n个KYYαβθ、KZZαβθ的平均值,单位:10-3μm2;
通过公式(15),计算不同的α、β对应的方差函数Wαβ,比较不同的α、β的Wαβ;其中,Wαβ最大值对应的渗流面ψ为最大差渗透率面,记录对应的倾角α,倾向β;Wαβ最小值对应的渗流面ψ为等渗透率面,记录对应的倾角α,倾向β。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610832079.3A CN106484958B (zh) | 2016-09-19 | 2016-09-19 | 一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型的*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610832079.3A CN106484958B (zh) | 2016-09-19 | 2016-09-19 | 一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型的*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106484958A true CN106484958A (zh) | 2017-03-08 |
CN106484958B CN106484958B (zh) | 2019-09-17 |
Family
ID=58268743
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610832079.3A Expired - Fee Related CN106484958B (zh) | 2016-09-19 | 2016-09-19 | 一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型的*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106484958B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107132586A (zh) * | 2017-04-19 | 2017-09-05 | 中国石油大学(华东) | 一种断陷盆地断裂带结构单元定量识别方法 |
CN108008117A (zh) * | 2017-12-21 | 2018-05-08 | 刘敬寿 | 裂缝发育模式与储层基质-裂缝渗透率敏感性预测方法 |
CN108061699A (zh) * | 2017-12-14 | 2018-05-22 | 刘敬寿 | 一种多期次、多尺度裂缝孔渗参数定量预测方法 |
CN109033538A (zh) * | 2018-06-30 | 2018-12-18 | 南京理工大学 | 一种基于实测结构面参数的裂隙岩体渗透张量的计算方法 |
CN109493397A (zh) * | 2018-09-30 | 2019-03-19 | 中国石油天然气股份有限公司 | 一种钻井数据处理方法、装置及*** |
CN109614634A (zh) * | 2018-10-25 | 2019-04-12 | 中国辐射防护研究院 | 一种水环境中放射性核素迁移扩散数值模拟方法及*** |
CN110020453A (zh) * | 2018-01-09 | 2019-07-16 | 姚君波 | 地层渗透率三维分解模拟器 |
US20220342113A1 (en) * | 2021-04-21 | 2022-10-27 | Baker Hughes Oilfield Operations Llc | Estimation of properties of a subterranean region using a synthetic physical model |
CN116956670A (zh) * | 2023-07-17 | 2023-10-27 | 长江大学 | 一种基于tpfa与mfd混合方法的投影嵌入式离散裂缝模型 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101942991A (zh) * | 2010-06-30 | 2011-01-12 | 中国石油大学(北京) | 裂缝各向异性油藏注水开发可预测物理模型建立方法 |
US20110087472A1 (en) * | 2009-04-30 | 2011-04-14 | Schlumberger Technology Corporation | Determining elastic and fluid flow properties of a fractured reservoir |
CN104200039A (zh) * | 2014-09-17 | 2014-12-10 | 中国石油大学(华东) | 一种构造裂缝产状定量预测方法 |
CN104500050A (zh) * | 2014-10-31 | 2015-04-08 | 中国石油大学(华东) | 一种裂缝性储层渗透率张量及各向异性定量预测方法 |
-
2016
- 2016-09-19 CN CN201610832079.3A patent/CN106484958B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110087472A1 (en) * | 2009-04-30 | 2011-04-14 | Schlumberger Technology Corporation | Determining elastic and fluid flow properties of a fractured reservoir |
CN101942991A (zh) * | 2010-06-30 | 2011-01-12 | 中国石油大学(北京) | 裂缝各向异性油藏注水开发可预测物理模型建立方法 |
CN104200039A (zh) * | 2014-09-17 | 2014-12-10 | 中国石油大学(华东) | 一种构造裂缝产状定量预测方法 |
CN104500050A (zh) * | 2014-10-31 | 2015-04-08 | 中国石油大学(华东) | 一种裂缝性储层渗透率张量及各向异性定量预测方法 |
Non-Patent Citations (2)
Title |
---|
ALIREZA JAFARI 等: "Estimation of equivalent fracture network permeability using fractal and statistical network properties", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 * |
王贵宾 等: "岩体节理三维网络模拟技术及渗透率张量分析", 《岩石力学与工程学报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107132586A (zh) * | 2017-04-19 | 2017-09-05 | 中国石油大学(华东) | 一种断陷盆地断裂带结构单元定量识别方法 |
CN107132586B (zh) * | 2017-04-19 | 2021-04-16 | 华北理工大学 | 一种断陷盆地断裂带结构单元定量识别方法 |
CN108061699B (zh) * | 2017-12-14 | 2020-01-21 | 中国石油大学(华东) | 一种多期次、多尺度裂缝孔渗参数定量预测方法 |
CN108061699A (zh) * | 2017-12-14 | 2018-05-22 | 刘敬寿 | 一种多期次、多尺度裂缝孔渗参数定量预测方法 |
CN108008117A (zh) * | 2017-12-21 | 2018-05-08 | 刘敬寿 | 裂缝发育模式与储层基质-裂缝渗透率敏感性预测方法 |
CN110020453B (zh) * | 2018-01-09 | 2023-03-10 | 姚君波 | 地层渗透率三维分解模拟器 |
CN110020453A (zh) * | 2018-01-09 | 2019-07-16 | 姚君波 | 地层渗透率三维分解模拟器 |
CN109033538A (zh) * | 2018-06-30 | 2018-12-18 | 南京理工大学 | 一种基于实测结构面参数的裂隙岩体渗透张量的计算方法 |
CN109033538B (zh) * | 2018-06-30 | 2022-07-22 | 南京理工大学 | 一种基于实测结构面参数的裂隙岩体渗透张量的计算方法 |
CN109493397A (zh) * | 2018-09-30 | 2019-03-19 | 中国石油天然气股份有限公司 | 一种钻井数据处理方法、装置及*** |
CN109614634A (zh) * | 2018-10-25 | 2019-04-12 | 中国辐射防护研究院 | 一种水环境中放射性核素迁移扩散数值模拟方法及*** |
CN109614634B (zh) * | 2018-10-25 | 2022-08-23 | 中国辐射防护研究院 | 一种水环境中放射性核素迁移扩散数值模拟方法及*** |
US20220342113A1 (en) * | 2021-04-21 | 2022-10-27 | Baker Hughes Oilfield Operations Llc | Estimation of properties of a subterranean region using a synthetic physical model |
CN116956670A (zh) * | 2023-07-17 | 2023-10-27 | 长江大学 | 一种基于tpfa与mfd混合方法的投影嵌入式离散裂缝模型 |
CN116956670B (zh) * | 2023-07-17 | 2024-01-23 | 长江大学 | 一种基于tpfa与mfd混合方法的投影嵌入式离散裂缝模型 |
Also Published As
Publication number | Publication date |
---|---|
CN106484958B (zh) | 2019-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106484958A (zh) | 一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型 | |
CN105672971B (zh) | 一种储层裂缝开启压力、开启次序及注水压力预测方法 | |
CN104500050B (zh) | 一种裂缝性储层渗透率张量及各向异性定量预测方法 | |
CN104747180B (zh) | 一种用于注水开发的缝洞型油藏分析方法及其应用 | |
Zhang et al. | Estimation of in situ stress along deep tunnels buried in complex geological conditions | |
CN106019405B (zh) | 储层裂缝建模方法及*** | |
Guo et al. | Simulation of three-dimensional tectonic stress fields and quantitative prediction of tectonic fracture within the Damintun Depression, Liaohe Basin, northeast China | |
CN104200039B (zh) | 一种构造裂缝产状定量预测方法 | |
CN104267429B (zh) | 确定地层压力的方法及装置 | |
CN106501856B (zh) | 亚地震断层定量预测方法 | |
CN104963657B (zh) | 基于高含水期监测数据约束的油藏数值模拟方法及装置 | |
CN106951660A (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
CN105822302A (zh) | 一种基于井地电位法的油水分布识别方法 | |
CN103914620B (zh) | 一种计算断层破裂带裂缝张开空间分布的方法 | |
CN107515423A (zh) | 一种高陡构造的三维地质模型构建方法 | |
CN115324560A (zh) | 利用地应力场模拟确定压裂诱发油气套管变形位置的方法 | |
CN109270578A (zh) | 滩浅海低品位油藏滚动开发方法 | |
Zhao et al. | Characteristics and application of present in situ stress field of a strike-slip fault: a 3D finite-element simulation study | |
Zhao et al. | Present in-situ stress distribution characteristics of strike-slip in SH Oilfield, Tarim Basin | |
Liu et al. | Geomechanical modeling and inversion Analysis of the in-situ stress field in deep marine shale formations: A case study of the Longmaxi Formation, Dingshan Area, China | |
Liu et al. | The Control Theory and Application for Well Pattern Optimization of Heterogeneous Sandstone Reservoirs | |
Han et al. | Development and calibrations of a coupled reservoir geomechanic model for Valhall field | |
Sosio et al. | Integrated geological, fluid flow and geomechanical model of a geothermal field | |
Luo et al. | Discrete element-based regression analysis of initial ground stress and application to an extra-long tunnel in China | |
Meng et al. | Interpretation of In-situ Stress at Baihetan project |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20190812 Address after: 266580 Qingdao Changjiang Road, Huangdao District, Shandong, No. 66 Applicant after: China Petroleum University (East China) Address before: 100083 Haidian District, Xueyuan Road, China University of Geosciences, (Beijing) Energy Institute Applicant before: China University of Geosciences (Beijing) |
|
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190917 Termination date: 20200919 |