油藏?cái)?shù)值模擬在河流相整裝油藏中的應(yīng)用

蘭麗麗

摘 要：【目的】利用油藏?cái)?shù)值模擬技術(shù)計(jì)算并模擬油藏開發(fā)歷程中的動(dòng)態(tài)變化，通過(guò)網(wǎng)格內(nèi)的流體流動(dòng)狀態(tài)及數(shù)值的計(jì)算來(lái)預(yù)測(cè)不同預(yù)案下的開發(fā)狀況，以預(yù)測(cè)在不同開采方案下的油氣生產(chǎn)情況。【方法】針對(duì)河流相整裝油藏類型，采取有針對(duì)性的參數(shù)調(diào)整方法并建立動(dòng)態(tài)模型，以精細(xì)油藏描述成果及三維油藏模型為基礎(chǔ)，以油田投產(chǎn)以來(lái)區(qū)塊全部油水井動(dòng)態(tài)資料為依據(jù)進(jìn)行研究?！窘Y(jié)果】對(duì)東營(yíng)坳陷東部區(qū)塊開展數(shù)值模擬，并對(duì)歷史動(dòng)態(tài)進(jìn)行精準(zhǔn)擬合，整體含水?dāng)M合率從85.3%提升至92.7%；總結(jié)了本區(qū)塊剩余油分布規(guī)律，提出4個(gè)有利開發(fā)區(qū)，增加可采儲(chǔ)量12.3萬(wàn)t?！窘Y(jié)論】油藏?cái)?shù)值模擬方案通過(guò)對(duì)區(qū)塊整體和單井的擬合，以及數(shù)據(jù)精準(zhǔn)描述地層開發(fā)過(guò)程及儲(chǔ)層動(dòng)態(tài)變化，找到剩余油的分布位置，為制定合理的開發(fā)調(diào)整方案提供可信的科學(xué)依據(jù)。

關(guān)鍵詞：河流相；整裝油藏；數(shù)值模擬；歷史擬合；剩余油

中圖分類號(hào)：TE347? ? 文獻(xiàn)標(biāo)志碼：A? ? ?文章編號(hào)：1003-5168（2024）03-0042-06

DOI：10.19968/j.cnki.hnkj.1003-5168.2024.03.009

Application of Reservoir Numerical Simulation in Fluvial Facies

Integrated Reservoirs

LAN Lili

（Shengli Oilfield Geophysical Exploration Research Institute， Dongying 257000，China）

Abstract： [Purposes] This paper aims to use reservoir numerical simulation technology to calculate and simulate the dynamic changes during the development process of the reservoir， predict the development status under different plans through the calculation of fluid flow status and numerical values within the grid， and predict the oil and gas production situation under different development plans. [Methods] Targeted parameter adjustment methods and dynamic model establishment are adopted for the type of fluvial facies integrated oil reservoirs. Based on the results of fine reservoir description and three-dimensional reservoir model， the dynamic data of all oil and water wells in the block since the oilfield was put into operation are used as the basis. [Findings] This article conducted numerical simulation on the eastern block of Dongying Depression， accurately fitting historical dynamics， and increasing the overall water content fitting rate from 85.3% to 92.7%. It also summarized the distribution pattern of remaining oil in this block， proposed four favorable development zones， and increased recoverable reserves by 123 000 tons. [Conclusions] Through the fitting of the whole block and single well， as well as the accurate description of the formation development process and the dynamic change of the reservoir， the numerical simulation scheme of the reservoir can find the distribution position of the remaining oil and provide a credible scientific basis for formulating a reasonable development adjustment scheme.

Keywords： fluvial facies; integrated oil reservoir; numerical simulation; historical fitting; remaining oil

0 引言

