收藏本站
《中国地质大学》 2009年
收藏 | 手机打开
二维码
手机客户端打开本文

低渗含裂缝砂岩油藏渗流规律及综合模拟研究

周锋德  
【摘要】: 宝浪油田地理位置处于新疆维吾尔自治区巴音郭楞蒙古自治州焉耆回族自治县境内。构造上位于焉耆盆地博湖坳陷北部凹陷宝浪苏木构造带。主要产油层为三工河组。沉积相以上辫状河三角洲前缘水下分流河道和下三角洲前缘水下分流河道为主。岩性包括砂砾岩、砾岩、粗砂岩、细砂和粉砂岩等。岩石成分成熟度和结构成熟度较低。 储层属辫状河三角洲上发育起来的粗粒低渗储层类型。单层砂体厚度从1m至6m不等,常见多层叠加。河道频繁迁移加剧了储层的非均质性。不同砂体内平均孔隙度从12%-14%,平均渗透率从2.972×10~(-3)μm~2-27.3×10~(-3)μm~2。储层属性和流体分布复杂,砂体厚度大(Ⅲ_1油组接近20m),同时,独特的高陡背斜构造及裂缝的发育,使该地区储层成为典型的非常规储层类型:低密度、低孔孔隙度、低渗透率、粗粒、高粘土矿物、裂缝和复杂构造。 开发13年以来,一些与地质有关的问题仍然影响该地区储层的有效开发。例如:砂体的三维分布;流动单元的分布;裂缝的特征及其对流体渗流的影响;含水率上升的控制因素,等。 本文主要利用地质统计学、随机模拟和数值模拟等方法,首次提出并开展了针对低渗含裂缝储层进行了综合表征和模拟的方法和实现;首次采用小尺度和大尺度数值模拟相结合的方法系统评价了不同裂缝及不同流动单元组合情况下的地质模型的开发特征及参数。 具体内容主要包括以下五个方面: 第一,流动单元综合研究:在地层划分、岩心分析、岩心测试、储层特征参数和岩石物理相的基础上,采用聚类方法划分了流动单元类型。通过多元回归方法建立了流动单元的预测模型。对比砂体的空间分布,分析了流动单元在平面和剖面上的分布规律。 第二,裂缝单元综合表征:通过岩心观察测定了岩心裂缝的密度、倾角、走向等参数;并通过测井和构造曲率方向分别建立了低角度裂缝和高角度裂缝的预测模型,预测并绘制了不同裂缝的裂缝密度平面分布图。 第三,对建立随机离散裂缝分布模型的方法和参数进行了研究,并对不同裂缝密度的随机离散分布模型进行了比较。 第四,通过对预测的低角度和高角度裂缝的平面分布图进行数字化,并对裂缝的计数密度进行了面积密度的转化,最后建立了研究区的低角度和高角度裂缝随机离散分布模型。并采用多界面(多级界面主要包括油组界面、小层界面和单砂体界面)及流动单元的约束方法建立了低渗基质的地质模型。 第五,对岩心规模及小尺度规模的不同流动单元及不同裂缝产状的模型进行了数值模拟和比较。并在储层地质模型的粗化基础上,对Ⅲ_1油组进行了数值模拟,在历史拟合的基础上,分析了剩余油分布和剩余可采储量的平面分布。 通过研究,取得的主要成果和认识如下: (1)通过FZI=4μm、FZI=2.72μm、FZI=1.85μm和FZI=1.36μm.可以把本地区流动单元类型可以分为5类,流动单元A、流动单元B、流动单元C、流动单元D和流动单元E。在流动单元基础上建立的储层孔渗测井解释模型预测的储层物性与岩心分析的物性的相关系数大于0.7,能满足精细储层建模的要求。 (2)研究区裂缝类型按成因主要分为两类:构造成因和构造-沉积复合成因。复合成因的层状裂缝最发育,断层有关的裂缝次之。将近60%裂缝属于复合成因,40%属于构造成因。裂缝按倾角主要分为两种类型,低角度裂缝和高角度裂缝。高角度裂缝走向以东北、南西和东南方向为主,平均为156°。其倾角分布范围为15°-85°之间,大多数在40°-75°之间分布,平均为54.2°。背斜的不同部位倾角不同,西北部大于东南部。低角度裂缝的倾角在0°-20°之间不等,常呈“薄饼状”,一些裂缝的宽度小于5cm。 (3)三孔隙度方法能够很好的预测该地区的低角度裂缝的密度,临界值为0.3。利用SPSS软件建立了低角度裂缝密度预测的多元回归方程,利用该方法预测的裂缝密度和岩性测量的裂缝之间的相关系数为0.716。在模型建立的基础上,编写了VB代码自动计算了除部分没有补偿中子、声波和密度测井的所有井,并分析了低角度裂缝的密度平面分布。高角度裂缝的预测模型通过分析裂缝发育指数和岩心裂缝密度的关系来建立。 (4)流动单元组合下的岩心尺度数值模拟结果是不同的流动单元储层单独开采下的采收率最高,其次是流动单元A和流动单元B组合,其他组合均低于这两种组合。流动单元组合下的小尺度模拟表明,在相对注入速度相等的情况下,井距对综合采出程度的影响较小;而注入速度对综合采出程度影响较大。 (5)四种不同倾角的裂缝模拟结果显示,在含水量达到98%时,流动单元A、流动单元B和流动单元C中采收率的关系是135°<45°<0°<90°<无裂缝模型。 (6)和没有裂缝的模型比较,具有裂缝的岩心当含水饱和度达到98%时,流动单元A和流动单元B的含水率减少,流动单元C的含水率增加。由于串通裂缝的层间流动影响,流动单元A的采收率增加,流动单元B和流动单元C的减少,总的采收率降低。 (7)离散裂缝等效模拟结果显示,裂缝的等效处理网格大小对基质渗透率较低的模型模拟结果影响较大。当垂向网格大小等于模型厚度时,模拟的含水率比离散裂缝模型的结果低20%。合理的垂向网格大小和基质渗透率成正比关系。 (8)多界面和流动单元共同控制储层地质建模方法能很好的控制砂体形态和砂体属性参数分布。流动单元A、B、C、D、E和泥岩的最大主变程分别为635m、429m、520m、558m、556m和783m。孔隙度在流动单元B、C、D和E中的最大主变程分别为483m、676m、516m和628m。孔隙度在流动单元B、C、D和E中的最小次变程分别为248m、200m、264m和232m。流动单元A、B、C、D和E的渗透率控制范围分别为52.5~162.6×10~(-3)μm~2、18.3~97×10~(-3)μm~2、6.6~51.6×10~(-3)μm~2、2.5~17.8×10~(-3)μm~2、1.4~16.8×10~(-3)μm~2。最后,建立了各种储层属性分布模型。 (9)运用随机模拟方法建立了不同裂缝参数条件下的裂缝分布模型,探讨了裂缝密度参数与裂缝模型结果参数的关系。当裂缝面积密度为0.05m~2/m~3、0.1m~2/m~3、0.5m~2/m~3、1m~2/m~3、2m~2/m~3和10m~2/m~3时,裂缝模型的平均孔隙度分别为0.0000、0.0002、0.0019、0.0028、0.0087和0.0461;平均渗透率分别为1×10~(-3)μm~2、10×10~(-3)μm~2、79×10~(-3)μm~2、119×10~(-3)μm~2、369×10~(-3)μm~2、1260×10~(-3)μm~2和1941×10~(-3)μm~2;裂缝形状系数平均值分别为3984m~(-2)、371m~(-2)、1476m~(-2)、3×10~6m~(-2)、8×10~5m~(-2)、和1×10~5m~(-2)。粗化的孔隙度和渗透率值与裂缝面密度成直线关系。 (10)利用数值模拟方法对(9)的模型在基质结合不同流动单元情况下的渗流特征进行了对比。发现不同模型的采收率和累计采油量与裂缝面积密度没有直线关系。采收率在基质为流动单元A、流动单元B和流动单元C时的拐点分别为2m~2/m~3、2m~2/m~3和1m~2/m~3;累计产油量在裂缝面密度大于0.5 m~2/m~3时递增,小于0.5 m~2/m~3时递减。注采井连线与裂缝走向平行时的采收率和累计产油量均小于相应的在注采井连线与裂缝走向垂直下的模拟结果。而含水饱和度达到98%的时间关系为:当基质为流动单元A时,垂直模型晚于平行模型,当基质为流动单元B和C时,垂直模型早于平行模型。 (11)运用随机模拟方法建立了Ⅲ_1油组的高角度裂缝和低角度裂缝离散分布模型。并利用Oda方法对该离散裂缝网络模型进行了粗化,得出了裂缝孔隙度、渗透率和裂缝形状指数。模型的平均孔隙度为0.4%,I、J和K方向的平均渗透率分别为11×10~(-3)μm~2、22×10~(-3)μm~2和13×10~(-3)μm~2。I、J和K方向的平均裂缝间距分别为2.78m、4.19m和0.58m。整个模型的裂缝形状指数为18836/m~2.在该模型的基础上,利用数值模拟方法对剩余油分布进行了研究。剩余油分布主要受注采井网不完善、次级断层、砂体分布和流动单元储层非均质性等因素影响控制和影响。
【学位授予单位】:中国地质大学
【学位级别】:博士
【学位授予年份】:2009
【分类号】:P618.13

