理工類畢業(yè)論文(論文)范例

1、 遼寧石油化工大學畢業(yè)設(shè)計 （論文） Graduation Project (Thesis) for Undergraduate of LSHU題題 目目TITLE學 院School專業(yè)班級Major&Class姓 名Name指導教師Supervisor20 年 月 日 論文獨創(chuàng)性聲明論文獨創(chuàng)性聲明本人所呈交的論文，是在指導教師指導下，獨立進行研究和開發(fā)工作所取得的成果。除文中已特別加以注明引用的內(nèi)容外，論文中不包含任何其他個人或集體已經(jīng)發(fā)表或撰寫過的研究成果。對本文的工作做出重要貢獻的個人和集體，均已在文中以明確方式標明并致謝。本聲明的法律結(jié)果由本人承擔。特此聲明。論文作者（簽名）：

2、 年 月 日遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙I摘 要針對稠油熱采過程中注入大量高溫蒸汽而造成井筒部分熱損失難以控制的問題，從熱工角度初步分析了高溫蒸汽在井筒部分的傳熱過程，建立了井筒及周圍地層部分的二維、瞬態(tài)傳熱數(shù)學模型，根據(jù)所求得的溫度場以及定量熱損失經(jīng)驗公式，探討了隔熱油管視導熱系數(shù)和蒸汽溫度對單位長度井筒熱損失的影響。研究表明：井筒單位長度熱損失隨著隔熱油管視導熱系數(shù)和蒸汽溫度的增加而增加，且隔熱油管視導熱系數(shù)對井筒單位長度熱損失的影響更大。高溫高壓蒸汽流經(jīng)井筒段后進入油層段，蒸汽通過射孔進入油層，并與油層進行傳質(zhì)傳熱，從而形成油層熱影響區(qū)域。本文利用有限元分析軟件 ANSYS，對遼

3、河某熱采井建立了井筒、地層和油層部分的二維、非穩(wěn)態(tài)數(shù)學模型，并模擬分析了油層和周圍地層溫度場，在此基礎(chǔ)之上探討了注汽時間對油層熱影響區(qū)域的影響。研究表明：油層熱影響區(qū)域先隨著注汽時間的增大而增大，而后基本不變，且注汽時間在 35 天最佳。關(guān)鍵詞：視導熱系數(shù)；注汽溫度；熱損失；熱影響區(qū)域；ANSYS遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙IIAbstractFor the purpose of reducing the heat loss in the mineshaft during the steam injection process, a two dimension, unsteady st

4、ate heat transfer model of the mineshaft and surrounding soil was established and then the temperature field was obtained. Base the Heat loss experience formula, the apparent heat conductivity and steam temperature to heat loss per length of the mineshaft influence were discussed. The results shown

5、that the heat loss per length of the mineshaft increased with the apparent heat conductivity and steam temperature increase, furthermore, the apparent heat conductivity affected larger than that of steam temperature.High temperature and high pressure steam flows through the wellbore segment backward

6、, such as the reservoir section, the steam through the perforation into the reservoir, and reservoir mass and heat transfer, and thus the formation of a reservoir heat-affected zone is build. Using the finite element analysis software ANSYS, a thermal recovery wellbore, formation and reservoir part

7、of the two non-steady-state mathematical model of Liaohe was build, and simulation Analysis of the temperature field of the reservoir and surrounding strata was made, and the relationgship between the time of steam injection and the reservoir heat-affected zone was discussed. The studys have shown t

8、hat: at first,reservoir heat-affected region increased with the increase of the steam injection time, and then kept unchang, the time of steam injection in 3 to 5 days was the best.Key words: Thermal conductivity; Steam temperature; Heat loss; Heat affected zone; ANSYS遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙I目錄1 1 緒論緒論.1

9、 11.1 目的意義 .11.2 國內(nèi)外的研究現(xiàn)狀 .21.2.1 井筒熱損失的研究 .21.2.2 油層熱影響區(qū)域的研究 .31.3 本文的研究內(nèi)容及創(chuàng)新點 .41.3.1 研究內(nèi)容 .41.3.2 創(chuàng)新點 .52 2 研究方法研究方法.6 62.1 有限元分析法簡介 .62.2 ANSYS 簡介.62.2.1 ANSYS 發(fā)展過程.72.2.2 ANSYS 技術(shù)特點.72.2.3 ANSYS 使用環(huán)境.82.2.4 ANSYS 程序功能.82.2.5 ANSYS 熱分析概述.92.2.6 ANSYS 求解步驟.93 3 井筒熱損失井筒熱損失.11113.1 問題描述 .113.2 數(shù)學模型

10、的建立 .113.2.1 能量平衡方程 .113.2.2 連續(xù)介質(zhì)的守恒方程 .13遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙II3.2.3 傳熱過程分析 .143.2.4 基本假設(shè) .163.2.5 控制方程 .163.2.6 邊界條件 .163.3 ANSYS 模擬.163.3.1 建立有限元模型 .163.3.2 施加載荷計算 .173.3.3 求解 .183.3.4 后處理 .183.4 模擬結(jié)果與分析 .183.4.1 隔熱油管視導熱系數(shù)對井筒熱損失的影響 .183.4.2 注汽溫度對單位長度井筒熱損失的影響 .204 4 油層熱影響區(qū)域油層熱影響區(qū)域.22224.1 問題描述 .224.

11、2 數(shù)學模型的建立 .224.2.1 連續(xù)介質(zhì)的守恒方程 .234.2.2 多孔介質(zhì)的守恒方程 .254.2.3 傳熱過程分析 .264.2.4 基本假設(shè) .274.2.5 邊界條件 .274.3 ANSYS 模擬.274.3.1 建立有限元模型 .274.3.2 施加載荷計算 .284.3.3 求解 .29遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙III4.3.4 后處理 .294.4 模擬結(jié)果與分析 .295 5 結(jié)論結(jié)論.3232參考文獻參考文獻.3333謝辭謝辭.3535遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙11 緒論1.1 目的意義稠油在世界石油資源中儲量豐富，據(jù)統(tǒng)計，其總量約為 10001