油藏?cái)?shù)值模擬方法是一種具象化的模擬計(jì)算工具，通過(guò)與三維油藏模型相結(jié)合，以滲流力學(xué)、流體工程力學(xué)等計(jì)算方法為基礎(chǔ)，研究地層中油、氣、水的滲流規(guī)律、驅(qū)油機(jī)理及剩余油分布，從而模擬出油藏整個(gè)開發(fā)階段的動(dòng)態(tài)變化過(guò)程。油藏?cái)?shù)值模擬方法的使用開始于20世紀(jì)30年代，基礎(chǔ)理論是基于達(dá)西滲流定律提出的。當(dāng)時(shí)的模型主要用于預(yù)測(cè)油藏的動(dòng)態(tài)變化、采收率和選擇合適的開采方式，直至70年代初才出現(xiàn)第一批商業(yè)化軟件。但由于數(shù)值模擬技術(shù)不完善、計(jì)算機(jī)費(fèi)用昂貴等缺點(diǎn)，油藏模擬僅用于少數(shù)專家的研究工作。隨著技術(shù)的逐步發(fā)展和完善，20世紀(jì)80年代初，油藏?cái)?shù)值模擬開始向大型多功能綜合性軟件系列發(fā)展。目前，隨著計(jì)算機(jī)應(yīng)用數(shù)學(xué)和油藏工程的發(fā)展，油藏?cái)?shù)值模擬方法在石油工業(yè)中得到迅速發(fā)展和廣泛應(yīng)用，并逐漸成為一種相對(duì)成熟的應(yīng)用技術(shù)［1-5］。

1 地質(zhì)概況

目前，華北地區(qū)東部老油田已進(jìn)入開發(fā)中后期，井網(wǎng)密集，采出程度高，綜合含水高，開發(fā)形勢(shì)日趨嚴(yán)峻。雖已進(jìn)入特高含水后期開發(fā)階段，但仍然是油田調(diào)整挖潛的主要陣地，對(duì)油田的穩(wěn)產(chǎn)具有重要意義。本次數(shù)值模擬研究以東營(yíng)坳陷G2區(qū)塊為例，該區(qū)塊是中間高、邊部低的牽引背斜構(gòu)造，構(gòu)造比較簡(jiǎn)單，地勢(shì)平緩，為整裝高滲油氣藏。該區(qū)自1986年投產(chǎn)，目前分Ng4、Ng5和Ng6三套層系開發(fā)，其中Ng4地質(zhì)儲(chǔ)量為665萬(wàn) t。目的層總共發(fā)育16個(gè)小層，主力層為42、45。主力層儲(chǔ)量占單元總儲(chǔ)量的57.5%。主力層45層整體厚度較小，砂體連片發(fā)育，42層北部層薄、南部砂體連片發(fā)育；非主力層44呈條帶狀分布，其他非主力小層油砂體呈土豆?fàn)盍闵⒎植肌?004年底開始注聚，注聚地層主要覆蓋Ng4和Ng5兩套層系。2011年中旬整體轉(zhuǎn)入后續(xù)水驅(qū)，累計(jì)注入溶液0.692 PV，注入聚合物23 589 t，提高采收率6.99%。

2 油藏?cái)?shù)值模擬

針對(duì)G2區(qū)塊的油藏構(gòu)造特征、儲(chǔ)層展布特征、流體滲流特征、歷史注采情況等動(dòng)靜態(tài)資料，前期建立了符合該區(qū)塊的三維地質(zhì)模型。本項(xiàng)目在三維油藏模型的基礎(chǔ)上，開展有針對(duì)性地適用于G2區(qū)塊的油藏?cái)?shù)值模擬技術(shù)研究。

2.1 建立油藏?cái)?shù)值模擬模型

2.1.1 網(wǎng)格類型的選擇。通過(guò)對(duì)正交網(wǎng)格、角點(diǎn)網(wǎng)格等網(wǎng)格類型的適用性分析，采用角點(diǎn)網(wǎng)格來(lái)建立本次數(shù)值模擬模型。角點(diǎn)網(wǎng)格相較于正交網(wǎng)格更加靈活，可以通過(guò)網(wǎng)格的扭轉(zhuǎn)與變形來(lái)更精準(zhǔn)地模擬斷層、尖滅、剝蝕、不規(guī)則邊界等地質(zhì)現(xiàn)象，是目前相對(duì)來(lái)說(shuō)更加成熟的算法，也是本次數(shù)值模擬的最佳選擇［6］。

2.1.2 數(shù)模模型的網(wǎng)格設(shè)計(jì)。

①建立網(wǎng)格的原則。數(shù)值模擬網(wǎng)格是在地質(zhì)模型網(wǎng)格的基礎(chǔ)上得到的。在建立網(wǎng)格時(shí)遵循了以下幾條原則：一是在條件允許的情況下盡量使網(wǎng)格正交，而且盡量使網(wǎng)格保持均勻；二是雖然角點(diǎn)網(wǎng)格已經(jīng)非常成熟，但在建立角點(diǎn)網(wǎng)格時(shí)不要過(guò)分扭曲網(wǎng)格；三是網(wǎng)格的走向應(yīng)該與油藏內(nèi)流體的流動(dòng)方向大體一致，其是網(wǎng)格劃分應(yīng)該遵守的最重要的一條原則，基于這一原則，考慮大部分水平井的走向，在主斷層設(shè)置網(wǎng)格控制線，使網(wǎng)格與主要大斷層的走向一致，以保證油藏內(nèi)大部分的水平井走向大致與網(wǎng)格方向一致；五是 DX/DY 基本接近于1；六是井與井之間應(yīng)不少于三個(gè)網(wǎng)格。

