浅议数值模拟收敛性调整方法.docxVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
浅议数值模拟收敛性调整方法   【摘要】本文以七区西馆下为例,对数值模拟中的不收敛性问题进行了分析,从地质模型网格划分、属性参数、pvt表、相对渗透率选取、模型初始化等角度进行分析并提出解决途径,改善了模型的收敛性,计算时间大幅缩短,提高了数值模拟工作效率。   【关键词】数值模拟;收敛性;调整方法   七区馆下段由天然能量开发转换为注水开发。在方案编写时采取了数值模拟研究。地质模型采用direct建立,网格步长20*20*0.9米,网格数量93*215*52=103.97万个,其中活网格数量24.59万个,采用eclipse油水两相黑油模型进行模拟,由于网格数量大,开发历史长,数模计算量大,模型建立初期,进行一次模拟时间长达78小时,因此,需要从检查模型的收敛性入手,提高模拟速度,提高模拟效果。   1.模型不收敛的原因分析   通过对数值模拟的基本原理和算法进行分析,模型不收敛的原因很多,网格参数,属性参数,流体pvt参数,岩石相渗曲线,毛管压力曲线,相渗曲线端点标定,初始化,井轨迹,垂直管流表都会造成模型不收敛,具体分析如下。   1.1网格部分:网格正交性差和网格尺寸相差太大是导致不收敛的主要原因之一。正交性差会给矩阵求解带来困难,而网格尺寸相差大会导致孔隙体积相差很大,大孔隙体积流到小孔隙体积常会造成不收敛。   1.2属性参数:不合理的插值计算会导致属性分布很差,主要是地质属性计算插值时由于异常点,导致属性分布不合理造成收敛性变差。   1.3流体pvt参数:流体pvt参数会有两种可能的问题,一是数据不合理导致了负总压缩系数,二是压力或气油比范围给的不够导致模型对pvt参数进行了外插。   1.4岩石相渗曲线和毛管压力曲线:eclipse不会对输入模型的相渗曲线和毛管压力曲线进行光滑,将会应用每一个输入饱和度和相渗值,所以要保证输入的参数是合理的。通常的问题有饱和度和相对渗透率的数据位数过多、饱和度值太接近,导致相渗曲线的倾角变化很大、饱和度有很小变化但相对渗透率发生了很大变化。   1.5模型初始化不合理:初始化最容易发生的问题是在初始时模型不稳定,流体在初始条件下就会发生流动,这也会导致模型不收敛。造成模型初始不稳定的主要有赋网格饱和度和压力值和初始含水饱和度场不匹配。   2.解决办法   2.1加强网格质量检查,提高网格正交性和分布范围。网格正交性差通常是在建角点网格时为描述断层或裂缝的走向而造成的。在此情况下,通过调整边界,使边界与主断层走向平行,这样一方面网格可以很好地描述断层,另一方面正交性也很好。重新定义网格划分,在平面上最好让网格大小能够较均匀,在没有井的地方网格可以很大,但最好能够从大到小均匀过渡,对局部大网格区域实施加密,纵向上对层厚再细分。对每层都在三维显示中检查。调整后,网格体积变化范围由50-850m3调整为120-520m3。同时在eclipse里用minpv关键字可以把小于设定孔隙体积的网格设为死网格,参与计算的活网格数量由24.59万个调整为22.18万个。   2.2光滑属性参数,截断异常属性值。对原模型进行修改,在尽量用地质模型的数据的前提下,利用光滑工具对异常值进行修正,并加一些控制点使属性合理分布。调整后x,y方向的渗透率相等或级差不大,在井连通网格的z方向渗透率大于30×10-3μm2。   2.3流体pvt参数。检查prt文件中的warning信息后发现,在油藏压力范围内有负总压缩系数的警告,通通过修改pvt表,把压力应该覆盖所有范围,包括注水后的压力上升后的压力,提高收敛性。在eclipse中通过extrapms关键字控制pvt插值,进一步提高收敛性。   2.4岩石相渗曲线和毛管压力曲线。对相渗曲线进行标准化,检查相渗曲线的导数,导数要光滑。同时保证相对渗透率数值有两位小数。在三维显示中检查标定完的pcw,通过给pcw调整为0.06mpa,控制毛管压力。在计算中,输出每个网格标定后的相渗曲线进行检查。重新进行调整。   2.5修改初始化方法。受地质认识和井的投产时间影响,测井二次解释的含油饱和度不一定符合实际因此,数模采用的模型网格的压力和饱和度值由直接赋值改为通过油水界面及参考压力来进行初始化计算。并进行毛管压力的端点标定,这样毛管压力会稳住每个网格的水,在初始条件下不会流动。本模型共划分了四个平衡分区,建立了四个解析水体,通过数模软件本身,实现模型初始化。   2.6控制计算参数。主要是调整tuning中的最大时间步。七区西馆下数值模拟最大时间步长由默认的365天调整为10天,最大线形迭代次数到70次,计算加快很多。  

文档评论(0)

软件开发 + 关注
官方认证
服务提供商

十余年的软件行业耕耘,可承接各类需求

认证主体深圳鼎云文化有限公司
IP属地陕西
统一社会信用代码/组织机构代码
91440300MA5G24KH9F

1亿VIP精品文档

相关文档