12、08噸，中國是世界上稠油資源豐富的國家之一，預(yù)測資源量約為 198 億噸。稠油由于粘度大，流動性差等特點，常規(guī)開采方法并不適用，基于稠油對溫度的敏感性特點，熱力采油已成為國內(nèi)外普遍認同的開采方法，即向地下油層內(nèi)注入高溫高壓蒸汽，一方面可以提高油層溫度、降低稠油粘度；另一方面可使油層壓力升高、增加驅(qū)油能力。常見的熱力采油方法有“蒸汽吞吐采油”和“蒸汽驅(qū)采油”兩種方式。 “蒸汽吞吐采油”是先向油井注入一定量的蒸汽，燜井一段時間，待蒸汽的熱能向油層擴散后，再開井生產(chǎn)的一種開采稠油的增產(chǎn)方法。蒸汽吞吐作業(yè)的過程分為三個階段，即注汽、燜井及回采?！罢羝?qū)采油”是稠油油藏經(jīng)過蒸汽吞吐采油之后，為進一步提高

13、采收率而采取的一項熱采方法，因為蒸汽吞吐采油只能采油各個油井附近油層中的原油，在油井與油井之間還留有大量的死油區(qū)。 “蒸汽驅(qū)采油”是由注入井連續(xù)不斷地往油層中注入高干度的蒸汽，蒸汽不斷地加熱油層，從而大大降低了地層原油的粘度。注入的蒸汽在地層中變?yōu)闊岬牧黧w，將原油驅(qū)趕到生產(chǎn)井的周圍，并被采到地面上來。在注蒸汽的過程中，由于蒸汽與地層之間存在溫差，蒸汽流動過程中又有摩擦阻力，因此注汽過程中存在著能量損失，注汽過程中的能量損失，特別是井筒中的能量損失，直接影響熱采效果。為了有效地減少這部分能量損失，提高熱采效果，需要對注汽井筒進行傳熱分析，進而探討隔熱油管視導熱系數(shù)和蒸汽溫度對單位長度井筒熱損失的

14、影響，這為熱采過程中熱量的高效、合理利用提供了一定的理論依據(jù)。蒸汽流經(jīng)井筒后進入目的層即油層，高溫高壓的蒸汽會通過射孔與油層進行傳質(zhì)傳熱，將熱量釋放給原油，使其溫度升高，粘度降低，同時增加油層壓力，達到驅(qū)油效果。但隨著蒸汽與油層熱量交換的進行，蒸汽自身溫度降低，逐漸冷凝成液體，阻止熱量繼續(xù)擴散，會形成熱影響區(qū)域。研究油層熱影響區(qū)域?qū)侠泶_定注汽時間和鉆井位置具有重要指導意遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙2義。1.2 國內(nèi)外的研究現(xiàn)狀1.2.1 井筒熱損失的研究自本世紀初以來，就有許多研究者致力于井筒熱損失的研究。井筒傳熱研究工作的先驅(qū) Ramey 以能量平衡為基礎(chǔ)，忽略了流體流動過程中動能

15、和摩擦的影響，最先提出了考慮時間和井深影響時的井筒液體溫度計算的理論模型。該模型作為井筒熱損研究的雛型，給出了井筒內(nèi)注入氣體和流體兩種工質(zhì)關(guān)于井深和時間的溫度解析表達式，為后來的井筒熱損研究奠定了基礎(chǔ)。但該模型中的時間函數(shù)要求注入流體的時間不得小于 7 天，所以此模型的精度在時間小于 7 天的時候受到限制。另外該模型僅適用于單相流，并在建立模型過程中假設(shè)井筒流體的物性參數(shù)不隨溫度變化，這與實際物理過程不相符1。國內(nèi)外許多研究者在 Ramey 模型的基礎(chǔ)上繼續(xù)做了多方面的研究，如 Hasan 和 Kabir 建立起井筒內(nèi)兩相流模型，改善了 Ramey 模型僅適用于單相流的不足，此模型可用于計算井

16、筒液體溫度的分布26。Pacheco 和 Farouq Ali 以流體流動的能量平衡方程為基礎(chǔ)，建立了能同時計算蒸汽壓力和蒸汽干度的數(shù)學模型7。該模型彌補了 Ramey 模型中沒有考慮水泥殼的影響的缺點。另外此模型結(jié)構(gòu)形式較復雜。但它可以用來較精確地計算蒸汽物性和熱損失隨深度和時間的變化。不過此模型還用到了 Ramey 模型中不完善的理論，如地層時間函數(shù) f ( )的計算，所以適用精度也受到一定限制。Satter 在使用 Ramey 模型時，考慮了注汽過程中的蒸汽相變影響，進而提出了一種預(yù)測冷凝流體干度隨深度及時間變化的計算方法，模型中提出的關(guān)系式8可用于計算飽和蒸汽注入時的井筒熱損失。Squ

17、ier 等人對熱水通過井筒注入時提出了一個完整的分析方法9；Huygen 和 Huitt10給出了在蒸汽注入期間的井筒熱損失理論和實驗結(jié)果，同時指出輻射熱損失的重要性；Willhite 提出了完整的計算總傳熱系數(shù)的方法11。Holst 和 Flock 以及 Earlougher 總結(jié)前面的一些成果，并考慮了蒸汽壓力的計算12,13；Eickmeier 等人采用有限差分模型計算注入初期過渡過程特性14。Hasan 等人又對井筒壓力降和熱損失計算進行了完善。但這些模型對于目前井筒注汽熱損失的計算尚有一些不足，遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙3首先它們在計算過程中，把某些參數(shù)看作是不隨溫度變化的

18、常量，這與實際并不完全相符；另外由于隔熱管的發(fā)展，上述模型中的總傳熱系數(shù)不能較真實地反映井筒現(xiàn)在結(jié)構(gòu)，給計算精度的提高帶來困難。我國對稠油熱采的相關(guān)研究比較晚。胡智勉、王彌康對注汽井筒總傳熱系數(shù)的計算方法提出了一些見解15,16。王彌康在 Ramey 模型的基礎(chǔ)上，從傳熱機理出發(fā)，給出了蒸汽井井筒熱傳遞定量計算的數(shù)學模型，該模型避免了 Ramey 公式僅適用于長時間注汽的不足，擴大了計算時間的適用范圍，從而改善了井筒傳熱計算的精度17。沈惠坊將地面管線與井筒進行了綜合考慮，探討了計算蒸汽壓力降、干度降及熱量損失的計算模型18。李景勤、陳艷華、王志國等人建立了計算井筒熱損數(shù)學模型，并提出了熱損合