②平面網(wǎng)格劃分。目前，數(shù)值模擬要求越來(lái)越精細(xì)，但計(jì)算機(jī)硬件是制約計(jì)算速度的瓶頸。為了建立現(xiàn)有的計(jì)算機(jī)條件下更精細(xì)的模型，對(duì)平面網(wǎng)格步長(zhǎng)選擇10 m、20 m、30 m建立三個(gè)模型進(jìn)行優(yōu)化選擇。通過(guò)對(duì)三種規(guī)模的模型進(jìn)行比較，在計(jì)算速度能接受的條件下，選擇了步長(zhǎng)為20 m的網(wǎng)格，使模型較為精細(xì)地體現(xiàn)油藏的非均質(zhì)性。

③縱向網(wǎng)格劃分?？v向上網(wǎng)格步長(zhǎng)的選擇，主要考慮兼顧計(jì)算速度和計(jì)算精度。在油層和油水過(guò)渡帶，縱向網(wǎng)格步長(zhǎng)基本保持0.5 m一個(gè)網(wǎng)格。最終確定網(wǎng)格規(guī)模為136×196×32個(gè)，總數(shù)為852 992 個(gè)，其中有效網(wǎng)格143 796個(gè)。

2.1.3 構(gòu)造模型和儲(chǔ)層模型。油藏?cái)?shù)值模擬模型中的構(gòu)造模型和儲(chǔ)層模型由前期建立的三維地質(zhì)模型粗化得到。模型粗化是使細(xì)網(wǎng)格的精細(xì)地質(zhì)模型“轉(zhuǎn)化”為粗網(wǎng)格模型的過(guò)程。在這一過(guò)程中，用一系列等效的粗網(wǎng)格去“替代”精細(xì)模型中的細(xì)網(wǎng)格，并使該等效網(wǎng)格模型能反映原模型的地質(zhì)特征及流動(dòng)響應(yīng)［7］。根據(jù)數(shù)模的要求，需要輸出粗化的孔滲模型、凈毛比模型。凈毛比模型、孔滲模型均為標(biāo)量，在粗化時(shí)主要應(yīng)用算術(shù)平均法和均方根法進(jìn)行粗化。其中凈毛比采用體積加權(quán)，孔滲采用體積和凈毛比加權(quán)，從而使粗化的模型能夠真實(shí)反映原模型的地質(zhì)特征及流動(dòng)響應(yīng)［8］。

2.1.4 流體模型。

①油水相對(duì)滲透率。根據(jù)區(qū)塊的2口取芯井開發(fā)試驗(yàn)資料，得出該區(qū)塊的相對(duì)滲透率曲線，如圖1所示，用于該油藏的模擬計(jì)算。由圖1可以計(jì)算得出：原始含油飽和度為0.665，流動(dòng)水飽和度為0.345，殘余油飽和度為0.203。

②高壓物性參數(shù)。根據(jù)油井投產(chǎn)初期實(shí)際生產(chǎn)資料，Ng1+2～Ng4小層地面脫氣原油密度為0.96 g/cm3，地面原油黏度為1 066 mPa·s，地下原油密度為0.89 g/cm3，地下原油黏度為38 mPa·s。

2.2 動(dòng)態(tài)模型

該模擬整理了148口油井和45口水井的生產(chǎn)、射孔和措施數(shù)據(jù)。從1986年8月至2019年3月每個(gè)月劃分一個(gè)時(shí)間階段，使生產(chǎn)作業(yè)史與實(shí)際作業(yè)情況保持一致。生產(chǎn)作業(yè)史包括射孔、補(bǔ)孔、改層、封堵、壓裂、酸化、封井等措施作業(yè)的具體時(shí)間和井深，結(jié)合井的注采動(dòng)態(tài)數(shù)據(jù)及產(chǎn)量，通過(guò)SCHEDULE模塊，建立油藏的動(dòng)態(tài)模型［9-10］。

2.3 歷史擬合