手机知网App
【相似文献】
中国期刊全文数据库 前10条
1 吕游;;浅析公众聚集场所火灾扑救中烟气的流动与控制[J];科技信息;2011年22期
2 刘晨;郭祖培;;SMSJ—44A型水模拟浇注系统实验台研制[J];实验技术与管理;1992年03期
3 ;[J];;年期
4 ;[J];;年期
5 ;[J];;年期
6 ;[J];;年期
7 ;[J];;年期
8 ;[J];;年期
9 ;[J];;年期
10 ;[J];;年期
中国重要会议论文全文数据库 前10条
1 程浩;郎兆新;;泡沫在多孔介质中的流动规律及其数值模拟[A];自然、工业与流动——第六届全国流体力学学术会议论文集[C];2001年
2 周锋德;姚光庆;魏忠元;刘斌;黄郑;;低渗含裂缝储层流动单元控制的地质建模[A];中国西部复杂油气藏地质与勘探技术研讨会论文集[C];2006年
3 霍云霈;杨新军;;太白山国内游客空间流动结构分析[A];中国地理学会2006年学术年会论文摘要集[C];2006年
4 李兆敏;董贤勇;杨宪;王渊;王德新;;泡沫流体冲砂洗井技术研究[A];科技、工程与经济社会协调发展——中国科协第五届青年学术年会论文集[C];2004年
5 杨正明;黄延章;尚根华;郝明强;;低渗透多孔介质中非线性渗流的流态与准数[A];中国力学学会学术大会'2005论文摘要集(下)[C];2005年
6 刘曰武;陈维良;;低渗及特低渗油藏井间渗流规律研究[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集(下)[C];2007年
7 李建国;;油底壳成形过程的三维有限元模拟[A];中国工程物理研究院科技年报(2002)[C];2002年
8 谷芳;刘春江;袁希钢;余国琮;;倾斜板上三维降膜流动的CFD模拟[A];第一届全国化学工程与生物化工年会论文摘要集(上)[C];2004年
9 刘春泽;夏惠芬;程林松;;黏弹性聚合物在不规则孔道中的驱油机理[A];中国力学学会学术大会'2005论文摘要集(下)[C];2005年
10 甘淑;何大明;冯彦;;纵向岭谷区流域景观变化与现代地表过程的关联分析[A];中国地理学会2006年学术年会论文摘要集[C];2006年
中国重要报纸全文数据库 前10条
1 包心鉴;流动人口研究的一个重要成果[N];中国人口报;2007年
2 通讯员  张卉桐 记者  张欣;吉林年内建成64座农村客运站[N];中国交通报;2006年
3 王迅;群采多空区下低贫损联合采矿技术通过鉴定[N];中国有色金属报;2006年
4 秦洪;资金回流二线蓝筹股[N];上海证券报;2007年
5 本报记者 梁懿娴;何时扭转“穷国向富国送钱”怪现状[N];国际商报;2008年
6 冯作文;开采人力资源的富矿[N];中国石油报;2004年
7 李曦辉;交通要素及其他[N];中华工商时报;2002年
8 袁之应无为县政府县长;多措并举促进农民工返乡创业[N];安徽日报;2007年
9 顾晶;三成居民有意五年内买车[N];南昌日报;2008年
10 记者 金衡;全国道路交通安全形势明显好转[N];人民公安报·交通安全周刊;2007年
中国博士学位论文全文数据库 前10条
1 周锋德;低渗含裂缝砂岩油藏渗流规律及综合模拟研究[D];中国地质大学;2009年
2 李峰;盘类件模锻过程金属变形模式及流动规律研究[D];哈尔滨工业大学;2007年
3 林日亿;蒸汽—氮气泡沫体系流动规律研究及应用[D];中国石油大学;2007年
4 王殿生;缝洞型介质流动机理实验与数值模拟研究[D];中国石油大学;2009年
5 黄善波;聚驱井采出液流动规律研究与应用[D];中国石油大学;2007年
6 穆文志;低渗透油层孔隙介质内微观流动机理研究[D];大庆石油学院;2009年
7 吴海进;高瓦斯低透气性煤层卸压增透理论与技术研究[D];中国矿业大学;2009年
8 徐依吉;超高压水射流理论与应用基础研究[D];西南石油学院;2004年
9 高昂;循环经济物质流特征与流动规律研究[D];西北大学;2010年
10 肖洋;横向流动条件下多孔水平动量射流掺混特性研究[D];河海大学;2005年
中国硕士学位论文全文数据库 前10条
1 刘颖;女性性服务者流动性的人类学研究[D];中央民族大学;2007年
2 李雪梅;煤炭资源流动规律及其替代机制研究[D];山东科技大学;2008年
3 徐永辉;泡沫酸特性研究及其应用[D];中国石油大学;2007年
4 丁娟;上海市徐家汇地铁商城群集流动规律研究及其疏散评估[D];同济大学;2007年
5 孙栋梁;云锡塘子凹松散岩体注浆加固试验研究[D];昆明理工大学;2006年
6 肖传发;板框式电化学反应器流动特性研究[D];郑州大学;2006年
7 姜海峰;非牛顿液体—气体多相流混输流动规律,水力计算方法及软件编制[D];大庆石油学院;2003年
8 王长军;基于动态最佳经济出矿品位的控制放矿技术研究[D];中南大学;2009年
9 丁顺利;大空间建筑火灾中烟气流动规律的研究[D];华北水利水电学院;2007年
10 李志;料仓中湿颗粒流动规律的数值仿真与试验研究[D];西安建筑科技大学;2003年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026