19、理分布的概念1921。倪學鋒提出了井筒內(nèi)參數(shù)計算模型22。劉文章提出了用物理模擬方法確定井筒總傳熱系數(shù)的方法23。1.2.2 油層熱影響區(qū)域的研究1968 年 Spillette Nielsen 通過建立一個二維模型，對油層注冷水和注熱水后溫度分布進行了模擬計算，把流場分布和溫度分布藕合在了一起進行迭代24。1970 年埃克米爾等人在研究井筒換熱和溫度分布時，建立了一套關(guān)于注液和產(chǎn)液期間，液體與井筒周圍地層間熱交換的有限差分模型。計算中將油管、環(huán)空液、套管、水泥環(huán)及地層情況全部考慮在內(nèi)。該模型可以有效地模擬實際井筒條件，計算短時間和長時間注液、產(chǎn)液時井筒內(nèi)溫度的變化，模型只是基于熱傳導的模式進

20、行分析的，沒有考慮注入層的對流傳熱和井筒中的對流傳熱。1970 年，Smith，R.C.，Steffensen，R.J.等人提出了利用溫度測井的計算機分析估計流動剖面25。他們通過改變影啊溫度特性的參數(shù)，對注水井的溫度測井進行了詳細的研究，研究結(jié)果表明溫度測井對地層的熱學特性以及井眼的幾何形狀具有敏感性:另外，指出了注入水的溫度和小流量的漏失對溫度測井的影響，并對溫度測井解釋提出了指導性意見。1972 年維脫霍特與梯克薩通過改變注入速率、注入深度和注入時間，來研究注入期間井下的溫度分布，并對注入層中的溫度分布進行了研究。1973 年 Steffenson.RJ，SmlthR.C.等人遼寧石油化

21、工大學畢業(yè)設(shè)計（論文）用紙4提出在油層中，由于流動的存在，勢必發(fā)生熱交換，要計算井眼的溫度必須確定井筒周圍油藏的大片區(qū)域內(nèi)的溫度分布，并提出了著名的 Joul 一 Thomson 效應(yīng)。1982 年 John Fagley，H.SeottFogler，C.BrentDavenport，R.S.Millhone 等人對關(guān)井前后井筒的熱傳遞進行了數(shù)值模擬26，提出關(guān)井前一段時間的熱注是一種解釋老井注入剖面的潛在的重要工具。1988 年袁益壯針對在油藏中水驅(qū)油的流動和傳熱問題提出了多孔介質(zhì)中完全可壓縮、可混溶驅(qū)動問題的差分方法。1994 年 FredericMaubeuge，MichelP.Dide

22、k，MartynB.Beardsen 等人提出了一個多層注入模型27。1995 年胡耀江等人建立了水飽和的一維多孔介質(zhì)在一端注入蒸汽后，熱液區(qū)溫度分布的數(shù)學模型，提出了運用拉氏變換及數(shù)值反演求解此問題的計算方法28。1996 年馮恩民等人在建立的溫度模型中，除了考慮熱傳導外，特別強調(diào)了水井中與油層內(nèi)傳質(zhì)傳熱，構(gòu)造了關(guān)于導熱系數(shù)、比熱等參變量的優(yōu)化識別模型29。還有部分學者利用解析法求解地層中的溫度分布。1.3 本文的研究內(nèi)容及創(chuàng)新點1.3.1 研究內(nèi)容本文研究內(nèi)容主要分為兩個部分，即井筒熱損失和油層熱影響區(qū)域。A井筒熱損失利用傳熱學、熱力學、流體力學等知識分析隔熱油管內(nèi)管至地層的傳熱過程，建立

23、連續(xù)介質(zhì)的守恒方程，并建立數(shù)學模型；利用有限元軟件 ANSYS 建立隔熱油管內(nèi)管至地層的有限元模型，并利用其熱分析模塊對模型進行熱分析，求得注汽過程中井筒的溫度場分布云圖；探討注汽溫度和油管視導熱系數(shù)對井筒熱損失的影響。B油層熱影響區(qū)域利用傳熱學、熱力學、流體力學等知識分析套管至油層的傳熱過程，建立連續(xù)介質(zhì)的守恒方程和多孔介質(zhì)的守恒方恒，并建立數(shù)學模型；利用有限元軟件 ANSYS 建立套管至油層的有限元模型，并利用其熱分析模塊對模遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙5型進行熱分析，求得不同時刻油層的溫度場分布云圖；探討注汽時間與油層熱影響區(qū)域之間的關(guān)系。1.3.2 創(chuàng)新點（1）在以往的井筒溫度

24、場及熱損失研究中，以數(shù)值解法為主，本文利用有限元軟件ANSYS 的熱分析模塊求得了井筒溫度場分布云圖，并結(jié)合熱損失計算經(jīng)驗公式，探討了注汽溫度和油管視導熱系數(shù)對熱損失的影響。與數(shù)值法相比，本文研究方法可以清楚、直觀地觀察到溫度場的分布情況。（2）利用有限元軟件 ANSYS 的熱分析模塊模擬分析了油層熱影響區(qū)域，并探討了注汽時間與油層熱影響區(qū)域之間的關(guān)系。遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙62 研究方法2.1 有限元分析法簡介有限元法（finite element method）是一種高效能、常用的計算方法。有限元法在早期是以變分原理為基礎(chǔ)發(fā)展起來的，所以它廣泛地應(yīng)用于以拉普拉斯方程和泊松方程

25、所描述的各類物理場中（這類場與泛函的極值問題有著緊密的聯(lián)系） 。自從 1969 年以來，某些學者在流體力學中應(yīng)用加權(quán)余數(shù)法中的迦遼金法（Galerkin）或最小二乘法等同樣獲得了有限元方程，因而有限元法可應(yīng)用于以任何微分方程所描述的各類物理場中，而不再要求這類物理場和泛函的極值問題有所聯(lián)系。其原理是將連續(xù)的求解域離散為一組單元的組合體，用在每個單元內(nèi)假設(shè)的近似函數(shù)來分片的表示求解域上待求的未知場函數(shù)，近似函數(shù)通常由未知場函數(shù)及其導數(shù)在單元各節(jié)點的數(shù)值插值函數(shù)來表達。從而使一個連續(xù)的無限自由度問題變成離散的有限自由度問題。有限元法作為目前工程中應(yīng)用較為廣泛的一種數(shù)值計算方法，以其獨有的計算優(yōu)勢得

