原创力文档- 1、本文档共45页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
油藏数值模拟大作业
1 差分方程组线性化
二维空间油水两相流体流动方程组
二维空间油水两相流体流动的差分方程组
油相：
水相：
其中：
令：
上式两边同乘以，则油水两相的差分方程为：
油相：
水相：
令：
上两式可化为：
差分方程线性化：
方程右端项处理：
油相：
水相：
利用IMPES方法，压力取隐式形式，饱和度和毛管压力取显示形式。
令
进一步可化简为：
其中：
油方程/Co2＋水方程/Cw2，得：
（1）
式中：
如果网格块上有井，需加上产量项，有两种情况：
1.定产量，直接在左端加上，注入为正，生产为负;
2.定井底流压，处理如下：
分别将产水量和产油量加到对应的油、水方程中，即：
2 求解差分方程
首先，对网格进行编号，见表2.1。
表2.1 网格编号
5 10 15 1 6 11 16 21 45 52 2 7 12 17 22 26 31 38 46 53 57 59 3 8 13 18 23 27 32 39 47 54 58 60 4 9 14 19 24 28 33 40 48 55 20 25 29 34 41 49 56 30 35 42 36 43 50 37 44 51 其中，18号井是注水井，11号井和46号井为生产井。
将网格编号后，就可计算出公式（1）的各个系数，并将其赋值到60*60的矩阵方程中，从而可计算出各网格的压力变换量，最终求得不同生产时间的油水两相压力分布和各相饱和度分布，以及生产井的生产数据。（其中矩阵方程的系数形状见图1）
图2.1 系数矩阵形状
利用C语言编制程序实现以上步骤，详细程序见附录。
3 求解结果
（1）分别解出100d时的各网格的油相压力及各网格的含水饱和度，见表3.1和表3.2。并根据表中数据分别绘制出地层压力分布和含水饱和度分布图，见图3.1和图3.2。
表3.1 100d时油相压力分布
5737.54 5731.54 5739.20 5793.84 5738.57 5737.52 5842.10 5902.66 5873.96 5848.08 5746.58 5844.99 5942.44 5850.47 6056.92 6081.26 6012.14 5953.05 5826.94 5849.02 5954.82 5957.56 5849.45 5850.15 5951.32 5936.06 6059.81 6162.36 6153.01 5957.53 5950.12 5853.40 5955.85 5956.64 5851.81 5852.36 5954.85 5958.71 6161.44 6183.26 6173.96 6060.40 6056.15 5955.60 5961.68 6162.08 6189.77 6184.64 6062.23 6058.05 5956.47 6204.12 6205.32 6104.27 6236.43 6106.21 6167.08 6257.25 6107.27 6167.54
表3.2 100d时含水饱和度分布
0.4646 0.4624 0.4686 0.4672 0.4703 0.0000 0.4671 0.4739 0.4597 0.4654 0.4651 0.4663 0.4622 0.4649 0.4665 0.4683 0.4679 0.4645 0.0000 0.4671 0.4622 0.4617 0.4643 0.4622 0.4624 0.4656 0.4650 0.4656 0.4652 0.4633 0.4615 0.4634 0.4631 0.
文档评论（0）