CN106291704B - 一种不同尺度裂缝面密度预测方法 - Google Patents
一种不同尺度裂缝面密度预测方法 Download PDFInfo
- Publication number
- CN106291704B CN106291704B CN201610833695.0A CN201610833695A CN106291704B CN 106291704 B CN106291704 B CN 106291704B CN 201610833695 A CN201610833695 A CN 201610833695A CN 106291704 B CN106291704 B CN 106291704B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- tomography
- fracture
- surface density
- 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
- 238000000034 method Methods 0.000 title claims abstract description 10
- 238000003325 tomography Methods 0.000 claims abstract description 42
- 238000011160 research Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000004141 dimensional analysis Methods 0.000 claims description 3
- 238000004379 similarity theory Methods 0.000 claims description 3
- 238000002591 computed tomography Methods 0.000 claims description 2
- 238000013467 fragmentation Methods 0.000 claims description 2
- 238000006062 fragmentation reaction Methods 0.000 claims description 2
- 238000011002 quantification Methods 0.000 claims description 2
- 238000012795 verification Methods 0.000 claims description 2
- 238000012546 transfer Methods 0.000 claims 1
- 238000011156 evaluation Methods 0.000 abstract description 5
- 238000004422 calculation algorithm Methods 0.000 abstract description 4
- 206010017076 Fracture Diseases 0.000 description 43
- 239000011435 rock Substances 0.000 description 5
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 235000003283 Pachira macrocarpa Nutrition 0.000 description 2
- 241001083492 Trapa Species 0.000 description 2
- 235000014364 Trapa natans Nutrition 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 235000009165 saligot Nutrition 0.000 description 2
- 102000016938 Catalase Human genes 0.000 description 1
- 108010053835 Catalase Proteins 0.000 description 1
- 206010019114 Hand fracture Diseases 0.000 description 1
- 241000124033 Salix Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/32—Transforming one recording into another or one representation into another
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及油气田勘探开发领域,尤其是一种不同尺度裂缝面密度预测方法。本发明在断裂带数据化后,通过分析断层、裂缝的相似性,选择合适的分形以及密度计算单元的边长,推导设计相应的算法,实现不同尺度的裂缝面密度表征。本发明对于不同尺度裂缝面密度预测等方面具有较高的实用价值,并且预测成本低廉、可操作性强,评价结果对预测油气勘探重点区域、裂缝优势分布区域以及后期裂缝孔渗参数的表征等多个方面有一定的指导意义。
Description
技术领域
本发明涉及油气田勘探开发领域,尤其是一种不同尺度裂缝面密度预测方法。
背景技术
在储层裂缝研究中,研究人员已经逐渐认识到不同尺度的裂缝在油气渗流中的作用不同,如何实现不同尺度的裂缝面密度预测已经成为制约裂缝预测的瓶颈。目前,裂缝的面密度预测主要集中在单一尺度的裂缝面密度预测,将不同尺度的裂缝面密度混淆一谈,往往影响裂缝参数的准确评价。
在断裂分析过程中,人们逐渐认识到断层、裂缝是同一应力场下作用下不同尺度的产物,在几何学、运动学以及动力学特征上具有统计意义的自相似性。岩石破碎过程具有自相似性,这种自相似性可以用分维定量表征,根据计算维数的定义,提出了不同尺度的裂缝面密度计算模型,开发相应的软件,实现了不同尺度的裂缝面密度定量预测。
发明内容
本发明旨在解决上述问题,提供了一种不同尺度裂缝面密度预测方法,它解决了裂缝的面密度预测主要集中在单一尺度的裂缝面密度预测,将不同尺度的裂缝面密度混淆一谈的问题。
本发明的技术方案为:一种不同尺度裂缝面密度预测方法,具体步骤如下:
第一步 通过地震解释、油田动态开发以及小层对比,获取断裂平面展布图
利用相关地震解释软件,通过精心地震解释获得工区的断裂平面展布图,通过油田动态开发以及小层对比确定工区次级断层的展布,形成研究区的断裂平面展布图。
第二步 工区断层-构造裂缝相似性分析
通过统计断层、裂缝的空间展布,裂缝、断层的规模、尺度分析,确定工区断层、裂缝的相似性,初步分析研究区断裂系统的相似性,论证基于相似性理论预测不同尺度的裂缝面密度的可行性。
第三步 断裂数据化以及点充填
在识别每条裂缝、断层的基础上,如图2所示,将每条断层沿着走向线依次取点,断层取点依次标记为(ai0、ai1、ai2...ain-1、ain),其中,i表示第i条断层或裂缝,n为该条断层的取点数目。设置断层的充填间距b,充填后得到的断层数据体为ΩD(X,Y);其中X,Y为充填点的坐标。
第四步 分形统计单元边长r以及裂缝面密度计算单元边长L确定
在断层分形评价中,断层在一定的尺度范围内存在自相似性,因此合理的设定分形统计单元边长r是断裂自相似性评价的关键,但诸多学者并未给出明确的网格单元边长的选取原则,这也影响断裂分形理论的推广影响;通过不断地变换分形统计单元的边长r,并分别用边长为r/3、r/4、r/5、r/6、r/7、r/8(不同ε)的栅格覆盖,得到工区断层体系单元边长r与信息维D、相关系数R2的分布关系。统计单元边长r越大,断裂信息维D总体呈增大趋势,相关系数R2越大,但不同统计单元间的差异也越小;因此选取信息维D、相关系数R2在趋于稳定后的最小值,在此边长下,一方面断裂具有较高的自相似性,另一方面,单元间同样有较大的差异性。
断层、裂缝的密度等参数有尺度效应,用不同的断层、裂缝密度计算单元边长L,得到的工区的裂缝参数有很大的差异,本发明专利通过不断地变换密度计算单元边长L,并求取对应的平均面密度ρLaver,拟合公式得到乘幂函数:
公式(1)中,ρLaver为不同单元边长L对应的工区平均面密度,单位:km/km2;L为不同的计算单元边长,单位:km;a、b为拟合系数,无量纲。
计算整个工区断裂的平均面密度ρaver:
公式(2)中,Lsum为工区内断层的总长度;S为工区面积。
公式(3)中Lopt为ρLaver=ρaver时,对应的边长,单位:km。合适的断裂面密度计算单元边长L介于1.8×Lopt~2.2×Lopt,从而确定断层、裂缝面密度计算单元边长L。
第五步 分形统计单元断裂参数C、D确定
在断裂构造的研究过程中,人们已逐渐认识到断裂系统在几何形态、构造演化以及成因动力上具有统计意义的自相似性;断裂分维D是断裂相似性的综合量化指标,其计算原理为:
公式(4)中,Pi是每个信息点落入第i个小栅格的概率;ε为栅格的边长,m。
I(r)=-Dln(r)+C (5)
公式(5)中,r为分形统计单元边长,m;D为断裂信息维,无量纲;C为断裂相似性拟合系数。
断裂插值充填后,在特定统计单元,统计落入其中的总点数Tsum,将统计单元划分为N(ε)个边长ε的栅格,计落入第i个栅格的点数为Ti,进而求取信息点落入不同栅格的概率Pi,I(ε)表示为:
通过不断地变换栅格的边长ε,得到不同ε对应的变量I(ε);对变量ln(ε)、I(ε)线性拟合后,得到断裂信息维以及对应的相关系数R2,通过编写程序,移动统计单元,计算得到统计单元内的信息维值D。在分形统计区间内,利用公式(5),通过拟合得到裂缝在分形统计区间内的信息维D以及截距C。
第六步 不同尺度的裂缝面密度计算模型
如图3所示,对于特定的面密度计算单元Ek,将其分为10×10网格,计算断层充填点在不同的小栅格内的概率δij:
公式(7)中,nij为第i排、第j列的小栅格内充填点的数目;nsum为计算单元Ek内充填点的总数目;
落入该计算单元Ek尺度为d的裂缝点数Nk表示为:
Nk=e-Dln(d)+C (8)
单元Ek内,不同的栅格内尺度为d的裂缝面密度ρs表示为:
公式(9)中,d为裂缝的尺度,m。
第七步 循环迭代,实现不同单元的不同尺度裂缝面密度预测,并进行模拟结果的对比验证。
利用上述算法,编制相应的程序,实现不同单元内栅格面密度的计算,依据岩心观测、薄片观测或者测井解释、地震反演,实现模拟结果验证。
本发明的有益效果是:本发明在断裂带数据化后,在确定断层、裂缝相似性的基础上,选择合适的分形以及密度计算单元的边长,推导设计相应的算法,实现不同尺度的裂缝面密度表征。本发明对于不同尺度裂缝面密度预测等方面具有较高的实用价值,并且预测成本低廉、可操作性强,评价结果对预测油气勘探重点区域、裂缝优势分布区域以及后期裂缝孔渗参数的表征等多个方面有一定的指导意义。
附图说明
图1为一种不同尺度裂缝面密度预测方法的流程图。
图2为断裂带数据化取点次序示意图。
图3为某一单元内不同栅格的断层充填点概率分布图。
图4金湖凹陷阜二段构造简图。
图5金湖凹陷阜二段断裂分布图。
图6金湖凹陷阜二段断层、裂缝走向展布相似性。
图7金湖凹陷阜二段断层长度与裂缝规模相似性。
图8金湖凹陷阜二段断层分形统计边长确定示意图。
图9金湖凹陷阜二段断裂密度统计边长确定示意图。
图10金湖凹陷阜二段断裂信息维D分布图。
图11金湖凹陷阜二段薄片尺度裂缝面密度分布图。
图12金湖凹陷阜二段岩心尺度裂缝面密度分布图。
具体实施方式
下面结合附图说明本发明的具体实施方式:
金湖凹陷在构造上位于苏北盆地东台坳陷的西部,是晚白垩世发育起来的断陷盆地。北起建湖隆起,南至天长凸起,西邻张八岭隆起,东为菱塘桥和柳堡两个低凸起,地理上跨越江苏、安徽两省,面积约5000km2,是苏北盆地中最大的一个沉积凹陷(图4)。在北东向建湖隆起、天长凸起的控制下,受北东向柳堡低凸起、菱塘桥低凸起及边界杨村断裂的制约,阜宁组沉积时期形成南断北超、南陡北缓的箕状凹陷,并呈北东向展布。
金湖凹陷结构为单断断超式即简单的半地堑,总体表现为东断西超,边界断层为杨村断层。吴堡事件后,北东向断裂发育,单一箕状凹陷在三垛期分为南北两个箕状断陷,北为北东向三河次凹为中心,东为石港断裂带,西为刘庄西北斜坡的东陡西缓的北东向展布的凹陷结构,南为以龙岗次凹为中心,东、西、北三面为斜坡的南陡北缓的箕状凹陷结构,凹陷内自东向西主要发育西部斜坡带、石港断裂带、汉涧斜坡带、卞闵杨构造带、唐湾构造带、宝应斜坡带等正向构造单元和三河次凹、汉涧次凹、龙岗次凹、汜水次凹等四个负向构造单元。金湖凹陷阜二段多尺度裂缝表征步骤为:
第1步 利用相关地震解释软件,通过精心地震解释获得工区的断裂平面展布图,通过油田动态开发以及小层对比确定工区次级断层的展布,形成金湖凹陷阜二段断裂平面展布图(图5)。
第2步 通过对金湖凹陷阜二段断层、裂缝统计,裂缝、断层的规模、尺度分析,确定工区断层、裂缝的相似性(图6、图7),初步判定基于相似性理论预测不同尺度的裂缝面密度是可行性。
第3步 如图2所示,将每条断层沿着走向线依次取点,断层取点依次标记为(ai0、ai1、ai2...ain-1、ain),其中,i表示第i条断层或裂缝,n为该条断层的取点数目。对断裂依次取点,完成对该区断裂数字化,统计各类断裂826条,数据化后的点数为22663个。依据区域小断裂的宽度,设定充填步长为10m,对断裂覆盖,在断裂带内充填的小格子数目为6756359个,充填后得到的断层数据体为ΩD(X,Y)。
第4步 利用公式(1)-公式(3),确定断层分形统计边长为10km(图8),确定断裂面密度计算单元边长为7km(图9),分形统计边长为10km。
第5步 利用公式(4)-(6),计算统计单元断裂分形统计参数C、D(图10)。
第6步 利用公式(7)-(9),输入不同的d值,得到薄片尺度的裂缝面密度分布(图11)、岩心尺度的裂缝面密度分布(图12)。
第7步 如表1所示,利用上述算法,编制相应的程序,实现不同单元内栅格面密度的计算,岩心裂缝线密度、孔隙度与预测结果对比表明,除在卞东地区预测出现较大的误差外,铜城、石港、桥河口、杨家坝以及闵桥地区吻合性较好;总体来看,模拟结果很好的体现了不同地区裂缝参数的差异性以及变化规律。
表1金湖凹陷裂缝数值模拟与岩心实测结果对比表
上面以举例方式对本发明进行了说明,但本发明不限于上述具体实施例,凡基于本发明所做的任何改动或变型均属于本发明要求保护的范围。
Claims (1)
1.一种不同尺度裂缝面密度预测方法,预测的步骤如下:
1)通过地震解释、油田动态开发以及小层对比,获取断裂平面展布图;
2)通过统计断层、裂缝的空间展布,裂缝、断层的规模、尺度分析,确定工区断层、裂缝的相似性,初步分析研究区断裂系统的相似性,论证基于相似性理论预测不同尺度的裂缝面密度的可行性;
3)在识别每条裂缝、断层的基础上,将每条断层沿着走向线依次取点,断层取点依次标记为ai0、ai1、ai2...ain-1、ain,其中,i表示第i条断层,n为该条断层的取点数目;设置断层的充填间距b,充填后得到的断层数据体为ΩD(X,Y);其中X,Y为充填点的坐标;
4)确定分形统计单元边长r以及裂缝面密度计算单元边长L,通过不断地变换分形统计单元的边长r,并分别用边长为r/3、r/4、r/5、r/6、r/7、r/8的栅格覆盖,得到工区断层体系单元边长r与信息维D、相关系数R2的分布关系,选取信息维D、相关系数R2在趋于稳定后的最小值作为断裂分形统计单元边长r;用不同的断层、裂缝密度计算单元边长L,得到的工区的裂缝参数有很大的差异,通过不断地变换断层、裂缝密度计算单元边长L,并求取对应的平均面密度ρLaver,拟合公式得到乘幂函数:
<mrow>
<msub>
<mi>&rho;</mi>
<msub>
<mi>L</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
</msub>
<mo>=</mo>
<mi>a</mi>
<mo>&CenterDot;</mo>
<msup>
<mi>L</mi>
<mi>b</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
上式中,ρLaver为不同单元边长L对应的工区平均面密度,单位:km/km2;L为不同的计算单元边长,单位:km;a、b为拟合系数,无量纲;同时计算整个工区断裂的平均面密度ρaver:
<mrow>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mi>u</mi>
<mi>m</mi>
</mrow>
</msub>
<mi>S</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
上式中,Lsum为工区内断层的总长度;S为工区面积;
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>o</mi>
<mi>p</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mroot>
<mrow>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>/</mo>
<mi>a</mi>
</mrow>
<mi>b</mi>
</mroot>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
上式中,Lopt为ρLaver=ρaver时,对应的边长,单位:km;合适的断层、裂缝面密度计算单元边长L介于1.8×Lopt~2.2×Lopt,从而确定断层、裂缝面密度计算单元边长L;
5)确定分形统计单元断裂参数C、D,断裂信息维D是断裂相似性的综合量化指标,其计算原理为:
<mrow>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mi>F</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<munder>
<mi>lim</mi>
<mrow>
<mi>&epsiv;</mi>
<mo>&RightArrow;</mo>
<mn>0</mn>
</mrow>
</munder>
<mfrac>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>&epsiv;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>&epsiv;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
上式中,Pi是每个信息点落入第i个小栅格的概率;ε为栅格的边长,m;
I(r)=-Dln(r)+C (5)
上式中,r为分形统计单元边长,m;D为断裂信息维,无量纲;C为断裂相似性拟合系数;统计落入其中的总点数Tsum,将统计单元划分为N(ε)个边长ε的栅格,计落入第i个栅格的点数为Ti,进而求取信息点落入不同栅格的概率Pi,I(ε)表示为:
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>&epsiv;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>&epsiv;</mi>
<mo>)</mo>
</mrow>
</mrow>
</munderover>
<mfrac>
<msub>
<mi>T</mi>
<mi>i</mi>
</msub>
<msub>
<mi>T</mi>
<mrow>
<mi>s</mi>
<mi>u</mi>
<mi>m</mi>
</mrow>
</msub>
</mfrac>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>T</mi>
<mrow>
<mi>s</mi>
<mi>u</mi>
<mi>m</mi>
</mrow>
</msub>
<msub>
<mi>T</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
通过不断地变换栅格的边长ε,得到不同ε对应的变量I(ε);对变量ln(ε)、I(ε)线性拟合后,得到断裂信息维以及对应的相关系数R2,计算得到统计单元内的信息维值D;在分形统计区间内,利用公式(5),通过拟合得到裂缝在分形统计区间内的信息维D以及截距C;
6)对于特定的面密度计算单元Ek,为了提高计算精度,将其分为10×10网格,计算断层充填点在不同的小栅格内的概率δij:
<mrow>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>n</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>n</mi>
<mrow>
<mi>s</mi>
<mi>u</mi>
<mi>m</mi>
</mrow>
</msub>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
公式(7)中,nij为第i排、第j列的小栅格内充填点的数目;nsum为计算单元Ek内充填点的总数目;
落入该计算单元Ek内尺度为d的裂缝点数Nk表示为:
Nk=e-Dln(d)+C (8)
单元Ek内,不同的栅格内尺度为d的裂缝面密度ρs表示为:
<mrow>
<msub>
<mi>&rho;</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>N</mi>
<mi>k</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
公式(9)中,d为裂缝的尺度,m,利用公式(7)-(9)计算不同尺度的裂缝面密度;
7)通过不同单元的循环迭代,实现不同单元的不同尺度裂缝面密度预测,并进行模拟结果的对比验证。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610833695.0A CN106291704B (zh) | 2016-09-20 | 2016-09-20 | 一种不同尺度裂缝面密度预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610833695.0A CN106291704B (zh) | 2016-09-20 | 2016-09-20 | 一种不同尺度裂缝面密度预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106291704A CN106291704A (zh) | 2017-01-04 |
CN106291704B true CN106291704B (zh) | 2018-02-27 |
Family
ID=57712232
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610833695.0A Expired - Fee Related CN106291704B (zh) | 2016-09-20 | 2016-09-20 | 一种不同尺度裂缝面密度预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106291704B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107784150B (zh) * | 2017-01-19 | 2020-11-24 | 中国石油大学(华东) | 一种裂缝封闭性预测方法 |
CN107065033B (zh) * | 2017-06-07 | 2019-06-25 | 中国石油大学(北京) | 一种致密砂岩储层微裂缝发育的识别方法及装置 |
CN108061699B (zh) * | 2017-12-14 | 2020-01-21 | 中国石油大学(华东) | 一种多期次、多尺度裂缝孔渗参数定量预测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102033242A (zh) * | 2010-10-22 | 2011-04-27 | 中国石油化工股份有限公司 | 一种深层倾斜裂缝储层地震振幅预测方法 |
CN104112057A (zh) * | 2013-04-19 | 2014-10-22 | 中国石油化工股份有限公司 | 一种大尺度裂缝性油藏数值模拟方法 |
CN104992067A (zh) * | 2015-07-14 | 2015-10-21 | 北京博达瑞恒科技有限公司 | 一种基于裂缝产状的裂缝分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2803068C (en) * | 2010-07-29 | 2016-10-11 | Exxonmobil Upstream Research Company | Method and system for reservoir modeling |
-
2016
- 2016-09-20 CN CN201610833695.0A patent/CN106291704B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102033242A (zh) * | 2010-10-22 | 2011-04-27 | 中国石油化工股份有限公司 | 一种深层倾斜裂缝储层地震振幅预测方法 |
CN104112057A (zh) * | 2013-04-19 | 2014-10-22 | 中国石油化工股份有限公司 | 一种大尺度裂缝性油藏数值模拟方法 |
CN104992067A (zh) * | 2015-07-14 | 2015-10-21 | 北京博达瑞恒科技有限公司 | 一种基于裂缝产状的裂缝分析方法 |
Non-Patent Citations (2)
Title |
---|
一种计算断裂信息维的方法及其应用;刘敬寿 等;《世界地质》;20150930;第34卷(第3期);第684-689页 * |
岩层厚度对碳酸盐岩构造裂缝面密度和分形分布的影响;孟庆峰 等;《高校地质学报》;20110930;第17卷(第3期);第462-468页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106291704A (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104500050B (zh) | 一种裂缝性储层渗透率张量及各向异性定量预测方法 | |
CN104459776B (zh) | 一种断裂分形特征优化计算方法 | |
CN106291704B (zh) | 一种不同尺度裂缝面密度预测方法 | |
CN108061699B (zh) | 一种多期次、多尺度裂缝孔渗参数定量预测方法 | |
CN104747180A (zh) | 一种用于注水开发的缝洞型油藏分析方法及其应用 | |
CN104866653A (zh) | 一种获取地下三维密度结构的方法 | |
CN103205972B (zh) | 一种分析基坑变形与坑外地面沉降关系的方法 | |
Yang et al. | Regional landslide zonation based on entropy method in Three Gorges area, China | |
CN104899928B (zh) | 基于稀疏钻孔点的三维地质建模方法 | |
CN106920176B (zh) | 一种矿集区尺度矿产资源估算方法及系统 | |
CN105321204B (zh) | 一种三维地质栅格体模型构建方法 | |
CN110097637A (zh) | 一种三维地质属性模型时空插值方法及系统 | |
Li et al. | Modelling the response of shallow groundwater levels to combined climate and water-diversion scenarios in Beijing-Tianjin-Hebei Plain, China | |
Xia et al. | Slope stability analysis based on group decision theory and fuzzy comprehensive evaluation | |
CN103065051A (zh) | 一种对岩体自动进行分级分区的方法 | |
CN107230252A (zh) | 基于曲流环形态数据库的点坝储层自动建模方法 | |
CN111006987A (zh) | 一种储层裂缝孔渗参数预测方法 | |
CN103778298A (zh) | 改进的模拟多孔介质中二维水流运动的多尺度有限元方法 | |
Wu et al. | An approach for deformation modulus mechanism of super-high arch dams | |
CN103605900A (zh) | 跨尺度复杂地质体地应力场识别方法及装置 | |
CN106980716A (zh) | 一种基于随机块体的地下洞室稳定分析方法 | |
CN106443782A (zh) | 一种断层、裂缝发育密度、均匀性以及组合样式评价方法 | |
CN110705168A (zh) | 构造应力场的模拟方法 | |
Yan et al. | Field study on deformation and stress characteristics of large open caisson during excavation in deep marine soft clay | |
CN104866682B (zh) | 一种基于地面形变进行页岩气勘探区构造应力场反演的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
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: 20180227 Termination date: 20200920 |