在研究區(qū)數(shù)值模擬過(guò)程中，充分應(yīng)用建模數(shù)模一體化技術(shù)，針對(duì)二元驅(qū)階段開發(fā)特點(diǎn)，通過(guò)相滲曲線、油水界面、黏濃關(guān)系、界面張力、非可及孔隙體積等參數(shù)的調(diào)整，完成了研究區(qū)的歷史擬合。

2.3.1 參數(shù)可調(diào)整范圍。在歷史擬合過(guò)程中，如果單以實(shí)際油藏動(dòng)態(tài)數(shù)據(jù)來(lái)確定模型參數(shù)，則參數(shù)調(diào)整自由度較大，導(dǎo)致矛盾較多。為了避免參數(shù)調(diào)整的隨意性，根據(jù)資料來(lái)源和質(zhì)量確定參數(shù)的調(diào)整次序和原則。

①盡量不做調(diào)整的參數(shù)。一是儲(chǔ)層孔隙度、有效厚度數(shù)據(jù)來(lái)自測(cè)井解釋資料，井點(diǎn)處參數(shù)一般不做調(diào)整，對(duì)井間插值的結(jié)果可作適當(dāng)修改；二是對(duì)構(gòu)造形態(tài)、儲(chǔ)層、厚度不作調(diào)整；三是對(duì)地層水、原油的PVT參數(shù)不作調(diào)整。

②可以適當(dāng)調(diào)整的參數(shù)：一是水區(qū)能量和巖石壓縮系數(shù)是調(diào)整油藏能量的主要參數(shù)；二是相對(duì)滲透率曲線和毛管壓力曲線可作為擬合含水和產(chǎn)水量的主要調(diào)整參數(shù)；三是滲透率參數(shù)主要來(lái)源于測(cè)井解釋結(jié)果，誤差較大，因此滲透率特別是方向滲透率可作一定的調(diào)整 ［11-13］。

2.3.2 油藏開發(fā)指標(biāo)的擬合。研究區(qū)歷史擬合包括整體地層壓力、綜合含水、單井含水率、產(chǎn)量等數(shù)據(jù)的擬合。在擬合的過(guò)程中，通過(guò)在參數(shù)合理的調(diào)整范圍以內(nèi)，對(duì)儲(chǔ)層的參數(shù)進(jìn)行調(diào)整，從而達(dá)到更高的擬合精度。

①主要修改參數(shù)。該模擬采取的是定液量生產(chǎn)，所以單井和區(qū)塊的計(jì)算液量和實(shí)際資料是一致的。為了提高單井產(chǎn)水量和含水率的擬合精度，首先整理區(qū)塊所有注采井的注采資料，繪制所有生產(chǎn)井與注水井的單井生產(chǎn)曲線，分析單井生產(chǎn)狀況與含水變化的原因，有目的地修改油藏參數(shù)［14］。主要根據(jù)各井含水情況不同，對(duì)以下參數(shù)進(jìn)行了調(diào)整：一是水區(qū)能量；二是在油水相滲曲線基礎(chǔ)上，對(duì)該曲線進(jìn)行插值和光滑處理，并進(jìn)行細(xì)微調(diào)整；三是根據(jù)各井二次解釋的有效厚度適當(dāng)修改毛管壓力，以調(diào)整各井區(qū)的油水過(guò)渡帶的大小；四是在模擬水驅(qū)油的過(guò)程中，對(duì)滲透率特別是垂向滲透率做了一定程度的調(diào)整；五是擬合底水錐進(jìn)時(shí)，對(duì)夾層的滲透性進(jìn)行了調(diào)整。

②區(qū)塊擬合結(jié)果。此次模擬對(duì)151口油井進(jìn)行了歷史擬合，擬合較好的井有132口，占87.4%，表明油藏的總體擬合效果較好，可以較精準(zhǔn)地表征整體區(qū)塊在開發(fā)過(guò)程中的地層變化、流體流動(dòng)特征等，也為后面剩余油分布規(guī)律研究和方案預(yù)測(cè)提供更精準(zhǔn)的參考，如圖2、圖3 所示。