26、到了廣泛的發(fā)展和應(yīng)用，ANSYS 軟件以其多物理場耦合分析功能而成為 CAE（工程設(shè)計中的計算機輔助工程）軟件的應(yīng)用主流，在熱分析工程應(yīng)用中得到了較為廣泛的應(yīng)用。2.2 ANSYS 簡介ANSYS 程序是融結(jié)構(gòu)、熱、流體、電磁、聲學于一體的大型通用有限元商用分析軟件，可廣泛應(yīng)用于核工業(yè)、鐵道、石油化工、航空航天、機械制造、能源、電子、造船、汽車交通、國防軍工、土木工程、生物醫(yī)學、輕工、地礦、水利、日用家電等各個領(lǐng)域。該軟件可在大多數(shù)計算機及操作系統(tǒng)中運行，從 PC 到工作站直至巨型計算機，ANSYS 文件在其所有的產(chǎn)品系列和工作平臺上均兼容；該軟件基于 Motif 的菜單系統(tǒng)使用戶能夠通過對話

27、框、下拉式菜單和子菜單進行數(shù)據(jù)輸入和功能選擇，此舉大大方便了用戶操作。它由世界上著名的有限元分析軟件公司美國 ANSYS 公司開發(fā)的，能與多數(shù) CAD 軟件配合使用，實現(xiàn)數(shù)據(jù)的共享和交換，是現(xiàn)代產(chǎn)品設(shè)計中的高級 CAD 工具之一。在 ANSYS 公司相繼收遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙7購 ICEM、CENTURY DYNAMICS、AAVID THERMAL、FLUENT 等世界著名有限元分析程序制造公司并將產(chǎn)品整合之后，ANSYS 實際上已成為世界上最通用和有效的商用有限元軟件。2.2.1 ANSYS 發(fā)展過程ANSYS 公司是由美國匹茲堡大學的 John Swanson 博士在 1

28、970 年創(chuàng)建的，其總部位于美國賓夕法尼亞州的匹茲堡，目前是世界 CAE 行業(yè)最大的公司。自創(chuàng)建以來，ANSYS 公司一直致力于設(shè)計分析軟件的開發(fā)，不斷吸取新的計算方法和計算技術(shù)，領(lǐng)導著世界有限元技術(shù)的發(fā)展。ANSYS 程序的最初版本與本文用的版本 ANSYS 10.0 相比有很大區(qū)別。最初版本僅僅提供了熱分析和線性分析功能，是一個批處理程序，而且只能在大型計算機上使用。20 世紀 70 年代初，隨著非線性、子結(jié)構(gòu)以及更多的單元類型的加入。ANSYS 程序發(fā)生了很大的變化，新技術(shù)的融入進一步滿足了用戶的需求。20 世紀 70 年代末。交互方式的加入是該軟件最為顯著的變化，此舉使得模型生成和結(jié)果

29、評價大為簡化。2.2.2 ANSYS 技術(shù)特點與其他有限元計算軟件相比，ANSYS 具有以下技術(shù)特征： 能實現(xiàn)多場及多場耦合功能。 集前后處理、分析求解及多場分析于一體。 獨一無二的優(yōu)化功能，唯一具有流場優(yōu)化功能的 CFD 軟件。 具有強大的非線性分析功能。 具備快速求解器。 最早采用并行計算技術(shù)的 FEA 軟件。 支持從微機、工作站到巨型機的所有硬件平臺。 從微機、工作站、大型機直至巨型機，與所有硬件平臺的全部數(shù)據(jù)文件兼容。 從微機、工作站、大型機直至巨型機，所有硬件平臺具有統(tǒng)一用戶界面。遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙8 可與大多數(shù)的 CAD 軟件集成并具備接口。 可進行智能網(wǎng)格劃分。

30、 具有多層次多框架的產(chǎn)品系列。 具備良好的用戶開發(fā)環(huán)境。2.2.3 ANSYS 使用環(huán)境ANSYS 程序是一個功能強大的有限元計算分析軟件包。它可運行于 PC、NT 工作站、UNIX 工作站以及巨型計算機等各類計算機及操作系統(tǒng)中，其數(shù)據(jù)文件在其所有的產(chǎn)品系列和工作平臺上均兼容。該軟件有多種不同版本，本文應(yīng)用版本為 ANSYS 10.0 版，其微機版本要求的操作系統(tǒng)為 Windows95/98/2000/XP，也可運行與 UNIX 系統(tǒng)中；微機版的基本硬件要求為：顯示分辨率為 1024768 像素，顯示內(nèi)存為 2MB 以上，硬盤大于 350MB，推薦使用 17 英寸顯示器。ANSYS 多物理場耦

31、合的功能，允許在同一模型上進行各種耦合計算，例如：熱-結(jié)構(gòu)耦合、熱-電耦合、磁-結(jié)構(gòu)耦合以及熱-電-磁-流體耦合，同時在 PC 上生成的模型可運行于工作站及巨型計算機上，所有這一切就保證了 ANSYS 用戶對多領(lǐng)域工程問題的求解。ANSYS 可與多種先進的 CAD（如 AutoCAD、Pro/Engineer、NASTRAN、Alogor、I-DEAS 等）軟件共享數(shù)據(jù)，利用 ANSYS 的數(shù)據(jù)接口，可以精確地將在 CAD 系統(tǒng)下生成的幾何模型數(shù)據(jù)傳輸?shù)?ANSYS，并通過必要的修補可準確地在模型上劃分網(wǎng)格并進行求解，這樣就可以節(jié)省用戶在創(chuàng)建模型的過程中所花費的大量時間，使用戶的工作效率大幅度

32、提高。2.2.4 ANSYS 程序功能ANSYS 程序主要包括 3 個部分30：前處理模塊、求解模塊和后處理模塊。前處理模塊提供了一個強大的實體建模及網(wǎng)格劃分工具，用戶可以方便地構(gòu)造有限元模型；求解模塊包括結(jié)構(gòu)分析（結(jié)構(gòu)線性分析、結(jié)構(gòu)非線性分析和結(jié)構(gòu)高度非線性分析） 、熱分析、流體動力學分析、電磁場分析、聲場分析、壓電分析以及多物理場的耦合分析，可模擬多種物理遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙9介質(zhì)的相互作用，具有靈敏度分析及優(yōu)化分析能力；后處理模塊可將計算結(jié)果以彩色等值線顯示、梯度顯示、矢量顯示、粒子流顯示、立體切片顯示、透明及半透明顯示等圖形方式顯示出來，也可將計算結(jié)果以圖表、曲線形式顯

