




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)
文檔簡介
1、1 1 什么是油藏數(shù)值模擬?什么是油藏數(shù)值模擬?2 2 油藏數(shù)值模擬軟件介紹油藏數(shù)值模擬軟件介紹3 3 油藏數(shù)值模擬能干些什么?油藏數(shù)值模擬能干些什么?4 4 油藏數(shù)值模擬研究步驟如何?油藏數(shù)值模擬研究步驟如何?5 CMG5 CMG操作簡介操作簡介油藏數(shù)值模擬的定義油藏模擬的概念及其所考慮到的因素油藏模擬的概念及其所考慮到的因素油藏模擬就是把油藏在三維的空間里分油藏模擬就是把油藏在三維的空間里分為許多離散的單元,并且模擬油藏及流為許多離散的單元,并且模擬油藏及流體在空間及一系列離散的時間步里的發(fā)體在空間及一系列離散的時間步里的發(fā)展變化。與物質(zhì)平衡方法一樣的是,系展變化。與物質(zhì)平衡方法一樣的是,
2、系統(tǒng)遵循物質(zhì)守恒原理。統(tǒng)遵循物質(zhì)守恒原理。建立油藏模擬軟件,一般包括:建立油藏模擬軟件,一般包括:(1 1)質(zhì)量守恒原理)質(zhì)量守恒原理(2 2)能量守恒原理)能量守恒原理(3 3)運(yùn)動方程(達(dá)西定律)運(yùn)動方程(達(dá)西定律)(4 4)狀態(tài)方程)狀態(tài)方程(5 5)輔助方程(如飽和度方程,毛管力)輔助方程(如飽和度方程,毛管力方程等方程等) 零維、一維、二維、三維零維、一維、二維、三維 單相、兩相、單相、兩相、 三相三相 單組分、兩組分、單組分、兩組分、N組分組分 單重介質(zhì)單重介質(zhì)(砂巖砂巖)、雙重介質(zhì)、雙重介質(zhì)(碳酸鹽巖碳酸鹽巖)分分 類類 零維、一維、二維、三維零維、一維、二維、三維 單相、兩相、
3、單相、兩相、 三相三相 單組分、兩組分、單組分、兩組分、N組分組分 單重介質(zhì)單重介質(zhì)(砂巖砂巖)、雙重介質(zhì)、雙重介質(zhì)(碳酸鹽巖碳酸鹽巖)按空間維數(shù)來分按空間維數(shù)來分按流體相數(shù)來分按流體相數(shù)來分按流體組分來分按流體組分來分按巖石類型來分按巖石類型來分按模型功能來分按模型功能來分 黑油模型黑油模型 狀態(tài)方程的組分模型狀態(tài)方程的組分模型 雙重介質(zhì)模型雙重介質(zhì)模型 熱采模型熱采模型 聚合物驅(qū)模型聚合物驅(qū)模型 三元復(fù)合驅(qū)模型三元復(fù)合驅(qū)模型 水平井模型水平井模型 1 1 什么是油藏數(shù)值模擬?什么是油藏數(shù)值模擬?2 2 油藏數(shù)值模擬軟件介紹油藏數(shù)值模擬軟件介紹3 3 油藏數(shù)值模擬能干些什么?油藏數(shù)值模擬能干
4、些什么?4 4 油藏數(shù)值模擬研究步驟如何?油藏數(shù)值模擬研究步驟如何?5 CMG5 CMG操作簡介操作簡介油藏數(shù)值模擬軟件對比油藏數(shù)值模擬軟件對比三大應(yīng)用軟件系列水驅(qū)水驅(qū)三次采油三次采油蒸汽驅(qū)蒸汽驅(qū)EclipseEclipseVIPVIPWorkBenchWorkBenchVIPVIPUtchemUtchem CMGCMGSSI-THERMSSI-THERM黑油模型軟件黑油模型軟件 序號序號項(xiàng)目項(xiàng)目VIPVIPCMGCMGECLECL1 1角點(diǎn)網(wǎng)格角點(diǎn)網(wǎng)格一般一般一般一般成熟成熟2 2水平井水平井水平井井筒無壓降水平井井筒無壓降(線源法)(線源法)水平井井筒有壓降水平井井筒有壓降(離散井筒法)(
5、離散井筒法)水平井井筒有壓降水平井井筒有壓降(離散井筒法)(離散井筒法)3 3裂縫描述裂縫描述較方便較方便較方便較方便較方便較方便4 4區(qū)域劃分區(qū)域劃分不允許交叉不允許交叉不允許交叉不允許交叉允許交叉允許交叉5 5流體追蹤流體追蹤APIAPI示蹤示蹤APIAPI和鹽水追蹤和鹽水追蹤6 6巖石性質(zhì)巖石性質(zhì)壓實(shí)壓實(shí)壓實(shí)壓實(shí)7 7相對滲透率相對滲透率方向相對方向相對滲透率滲透率方向相對方向相對滲透率滲透率8 8前處理前處理方便方便一般一般方便方便9 9后處理后處理較強(qiáng)較強(qiáng)一般一般強(qiáng)強(qiáng)1010三維圖形三維圖形強(qiáng)強(qiáng)較強(qiáng)較強(qiáng)強(qiáng)強(qiáng)總體評價總體評價優(yōu)優(yōu) 黑油模型的特殊功能黑油模型的特殊功能角點(diǎn)網(wǎng)格技術(shù)角點(diǎn)網(wǎng)格
6、技術(shù)各種類型的局部網(wǎng)格加密技術(shù)各種類型的局部網(wǎng)格加密技術(shù)斜斷層處理方法斜斷層處理方法不同巖石類型劃分方法不同巖石類型劃分方法不同類型不同類型PVTPVT參數(shù)劃分方法參數(shù)劃分方法飽和度標(biāo)定技術(shù)飽和度標(biāo)定技術(shù)儲層物性時變處理技術(shù)儲層物性時變處理技術(shù)天然裂縫和人工壓裂裂縫模型描述技術(shù)。天然裂縫和人工壓裂裂縫模型描述技術(shù)。 八八項(xiàng)項(xiàng)建建模模技技術(shù)術(shù)三次采油軟件三次采油軟件 序號序號項(xiàng)目項(xiàng)目VIPVIPCMGCMGECLECL1 1聚液粘度依賴聚液粘度依賴聚合物濃度、鹽聚合物濃度、鹽濃度剪切速度濃度剪切速度聚合物濃度聚合物濃度聚合物濃度,鹽聚合物濃度,鹽濃度濃度2 2吸附依賴吸附依賴聚合物濃度、鹽聚合物
7、濃度、鹽濃度、滲透率濃度、滲透率聚合物濃度聚合物濃度滲透率滲透率聚合物濃度聚合物濃度3 3剪切降粘剪切降粘可逆可逆可逆可逆4 4滲透率減小依賴滲透率減小依賴聚合物濃度、鹽聚合物濃度、鹽濃度、滲透率濃度、滲透率滲透率滲透率聚合物濃度聚合物濃度5 5聚合物組分?jǐn)U散聚合物組分?jǐn)U散有有有有6 6復(fù)合驅(qū)復(fù)合驅(qū)聚合物、表活劑聚合物、表活劑7 7陽離子交換陽離子交換有有總體評價總體評價優(yōu)優(yōu)CMGCMG軟件界面軟件界面CMG公司公司ECLIPSEECLIPSE軟件界面軟件界面Schlumberger公司公司VIPVIP軟件界面軟件界面Landmark公司公司1 1 什么是油藏數(shù)值模擬?什么是油藏數(shù)值模擬?2
8、2 油藏數(shù)值模擬軟件介紹油藏數(shù)值模擬軟件介紹3 3 油藏數(shù)值模擬能干些什么?油藏數(shù)值模擬能干些什么?4 4 油藏數(shù)值模擬研究步驟如何?油藏數(shù)值模擬研究步驟如何?5 CMG5 CMG操作簡介操作簡介油藏模擬的目的油藏模擬的目的 1.1.初期開發(fā)方案的模擬初期開發(fā)方案的模擬 2.2.對已開發(fā)油田歷史模擬對已開發(fā)油田歷史模擬 3.3.動態(tài)預(yù)測動態(tài)預(yù)測1). 1). 評價開發(fā)方式;如:枯竭開采、注評價開發(fā)方式;如:枯竭開采、注水開發(fā)等。水開發(fā)等。 2).2).選擇合理井網(wǎng)、開發(fā)層系、確定井位;選擇合理井網(wǎng)、開發(fā)層系、確定井位; 3).3).選擇合理的注采方式、注采比;選擇合理的注采方式、注采比;4).
9、4).對油藏和流體性質(zhì)敏感性研究。對油藏和流體性質(zhì)敏感性研究。1).1).核實(shí)地質(zhì)儲量,確定基本的驅(qū)替機(jī)理核實(shí)地質(zhì)儲量,確定基本的驅(qū)替機(jī)理(如:是天然驅(qū),還是注水開發(fā));(如:是天然驅(qū),還是注水開發(fā)); 2).2).確定產(chǎn)液量和生產(chǎn)周期;確定產(chǎn)液量和生產(chǎn)周期;3).3).確定油藏和流體特性;確定油藏和流體特性;4).4).提出問題、潛力所在區(qū)域提出問題、潛力所在區(qū)域。1).1).開發(fā)指標(biāo)預(yù)測及經(jīng)濟(jì)評價開發(fā)指標(biāo)預(yù)測及經(jīng)濟(jì)評價2).2).評價提高采收率的方法(如:一次采評價提高采收率的方法(如:一次采油、注水、注氣、化學(xué)驅(qū)等)油、注水、注氣、化學(xué)驅(qū)等)3).3).剩余油飽和度分布規(guī)律的研究剩余油飽
10、和度分布規(guī)律的研究,再現(xiàn),再現(xiàn)生產(chǎn)歷史動態(tài)生產(chǎn)歷史動態(tài)4).4).潛力評價和提高采收率的方向潛力評價和提高采收率的方向5).5).專題和機(jī)理問題的研究專題和機(jī)理問題的研究油藏描述油藏描述 油藏描述是油田開發(fā)的基礎(chǔ),是一項(xiàng)系統(tǒng)工程,由多學(xué)科各種方法聯(lián)合研究的結(jié)油藏描述是油田開發(fā)的基礎(chǔ),是一項(xiàng)系統(tǒng)工程,由多學(xué)科各種方法聯(lián)合研究的結(jié)果。油藏數(shù)值模擬作為一種方法,在油藏描述中起了一定的作用。果。油藏數(shù)值模擬作為一種方法,在油藏描述中起了一定的作用。 1) 1) 孔隙結(jié)構(gòu)研究孔隙結(jié)構(gòu)研究 -10m-10m級級 利用利用CT CT 、核磁共振、圖象分析儀、微觀驅(qū)油機(jī)理、毛管壓力實(shí)驗(yàn)等進(jìn)行研究、核磁共振、圖
11、象分析儀、微觀驅(qū)油機(jī)理、毛管壓力實(shí)驗(yàn)等進(jìn)行研究 2)2)巖石物理性質(zhì)研究巖石物理性質(zhì)研究 -cm -cm 級級 利用油層物理實(shí)驗(yàn)室進(jìn)行利用油層物理實(shí)驗(yàn)室進(jìn)行、K K、S S、k kr r等研究。等研究。 3)3)測井研究測井研究 -m-m級級 利用測井方法進(jìn)行利用測井方法進(jìn)行、K K、S S等研究。等研究。 4)4)井間地震研究井間地震研究 1-100m1-100m級級 利用井間地震進(jìn)行油層連續(xù)性、油層厚度等研究。利用井間地震進(jìn)行油層連續(xù)性、油層厚度等研究。 5) 5) 井間監(jiān)測井間監(jiān)測 m-1000mm-1000m 利用試井和示蹤劑等方法來研究主力油層的利用試井和示蹤劑等方法來研究主力油層的
12、、K K、S S等。等。 6)6)數(shù)值模擬網(wǎng)格數(shù)值模擬網(wǎng)格 10m-1000m10m-1000m 將上述研究的將上述研究的、K K、S S等參數(shù)填入數(shù)值模擬網(wǎng)格。等參數(shù)填入數(shù)值模擬網(wǎng)格。 7)7)油藏數(shù)值模擬油藏數(shù)值模擬 1000m-10000m1000m-10000m 利用油藏數(shù)值方法研究油藏靜態(tài)參數(shù)和利用油藏數(shù)值方法研究油藏靜態(tài)參數(shù)和P P、S S分布。分布。 8)8)小井距生產(chǎn)試驗(yàn)區(qū)小井距生產(chǎn)試驗(yàn)區(qū) 如大慶、孤東等小井距研究彈性驅(qū)、水驅(qū)和各種提高采收率的開發(fā)規(guī)律如大慶、孤東等小井距研究彈性驅(qū)、水驅(qū)和各種提高采收率的開發(fā)規(guī)律。驅(qū)油機(jī)理研究驅(qū)油機(jī)理研究層內(nèi)油水運(yùn)動機(jī)理研究油水運(yùn)動的作用力油水
13、運(yùn)動的作用力Pghkcos驅(qū)動力驅(qū)動力重重 力力毛管力毛管力控制地質(zhì)儲量、井間連通狀控制地質(zhì)儲量、井間連通狀況水體大小、驅(qū)動能量況水體大小、驅(qū)動能量 、 產(chǎn)能產(chǎn)能, ,為方案研究提供依據(jù)為方案研究提供依據(jù)層系、井網(wǎng)、井距、層系、井網(wǎng)、井距、開發(fā)方式等優(yōu)化開發(fā)方式等優(yōu)化敏感性分析方案指標(biāo)預(yù)測敏感性分析方案指標(biāo)預(yù)測驅(qū)替機(jī)理、合理注入條件驅(qū)替機(jī)理、合理注入條件及參數(shù)及參數(shù) 、動態(tài)預(yù)測、動態(tài)預(yù)測剩余油分布研究剩余油分布研究優(yōu)化調(diào)整部署優(yōu)化調(diào)整部署早期評價早期評價模擬模擬新油田開新油田開發(fā)方案模發(fā)方案模擬擬老油田調(diào)老油田調(diào)整方案模整方案模擬擬提高采收提高采收率模擬率模擬開發(fā)前期開發(fā)前期開發(fā)初期開發(fā)初期開
14、發(fā)中期開發(fā)中期開發(fā)后期開發(fā)后期油藏動態(tài)預(yù)測油藏動態(tài)預(yù)測1 1 什么是油藏數(shù)值模擬?什么是油藏數(shù)值模擬?2 2 油藏數(shù)值模擬軟件介紹油藏數(shù)值模擬軟件介紹3 3 油藏數(shù)值模擬能干些什么?油藏數(shù)值模擬能干些什么?4 4 油藏數(shù)值模擬研究步驟如何?油藏數(shù)值模擬研究步驟如何?5 CMG5 CMG操作簡介操作簡介l油藏模擬計劃的擬訂l數(shù)值模擬軟件選取l數(shù)值模擬研究所需數(shù)據(jù)準(zhǔn)備l建立數(shù)值模擬數(shù)據(jù)文件l歷史擬合l動態(tài)預(yù)測l開發(fā)報告編寫每一步每一步的目的?的目的?1 1 油藏模擬計劃的擬訂油藏模擬計劃的擬訂2 2 模擬器的選取模擬器的選取天然裂縫油氣藏的模擬應(yīng)用雙重介質(zhì)的裂縫模型;天然裂縫油氣藏的模擬應(yīng)用雙重介
15、質(zhì)的裂縫模型;凝析油氣田開發(fā)以及注氣油田開采應(yīng)用組份模型;凝析油氣田開發(fā)以及注氣油田開采應(yīng)用組份模型;視油、氣田含富氣或平氣大小可采用組份或黑油模視油、氣田含富氣或平氣大小可采用組份或黑油模型;型;稠油開發(fā)用熱采注蒸汽模型;稠油開發(fā)用熱采注蒸汽模型;注各種化學(xué)劑的三次采油應(yīng)用化學(xué)模型;注各種化學(xué)劑的三次采油應(yīng)用化學(xué)模型;而對一般油藏,模擬可選用常規(guī)油氣田開發(fā)的黑油而對一般油藏,模擬可選用常規(guī)油氣田開發(fā)的黑油模型。模型。3 數(shù)據(jù)準(zhǔn)備 1.1.表格數(shù)據(jù)表格數(shù)據(jù) 主要指巖石物性、流體性質(zhì)。主要指巖石物性、流體性質(zhì)。 油氣油氣PVTPVT數(shù)據(jù)表(高壓物性分析);數(shù)據(jù)表(高壓物性分析); 水及巖石水及巖
16、石PVTPVT性質(zhì)(高壓物性分析);性質(zhì)(高壓物性分析); 油水相對滲透曲線;油水相對滲透曲線; 毛管壓力曲線(巖心壓汞實(shí)驗(yàn))等;毛管壓力曲線(巖心壓汞實(shí)驗(yàn))等; 井筒流動數(shù)據(jù)。井筒流動數(shù)據(jù)。 2.2.網(wǎng)格數(shù)據(jù)網(wǎng)格數(shù)據(jù) 油藏頂面海拔深度;油藏頂面海拔深度; 砂層厚度(有效);砂層厚度(有效); 孔、滲、飽參數(shù)巖石類型等。孔、滲、飽參數(shù)巖石類型等。 3. 3.動態(tài)數(shù)據(jù)動態(tài)數(shù)據(jù) 完井?dāng)?shù)據(jù):射孔、補(bǔ)孔、壓裂、堵水、解堵日期、層位、井指數(shù)等;完井?dāng)?shù)據(jù):射孔、補(bǔ)孔、壓裂、堵水、解堵日期、層位、井指數(shù)等; 生產(chǎn)數(shù)據(jù):平均日產(chǎn)油、日產(chǎn)水、日產(chǎn)氣、平均油氣比和含水比等;生產(chǎn)數(shù)據(jù):平均日產(chǎn)油、日產(chǎn)水、日產(chǎn)氣、平
17、均油氣比和含水比等; 壓力數(shù)據(jù):井底流壓、網(wǎng)格壓力等。壓力數(shù)據(jù):井底流壓、網(wǎng)格壓力等。 動態(tài)監(jiān)測資料(分層測試、吸水、產(chǎn)液剖面等)動態(tài)監(jiān)測資料(分層測試、吸水、產(chǎn)液剖面等) 4.4.其它數(shù)據(jù)其它數(shù)據(jù) 主要包括算法選擇、輸入輸出控制、油水井約束界限、油井定壓定產(chǎn)等主要包括算法選擇、輸入輸出控制、油水井約束界限、油井定壓定產(chǎn)等參數(shù)。參數(shù)。各砂巖組個數(shù)及名稱各砂巖組個數(shù)及名稱小層名稱小層名稱沉積類型沉積類型連通情況連通情況油層分層數(shù)據(jù)油層分層數(shù)據(jù)斷層斷層數(shù)據(jù)數(shù)據(jù)斷層方向斷層方向斷點(diǎn)深度斷點(diǎn)深度斷距斷距密封情況密封情況小層平面圖小層平面圖頂部深度等值圖頂部深度等值圖有效厚度等值圖有效厚度等值圖滲透率等
18、值圖滲透率等值圖孔隙度等值圖孔隙度等值圖地質(zhì)構(gòu)架地質(zhì)構(gòu)架屬性數(shù)據(jù)屬性數(shù)據(jù)小層靜態(tài)數(shù)據(jù)小層靜態(tài)數(shù)據(jù) 油層流體性質(zhì)數(shù)據(jù)油層流體性質(zhì)數(shù)據(jù)油、氣、水的高壓物性資料油、氣、水的高壓物性資料巖石的壓縮系數(shù)巖石的壓縮系數(shù)注意:輸入時最好將注意:輸入時最好將PVT數(shù)據(jù)繪成曲線,以保證輸入模型數(shù)據(jù)繪成曲線,以保證輸入模型中的中的PVT曲線光滑曲線光滑包括:原油的包括:原油的Bo-pBo-p曲線曲線,Uo-pUo-p曲線曲線,Rso-pRso-p曲線,曲線,UoUo-T-T曲曲線,此外還有原油的地下密度,地面密度,飽和壓力,線,此外還有原油的地下密度,地面密度,飽和壓力,地面粘度,壓縮系數(shù)等。地面粘度,壓縮系數(shù)等
19、。 氣體的氣體的Bg-pBg-p曲線,曲線, Ug-pUg-p曲線,地下地面密度及壓曲線,地下地面密度及壓縮系數(shù)地層水的粘度,礦化度,體積系數(shù)和壓縮系數(shù)??s系數(shù)地層水的粘度,礦化度,體積系數(shù)和壓縮系數(shù)。地層壓力原油體積系數(shù)Pb 巖石和流體的流動性質(zhì)數(shù)據(jù)巖石和流體的流動性質(zhì)數(shù)據(jù) 油水相滲曲線油水相滲曲線 油氣相滲曲線油氣相滲曲線 毛管力曲線毛管力曲線注意:油氣相滲曲線資料少,可借用注意:油氣相滲曲線資料少,可借用 相滲曲線和毛管力曲線飽和度端點(diǎn)值應(yīng)匹配相滲曲線和毛管力曲線飽和度端點(diǎn)值應(yīng)匹配 毛管力曲線對純油區(qū)可以忽略,但對過渡帶必須考慮毛管力曲線對純油區(qū)可以忽略,但對過渡帶必須考慮 初始(初始(
20、t=0t=0)的壓力、含油飽和度分布)的壓力、含油飽和度分布. .初始壓力,壓力平衡化算法初始壓力,壓力平衡化算法. .初始含油飽和度,直接賦值初始含油飽和度,直接賦值 初始數(shù)據(jù)初始數(shù)據(jù)輸入各生產(chǎn)階段每口井的生產(chǎn)動態(tài)數(shù)據(jù)輸入各生產(chǎn)階段每口井的生產(chǎn)動態(tài)數(shù)據(jù) 生產(chǎn)動態(tài)數(shù)據(jù)指各生產(chǎn)井及注水井在各個開發(fā)時刻的日生產(chǎn)動態(tài)數(shù)據(jù)指各生產(chǎn)井及注水井在各個開發(fā)時刻的日油量,日水量,含水率,油氣比,井底壓力,井口壓力,日油量,日水量,含水率,油氣比,井底壓力,井口壓力,日注量等數(shù)據(jù),并將這些數(shù)據(jù)按一定時間步長整理,做成一系注量等數(shù)據(jù),并將這些數(shù)據(jù)按一定時間步長整理,做成一系列時間卡,按定產(chǎn)條件進(jìn)行輸入,若有的井與
21、其它層位進(jìn)行列時間卡,按定產(chǎn)條件進(jìn)行輸入,若有的井與其它層位進(jìn)行合采合注的,應(yīng)視具體情況進(jìn)行劈產(chǎn)劈注。合采合注的,應(yīng)視具體情況進(jìn)行劈產(chǎn)劈注。 對各井進(jìn)行的措施也應(yīng)給以考慮,各井射孔位置、日期,對各井進(jìn)行的措施也應(yīng)給以考慮,各井射孔位置、日期,補(bǔ)孔記錄,壓裂堵水記錄,大修記錄等各項(xiàng)措施記錄。補(bǔ)孔記錄,壓裂堵水記錄,大修記錄等各項(xiàng)措施記錄。 測壓資料測壓資料 出油剖面、吸水剖面出油剖面、吸水剖面 見水日期、層位、來水方向等見水日期、層位、來水方向等其它有關(guān)數(shù)據(jù)其它有關(guān)數(shù)據(jù)供歷史擬合時參考供歷史擬合時參考1 1 什么是油藏數(shù)值模擬?什么是油藏數(shù)值模擬?2 2 油藏數(shù)值模擬軟件介紹油藏數(shù)值模擬軟件介紹
22、3 3 油藏數(shù)值模擬能干些什么?油藏數(shù)值模擬能干些什么?4 4 油藏數(shù)值模擬研究步驟如何?油藏數(shù)值模擬研究步驟如何?5 CMG5 CMG操作簡介操作簡介l數(shù)值模擬數(shù)據(jù)文件(前處理建模)l歷史擬合l結(jié)果分析(后處理)l動態(tài)預(yù)測l開發(fā)報告編寫每一步每一步的目的?的目的?CMG主界面CMGCMG所包括所包括模塊模塊 數(shù)據(jù)文件建數(shù)據(jù)文件建立模塊立模塊黑油模型模黑油模型模塊塊組分模型模組分模型模塊塊三維立體顯三維立體顯示模塊示模塊自定義平面自定義平面圖顯示模塊圖顯示模塊自定義表表自定義表表輸出模塊輸出模塊熱采和化學(xué)熱采和化學(xué)驅(qū)模塊驅(qū)模塊相態(tài)及流體性質(zhì)相態(tài)及流體性質(zhì)計算模塊計算模塊轉(zhuǎn)換模塊轉(zhuǎn)換模塊Sect
23、ionContent1I/O Control Section名字、輸入輸出單位、各控制參數(shù)等。名字、輸入輸出單位、各控制參數(shù)等。2Reservoir Description網(wǎng)格數(shù)目、井口半徑、深度、孔隙度、網(wǎng)格數(shù)目、井口半徑、深度、孔隙度、ijkijk方向滲透率、有效厚度等。方向滲透率、有效厚度等。3Component Properties選擇模式、選擇模式、PVTPVT參數(shù)(油氣水三相密度、參數(shù)(油氣水三相密度、粘度、體積系數(shù))、油氣水巖石的壓縮粘度、體積系數(shù))、油氣水巖石的壓縮系數(shù)等系數(shù)等4Rock-Fluid Property油水、氣水相對滲透率曲線油水、氣水相對滲透率曲線5Initia
24、l Condition飽和度、泡點(diǎn)壓力、參考壓力、參考深飽和度、泡點(diǎn)壓力、參考壓力、參考深度、油水、氣水界面深度度、油水、氣水界面深度。6Numerical Control數(shù)值模式選擇、最大最小生產(chǎn)時間。數(shù)值模式選擇、最大最小生產(chǎn)時間。7Well Data生產(chǎn)時間、井的類別、生產(chǎn)流體、最大生產(chǎn)時間、井的類別、生產(chǎn)流體、最大產(chǎn)量、最小井底流壓、表皮系數(shù)、射孔產(chǎn)量、最小井底流壓、表皮系數(shù)、射孔位置、射孔時間。位置、射孔時間。網(wǎng)格及模擬層的劃分網(wǎng)格及模擬層的劃分模擬區(qū)域的選擇模擬區(qū)域的選擇 明確數(shù)模所研究的區(qū)域,確定模擬的明確數(shù)模所研究的區(qū)域,確定模擬的范圍,井?dāng)?shù)、層范圍,井?dāng)?shù)、層位和數(shù)模的時間位和
25、數(shù)模的時間長短長短模擬層的劃分模擬層的劃分 縱向上進(jìn)行網(wǎng)格劃分,也就是縱向上進(jìn)行網(wǎng)格劃分,也就是確定數(shù)模的層位和層數(shù)確定數(shù)模的層位和層數(shù),數(shù)模時所考慮的層位與生產(chǎn)層位有較大區(qū)別,在劃分時數(shù)模時所考慮的層位與生產(chǎn)層位有較大區(qū)別,在劃分時應(yīng)考慮以下幾點(diǎn):應(yīng)考慮以下幾點(diǎn):在能夠說明問題的前提下,在能夠說明問題的前提下,盡量少用縱向網(wǎng)格盡量少用縱向網(wǎng)格,因?yàn)榭v,因?yàn)榭v向?qū)訑?shù)對網(wǎng)格節(jié)點(diǎn)數(shù)影響較大向?qū)訑?shù)對網(wǎng)格節(jié)點(diǎn)數(shù)影響較大對層間對層間有較厚隔層有較厚隔層,分隔狀況良好的生產(chǎn)層位,應(yīng)適當(dāng),分隔狀況良好的生產(chǎn)層位,應(yīng)適當(dāng)進(jìn)行分開劃分,以準(zhǔn)確描述其層內(nèi)流體分布;對有一定進(jìn)行分開劃分,以準(zhǔn)確描述其層內(nèi)流體分布;對
26、有一定隔層,但隔層,但分隔不明顯分隔不明顯,且大面積連通的層位可考慮將兩,且大面積連通的層位可考慮將兩層合并處理;對層合并處理;對分布較小的小砂體分布較小的小砂體,且與其它層位儲量,且與其它層位儲量差別比較大的小層,可考慮將其平加到其相鄰層位上去,差別比較大的小層,可考慮將其平加到其相鄰層位上去,這樣對結(jié)果影響不大這樣對結(jié)果影響不大對對層位較厚的大層層位較厚的大層,由于韻律性不同,層內(nèi)流體分布不,由于韻律性不同,層內(nèi)流體分布不同,應(yīng)將其分開。上下層之間的聯(lián)系靠垂向滲透率來控同,應(yīng)將其分開。上下層之間的聯(lián)系靠垂向滲透率來控制,一般取平面滲透率的制,一般取平面滲透率的1/51/5到到1/10.1/
27、10. 網(wǎng)格劃分包括兩個方面:網(wǎng)格劃分包括兩個方面:網(wǎng)格方向及網(wǎng)格尺寸網(wǎng)格方向及網(wǎng)格尺寸。因?yàn)橐话隳P投伎紤]的封閉邊界,所以因?yàn)橐话隳P投伎紤]的封閉邊界,所以網(wǎng)格邊界應(yīng)盡量網(wǎng)格邊界應(yīng)盡量與封閉的天然邊界一致與封閉的天然邊界一致,如斷層、尖滅、油水邊界、井,如斷層、尖滅、油水邊界、井排方向等,這樣做可以減少無效網(wǎng)格的數(shù)目,同時也可排方向等,這樣做可以減少無效網(wǎng)格的數(shù)目,同時也可以較精確地描述邊界形態(tài)。以較精確地描述邊界形態(tài)。網(wǎng)格的定向應(yīng)盡量網(wǎng)格的定向應(yīng)盡量包含全部井位及將來可能進(jìn)行的加密包含全部井位及將來可能進(jìn)行的加密井井,須取全取準(zhǔn)。,須取全取準(zhǔn)。網(wǎng)格方向應(yīng)注意考慮以下幾點(diǎn):網(wǎng)格方向應(yīng)注意考慮
28、以下幾點(diǎn):網(wǎng)格取向還應(yīng)考慮:網(wǎng)格取向還應(yīng)考慮:流體流動的方向和油藏內(nèi)天然勢流體流動的方向和油藏內(nèi)天然勢梯度的方向梯度的方向,對于不同的網(wǎng)格系統(tǒng),其動態(tài)可能會不,對于不同的網(wǎng)格系統(tǒng),其動態(tài)可能會不同,尤其在流度比不利的情況下,更應(yīng)注意網(wǎng)格所造同,尤其在流度比不利的情況下,更應(yīng)注意網(wǎng)格所造成的影響,也可用數(shù)學(xué)上的方法來解決。成的影響,也可用數(shù)學(xué)上的方法來解決。應(yīng)考慮應(yīng)考慮網(wǎng)格方向與滲透率主軸相平行網(wǎng)格方向與滲透率主軸相平行應(yīng)應(yīng)盡量減少無效網(wǎng)格數(shù)目盡量減少無效網(wǎng)格數(shù)目網(wǎng)格尺寸應(yīng)考慮:網(wǎng)格尺寸應(yīng)考慮:考慮考慮機(jī)器能力,選取網(wǎng)格尺寸機(jī)器能力,選取網(wǎng)格尺寸。網(wǎng)格越小,節(jié)點(diǎn)數(shù)越。網(wǎng)格越小,節(jié)點(diǎn)數(shù)越多,多,C
29、PUCPU所需量越高,費(fèi)用、時間也就越多,同時還所需量越高,費(fèi)用、時間也就越多,同時還應(yīng)考慮模型節(jié)點(diǎn)所能承受的能力應(yīng)考慮模型節(jié)點(diǎn)所能承受的能力精細(xì)油藏模擬時,精細(xì)油藏模擬時,應(yīng)使用足夠多的網(wǎng)格,應(yīng)使用足夠多的網(wǎng)格,使其能準(zhǔn)確使其能準(zhǔn)確反映油藏結(jié)構(gòu)和參數(shù)在空間中的變化規(guī)律,不能以大反映油藏結(jié)構(gòu)和參數(shù)在空間中的變化規(guī)律,不能以大網(wǎng)格掩蓋了其間的變化,如小尖滅、小構(gòu)造和小砂體網(wǎng)格掩蓋了其間的變化,如小尖滅、小構(gòu)造和小砂體在在研究流體運(yùn)動規(guī)律時,必須使用較多的網(wǎng)格研究流體運(yùn)動規(guī)律時,必須使用較多的網(wǎng)格,以便能夠,以便能夠控制和跟蹤流體界面的流動。若采用的網(wǎng)格太粗,就會使控制和跟蹤流體界面的流動。若采用
30、的網(wǎng)格太粗,就會使網(wǎng)格飽和度變化平均化,不能刻劃出流體的變化趨勢網(wǎng)格飽和度變化平均化,不能刻劃出流體的變化趨勢 避免尺寸大小突變避免尺寸大小突變,相鄰網(wǎng)格尺寸比小于,相鄰網(wǎng)格尺寸比小于2-32-3,否則會帶來,否則會帶來很大的截斷誤差很大的截斷誤差 網(wǎng)格尺寸應(yīng)與井位相適應(yīng)網(wǎng)格尺寸應(yīng)與井位相適應(yīng),保證一個網(wǎng)格內(nèi)只能有一口井。,保證一個網(wǎng)格內(nèi)只能有一口井。兩口井之間至少有一個空網(wǎng)格兩口井之間至少有一個空網(wǎng)格 井井密集區(qū)和主要模擬區(qū)應(yīng)該適當(dāng)增加網(wǎng)格數(shù)目密集區(qū)和主要模擬區(qū)應(yīng)該適當(dāng)增加網(wǎng)格數(shù)目,井周圍通,井周圍通常用密網(wǎng)格,其他地區(qū)可用稀網(wǎng)格常用密網(wǎng)格,其他地區(qū)可用稀網(wǎng)格模擬層的劃分模擬小模擬小層層油砂
31、體油砂體地質(zhì)儲量地質(zhì)儲量(104t)1E2d2324.22E2d23310.03E2d23415.74E2d2364.55E2d23731.6合計合計 66縱向縱向江蘇油田真江蘇油田真1212斷塊斷塊E E2 2d d2 23 3網(wǎng)格大?。壕W(wǎng)格大?。?52525m25m,網(wǎng)格個數(shù):,網(wǎng)格個數(shù):66 66 2222平面平面三維圖三維圖l數(shù)值模擬數(shù)據(jù)文件(前處理建模)l歷史擬合l結(jié)果分析(后處理)l動態(tài)預(yù)測l開發(fā)報告編寫每一步每一步的目的?的目的?歷史擬合歷史擬合歷史擬合目的歷史擬合目的 歷史擬合的目的是比較客觀地認(rèn)識油田的過去歷史擬合的目的是比較客觀地認(rèn)識油田的過去和現(xiàn)狀,為開發(fā)動態(tài)預(yù)測打下基礎(chǔ)
32、。和現(xiàn)狀,為開發(fā)動態(tài)預(yù)測打下基礎(chǔ)。歷史擬合方法歷史擬合方法 直接求解直接求解 試湊法試湊法 自動歷史擬合自動歷史擬合歷史擬合指標(biāo)歷史擬合指標(biāo) 歷史擬合包括全油藏的擬合和單井的擬合,擬合歷史擬合包括全油藏的擬合和單井的擬合,擬合的指標(biāo)有:的指標(biāo)有: 儲量儲量、壓力、壓力、產(chǎn)量產(chǎn)量、含水率等含水率等參數(shù)調(diào)整原則參數(shù)調(diào)整原則 不同參數(shù)的不同組合會得到相同的計算結(jié)不同參數(shù)的不同組合會得到相同的計算結(jié)果。為了避免參數(shù)修改的任意性,在歷史擬合果。為了避免參數(shù)修改的任意性,在歷史擬合開始前必須開始前必須確定各參數(shù)的可調(diào)范圍確定各參數(shù)的可調(diào)范圍,判斷參數(shù)判斷參數(shù)來源是否可靠來源是否可靠,確定的參數(shù)一般不修改,
33、或只,確定的參數(shù)一般不修改,或只在較小范圍內(nèi)修改;不確定的參數(shù)允許修改,在較小范圍內(nèi)修改;不確定的參數(shù)允許修改,可在較大范圍內(nèi)修改。可在較大范圍內(nèi)修改。 具體應(yīng)用時應(yīng)根據(jù)實(shí)際情況而定具體應(yīng)用時應(yīng)根據(jù)實(shí)際情況而定2.2.可調(diào)參數(shù)及調(diào)參范圍可調(diào)參數(shù)及調(diào)參范圍 滲透率滲透率孔隙度孔隙度有效厚度有效厚度 油層測井解釋的有效厚度與取芯資料對比,油層測井解釋的有效厚度與取芯資料對比,一般偏高一般偏高3030左右,主要是鈣質(zhì)層和泥質(zhì)夾層左右,主要是鈣質(zhì)層和泥質(zhì)夾層沒有完全扣除,因此可調(diào)范圍為沒有完全扣除,因此可調(diào)范圍為-30%-30%0 0。 油層孔隙度變化范圍不大。因此孔隙度視油層孔隙度變化范圍不大。因此
34、孔隙度視為確定參數(shù),不做修改,或允許改動范圍在為確定參數(shù),不做修改,或允許改動范圍在3 3。巖石和流體的壓縮系數(shù)巖石和流體的壓縮系數(shù) 流體壓縮系數(shù)流體壓縮系數(shù)是實(shí)驗(yàn)測定的,變化范是實(shí)驗(yàn)測定的,變化范圍較小,認(rèn)為是確定的。而圍較小,認(rèn)為是確定的。而巖石壓縮系巖石壓縮系數(shù)數(shù)雖然也是實(shí)驗(yàn)測定的,但受巖石內(nèi)飽雖然也是實(shí)驗(yàn)測定的,但受巖石內(nèi)飽和流體和應(yīng)力狀態(tài)等的影響,有一定的和流體和應(yīng)力狀態(tài)等的影響,有一定的變化范圍,而且與有效厚度相連的非有變化范圍,而且與有效厚度相連的非有效部分也有一定彈性作用,考慮這部分效部分也有一定彈性作用,考慮這部分的影響,允許巖石的壓縮系數(shù)可以的影響,允許巖石的壓縮系數(shù)可以擴(kuò)
35、大擴(kuò)大一倍一倍。油氣油氣PVT性質(zhì)性質(zhì) 一般認(rèn)為是確定參數(shù),在具體情況下也允一般認(rèn)為是確定參數(shù),在具體情況下也允許作適當(dāng)修改,應(yīng)根據(jù)具體情況而定。許作適當(dāng)修改,應(yīng)根據(jù)具體情況而定。初始流體飽和度初始流體飽和度 一般認(rèn)為是確定參數(shù),必要時允許小范一般認(rèn)為是確定參數(shù),必要時允許小范圍修改。圍修改。相對滲透率曲線相對滲透率曲線歷史擬合的一般操作方法歷史擬合的一般操作方法 數(shù)值模擬的數(shù)值模擬的第一步是對油藏的地質(zhì)儲量作出擬合第一步是對油藏的地質(zhì)儲量作出擬合,因,因?yàn)榈刭|(zhì)儲量是一個比較敏感的參數(shù),一般上報儲量會并確為地質(zhì)儲量是一個比較敏感的參數(shù),一般上報儲量會并確定下來,以后的各種開發(fā)指標(biāo)均以此地質(zhì)儲量
36、為準(zhǔn)。當(dāng)?shù)囟ㄏ聛?,以后的各種開發(fā)指標(biāo)均以此地質(zhì)儲量為準(zhǔn)。當(dāng)?shù)刭|(zhì)儲量有差別時,質(zhì)儲量有差別時,應(yīng)先對輸入的參數(shù)做檢查,看是否存在應(yīng)先對輸入的參數(shù)做檢查,看是否存在較大誤差。較大誤差。地質(zhì)儲量地質(zhì)儲量相關(guān)參數(shù):飽和度、有效厚度、孔隙度、壓縮系數(shù)等相關(guān)參數(shù):飽和度、有效厚度、孔隙度、壓縮系數(shù)等模擬層模擬層油砂體油砂體實(shí)際儲量實(shí)際儲量(104t)計算儲量計算儲量(104t)誤差()誤差()1E2d2324.23.9256.52E2d23310.09.6343.73E2d23415.716.2143.34E2d2364.54.1208.45E2d23731.632.5733.1合計合計 6666.466
37、0.7儲量壓力擬合壓力擬合 從上圖中可以看出,從整個變化趨勢上看,擬合壓力從上圖中可以看出,從整個變化趨勢上看,擬合壓力低于觀察壓力,但兩條曲線幾乎平行,這是典型的油藏滲低于觀察壓力,但兩條曲線幾乎平行,這是典型的油藏滲透率給低了的情況,若將油藏透率給低了的情況,若將油藏滲透率增加為原來的滲透率增加為原來的1.51.5,則可以實(shí)現(xiàn)最終擬合。則可以實(shí)現(xiàn)最終擬合。 從上圖中可以看出,計算壓力值低于觀察值且繼續(xù)發(fā)從上圖中可以看出,計算壓力值低于觀察值且繼續(xù)發(fā)散,這是典型的由于遠(yuǎn)離井區(qū)的滲透率和孔隙度給定的過散,這是典型的由于遠(yuǎn)離井區(qū)的滲透率和孔隙度給定的過低造成的結(jié)果,在此情況下可保持低造成的結(jié)果,
38、在此情況下可保持油藏滲透率不變油藏滲透率不變,將,將水水層的滲透率提高到原來的層的滲透率提高到原來的1.51.5倍倍,就可以實(shí)現(xiàn)擬合。,就可以實(shí)現(xiàn)擬合。含水率擬合含水率擬合 例如位于油區(qū)的油井開井例如位于油區(qū)的油井開井初期見水初期見水,而實(shí)際情況并,而實(shí)際情況并沒有見水,有時出現(xiàn)沒有見水,有時出現(xiàn)見水滯后見水滯后,其主要原因可能是:,其主要原因可能是:相滲關(guān)系的端點(diǎn)值與初始流體飽和度不匹配,需要相滲關(guān)系的端點(diǎn)值與初始流體飽和度不匹配,需要進(jìn)行調(diào)整和修正進(jìn)行調(diào)整和修正。見水過早見水過早見水滯后見水滯后水相相對滲透率曲線調(diào)整水相相對滲透率曲線調(diào)整 調(diào)整滲透率的分布或原油粘度后,壓力分布和調(diào)整滲透率
39、的分布或原油粘度后,壓力分布和單井壓力都發(fā)生變化,當(dāng)調(diào)整相對滲透率曲線單井壓力都發(fā)生變化,當(dāng)調(diào)整相對滲透率曲線時,由于該相流量也隨之發(fā)生相應(yīng)的改變,必時,由于該相流量也隨之發(fā)生相應(yīng)的改變,必然導(dǎo)致壓力變化,例如當(dāng)水的相對滲透率曲線然導(dǎo)致壓力變化,例如當(dāng)水的相對滲透率曲線下移時,由于水相的流量減少,必然導(dǎo)致流動下移時,由于水相的流量減少,必然導(dǎo)致流動時壓力值減少,所以在時壓力值減少,所以在含水率擬合好后,必須含水率擬合好后,必須再對壓力進(jìn)行必要的調(diào)整。再對壓力進(jìn)行必要的調(diào)整。 壓力、含水率、氣油比等都已經(jīng)擬合好后就可壓力、含水率、氣油比等都已經(jīng)擬合好后就可以進(jìn)行動態(tài)預(yù)測了。以進(jìn)行動態(tài)預(yù)測了。指標(biāo)
40、指標(biāo)實(shí)際值實(shí)際值計算值計算值相對誤差相對誤差累產(chǎn)油累產(chǎn)油(104t)20.92921.1591.1%累注水累注水(104m3)64.158564.90501.2%綜合含水綜合含水(%)93.091.122.0%采出程度采出程度(%)31.7131.830.4%地層壓力地層壓力(MPa) 21.08開發(fā)指標(biāo)區(qū)塊日產(chǎn)油020406080100120140051015202530時間(a )日產(chǎn)油(m 3 / d)實(shí)際值計算值0510152025051015202530時間(a )累產(chǎn)油(1 0 4 t)實(shí)際值計算值區(qū)塊累產(chǎn)油050100150200250300051015202530時間(a )日
41、注水(m 3 / d)實(shí)際值計算值區(qū)塊日注水區(qū)塊累注水010203040506070051015202530時間(a )累注水(1 0 4 m3)計算值實(shí)際值區(qū)塊含水率0102030405060708090100051015202530時間(a )含水率()實(shí)際值計算值01020304050607080901000.05.010.015.020.025.030.0時間(a )含水率()計算值實(shí)際值單井含水真6l數(shù)值模擬數(shù)據(jù)文件(前處理建模)l歷史擬合l結(jié)果分析(后處理)l動態(tài)預(yù)測l開發(fā)報告編寫每一步每一步的目的?的目的? 只有通過對擬合結(jié)果的分析,才能確定只有通過對擬合結(jié)果的分析,才能確定出需要進(jìn)行調(diào)整的方案。其分析的主要內(nèi)容出需要進(jìn)行調(diào)整的方案。其分析的主要內(nèi)容是,根據(jù)是,根據(jù)剩余油飽和度剩余油飽和度分布,分析剩余油的分布,分析剩余油的控制狀況及所在位置,由控制狀況及所在位置,由層間開發(fā)指標(biāo)層間開發(fā)指標(biāo)分析分析各層動用狀況,對層內(nèi)剩余油飽和度進(jìn)行研各層動用狀況,對層內(nèi)剩余油飽和度進(jìn)行研究,尋找水淹規(guī)律,由究,尋找水淹規(guī)律,由剩余儲量分布,水淹剩余儲量分布,水淹分布分布提供新的井點(diǎn),提出各種提高采收率方提供新的井點(diǎn),提出各種提
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 蘇教版數(shù)學(xué)三年級上冊單元測試卷-第五單元-解決問題的策略(含答案)-
- 2025年湖南吉利汽車職業(yè)技術(shù)學(xué)院單招職業(yè)技能測試題庫及答案一套
- 2025年河南物流職業(yè)學(xué)院單招職業(yè)適應(yīng)性測試題庫匯編
- 2025年河北青年管理干部學(xué)院單招職業(yè)傾向性測試題庫參考答案
- 2024年H-系列卷材涂料項(xiàng)目資金申請報告代可行性研究報告
- 2024年體育公園項(xiàng)目投資申請報告代可行性研究報告
- 2025年河北女子職業(yè)技術(shù)學(xué)院單招職業(yè)適應(yīng)性測試題庫一套
- 物理-云南省師范大學(xué)附屬中學(xué)2025屆高三下學(xué)期開學(xué)考試試題和答案
- 第22課《智取生辰綱》教學(xué)設(shè)計 2024-2025學(xué)年統(tǒng)編版語文九年級上冊
- 汽車旅館裝修合同解除
- 人教版八年級數(shù)學(xué)下冊《第十六章二次根式》專題復(fù)習(xí)附帶答案
- 2024屆武漢武昌區(qū)五校聯(lián)考數(shù)學(xué)九年級第一學(xué)期期末經(jīng)典試題含解析
- 高考復(fù)習(xí)概率中的遞推數(shù)列問題課件
- 生物工程設(shè)備課件
- 詐騙控告書模板
- 善借者贏天下(2017甘肅慶陽中考議論文閱讀試題含答案)
- 新聞采訪與寫作課件第十章采訪的實(shí)施現(xiàn)場觀察
- 八年級數(shù)學(xué)下冊《三角形的證明》單元測試卷(附答案解析)
- 國內(nèi)公務(wù)接待清單
- 《調(diào)整心態(tài)迎接中考》主題班會
- 領(lǐng)導(dǎo)科學(xué)與領(lǐng)導(dǎo)藝術(shù)
評論
0/150
提交評論