③單井動(dòng)態(tài)指標(biāo)擬合。區(qū)塊動(dòng)態(tài)指標(biāo)擬合反映的是整個(gè)開發(fā)變化規(guī)律，由于儲(chǔ)層空間的非均質(zhì)性，即使區(qū)塊指標(biāo)擬合良好，也會(huì)存在局部擬合上的較大誤差。因此，必須在區(qū)塊擬合的基礎(chǔ)上進(jìn)行單井動(dòng)態(tài)歷史擬合，以期從宏觀到微觀把握油藏動(dòng)態(tài)規(guī)律，提高后期預(yù)測(cè)分析的精度［15］。單井含水率的擬合主要是通過(guò)對(duì)相對(duì)滲透率曲線端點(diǎn)標(biāo)定與修正來(lái)達(dá)到單井含水率擬合的目的。由于油井計(jì)算控制條件是定液生產(chǎn)，在含水?dāng)M合完成后，日產(chǎn)油指標(biāo)自然得到較好的反映，具體如圖4、圖5所示。

3 剩余油分布規(guī)律

研究區(qū)Ng4單元自1986年投產(chǎn)以來(lái)，經(jīng)歷了天然能量、水驅(qū)、聚驅(qū)、后續(xù)水驅(qū)等開發(fā)階段，綜合含水率高達(dá)96.5%，采出程度為30.3%。隨著儲(chǔ)層及動(dòng)態(tài)認(rèn)識(shí)的不斷深入，暴露出構(gòu)造變化較大、注采流線固定、平面層間矛盾突出等問題， 因此有必要對(duì)研究區(qū)精細(xì)地質(zhì)和剩余油分布進(jìn)行研究，以進(jìn)一步提高采收率。

3.1 平面剩余油分布特征

3.1.1 部分構(gòu)造邊部剩余油相對(duì)富集。A井：南部剩余油飽和度達(dá)0.64，儲(chǔ)量豐度達(dá)45×104 t/km2； B井：東北部剩余油飽和度達(dá)0.63，儲(chǔ)量豐度達(dá)63×104 t/km2，具體如圖6、圖7所示。

3.1.2 油藏內(nèi)部無(wú)井控制處剩余油相對(duì)富集。C井：東部剩余油飽和度達(dá)0.64，儲(chǔ)量豐度達(dá)52×104? t/km2；D井：南部的剩余油飽和度達(dá)0.62，儲(chǔ)量豐度達(dá)44×104? t/km2，具體如圖8、圖9所示。

3.1.3 無(wú)井控制的獨(dú)立砂體剩余油相對(duì)富集。研究區(qū)部分小層砂體較零散，獨(dú)立砂體發(fā)育，其中部分砂體無(wú)生產(chǎn)井，剩余油比較富集，如Ng44的F井西北部，具體如圖10、圖11所示。

3.2 縱向剩余油分布特征

本次數(shù)模Ng4單元總共10個(gè)小層，主力層為42、45小層，主力層儲(chǔ)量占單元總儲(chǔ)量的72.1%；非主力層32、35、41、43、44儲(chǔ)量占總儲(chǔ)量的26.7%；低儲(chǔ)量層為31、33、34，儲(chǔ)量占總儲(chǔ)量的1.2%。

3.2.1 低儲(chǔ)量層基本處于未開發(fā)狀態(tài)。低儲(chǔ)量油層31、33、34剩余油呈土豆?fàn)盍闵⒎植?，含油面積小，儲(chǔ)量低。某井1992年1月鉆開34層，1992年3、4月與42層合采，產(chǎn)量低，時(shí)間短。 31、33層無(wú)油水井，完全處于未開發(fā)狀態(tài)。這3層由于含油面積小儲(chǔ)量低，在后期的開發(fā)方案調(diào)整中也基本不予考慮，具體如圖12、圖13所示。

3.2.2 非主力油層開發(fā)效果較差，但部分層還有一定的開發(fā)潛力。非主力油層由于砂體不連片，呈土豆?fàn)盍闵⒎植迹W(wǎng)很難完善，注采對(duì)應(yīng)率也很低，因此開發(fā)效果較差。如Ng32、Ng35、Ng41層地質(zhì)儲(chǔ)量、剩余儲(chǔ)量及采出程度都很低。但Ng43、Ng44砂體相對(duì)較好，尤其Ng44呈條帶分布，儲(chǔ)量相對(duì)較多，剩余儲(chǔ)量較豐富，具有一定的開發(fā)潛力，具體如圖14、圖15所示。

3.2.3 主力油層開發(fā)效果較理想。主力層Ng42、Ng45 的地質(zhì)儲(chǔ)量、采出程度及剩余儲(chǔ)量較高，開發(fā)效果較理想。主要原因：一是儲(chǔ)層發(fā)育好，砂體連片分布，油藏連通性好；二是注采井網(wǎng)相對(duì)較完善，注采對(duì)應(yīng)率高，具體如圖16、圖17所示。