33、示或輸出。ANSYS 程序提供了 100 種以上的單元類型，用來模擬實際工程中的各種結(jié)構(gòu)和材料。2.2.5 ANSYS 熱分析概述在實際生產(chǎn)過程中，常常會遇到多種多樣的熱量傳遞問題：如計算某個系統(tǒng)或部件的溫度分布、熱量的獲取或損失、熱梯度、熱流密度、熱應(yīng)力、相變等。所涉及的部門包括：能源、化工、冶金、建筑、電子、航空航天、農(nóng)業(yè)、制冷、船舶等。以機械加工為例，往往需要估算和控制工件溫度場；分析不同工作條件下，不同材料及幾何形狀對溫度場變化的影響以及防止加工過程工件中缺陷的產(chǎn)生。ANSYS 作為新穎的有限元分析軟件在熱分析問題方面具有強大的功能，而且界面友好，易于掌握。用戶可以隨心所欲地選擇圖形用

34、戶界面方式（GUI）或命令流方式進行計算。ANSYS 的熱分析模塊一般包括ANSYS/Muliphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED 五種，其中 ANSYS/FLOTRAN 不含相變熱分析。ANSYS 熱分析基于能量守恒原理的熱平衡方程，用有限元法計算物體內(nèi)部各節(jié)點的溫度，并導出其他熱物理參數(shù)。運用ANSYS 軟件可進行熱傳導、熱對流、熱輻射、相變、熱應(yīng)力及接觸熱阻等問題的分析求解。2.2.6 ANSYS 求解步驟ANSYS 求解模型的過程中的主要步驟：定義單元類型及選項，根據(jù)所分析的問題在單元庫中選擇分析問題

35、的單元類型；定義單元是常數(shù)，單元是常數(shù)是與特定單元有關(guān)的量；定義材料屬性，在材料屬欄中選擇性的輸入熱導率、泊松比、密度等所需參數(shù)；建立幾何模型；遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙10定義網(wǎng)格并劃分參數(shù)；網(wǎng)格劃分以創(chuàng)建的模型。ANSYS 在后處理器中可以查看結(jié)果：顯示等值線；顯示分析結(jié)果數(shù)據(jù)；對結(jié)果數(shù)據(jù)進行計算和路徑操作；繪制云圖。結(jié)合本文研究內(nèi)容及 ANSYS 10.0 軟件的 ANSYS/Thermal 熱分析模塊功能，ANSYS 10.0 適用于本文研究內(nèi)容。遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙113 井筒熱損失3.1 問題描述以遼河油田某鉆井為研究對象，注汽時間 T=5d，油套環(huán)空導熱

36、系數(shù) yt為0.023W/mK，套管導熱系數(shù) tg為 43.2 W/mK，水泥環(huán)導熱系數(shù) cem為 0.81 W/mK，地層導熱系數(shù) cm為 1.745 W/mK，地層熱擴散系數(shù) c 為 0.027 m2/h，地層初始溫度 Tcm=40oC，蒸汽與隔熱油管內(nèi)壁對流換熱系數(shù) 為 10 kw/m2oC。3.2 數(shù)學模型的建立隔熱油管內(nèi)管內(nèi)半徑 r1為 0.031m，外半徑 r2為 0.0365m；隔熱油管外管內(nèi)半徑 r3為0.05015m，外半徑 r4為 0.05715m；套管內(nèi)半徑 r5為 0.0809m，外半徑 r6為 0.0889m；水泥環(huán)外半徑 r7為 0.1239m；地層半徑 r8為 0

