CN106484958B - 一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型的系统 - Google Patents
一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型的系统 Download PDFInfo
- Publication number
- CN106484958B CN106484958B CN201610832079.3A CN201610832079A CN106484958B CN 106484958 B CN106484958 B CN 106484958B CN 201610832079 A CN201610832079 A CN 201610832079A CN 106484958 B CN106484958 B CN 106484958B
- Authority
- CN
- China
- Prior art keywords
- crack
- permeability
- formula
- seepage flow
- fracture
- 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.)
- Expired - Fee Related
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)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (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 W B,1994;Chen Z,2000;秦积舜;2002;]给出了不同条件下的测试数值。
第四步建立真三维裂缝渗透率张量计算几何模型
如图2所示,设裂缝面Φ的倾角为δ,倾向为ω,以裂缝为参照物建立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αβθ、KYYαβθ、KZZαβθ。
利用公式(1)-(14),设置α=0°,β为[0-360°)的任意数值;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ、KZZαβθ,得到水平面的渗透率变化。
利用公式(1)-(14),设置α=90°,β=0°;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ、KZZαβθ,得到东西向剖面的渗透率变化。
利用公式(1)-(14),设置α=90°,β=90°;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ、KZZαβθ,得到南北向剖面的渗透率变化。
利用公式(1)-(14),设置α范围为[0-90°],β的范围为[0-360°),步长,⊿sp;计算不同的α、β对应的平面内——即不同的渗流面ψ内,KYYαβθ、KZZαβθ,每个渗流面ψ内计算得到的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αβ最小值对应的渗流面ψ为等渗透率面,记录对应的倾角α,倾向β。所述的KXXαβθ为X方向渗透率主值,KYYαβθ为Y方向渗透率主值,KZZαβθ为Z方向渗透率主值。
本发明的有益效果是:本发明专利基于岩心统计、薄片观测以及测井解释,通过现今地应力分析,确定裂缝地下开度,利用古地磁定向、岩层产状法或地层倾角法确定裂缝的三维展布,进而建立真三维裂缝渗透率计算模型,推导可编程算法,编制相应的程序,实现准确的计算三维等渗透率面、最大差渗透率面、不同方向的渗透率剖面等。本发明由严格的数学算法推导组成,对相应的地质信息数字化后,可以利用计算机编程语言开发相应的计算程序,实现井筒裂缝真三维空间渗透率张量定量计算。本发明对于真三维空间裂缝渗透率张量及井网设计、确定注水井与采油井的空间位置关系、寻找最优钻井轨迹等油田开发方案的部署具有较高的实用价值,并且预测成本低廉、可操作性强,预测结果可靠。
附图说明
图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)推导的算法编制相应的程序,通过循环迭代变量α、β、θ,其中的α范围为[0-90°],β的范围为[0-360°),θ的范围为[0-90°),根据计算精度需求,设置循环迭代步长⊿sp=1°,将α、β、θ逐次累加⊿sp,利用公式(12)-(14),计算坐标系O-XYZ不同的α、β、θ对应的KXXαβθ、KYYαβθ、KZZαβθ。
利用公式(1)-(14),设置α=0°,β=0°;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ、KZZαβθ,得到水平面的渗透率变化(图6)。
利用公式(1)-(14),设置α=90°,β=0°;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ、KZZαβθ,得到东西向剖面的渗透率变化(图7)。
利用公式(1)-(14),设置α=90°,β=90°;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ、KZZαβθ,得到南北向剖面的渗透率变化(图8)。
利用公式(1)-(14),设置α范围为[0-90°],β的范围为[0-360°),计算不同的渗流面ψ内计算得到的KYYαβθ、KZZαβθ数目为180,180个KYYαβθ、KZZαβθ的平均值为Kaver;计算不同的α、β对应的方差函数Wαβ;通过公式(15),计算不同的α、β对应的方差函数Wαβ,比较不同的α、β的Wαβ。其中,Wαβ最大值对应的渗流面ψ为最大差渗透率面(图9);Wαβ最小值对应的渗流面ψ为等渗透率面(图10)。
上面以举例方式对本发明进行了说明,但本发明不限于上述具体实施例,凡基于本发明所做的任何改动或变型均属于本发明要求保护的范围。
Claims (1)
1.一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型的系统,步骤如下:
1)利用岩心观测、测井资料解释对裂缝的开度测量,确定裂缝的组系特征;依据成像资料、测井或者岩心观察,求取裂缝的线密度;
2)在岩心裂缝观测中,对非定向岩心采用古地磁定向、岩层产状法或地层倾角法确定每条裂缝的产状;在成像资料、声波测井解释裂缝时,利用地层倾角测井或者成像测井资料对裂缝定向;
3)利用水力压裂法、波速法、声发射法、应力解除法、光弹性应力法对工区现今地应力方位进行了判断,并通过测井、压裂资料以及物理实验计算了关键井现今地应力的数值,利用地震反演、测井解释、岩石力学实验确定岩石的力学属性;在此基础上,利用软件建立有限元模型,采用多井多层段实测地应力的大小、方向为约束条件,对研究区现今地应力场进行了数值模拟;利用模拟结果,确定井筒裂缝发育处的现今地应力大小、方向;利用公式(1)计算裂缝在地下的开度:
公式(1)中,b为裂缝的现今地下实际开度,单位:m;b0为裂缝的原始开度,单位:m;σ′n为有效正应力,单位:MPa;bres代表裂缝面承受最大有效正应力时的裂缝开度,单位:m;σnref为使裂缝开度减小90%的有效正应力,单位:MPa;
4)如图2所示,建立真三维裂缝渗透率张量计算几何模型,设裂缝面Φ的倾角为δ,倾向为ω,以裂缝为参照物建立O-ABC坐标系,O-ABC坐标系中的三个坐标轴分别对应于裂缝面的法线方向、裂缝走向方向、裂缝面内垂直于裂缝走向线的方向;设任意渗流面ψ的倾角为α,倾向为β,其单位法向矢量为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坐标系中的表达式表示为:
公式(7)中δ为裂缝面Φ的倾角,公式(8)中e、f、g的含义如公式(6)所示;
在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)利用公式(1)-(14)推导的算法编制相应的程序,通过循环迭代变量α、β、θ,所述的α为渗流面ψ的倾角,β为渗流面ψ的倾向,θ为OY轴与OP轴的夹角,其中的α范围为[0-90°],β的范围为[0-360°),θ的范围为[0-90°),根据计算精度需求,设置循环迭代步长⊿sp,将α、β、θ逐次累加⊿sp,利用公式(12)-(14),计算坐标系O-XYZ不同的α、β、θ对应的KXXαβθ、KYYαβθ、KZZαβθ;
利用公式(1)-(14),设置α=0°,β为[0-360°)的任意数值;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ、KZZαβθ,得到水平面的渗透率变化;
利用公式(1)-(14),设置α=90°,β=0°;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ、KZZαβθ,得到东西向剖面的渗透率变化;
利用公式(1)-(14),设置α=90°,β=90°;调节θ的大小,步长⊿sp,范围为[0-90°),输出KYYαβθ、KZZαβθ,得到南北向剖面的渗透率变化;
利用公式(1)-(14),设置α范围为[0-90°],β的范围为[0-360°),步长,⊿sp;计算不同的α、β对应的平面内——即不同的渗流面ψ内,KYYαβθ、KZZαβθ,每个渗流面ψ内计算得到的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;所述的KXXαβθ为X方向渗透率主值,KYYαβθ为Y方向渗透率主值,KZZαβθ为Z方向渗透率主值;
通过公式(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 CN106484958A (zh) | 2017-03-08 |
CN106484958B true 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) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107132586B (zh) * | 2017-04-19 | 2021-04-16 | 华北理工大学 | 一种断陷盆地断裂带结构单元定量识别方法 |
CN108061699B (zh) * | 2017-12-14 | 2020-01-21 | 中国石油大学(华东) | 一种多期次、多尺度裂缝孔渗参数定量预测方法 |
CN108008117B (zh) * | 2017-12-21 | 2020-04-03 | 中国石油大学(华东) | 裂缝发育模式与储层基质-裂缝渗透率敏感性预测方法 |
CN110020453B (zh) * | 2018-01-09 | 2023-03-10 | 姚君波 | 地层渗透率三维分解模拟器 |
CN109033538B (zh) * | 2018-06-30 | 2022-07-22 | 南京理工大学 | 一种基于实测结构面参数的裂隙岩体渗透张量的计算方法 |
CN109493397B (zh) * | 2018-09-30 | 2022-02-01 | 中国石油天然气股份有限公司 | 一种钻井数据处理方法、装置及系统 |
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 |
CN116956670B (zh) * | 2023-07-17 | 2024-01-23 | 长江大学 | 一种基于tpfa与mfd混合方法的投影嵌入式离散裂缝模型 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101942991A (zh) * | 2010-06-30 | 2011-01-12 | 中国石油大学(北京) | 裂缝各向异性油藏注水开发可预测物理模型建立方法 |
CN104200039A (zh) * | 2014-09-17 | 2014-12-10 | 中国石油大学(华东) | 一种构造裂缝产状定量预测方法 |
CN104500050A (zh) * | 2014-10-31 | 2015-04-08 | 中国石油大学(华东) | 一种裂缝性储层渗透率张量及各向异性定量预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8781806B2 (en) * | 2009-04-30 | 2014-07-15 | Schlumberger Technology Corporation | Determining elastic and fluid flow properties of a fractured reservoir |
-
2016
- 2016-09-19 CN CN201610832079.3A patent/CN106484958B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
---|
Estimation of equivalent fracture network permeability using fractal and statistical network properties;Alireza Jafari 等;《journal of Petroleum Science and Engineering》;20120621;第110-123页 |
岩体节理三维网络模拟技术及渗透率张量分析;王贵宾 等;《岩石力学与工程学报》;20041130;第23卷(第21期);第3591-3594页 |
Also Published As
Publication number | Publication date |
---|---|
CN106484958A (zh) | 2017-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106484958B (zh) | 一种基于井筒裂缝观测的三维裂缝渗透率张量计算模型的系统 | |
CN104977618B (zh) | 一种评价页岩气储层及寻找甜点区的方法 | |
US10101498B2 (en) | Well survivability in multidimensional geomechanical space | |
CN105672971B (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 | |
CN106951660A (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
CN104863574B (zh) | 一种适用于致密砂岩储层的流体识别方法 | |
CN105510993A (zh) | 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 | |
Liu et al. | Asymmetric propagation mechanism of hydraulic fracture networks in continental reservoirs | |
Lai et al. | Earth stress and reservoir quality evaluation in high and steep structure: The Lower Cretaceous in the Kuqa Depression, Tarim Basin, China | |
CN107748399A (zh) | 利用重力界面反演识别山前带深部构造层方法 | |
Feng et al. | 3D numerical simulation of heterogeneous in situ stress field in low-permeability reservoirs | |
Liu et al. | Quantitative multiparameter prediction of fractured tight sandstone reservoirs: a case study of the Yanchang Formation of the Ordos Basin, Central China | |
Zeng et al. | Unreliable determination of in situ stress orientation by borehole breakouts in fractured tight reservoirs: A case study of the upper Eocene Hetaoyuan Formation in the Anpeng field, Nanxiang Basin, China | |
Feng et al. | Near-surface stress measurements in the Longmenshan fault belt after the 2008 Wenchuan Ms8. 0 earthquake | |
CN111352172A (zh) | 一种用井震联合法获取铀异常在砂体中空间分布位置的方法 | |
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 | |
Zhao et al. | Present in-situ stress distribution characteristics of strike-slip in SH Oilfield, Tarim Basin | |
Louie et al. | Advanced seismic imaging for geothermal development | |
Feng et al. | Development characteristics and quantitative prediction of multiperiod fractures in superdeep thrust-fold belt | |
Guo et al. | A new method for identifying fractures in tight sandstone of a gentle structural area using well logs | |
Ou et al. | Three-dimensional discrete network modeling of structural fractures based on the geometric restoration of structure surface: Methodology and its application | |
Zeng et al. | Variation in the orientation of the maximum horizontal stress in thick channel-fill sandstones with low-permeability: A case of the Bonan Oilfield in the Bohai Bay Basin, eastern China | |
Du et al. | Integrated shale gas reservoir modeling | |
Liang et al. | Stepwise northward compression in the northeastern Tibetan Plateau: Insights from the chronology of the Baima Basin |
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 |
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) |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190917 Termination date: 20200919 |
|
CF01 | Termination of patent right due to non-payment of annual fee |