3.2.4 剩余油主要分布在主力油層。研究區(qū)Ng3-4小層儲(chǔ)量統(tǒng)計(jì)如圖18所示。由圖18可知，主力油層Ng42、Ng45的地質(zhì)儲(chǔ)量和累積產(chǎn)油量都比其他層高，開發(fā)效果最好，同時(shí)剩余儲(chǔ)量也最多，分別占總儲(chǔ)量的27.88%和38.12%，是下一步挖潛的重點(diǎn)。

4 結(jié)論

①在特高含水開發(fā)后期，儲(chǔ)層特征、地層含油性等較開發(fā)初期有很大的變化，通過(guò)油藏?cái)?shù)值模擬來(lái)演化開發(fā)過(guò)程中的動(dòng)態(tài)變化已成為老區(qū)油藏描述中不可或缺的工具，其研究成果可直接應(yīng)用于油田的實(shí)際生產(chǎn)，并已經(jīng)取得較好的應(yīng)用效果。

②對(duì)于不同類型油藏的實(shí)際情況，可制定有針對(duì)性的開發(fā)方案。針對(duì)不同沉積相油藏、不同流體特征油藏等，可通過(guò)調(diào)整相應(yīng)參數(shù)來(lái)進(jìn)行精準(zhǔn)模擬，最終確定剩余油分布規(guī)律，為下一步預(yù)測(cè)與挖潛提供依據(jù)。

參考文獻(xiàn)：

［1］曹國(guó)華，陳鵬，李彩霞.油藏?cái)?shù)值模擬應(yīng)用及未來(lái)發(fā)展趨勢(shì)［J］.內(nèi)蒙古師范大學(xué)學(xué)報(bào)（自然科學(xué)漢文版），2019，48（5）：465-470.

［2］李清清.多層合采井細(xì)分采油潛力研究［D］.大慶：東北石油大學(xué)，2020.

［3］陳一鶴.×油田A區(qū)塊剩余油分布數(shù)值模擬研究［J］.內(nèi)蒙古石油化工，2015，41（13）：114-118.

［4］林博. 勝坨油田河流相儲(chǔ)層建筑結(jié)構(gòu)分析與剩余油分布研究［D］.青島：中國(guó)石油大學(xué)，2007.

［5］馬躍.松遼盆地乾153-59井區(qū)三維地質(zhì)建模和數(shù)值模擬［D］.北京：中國(guó)地質(zhì)大學(xué)（北京），2014.

［6］周為喜，陳玉華，楊永國(guó)，等.基于角點(diǎn)網(wǎng)格的煤層三維建模與可視化研究［J］.煤田地質(zhì)與勘探，2016，44（5）：53-57.

［7］朱玉雙，羅江華，張淑娟，等.等效粗化技術(shù)在任丘碳酸鹽巖油藏中的應(yīng)用：以任丘潛山油藏為例［J］.西北大學(xué)學(xué)報(bào)（自然科學(xué)版），2017，47（3）：455-460.

［8］趙輝，沈煥文，王碧濤，等.儲(chǔ)層三維地質(zhì)建模技術(shù)在吳420油藏的應(yīng)用［J］.內(nèi)蒙古石油化工，2011，37（8）：295-297.

［9］徐強(qiáng)，宋建國(guó).墾71斷塊東二段油藏建模研究［J］.內(nèi)蒙古石油化工，2007（4）：89-91.

［10］董文斌，趙中偉.VIP（油藏?cái)?shù)值模擬）在文79南塊的應(yīng)用［J］.內(nèi)蒙古石油化工，2005（5）：150-151.

［11］張丁涌.史深100區(qū)塊地應(yīng)力及整體壓裂研究與應(yīng)用［D］.成都：西南石油學(xué)院，2004.

［12］吳旭光. W10塊油藏剩余油分布與水力壓裂技術(shù)研究與應(yīng)用［D］.成都：西南石油大學(xué)，2006.

［13］孟強(qiáng).超稠油水平井蒸汽吞吐注采參數(shù)優(yōu)化研究［D］.青島：中國(guó)石油大學(xué)，2008.

［14］葛永濤.真武油田真35斷塊聚合物驅(qū)優(yōu)化方案設(shè)計(jì)［D］.青島：中國(guó)石油大學(xué)，2008.

［15］劉俊山.喬莊油田剩余油分布規(guī)律研究［J］.內(nèi)江科技，2012，33（11）：117，116.