37、.2478m，詳見圖 3-1 井筒及周圍地層截面示意圖。圖圖 3-1 井筒及其周圍地層截面示意圖井筒及其周圍地層截面示意圖3.2.1 能量平衡方程依據(jù)熱平衡原理，選取井筒內(nèi)的一微元流體單元體為研究對象，如圖 3-2 所示。r8r7r6r5r2r1r4r3蒸 汽隔熱層環(huán) 空 水泥環(huán)地層 遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙12圖圖 3-2 流體單元能量平衡模型流體單元能量平衡模型當蒸汽縱向流入或流出該單元體的控制面 1、2 時，攜帶的能量可表示為： （3-2 1i ).2(2，iiioigzvhmE1）式中：蒸汽攜帶的能量，W；iE 蒸汽流量，kg/s；om蒸汽焓值，J/kg ；ih蒸汽流速，m

38、/s ；iv 高度，m ；iz重力加速度，9.81m/s2。g該流體單元的能量平衡表達式為： （3-)2()2(1211022220zgvhmzgvhmQ2）E1E2123Q遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙13式中：單位時間的徑向散熱量，W。Q蒸汽的物性參數(shù)在井筒截面上是不變的，但不同截面處的物性參數(shù)是隨溫度變化的。由式（3-2）看出蒸汽沿井筒縱向的焓值在減少，蒸汽的能量是按遞減規(guī)律傳至地層的。3.2.2 連續(xù)介質(zhì)的守恒方程A井筒內(nèi)流體的能量平衡方程井筒油管內(nèi)高溫蒸汽在注入過程中遵循熱力學第一規(guī)律，其能量平衡方程為： （3-dzdzzdzzzzzQvgHvgzH2dz)-(z22-dz-z

39、-23）式中：、z 及 z-dz 處的焓，kJ；zHdzzH- vz、 vz-dz z 及 z-dz 處的速度，m/h；、z 及 z-dz 處的流體密度，kg/m3 ；zdz-zQ dz蒸汽散熱量，kJ 。公式（3-3）表示包括蒸汽焓、勢能和動能在 z 處的蒸汽能量與 z + dz 處的蒸汽能量及蒸汽散熱量的平衡方程。B井筒內(nèi)流體的連續(xù)性方程高溫蒸汽在油管內(nèi)流動時，質(zhì)量流量保持不變，因此滿足流體的連續(xù)性方程，即 （3-0mvAvAdzzdzzdzzzzz4）式中：、z 及 z+dz 處的橫截面積，m2；zAdzzAvz、 vz+dz z 及 z+dz 處的速度，m/h；、z 及 z+dz 處的

40、流體密度，kg/m3 。zdzzC井筒內(nèi)流體的動量平衡方程把蒸汽在油管內(nèi)的流動看作是兩相流的流動，為描述其流動特性，通常從兩相流壓力遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙14降的研究出發(fā)。管內(nèi)穩(wěn)態(tài)均質(zhì)兩相流滿足動量平衡方程，即 （3-vdvdzdvfgdzdPmmm-2-125）式中：段內(nèi)的壓力降，MPa；dPdz蒸汽密度 kg/m3；m摩阻系數(shù)；f井筒內(nèi)徑，m。1d表示段內(nèi)的勢能變化， 表示段內(nèi)的摩擦損失，表示段內(nèi)gdzmdzdzvdvmdz的動能變化。也就是說管內(nèi)兩相流的壓力降是摩擦損失、勢能變化和動能變化的綜合結(jié)果。D油管內(nèi)濕蒸汽的焓由于濕蒸汽為汽水兩相混合物，所以它的焓可用相同溫度或壓力

41、下的飽和水、飽和蒸汽及濕蒸汽的干度來表示，即 （3-)-1 (xhxhhlg6）式中：蒸汽干度；x濕蒸汽的焓，kJ/kg；h飽和蒸汽的焓，kJ/kg；gh飽和水的焓，kJ/kg。lh3.2.3 傳熱過程分析 根據(jù)熱力學31和傳熱學32知識對井筒傳熱進行分析，由于水泥環(huán)外側(cè)的地層體積和熱容量比水泥環(huán)內(nèi)側(cè)的井筒大得多，水泥環(huán)與地層交界處的溫度 隨注入時間的推移緩慢上kt升。表征熱量在地層中擴散速率的參數(shù)熱擴散系數(shù) 的數(shù)量級僅為 10-7 2m /s ，因而dzdvfm122遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙15與未受熱影響的遠處地層溫度的差值不大；相反，隔熱管內(nèi)蒸汽溫度與的差別較ktcitftk

42、t大，有一個數(shù)量級的差別。因而，可以假設(shè)水泥環(huán)內(nèi)側(cè)（簡稱井筒內(nèi)）熱量傳遞是穩(wěn)態(tài)的，而水泥層外側(cè)（簡稱地層內(nèi)）熱量傳遞則按非穩(wěn)態(tài)處理。A井筒內(nèi)穩(wěn)態(tài)傳熱井筒內(nèi)傳熱過程主要包括：蒸汽與隔熱油管內(nèi)壁之間的對流換熱、隔熱油管多孔隔熱材料中伴有熱輻射的導熱、油套管環(huán)空中對流、輻射同時作用的換熱以及套管、水泥環(huán)和地層的導熱等環(huán)節(jié)。對于注汽井一微元段 dz 的徑向熱損失，可表示為： （3-)-(22kfttdzkrdQ7）式中：段內(nèi)井筒的徑向熱損失，W；dQdz隔熱油管外半徑，m；2r以隔熱油管內(nèi)管外表面為基準面積的總傳熱系數(shù)，W/（m2k） 。k當井筒內(nèi)隔熱管抽真空時， 可表示為：k （3-17654321

43、111RRRRRRRk8）式中：蒸汽凝結(jié)換熱熱阻；1R隔熱油管熱阻；2R隔熱管內(nèi)隔熱層換熱熱阻；3R環(huán)空輻射換熱熱阻；4R環(huán)空對流熱熱阻；5R套管壁導熱熱阻；6R水泥環(huán)導熱熱阻；7R遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙16B地層內(nèi)非穩(wěn)態(tài)傳熱在井筒穩(wěn)態(tài)傳熱的同一微元段 dz 內(nèi)，地層非穩(wěn)態(tài)導熱方程如下： （3-eeetartrrt11229）初始條件時：0 （3-ZatttTbcie10）邊界條件時：0 （3-0rte11） （3-772-rrrtdzrdQee12）式中：地層內(nèi)某處溫度， oC；et， )(rfte地層導熱系數(shù)，W/（mK） ；e地層擴散系數(shù)，。/smc/2p，ee3.2.4

44、基本假設(shè)（1）從油管內(nèi)的蒸汽到地層的熱量傳遞為一維穩(wěn)態(tài)傳熱，地層為非穩(wěn)態(tài)導熱；（2）視油管內(nèi)的蒸汽流為一維均質(zhì)兩相流，且無相變及泄露現(xiàn)象；（3）忽略地層溫度梯度的影響。3.2.5 控制方程根據(jù)傳熱基本理論，二維瞬態(tài)傳熱問題的控制方程為： （3-TyTxT1222213）遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙17式中：x，y 分別指的是 x 和 y 方向； 是時間，s； 是熱擴散率，m2/s。3.2.6 邊界條件蒸汽與油管內(nèi)壁之間屬于對流換熱，所以油管內(nèi)壁為第三類邊界，應(yīng)滿足： （3-）（fTTanTk14）在不考慮地溫梯度的情況下，地層溫度為定值，因此地層邊界處為第一類邊界條件，其他邊界均為絕熱

45、邊界。3.3 ANSYS 模擬利用 ANSYS 10.0 熱分析軟件，做出在不同隔熱油管視導熱系數(shù)和注汽溫度時的井筒溫度場分布情況。3.3.1 建立有限元模型（1）打開 ANSYS 10.0 有限元分析軟件，過濾圖形用戶界面；（2）確定工作文件名“heat loss in well” ，單位制默認為國際單位制；（3）進入 PREP7 前處理器，選擇 plane55 熱分析單元，并設(shè)置成軸對稱單元；（4）將已知物性參數(shù)定義為材料熱物性參數(shù)；（5）根據(jù)圖 3-1 建立有限元模型，如圖 3-3；（6）采用自由網(wǎng)格劃分方法對模型進行網(wǎng)格劃分，結(jié)果如圖 3-4；3.3.2 施加載荷計算A定義分析類型定義

46、新的分析類型為瞬態(tài)分析“Transient” 。B設(shè)置瞬態(tài)熱分析的初始條件遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙18由于初始溫度場是不均勻的且又是未知的，所以首先要進行穩(wěn)態(tài)熱分析，為穩(wěn)態(tài)熱分析建立初始條件：在井筒壁處施加對流載荷，在地層處施加原始地層溫度場載荷，其他邊界為絕熱邊界。圖圖 3-3 井筒有限元模型井筒有限元模型圖圖 3-4 網(wǎng)格劃分結(jié)果網(wǎng)格劃分結(jié)果3.3.3 求解采用 GUI 命令：MainSolutionLoad Step OptsTime FrequencTime-Time Step，設(shè)置計算終止時間為 0.01s，時間步長為 0.01s，并將時間積分效果設(shè)定為 Off，對模型進

47、行穩(wěn)態(tài)熱分析。穩(wěn)態(tài)熱分析結(jié)束后，將時間積分效果設(shè)定為 On，并刪除所有節(jié)點的溫度載荷，計算終遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙19止時間設(shè)定為 5 天，時間步長設(shè)定為 1 小時，對模型進行瞬態(tài)熱分析。3.3.4 后處理對于瞬態(tài)熱分析，ANSYS 提供了兩種后處理方式，即 POST 1 和 POST 26。本文采用POST 1 進行后處理，利用 GUI 命令：Main MenuGeneral PostprocPlot Results 顯示溫度場彩色云圖，并拾取套管內(nèi)徑邊緣處溫度。 3.4 模擬結(jié)果與分析取水泥環(huán)與地層之間截面的平均溫度作為溫度參數(shù)，并結(jié)合文獻33中給出的熱損失經(jīng)驗公式，分別算出

48、在不同的隔熱油管視導熱系數(shù)及注汽溫度下井筒熱損失情況。3.4.1 隔熱油管視導熱系數(shù)對井筒熱損失的影響現(xiàn)在使用的隔熱油管多為抽真空隔熱油管，結(jié)合預(yù)應(yīng)力隔熱油管標準（SY/T 5324-94）中提出的視導熱系數(shù)概念，把隔熱管的隔熱層內(nèi)以導熱、對流和輻射三種方式傳遞的熱量視為以純導熱方式傳遞的熱量。其中，視導熱系數(shù)數(shù)值劃分為五個范圍等級，如表 3-1 所示，因此，分別選取視導熱系數(shù)為 0.07、0.05、0.03、0.01 和 0.005 W/mK。 表表 3-1 隔熱油管等級分類隔熱油管等級分類 單位：W/mK隔熱性能等級ABCDE視導熱系數(shù)0.060.080.040.060.020.040.0

49、060.020.0020.006在實際的稠油熱采過程中，注汽溫度一般為 300350之間，基于此，本文以蒸汽溫度為 312情況下進行模擬，圖 3-5 為在視導熱系數(shù)為 0.07 時的溫度分布云圖，結(jié)合熱損失經(jīng)驗公式，得到單位長度井筒熱損失與視導熱系數(shù)之間的關(guān)系，如表 3-2 和圖 3-6 所示。遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙20圖圖 3-5 視導熱系數(shù)為視導熱系數(shù)為 0.07 W/mK 時的溫度分布云圖時的溫度分布云圖表表 3-2 不同視導熱系數(shù)下的計算結(jié)果不同視導熱系數(shù)下的計算結(jié)果序號視導熱系數(shù)（W/mK）套管溫度（oC）熱損失（W/m）123450.0050.010.030.050.

50、0741.58842.84846.04647.79748.9022.7154.87010.33813.33215.222遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙21圖圖 3-6 隔熱油管視導熱系數(shù)對井筒熱損失的影響隔熱油管視導熱系數(shù)對井筒熱損失的影響根據(jù)圖 3-6 分析可知，井筒單位長度熱損失整體上隨著隔熱油管視導熱系數(shù)的增大而增大。其中，對于 A 級的隔熱油管，單位長度井筒熱損失約達 15W/m 左右，而對于 E 級的隔熱油管，單位長度井筒熱損失僅約為 3.5W/m，兩者單位長度井筒熱損失相差約11.5W/m，因此，選用視導熱系數(shù)高級別的隔熱油管對于降低井筒單位長度的熱損失具有明顯的影響。此外，

51、當視導熱系數(shù)從 0.07W/mK 提高至 0.03 W/mK 后，單位長度熱損失的變化速率隨著視導熱系數(shù)的減小而出現(xiàn)明顯的降低。綜上所述，在實際的應(yīng)用過程中，建議選用隔熱油管視導熱系數(shù)在 0.030.005W/mK。3.4.2 注汽溫度對單位長度井筒熱損失的影響上述視導熱系數(shù)的建議范圍內(nèi)，在視導熱系數(shù)為 0.03 W/mK 情況下，分別研究以注汽溫度分別為 312、324、333、346 和 350oC 五種情況下，注汽溫度對單位長度井筒熱損失的影響，詳細情況如表 3-3 和圖 3-7 所示，圖 3-8 為當注汽溫度為 324oC 時的溫度分布云圖。表表 3-3 不同注汽溫度下的計算結(jié)果不同注

52、汽溫度下的計算結(jié)果序號注汽溫度（oC）套管溫度（oC）熱損失（W/m）1231232448.90249.29715.22215.897遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙2234533334635049.59249.98750.08216.40117.07717.239根據(jù)上圖分析所知，單位長度井筒的熱損失隨著注汽溫度的升高幾乎呈線性增大的趨勢。在注汽溫度為 312oC 情況下，單位長度井筒熱損失約為 15.2W/m，在注汽溫度在350oC 情況下，單位長度井筒熱損失約為 17.2W/m，由此可見，注汽溫度的提高對于井筒單位長度熱損失的影響較小，每提高 50oC，井筒單位長度熱損失相差約 2W

53、/m。圖圖 3-7 注汽溫度對井筒熱損失的影響注汽溫度對井筒熱損失的影響遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙23圖圖 3-8 注汽溫度為注汽溫度為 324oC 時的溫度分布云圖時的溫度分布云圖遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙244 油層熱影響區(qū)域4.1 問題描述高溫高壓的蒸汽流經(jīng)井筒后進入封隔器以下的油層區(qū)，通過射孔與油層進行傳質(zhì)傳熱，傳熱過程如圖 4-1 所示。圖圖 4-1 油層傳熱過程示意圖油層傳熱過程示意圖以遼河油田某鉆井為研究對象，注汽時間 T=10d，套管導熱系數(shù) tg為 43.2 W/mK，水泥環(huán)導熱系數(shù) cem為 1.2 W/mK，油層綜合導熱系數(shù) cm為 2.53 W/mK

54、，地層導熱系數(shù)為1.968 W/mK，地層熱擴散系數(shù) c 為 0.027 m2/h，地層初始溫度 Tcm=40oC，蒸汽與隔熱油管內(nèi)壁對流換熱系數(shù) 為 10 kw/m2oC。4.2 數(shù)學模型的建立套管內(nèi)半徑 r1為 0.0809m，外半徑 r2為 0.0889m；水泥環(huán)外半徑 r3為 0.1239m；地層遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙25和油層半徑 r4為 20m，詳見圖 4-2 油層截面示意圖。蒸 汽水 泥 環(huán)地 層油 層套 管r1r2r30.5mm1m圖圖 4-2 油層截面示意圖油層截面示意圖4.2.1 連續(xù)介質(zhì)的守恒方程求解與流體中動量、熱量和質(zhì)量傳遞有關(guān)的問題時，首先應(yīng)根據(jù)基本的

55、守恒定律對該物理系統(tǒng)建立一組方程。這些方程包括:能量平衡方程；連續(xù)性方程；焓方程。A連續(xù)介質(zhì)的能量平衡方程井筒油管內(nèi)高溫蒸汽在注入過程中遵循熱力學第一規(guī)律，其能量平衡方程為： （4-dzdzzdzzzzzQvgHvgzH2dz)-(z22-dz-z-21）遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙26式中：Hz、Hz-dzz 及 z-dz 處的焓，kJ； vz、 vz-dz z 及 z-dz 處的速度，m/h；、z 及 z-dz 處的流體密度，kg/m3 ；zdz-zQ dz蒸汽散熱量，kJ 。上式表示包括蒸汽焓、勢能和動能在 z 處的蒸汽能量與 z + dz 處的蒸汽能量及蒸汽散熱量的平衡方程。B

56、連續(xù)介質(zhì)的連續(xù)性方程高溫蒸汽在油管內(nèi)流動時，質(zhì)量流量保持不變，因此滿足流體的連續(xù)性方程，即 （4-0mvAvAdzzdzzdzzzzz2）式中：、z 及 z+dz 處的橫截面積，m2；zAdzzAvz、 vz+dz z 及 z+dz 處的速度，m/h；、z 及 z+dz 處的流體密度，kg/m3 。zdzzC連續(xù)介質(zhì)的動量平衡方程把蒸汽在油管內(nèi)的流動看作是兩相流的流動，為描述其流動特性，通常從兩相流壓力降的研究出發(fā)。管內(nèi)穩(wěn)態(tài)均質(zhì)兩相流滿足動量平衡方程，即 （4-vdvdzdvfgdzdPmmm-2-123）式中：段內(nèi)的壓力降，MPa；dPdz蒸汽密度 kg/m3；m摩阻系數(shù)；f井筒內(nèi)徑，m。1

57、d表示段內(nèi)的勢能變化， 表示段內(nèi)的摩擦損失，表示段內(nèi)gdzmdzdzvdvmdzdzdvfm122遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙27的動能變化。D油管內(nèi)濕蒸汽的焓由于濕蒸汽為汽水兩相混合物，所以它的焓可用相同溫度或壓力下的飽和水、飽和蒸汽及濕蒸汽的干度來表示，即 （4-)-1 (xhxhhlg4）式中：蒸汽干度；x濕蒸汽的焓，kJ/kg；h飽和蒸汽的焓，kJ/kg；gh飽和水的焓，kJ/kg。lh4.2.2 多孔介質(zhì)的守恒方程A多孔介質(zhì)的質(zhì)量守恒方程多孔介質(zhì)的質(zhì)量守恒方程為： （4-01vt5）式中：密度，kg/m3；介質(zhì)的孔隙度；滲流速度，m/d。vB多孔介質(zhì)的動力學方程a. 單相流在

58、只有飽和牛頓流體的空隙介質(zhì)中，達西定律如下： （4-)-(-Dgpkv6）遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙28式中：流體的滲流速度，m/d；v滲透率，；k2m粘度，Pas；流體的密度，kg/m3；壓力，MPa；pD垂直高度，m；G重力加速度，m/s2。b. 多相流多相流的滲流微分方程為： (4-7)()-(11llllrlstqDgpkk式中:sll 相流體的飽和度；ql流體流入或流出的流體質(zhì)量流量，kg/s。4.2.3 傳熱過程分析A井筒內(nèi)穩(wěn)態(tài)傳熱井筒內(nèi)傳熱過程主要包括：蒸汽與套管內(nèi)壁之間的對流換熱、套管與水泥環(huán)和水泥環(huán)與油層的導熱等環(huán)節(jié)。對于注汽井一微元段 dz 的徑向熱流量，可表示為

59、： （4-)-(22kfttdzkrdQ8）式中：段內(nèi)井筒的徑向熱流量，W；dQdz套管外半徑，m；2r以套管外表面為基準面積的總傳熱系數(shù)，W/（m2k） 。k的表達式為：k遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙29 （4-1321RRRk9）式中：蒸汽與套管內(nèi)壁對流換熱熱阻；1R套管壁導熱熱阻；2R水泥環(huán)導熱熱阻。3RB油層傳熱分析 井筒蒸汽與套管通過對流換熱向井筒周圍散熱，同時產(chǎn)生從射孔到油層的傳質(zhì)熱傳導。由于在有限元中無法解決傳質(zhì)問題，因此一般簡化為具有一定半徑的加熱區(qū)模型。非穩(wěn)態(tài)無內(nèi)熱源軸對稱溫度場的微分方程為： （4-0c-),(p2222tTrTrTrxTrktrxTD10）應(yīng)用 G

60、alerkin 法可得到軸對稱溫度場整體區(qū)域“變分”計算的基本方程為: (4-11)TlDlplllPdSnTrkWdxdrtTrWcrTrWxTxWkrTJ0-利用溫度場變分方程對模型進行剖分離散、單元分析及最終總體合成，即可完成模型溫度場的求解。4.2.4 基本假設(shè)沿井深方向上的熱能不發(fā)生傳遞；視蒸汽流為一維均質(zhì)兩相流，且無相變；忽略油層溫度梯度的影響；以井眼中心為中心的同心圓上的溫度是相同的。4.2.5 邊界條件蒸汽與套管內(nèi)壁之間屬于對流換熱，所以套管內(nèi)壁為第三類邊界，應(yīng)滿足： （4-）（fTTanTk遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙3012）加熱區(qū)為第一類邊界條件，由于不受油層溫度梯度的影