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

PAGE遼寧石油化工大學畢業(yè)設(shè)計（論文）Graduation

Project(Thesis)forUndergraduateofLSHU題目TITLE學院School專業(yè)班級Major&Class姓名Name指導教師Supervisor20年月日論文獨創(chuàng)性聲明本人所呈交的論文，是在指導教師指導下，獨立進行研究和開發(fā)工作所取得的成果。除文中已特別加以注明引用的內(nèi)容外，論文中不包含任何其他個人或集體已經(jīng)發(fā)表或撰寫過的研究成果。對本文的工作做出重要貢獻的個人和集體，均已在文中以明確方式標明并致謝。本聲明的法律結(jié)果由本人承擔。特此聲明。論文作者（簽名）：年月日遼寧石油化工大學畢業(yè)設(shè)計（論文）用紙PAGEIII摘要針對稠油熱采過程中注入大量高溫蒸汽而造成井筒部分熱損失難以控制的問題，從熱工角度初步分析了高溫蒸汽在井筒部分的傳熱過程，建立了井筒及周圍地層部分的二維、瞬態(tài)傳熱數(shù)學模型，根據(jù)所求得的溫度場以及定量熱損失經(jīng)驗公式，探討了隔熱油管視導熱系數(shù)和蒸汽溫度對單位長度井筒熱損失的影響。研究表明：井筒單位長度熱損失隨著隔熱油管視導熱系數(shù)和蒸汽溫度的增加而增加，且隔熱油管視導熱系數(shù)對井筒單位長度熱損失的影響更大。高溫高壓蒸汽流經(jīng)井筒段后進入油層段，蒸汽通過射孔進入油層，并與油層進行傳質(zhì)傳熱，從而形成油層熱影響區(qū)域。本文利用有限元分析軟件ANSYS，對遼河某熱采井建立了井筒、地層和油層部分的二維、非穩(wěn)態(tài)數(shù)學模型，并模擬分析了油層和周圍地層溫度場，在此基礎(chǔ)之上探討了注汽時間對油層熱影響區(qū)域的影響。研究表明：油層熱影響區(qū)域先隨著注汽時間的增大而增大，而后基本不變，且注汽時間在3~5天最佳。關(guān)鍵詞：視導熱系數(shù)；注汽溫度；熱損失；熱影響區(qū)域；ANSYSAbstractForthepurposeofreducingtheheatlossinthemineshaftduringthesteaminjectionprocess,atwodimension,unsteadystateheattransfermodelofthemineshaftandsurroundingsoilwasestablishedandthenthetemperaturefieldwasobtained.BasetheHeatlossexperienceformula,theapparentheatconductivityandsteamtemperaturetoheatlossperlengthofthemineshaftinfluencewerediscussed.Theresultsshownthattheheatlossperlengthofthemineshaftincreasedwiththeapparentheatconductivityandsteamtemperatureincrease,furthermore,theapparentheatconductivityaffectedlargerthanthatofsteamtemperature.Hightemperatureandhighpressuresteamflowsthroughthewellboresegmentbackward,suchasthereservoirsection,thesteamthroughtheperforationintothereservoir,andreservoirmassandheattransfer,andthustheformationofareservoirheat-affectedzoneisbuild.UsingthefiniteelementanalysissoftwareANSYS,athermalrecoverywellbore,formationandreservoirpartofthetwonon-steady-statemathematicalmodelofLiaohewasbuild,andsimulationAnalysisofthetemperaturefieldofthereservoirandsurroundingstratawasmade,andtherelationshipbetweenthetimeofsteaminjectionandthereservoirheat-affectedzonewasdiscussed.Thestudieshaveshownthat:atfirst,reservoirheat-affectedregionincreasedwiththeincreaseofthesteaminjectiontime,andthenkeptinstall,thetimeofsteaminjectionin3to5dayswasthebest.Keywords:Thermalconductivity;Steamtemperature;Heatloss;Heataffectedzone;ANSYS目錄1緒論 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研究方法 62.1有限元分析法簡介 62.2ANSYS簡介 62.2.1ANSYS發(fā)展過程 72.2.2ANSYS技術(shù)特點 72.2.3ANSYS使用環(huán)境 82.2.4ANSYS程序功能 82.2.5ANSYS熱分析概述 92.2.6ANSYS求解步驟 93井筒熱損失 113.1問題描述 113.2數(shù)學模型的建立 113.2.1能量平衡方程 113.2.2連續(xù)介質(zhì)的守恒方程 133.2.3傳熱過程分析 143.2.4基本假設(shè) 163.2.5控制方程 163.2.6邊界條件 163.3ANSYS模擬 163.3.1建立有限元模型 163.3.2施加載荷計算 173.3.3求解 183.3.4后處理 183.4模擬結(jié)果與分析 183.4.1隔熱油管視導熱系數(shù)對井筒熱損失的影響 183.4.2注汽溫度對單位長度井筒熱損失的影響 204油層熱影響區(qū)域 224.1問題描述 224.2數(shù)學模型的建立 224.2.1連續(xù)介質(zhì)的守恒方程 234.2.2多孔介質(zhì)的守恒方程 254.2.3傳熱過程分析 264.2.4基本假設(shè) 274.2.5邊界條件 274.3ANSYS模擬 274.3.1建立有限元模型 274.3.2施加載荷計算 284.3.3求解 294.3.4后處理 294.4模擬結(jié)果與分析 295結(jié)論 32參考文獻 33謝辭 35 PAGE391緒論1.1目的意義稠油在世界石油資源中儲量豐富，據(jù)統(tǒng)計，其總量約為1000×108噸，中國是世界上稠油資源豐富的國家之一，預測資源量約為198億噸。稠油由于粘度大，流動性差等特點，常規(guī)開采方法并不適用，基于稠油對溫度的敏感性特點，熱力采油已成為國內(nèi)外普遍認同的開采方法，即向地下油層內(nèi)注入高溫高壓蒸汽，一方面可以提高油層溫度、降低稠油粘度；另一方面可使油層壓力升高、增加驅(qū)油能力。常見的熱力采油方法有“蒸汽吞吐采油”和“蒸汽驅(qū)采油”兩種方式?！罢羝掏虏捎汀笔窍认蛴途⑷胍欢康恼羝瑺F井一段時間，待蒸汽的熱能向油層擴散后，再開井生產(chǎn)的一種開采稠油的增產(chǎn)方法。蒸汽吞吐作業(yè)的過程分為三個階段，即注汽、燜井及回采?！罢羝?qū)采油”是稠油油藏經(jīng)過蒸汽吞吐采油之后，為進一步提高采收率而采取的一項熱采方法，因為蒸汽吞吐采油只能采油各個油井附近油層中的原油，在油井與油井之間還留有大量的死油區(qū)?！罢羝?qū)采油”是由注入井連續(xù)不斷地往油層中注入高干度的蒸汽，蒸汽不斷地加熱油層，從而大大降低了地層原油的粘度。注入的蒸汽在地層中變?yōu)闊岬牧黧w，將原油驅(qū)趕到生產(chǎn)井的周圍，并被采到地面上來。在注蒸汽的過程中，由于蒸汽與地層之間存在溫差，蒸汽流動過程中又有摩擦阻力，因此注汽過程中存在著能量損失，注汽過程中的能量損失，特別是井筒中的能量損失，直接影響熱采效果。為了有效地減少這部分能量損失，提高熱采效果，需要對注汽井筒進行傳熱分析，進而探討隔熱油管視導熱系數(shù)和蒸汽溫度對單位長度井筒熱損失的影響，這為熱采過程中熱量的高效、合理利用提供了一定的理論依據(jù)。蒸汽流經(jīng)井筒后進入目的層即油層，高溫高壓的蒸汽會通過射孔與油層進行傳質(zhì)傳熱，將熱量釋放給原油，使其溫度升高，粘度降低，同時增加油層壓力，達到驅(qū)油效果。但隨著蒸汽與油層熱量交換的進行，蒸汽自身溫度降低，逐漸冷凝成液體，阻止熱量繼續(xù)擴散，會形成熱影響區(qū)域。研究油層熱影響區(qū)域?qū)侠泶_定注汽時間和鉆井位置具有重要指導意義。1.2國內(nèi)外的研究現(xiàn)狀1.2.1井筒熱損失的研究自本世紀初以來，就有許多研究者致力于井筒熱損失的研究。井筒傳熱研究工作的先驅(qū)Ramey以能量平衡為基礎(chǔ)，忽略了流體流動過程中動能和摩擦的影響，最先提出了考慮時間和井深影響時的井筒液體溫度計算的理論模型。該模型作為井筒熱損研究的雛型，給出了井筒內(nèi)注入氣體和流體兩種工質(zhì)關(guān)于井深和時間的溫度解析表達式，為后來的井筒熱損研究奠定了基礎(chǔ)。但該模型中的時間函數(shù)要求注入流體的時間不得小于7天，所以此模型的精度在時間小于7天的時候受到限制。另外該模型僅適用于單相流，并在建立模型過程中假設(shè)井筒流體的物性參數(shù)不隨溫度變化，這與實際物理過程不相符[1]。國內(nèi)外許多研究者在Ramey模型的基礎(chǔ)上繼續(xù)做了多方面的研究，如Hasan和Kabir建立起井筒內(nèi)兩相流模型，改善了Ramey模型僅適用于單相流的不足，此模型可用于計算井筒液體溫度的分布[2~6]。Pacheco和FarouqAli以流體流動的能量平衡方程為基礎(chǔ)，建立了能同時計算蒸汽壓力和蒸汽干度的數(shù)學模型[7]。該模型彌補了Ramey模型中沒有考慮水泥殼的影響的缺點。另外此模型結(jié)構(gòu)形式較復雜。但它可以用來較精確地計算蒸汽物性和熱損失隨深度和時間的變化。不過此模型還用到了Ramey模型中不完善的理論，如地層時間函數(shù)f(τ)的計算，所以適用精度也受到一定限制。Satter在使用Ramey模型時，考慮了注汽過程中的蒸汽相變影響，進而提出了一種預測冷凝流體干度隨深度及時間變化的計算方法，模型中提出的關(guān)系式[8]可用于計算飽和蒸汽注入時的井筒熱損失。Squier等人對熱水通過井筒注入時提出了一個完整的分析方法[9]；Huygen和Huitt[10]給出了在蒸汽注入期間的井筒熱損失理論和實驗結(jié)果，同時指出輻射熱損失的重要性；Willhite提出了完整的計算總傳熱系數(shù)的方法[11]。Holst和Flock以及Earlougher總結(jié)前面的一些成果，并考慮了蒸汽壓力的計算[12,13]；Eickmeier等人采用有限差分模型計算注入初期過渡過程特性[14]。Hasan等人又對井筒壓力降和熱損失計算進行了完善。但這些模型對于目前井筒注汽熱損失的計算尚有一些不足，首先它們在計算過程中，把某些參數(shù)看作是不隨溫度變化的常量，這與實際并不完全相符；另外由于隔熱管的發(fā)展，上述模型中的總傳熱系數(shù)不能較真實地反映井筒現(xiàn)在結(jié)構(gòu)，給計算精度的提高帶來困難。我國對稠油熱采的相關(guān)研究比較晚。胡智勉、王彌康對注汽井筒總傳熱系數(shù)的計算方法提出了一些見解[15,16]。王彌康在Ramey模型的基礎(chǔ)上，從傳熱機理出發(fā)，給出了蒸汽井井筒熱傳遞定量計算的數(shù)學模型，該模型避免了Ramey公式僅適用于長時間注汽的不足，擴大了計算時間的適用范圍，從而改善了井筒傳熱計算的精度[17]。沈惠坊將地面管線與井筒進行了綜合考慮，探討了計算蒸汽壓力降、干度降及熱量損失的計算模型[18]。李景勤、陳艷華、王志國等人建立了計算井筒熱損數(shù)學模型，并提出了熱損合理分布的概念[19~21]。倪學鋒提出了井筒內(nèi)參數(shù)計算模型[22]。劉文章提出了用物理模擬方法確定井筒總傳熱系數(shù)的方法[23]。1.2.2油層熱影響區(qū)域的研究1968年SpilletteNielsen通過建立一個二維模型，對油層注冷水和注熱水后溫度分布進行了模擬計算，把流場分布和溫度分布藕合在了一起進行迭代[24]。1970年埃克米爾等人在研究井筒換熱和溫度分布時，建立了一套關(guān)于注液和產(chǎn)液期間，液體與井筒周圍地層間熱交換的有限差分模型。計算中將油管、環(huán)空液、套管、水泥環(huán)及地層情況全部考慮在內(nèi)。該模型可以有效地模擬實際井筒條件，計算短時間和長時間注液、產(chǎn)液時井筒內(nèi)溫度的變化，模型只是基于熱傳導的模式進行分析的，沒有考慮注入層的對流傳熱和井筒中的對流傳熱。1970年，Smith，R.C.，Steffensen，R.J.等人提出了利用溫度測井的計算機分析估計流動剖面[25]。他們通過改變影啊溫度特性的參數(shù)，對注水井的溫度測井進行了詳細的研究，研究結(jié)果表明溫度測井對地層的熱學特性以及井眼的幾何形狀具有敏感性:另外，指出了注入水的溫度和小流量的漏失對溫度測井的影響，并對溫度測井解釋提出了指導性意見。1972年維脫霍特與梯克薩通過改變注入速率、注入深度和注入時間，來研究注入期間井下的溫度分布，并對注入層中的溫度分布進行了研究。1973年Steffenson.RJ，SmlthR.C.等人提出在油層中，由于流動的存在，勢必發(fā)生熱交換，要計算井眼的溫度必須確定井筒周圍油藏的大片區(qū)域內(nèi)的溫度分布，并提出了著名的Joul一Thomson效應。1982年JohnFagley，H.SeottFogler，C.BrentDavenport，R.S.Millhone等人對關(guān)井前后井筒的熱傳遞進行了數(shù)值模擬[26]，提出關(guān)井前一段時間的熱注是一種解釋老井注入剖面的潛在的重要工具。1988年袁益壯針對在油藏中水驅(qū)油的流動和傳熱問題提出了多孔介質(zhì)中完全可壓縮、可混溶驅(qū)動問題的差分方法。1994年FredericMaubeuge，MichelP.Didek，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.井筒熱損失=1\*GB3①利用傳熱學、熱力學、流體力學等知識分析隔熱油管內(nèi)管至地層的傳熱過程，建立連續(xù)介質(zhì)的守恒方程，并建立數(shù)學模型；=2\*GB3②利用有限元軟件ANSYS建立隔熱油管內(nèi)管至地層的有限元模型，并利用其熱分析模塊對模型進行熱分析，求得注汽過程中井筒的溫度場分布云圖；=3\*GB3③探討注汽溫度和油管視導熱系數(shù)對井筒熱損失的影響。B.油層熱影響區(qū)域=1\*GB3①利用傳熱學、熱力學、流體力學等知識分析套管至油層的傳熱過程，建立連續(xù)介質(zhì)的守恒方程和多孔介質(zhì)的守恒方恒，并建立數(shù)學模型；=2\*GB3②利用有限元軟件ANSYS建立套管至油層的有限元模型，并利用其熱分析模塊對模型進行熱分析，求得不同時刻油層的溫度場分布云圖；=3\*GB3③探討注汽時間與油層熱影響區(qū)域之間的關(guān)系。1.3.2創(chuàng)新點（1）在以往的井筒溫度場及熱損失研究中，以數(shù)值解法為主，本文利用有限元軟件ANSYS的熱分析模塊求得了井筒溫度場分布云圖，并結(jié)合熱損失計算經(jīng)驗公式，探討了注汽溫度和油管視導熱系數(shù)對熱損失的影響。與數(shù)值法相比，本文研究方法可以清楚、直觀地觀察到溫度場的分布情況。（2）利用有限元軟件ANSYS的熱分析模塊模擬分析了油層熱影響區(qū)域，并探討了注汽時間與油層熱影響區(qū)域之間的關(guān)系。2研究方法2.1有限元分析法簡介有限元法（finiteelementmethod）是一種高效能、常用的計算方法。有限元法在早期是以變分原理為基礎(chǔ)發(fā)展起來的，所以它廣泛地應用于以拉普拉斯方程和泊松方程所描述的各類物理場中（這類場與泛函的極值問題有著緊密的聯(lián)系）。自從1969年以來，某些學者在流體力學中應用加權(quán)余數(shù)法中的迦遼金法（Galerkin）或最小二乘法等同樣獲得了有限元方程，因而有限元法可應用于以任何微分方程所描述的各類物理場中，而不再要求這類物理場和泛函的極值問題有所聯(lián)系。其原理是將連續(xù)的求解域離散為一組單元的組合體，用在每個單元內(nèi)假設(shè)的近似函數(shù)來分片的表示求解域上待求的未知場函數(shù)，近似函數(shù)通常由未知場函數(shù)及其導數(shù)在單元各節(jié)點的數(shù)值插值函數(shù)來表達。從而使一個連續(xù)的無限自由度問題變成離散的有限自由度問題。有限元法作為目前工程中應用較為廣泛的一種數(shù)值計算方法，以其獨有的計算優(yōu)勢得到了廣泛的發(fā)展和應用，ANSYS軟件以其多物理場耦合分析功能而成為CAE（工程設(shè)計中的計算機輔助工程）軟件的應用主流，在熱分析工程應用中得到了較為廣泛的應用。2.2ANSYS簡介ANSYS程序是融結(jié)構(gòu)、熱、流體、電磁、聲學于一體的大型通用有限元商用分析軟件，可廣泛應用于核工業(yè)、鐵道、石油化工、航空航天、機械制造、能源、電子、造船、汽車交通、國防軍工、土木工程、生物醫(yī)學、輕工、地礦、水利、日用家電等各個領(lǐng)域。該軟件可在大多數(shù)計算機及操作系統(tǒng)中運行，從PC到工作站直至巨型計算機，ANSYS文件在其所有的產(chǎn)品系列和工作平臺上均兼容；該軟件基于Motif的菜單系統(tǒng)使用戶能夠通過對話框、下拉式菜單和子菜單進行數(shù)據(jù)輸入和功能選擇，此舉大大方便了用戶操作。它由世界上著名的有限元分析軟件公司美國ANSYS公司開發(fā)的，能與多數(shù)CAD軟件配合使用，實現(xiàn)數(shù)據(jù)的共享和交換，是現(xiàn)代產(chǎn)品設(shè)計中的高級CAD工具之一。在ANSYS公司相繼收購ICEM、CENTURYDYNAMICS、AAVIDTHERMAL、FLUENT等世界著名有限元分析程序制造公司并將產(chǎn)品整合之后，ANSYS實際上已成為世界上最通用和有效的商用有限元軟件。2.2.1ANSYS發(fā)展過程ANSYS公司是由美國匹茲堡大學的JohnSwanson博士在1970年創(chuàng)建的，其總部位于美國賓夕法尼亞州的匹茲堡，目前是世界CAE行業(yè)最大的公司。自創(chuàng)建以來，ANSYS公司一直致力于設(shè)計分析軟件的開發(fā)，不斷吸取新的計算方法和計算技術(shù)，領(lǐng)導著世界有限元技術(shù)的發(fā)展。ANSYS程序的最初版本與本文用的版本ANSYS10.0相比有很大區(qū)別。最初版本僅僅提供了熱分析和線性分析功能，是一個批處理程序，而且只能在大型計算機上使用。20世紀70年代初，隨著非線性、子結(jié)構(gòu)以及更多的單元類型的加入。ANSYS程序發(fā)生了很大的變化，新技術(shù)的融入進一步滿足了用戶的需求。20世紀70年代末。交互方式的加入是該軟件最為顯著的變化，此舉使得模型生成和結(jié)果評價大為簡化。2.2.2ANSYS技術(shù)特點與其他有限元計算軟件相比，ANSYS具有以下技術(shù)特征：●能實現(xiàn)多場及多場耦合功能?！窦昂筇幚?、分析求解及多場分析于一體?！癃氁粺o二的優(yōu)化功能，唯一具有流場優(yōu)化功能的CFD軟件?！窬哂袕姶蟮姆蔷€性分析功能?！窬邆淇焖偾蠼馄鳌！褡钤绮捎貌⑿杏嬎慵夹g(shù)的FEA軟件?！裰С謴奈C、工作站到巨型機的所有硬件平臺?！駨奈C、工作站、大型機直至巨型機，與所有硬件平臺的全部數(shù)據(jù)文件兼容?！駨奈C、工作站、大型機直至巨型機，所有硬件平臺具有統(tǒng)一用戶界面。●可與大多數(shù)的CAD軟件集成并具備接口。●可進行智能網(wǎng)格劃分。●具有多層次多框架的產(chǎn)品系列?！窬邆淞己玫挠脩糸_發(fā)環(huán)境。2.2.3ANSYS使用環(huán)境ANSYS程序是一個功能強大的有限元計算分析軟件包。它可運行于PC、NT工作站、UNIX工作站以及巨型計算機等各類計算機及操作系統(tǒng)中，其數(shù)據(jù)文件在其所有的產(chǎn)品系列和工作平臺上均兼容。該軟件有多種不同版本，本文應用版本為ANSYS10.0版，其微機版本要求的操作系統(tǒng)為Windows95/98/2000/XP，也可運行與UNIX系統(tǒng)中；微機版的基本硬件要求為：顯示分辨率為1024×768像素，顯示內(nèi)存為2MB以上，硬盤大于350MB，推薦使用17英寸顯示器。ANSYS多物理場耦合的功能，允許在同一模型上進行各種耦合計算，例如：熱-結(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ù)紸NSYS，并通過必要的修補可準確地在模型上劃分網(wǎng)格并進行求解，這樣就可以節(jié)省用戶在創(chuàng)建模型的過程中所花費的大量時間，使用戶的工作效率大幅度提高。2.2.4ANSYS程序功能ANSYS程序主要包括3個部分[30]：前處理模塊、求解模塊和后處理模塊。前處理模塊提供了一個強大的實體建模及網(wǎng)格劃分工具，用戶可以方便地構(gòu)造有限元模型；求解模塊包括結(jié)構(gòu)分析（結(jié)構(gòu)線性分析、結(jié)構(gòu)非線性分析和結(jié)構(gòu)高度非線性分析）、熱分析、流體動力學分析、電磁場分析、聲場分析、壓電分析以及多物理場的耦合分析，可模擬多種物理介質(zhì)的相互作用，具有靈敏度分析及優(yōu)化分析能力；后處理模塊可將計算結(jié)果以彩色等值線顯示、梯度顯示、矢量顯示、粒子流顯示、立體切片顯示、透明及半透明顯示等圖形方式顯示出來，也可將計算結(jié)果以圖表、曲線形式顯示或輸出。ANSYS程序提供了100種以上的單元類型，用來模擬實際工程中的各種結(jié)構(gòu)和材料。2.2.5ANSYS熱分析概述在實際生產(chǎn)過程中，常常會遇到多種多樣的熱量傳遞問題：如計算某個系統(tǒng)或部件的溫度分布、熱量的獲取或損失、熱梯度、熱流密度、熱應力、相變等。所涉及的部門包括：能源、化工、冶金、建筑、電子、航空航天、農(nóng)業(yè)、制冷、船舶等。以機械加工為例，往往需要估算和控制工件溫度場；分析不同工作條件下，不同材料及幾何形狀對溫度場變化的影響以及防止加工過程工件中缺陷的產(chǎn)生。ANSYS作為新穎的有限元分析軟件在熱分析問題方面具有強大的功能，而且界面友好，易于掌握。用戶可以隨心所欲地選擇圖形用戶界面方式（GUI）或命令流方式進行計算。ANSYS的熱分析模塊一般包括ANSYS/Muliphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五種，其中ANSYS/FLOTRAN不含相變熱分析。ANSYS熱分析基于能量守恒原理的熱平衡方程，用有限元法計算物體內(nèi)部各節(jié)點的溫度，并導出其他熱物理參數(shù)。運用ANSYS軟件可進行熱傳導、熱對流、熱輻射、相變、熱應力及接觸熱阻等問題的分析求解。2.2.6ANSYS求解步驟ANSYS求解模型的過程中的主要步驟：=1\*GB3①定義單元類型及選項，根據(jù)所分析的問題在單元庫中選擇分析問題的單元類型；=2\*GB3②定義單元是常數(shù)，單元是常數(shù)是與特定單元有關(guān)的量；=3\*GB3③定義材料屬性，在材料屬欄中選擇性的輸入熱導率、泊松比、密度等所需參數(shù)；=4\*GB3④建立幾何模型；=5\*GB3⑤定義網(wǎng)格并劃分參數(shù)；=6\*GB3⑥網(wǎng)格劃分以創(chuàng)建的模型。ANSYS在后處理器中可以查看結(jié)果：=1\*GB3①顯示等值線；=2\*GB3②顯示分析結(jié)果數(shù)據(jù)；=3\*GB3③對結(jié)果數(shù)據(jù)進行計算和路徑操作；=4\*GB3④繪制云圖。結(jié)合本文研究內(nèi)容及ANSYS10.0軟件的ANSYS/Thermal熱分析模塊功能，ANSYS10.0適用于本文研究內(nèi)容。3井筒熱損失3.1問題描述以遼河油田某鉆井為研究對象，注汽時間T=5d，油套環(huán)空導熱系數(shù)λyt為0.023W/m?K，套管導熱系數(shù)λtg為43.2W/m?K，水泥環(huán)導熱系數(shù)λcem為0.81W/m?K，地層導熱系數(shù)λcm為1.745W/m?K，地層熱擴散系數(shù)c為0.027m2/h，地層初始溫度Tcm=40oC，蒸汽與隔熱油管內(nèi)壁對流換熱系數(shù)α為10kw/m2?oC3.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.2478m，詳見圖3-1井筒及周圍地層截面示意圖。r8r7r8r7r6r5r2r1r4r3蒸汽隔熱層環(huán)空水泥環(huán)地層3.2.1能量平衡方程依據(jù)熱平衡原理，選取井筒內(nèi)的一微元流體單元體為研究對象，如圖3-2所示。E1E2E1E2123Q當蒸汽縱向流入或流出該單元體的控制面1、2時，攜帶的能量可表示為：（3-1）式中：—蒸汽攜帶的能量，W；—蒸汽流量，kg/s；—蒸汽焓值，J/kg；—蒸汽流速，m/s；—高度，m；—重力加速度，9.81m/s2。該流體單元的能量平衡表達式為：（3-2）式中：—單位時間的徑向散熱量，W。蒸汽的物性參數(shù)在井筒截面上是不變的，但不同截面處的物性參數(shù)是隨溫度變化的。由式（3-2）看出蒸汽沿井筒縱向的焓值在減少，蒸汽的能量是按遞減規(guī)律傳至地層的。3.2.2連續(xù)介質(zhì)的守恒方程A.井筒內(nèi)流體的能量平衡方程井筒油管內(nèi)高溫蒸汽在注入過程中遵循熱力學第一規(guī)律，其能量平衡方程為：（3-3）式中：、—z及z-dz處的焓，kJ；vz、vz-dz—z及z-dz處的速度，m/h；、—z及z-dz處的流體密度，kg/m3；Qdz—蒸汽散熱量，kJ。公式（3-3）表示包括蒸汽焓、勢能和動能在z處的蒸汽能量與z+dz處的蒸汽能量及蒸汽散熱量的平衡方程。B.井筒內(nèi)流體的連續(xù)性方程高溫蒸汽在油管內(nèi)流動時，質(zhì)量流量保持不變，因此滿足流體的連續(xù)性方程，即（3-4）式中：、—z及z+dz處的橫截面積，m2；vz、vz+dz—z及z+dz處的速度，m/h；、—z及z+dz處的流體密度，kg/m3。C.井筒內(nèi)流體的動量平衡方程把蒸汽在油管內(nèi)的流動看作是兩相流的流動，為描述其流動特性，通常從兩相流壓力降的研究出發(fā)。管內(nèi)穩(wěn)態(tài)均質(zhì)兩相流滿足動量平衡方程，即（3-5）式中：—段內(nèi)的壓力降，MPa；—蒸汽密度kg/m3；—摩阻系數(shù)；—井筒內(nèi)徑，m。表示段內(nèi)的勢能變化，表示段內(nèi)的摩擦損失，表示段內(nèi)的動能變化。也就是說管內(nèi)兩相流的壓力降是摩擦損失、勢能變化和動能變化的綜合結(jié)果。D.油管內(nèi)濕蒸汽的焓由于濕蒸汽為汽水兩相混合物，所以它的焓可用相同溫度或壓力下的飽和水、飽和蒸汽及濕蒸汽的干度來表示，即（3-6）式中：—蒸汽干度；—濕蒸汽的焓，kJ/kg；—飽和蒸汽的焓，kJ/kg；—飽和水的焓，kJ/kg。3.2.3傳熱過程分析 根據(jù)熱力學[31]和傳熱學[32]知識對井筒傳熱進行分析，由于水泥環(huán)外側(cè)的地層體積和熱容量比水泥環(huán)內(nèi)側(cè)的井筒大得多，水泥環(huán)與地層交界處的溫度隨注入時間的推移緩慢上升。表征熱量在地層中擴散速率的參數(shù)——熱擴散系數(shù)α的數(shù)量級僅為10-72m/s，因而與未受熱影響的遠處地層溫度的差值不大；相反，隔熱管內(nèi)蒸汽溫度與的差別較大，有一個數(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-7）式中：—段內(nèi)井筒的徑向熱損失，W；—隔熱油管外半徑，m；—以隔熱油管內(nèi)管外表面為基準面積的總傳熱系數(shù)，W/（m2?k）。當井筒內(nèi)隔熱管抽真空時，可表示為：（3-8）式中：—蒸汽凝結(jié)換熱熱阻；—隔熱油管熱阻；—隔熱管內(nèi)隔熱層換熱熱阻；—環(huán)空輻射換熱熱阻；—環(huán)空對流熱熱阻；—套管壁導熱熱阻；—水泥環(huán)導熱熱阻。B.地層內(nèi)非穩(wěn)態(tài)傳熱在井筒穩(wěn)態(tài)傳熱的同一微元段dz內(nèi)，地層非穩(wěn)態(tài)導熱方程如下：（3-9）初始條件時：（3-10）邊界條件時：（3-11）（3-12）式中：—地層內(nèi)某處溫度，oC；—地層導熱系數(shù)，W/（m?K）；—地層擴散系數(shù)，。3.2.4基本假設(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-13）式中：x，y分別指的是x和y方向；τ是時間，s；α是熱擴散率，m2/s。3.2.6邊界條件蒸汽與油管內(nèi)壁之間屬于對流換熱，所以油管內(nèi)壁為第三類邊界，應滿足：（3-14）在不考慮地溫梯度的情況下，地層溫度為定值，因此地層邊界處為第一類邊界條件，其他邊界均為絕熱邊界。3.3ANSYS模擬利用ANSYS10.0熱分析軟件，做出在不同隔熱油管視導熱系數(shù)和注汽溫度時的井筒溫度場分布情況。3.3.1建立有限元模型（1）打開ANSYS10.0有限元分析軟件，過濾圖形用戶界面；（2）確定工作文件名“heatlossinwell”，單位制默認為國際單位制；（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.定義分析類型定義新的分析類型為瞬態(tài)分析“Transient”。B.設(shè)置瞬態(tài)熱分析的初始條件由于初始溫度場是不均勻的且又是未知的，所以首先要進行穩(wěn)態(tài)熱分析，為穩(wěn)態(tài)熱分析建立初始條件：在井筒壁處施加對流載荷，在地層處施加原始地層溫度場載荷，其他邊界為絕熱邊界。圖3-3井筒有限元模型圖3-4網(wǎng)格劃分結(jié)果3.3.3求解采用GUI命令：Main｜Solution｜LoadStepOpts｜TimeFrequenc｜Time-TimeStep，設(shè)置計算終止時間為0.01s，時間步長為0.01s，并將時間積分效果設(shè)定為Off，對模型進行穩(wěn)態(tài)熱分析。穩(wěn)態(tài)熱分析結(jié)束后，將時間積分效果設(shè)定為On，并刪除所有節(jié)點的溫度載荷，計算終止時間設(shè)定為5天，時間步長設(shè)定為1小時，對模型進行瞬態(tài)熱分析。3.3.4后處理對于瞬態(tài)熱分析，ANSYS提供了兩種后處理方式，即POST1和POST26。本文采用POST1進行后處理，利用GUI命令：MainMenu｜GeneralPostproc｜PlotResults顯示溫度場彩色云圖，并拾取套管內(nèi)徑邊緣處溫度。3.4模擬結(jié)果與分析取水泥環(huán)與地層之間截面的平均溫度作為溫度參數(shù)，并結(jié)合文獻[33]中給出的熱損失經(jīng)驗公式，分別算出在不同的隔熱油管視導熱系數(shù)及注汽溫度下井筒熱損失情況。3.4.1隔熱油管視導熱系數(shù)對井筒熱損失的影響現(xiàn)在使用的隔熱油管多為抽真空隔熱油管，結(jié)合預應力隔熱油管標準（SY/T5324-94）中提出的視導熱系數(shù)概念，把隔熱管的隔熱層內(nèi)以導熱、對流和輻射三種方式傳遞的熱量視為以純導熱方式傳遞的熱量。其中，視導熱系數(shù)數(shù)值劃分為五個范圍等級，如表3-1所示，因此，分別選取視導熱系數(shù)為0.07、0.05、0.03、0.01和0.005W/m?K。表3-1隔熱油管等級分類單位：W/m?K隔熱性能等級ABCDE視導熱系數(shù)0.06≤λ<0.080.04≤λ<0.060.02≤λ<0.040.006≤λ<0.020.002≤λ<0.006在實際的稠油熱采過程中，注汽溫度一般為300~350℃之間，基于此，本文以蒸汽溫度為312℃情況下進行模擬，圖3-5為在視導熱系數(shù)為0.07時的溫度分布云圖，結(jié)合熱損失經(jīng)驗公式，得到單位長度井筒熱損失與視導熱系數(shù)之間的關(guān)系，如表3-2和圖圖3-5視導熱系數(shù)為0.07W/m?K時的溫度分布云圖表3-2不同視導熱系數(shù)下的計算結(jié)果序號視導熱系數(shù)（W/m?K）套管溫度（oC）熱損失（W/m）123450.0050.010.030.050.0741.58842.84846.04647.79748.9022.7154.87010.33813.33215.222圖3-6隔熱油管視導熱系數(shù)對井筒熱損失的影響根據(jù)圖3-6分析可知，井筒單位長度熱損失整體上隨著隔熱油管視導熱系數(shù)的增大而增大。其中，對于A級的隔熱油管，單位長度井筒熱損失約達15W/m左右，而對于E級的隔熱油管，單位長度井筒熱損失僅約為3.5W/m，兩者單位長度井筒熱損失相差約11.5W/m，因此，選用視導熱系數(shù)高級別的隔熱油管對于降低井筒單位長度的熱損失具有明顯的影響。此外，當視導熱系數(shù)從0.07W/m?K提高至0.03W/m?K后，單位長度熱損失的變化速率隨著視導熱系數(shù)的減小而出現(xiàn)明顯的降低。綜上所述，在實際的應用過程中，建議選用隔熱油管視導熱系數(shù)在0.03~0.005W/m?K。3.4.2注汽溫度對單位長度井筒熱損失的影響上述視導熱系數(shù)的建議范圍內(nèi)，在視導熱系數(shù)為0.03W/m?K情況下，分別研究以注汽溫度分別為312、324、333、346和350oC五種情況下，注汽溫度對單位長度井筒熱損失的影響，詳細情況如表3-3和圖3-7所示，圖3-8為當注汽溫度為324oC時的溫度分布云圖。表3-3不同注汽溫度下的計算結(jié)果序號注汽溫度（oC）套管溫度（oC）熱損失（W/m）1234531232433334635048.90249.29749.59249.98750.08215.22215.89716.40117.07717.239根據(jù)上圖分析所知，單位長度井筒的熱損失隨著注汽溫度的升高幾乎呈線性增大的趨勢。在注汽溫度為312oC情況下，單位長度井筒熱損失約為15.2W/m，在注汽溫度在350oC情況下，單位長度井筒熱損失約為17.2W/m，由此可見，注汽溫度的提高對于井筒單位長度熱損失的影響較小，每提高50oC，井筒單位長度熱損失相差約2W/m。圖3-7注汽溫度對井筒熱損失的影響圖3-8注汽溫度為324oC時的溫度分布云圖4油層熱影響區(qū)域4.1問題描述高溫高壓的蒸汽流經(jīng)井筒后進入封隔器以下的油層區(qū)，通過射孔與油層進行傳質(zhì)傳熱，傳熱過程如圖4-1所示。圖4-1油層傳熱過程示意圖以遼河油田某鉆井為研究對象，注汽時間T=10d，套管導熱系數(shù)λtg為43.2W/m?K，水泥環(huán)導熱系數(shù)λcem為1.2W/m?K，油層綜合導熱系數(shù)λcm為2.53W/m?K，地層導熱系數(shù)為1.968W/m?K，地層熱擴散系數(shù)c為0.027m2/h，地層初始溫度Tcm=40oC，蒸汽與隔熱油管內(nèi)壁對流換熱系數(shù)α為10kw/m2?oC4.2數(shù)學模型的建立套管內(nèi)半徑r1為0.0809m，外半徑r2為0.0889m；水泥環(huán)外半徑r3為0.1239m；地層和油層半徑r4為20m，詳見圖4-2油層蒸汽蒸汽水泥環(huán)地層油層套管r1r2r30.5mm1m圖4-2油層截面示意圖4.2.1連續(xù)介質(zhì)的守恒方程求解與流體中動量、熱量和質(zhì)量傳遞有關(guān)的問題時，首先應根據(jù)基本的守恒定律對該物理系統(tǒng)建立一組方程。這些方程包括:=1\*GB3①能量平衡方程；=2\*GB3②連續(xù)性方程；=3\*GB3③焓方程。A.連續(xù)介質(zhì)的能量平衡方程井筒油管內(nèi)高溫蒸汽在注入過程中遵循熱力學第一規(guī)律，其能量平衡方程為：（4-1）式中：Hz、Hz-dz—z及z-dz處的焓，kJ；vz、vz-dz—z及z-dz處的速度，m/h；、—z及z-dz處的流體密度，kg/m3；Qdz—蒸汽散熱量，kJ。上式表示包括蒸汽焓、勢能和動能在z處的蒸汽能量與z+dz處的蒸汽能量及蒸汽散熱量的平衡方程。B.連續(xù)介質(zhì)的連續(xù)性方程高溫蒸汽在油管內(nèi)流動時，質(zhì)量流量保持不變，因此滿足流體的連續(xù)性方程，即（4-2）式中：、—z及z+dz處的橫截面積，m2；vz、vz+dz—z及z+dz處的速度，m/h；、—z及z+dz處的流體密度，kg/m3。C.連續(xù)介質(zhì)的動量平衡方程把蒸汽在油管內(nèi)的流動看作是兩相流的流動，為描述其流動特性，通常從兩相流壓力降的研究出發(fā)。管內(nèi)穩(wěn)態(tài)均質(zhì)兩相流滿足動量平衡方程，即（4-3）式中：—段內(nèi)的壓力降，MPa；—蒸汽密度kg/m3；—摩阻系數(shù)；—井筒內(nèi)徑，m。表示段內(nèi)的勢能變化，表示段內(nèi)的摩擦損失，表示段內(nèi)的動能變化。D.油管內(nèi)濕蒸汽的焓由于濕蒸汽為汽水兩相混合物，所以它的焓可用相同溫度或壓力下的飽和水、飽和蒸汽及濕蒸汽的干度來表示，即（4-4）式中：—蒸汽干度；—濕蒸汽的焓，kJ/kg；—飽和蒸汽的焓，kJ/kg；—飽和水的焓，kJ/kg。4.2.2多孔介質(zhì)的守恒方程A.多孔介質(zhì)的質(zhì)量守恒方程多孔介質(zhì)的質(zhì)量守恒方程為：（4-5）式中：—密度，kg/m3；—介質(zhì)的孔隙度；—滲流速度，m/d。B.多孔介質(zhì)的動力學方程a.單相流在只有飽和牛頓流體的空隙介質(zhì)中，達西定律如下：（4-6）式中：—流體的滲流速度，m/d；—滲透率，；—粘度，Pa?s；—流體的密度，kg/m3；—壓力，MPa；D—垂直高度，m；G—重力加速度，m/s2。b.多相流多相流的滲流微分方程為：(4-7)式中：sl—l相流體的飽和度；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的徑向熱流量，可表示為：（4-8）式中：—段內(nèi)井筒的徑向熱流量，W；—套管外半徑，m；—以套管外表面為基準面積的總傳熱系數(shù)，W/（m2?k）。的表達式為：（4-9）式中：—蒸汽與套管內(nèi)壁對流換熱熱阻；—套管壁導熱熱阻；—水泥環(huán)導熱熱阻。B.油層傳熱分析井筒蒸汽與套管通過對流換熱向井筒周圍散熱，同時產(chǎn)生從射孔到油層的傳質(zhì)熱傳導。由于在有限元中無法解決傳質(zhì)問題，因此一般簡化為具有一定半徑的加熱區(qū)模型。非穩(wěn)態(tài)無內(nèi)熱源軸對稱溫度場的微分方程為：（4-10）應用Galerkin法可得到軸對稱溫度場整體區(qū)域“變分”計算的基本方程為:(4-11)利用溫度場變分方程對模型進行剖分離散、單元分析及最終總體合成，即可完成模型溫度場的求解。4.2.4基本假設(shè)=1\*GB3①沿井深方向上的熱能不發(fā)生傳遞；=2\*GB3②視蒸汽流為一維均質(zhì)兩相流，且無相變；=3\*GB3③忽略油層溫度梯度的影響；=4\*GB3④以井眼中心為中心的同心圓上的溫度是相同的。4.2.5邊界條件蒸汽與套管內(nèi)壁之間屬于對流換熱，所以套管內(nèi)壁為第三類邊界，應滿足：（4-12）加熱區(qū)為第一類邊界條件，由于不受油層溫度梯度的影響，所以油層為第一類邊界條件，其他為絕熱邊界。4.3ANSYS模擬4.3.1建立有限元模型=1\*GB3①打開ANSYS10.0有限元分析軟件，過濾圖形用戶界面；=2\*GB3②確定工作文件名“youcengreyingxiangquyu”，單位制默認為國際單位制；=3\*GB3③進入PREP7前處理器，選擇plane55熱分析單元，并設(shè)置成軸對稱單元；=4\*GB3④將已知物性參數(shù)定義為材料熱物性參數(shù)；=5\*GB3⑤根據(jù)問題描述建立有限元模型，結(jié)果如圖4-3所示；圖4-3油層有限元模型（局部放大圖）=6\*GB3⑥采用自由網(wǎng)格劃分方法對模型進行網(wǎng)格劃分，結(jié)果如圖4-4所示。圖4-4網(wǎng)格劃分結(jié)果（局部放大圖）4.3.2施加載荷計算A.定義分析類型定義新的分析類型為瞬態(tài)分析“transient”。B.設(shè)置瞬態(tài)熱分析的初始條件由于初始溫度場是不均勻的且又是未知的，所以首先要進行穩(wěn)態(tài)熱分析，為穩(wěn)態(tài)熱分析建立初始條件：在套管內(nèi)壁處施加對流載荷，在油層和地層處施加原始溫度場載荷，其他邊界為絕熱邊界。4.3.3求解采用GUI命令：Main｜Solution｜LoadStepOpts｜TimeFrequenc｜Time-TimeStep，設(shè)置計算終止時間為0.01s，時間步長為0.01s，并將時間積分效果設(shè)定為Off，對模型進行穩(wěn)態(tài)熱分析。穩(wěn)態(tài)熱分析結(jié)束后，將時間積分效果設(shè)定為On，并刪除所有節(jié)點的溫度載荷，計算終止時間設(shè)定為10天，時間步長設(shè)定為1小時，對模型進行瞬態(tài)熱分析。4.3.4后處理采用POST1進行后處理，利用GUI命令：MainMenu｜GeneralPostproc｜PlotResults顯示溫度場彩色云圖。4.4模擬結(jié)果與分析模擬油層溫度場分布云圖如圖4-5至圖4-14所示。圖4-5注汽第1天溫度分布云圖圖4-6注汽第2天溫度分布云圖圖4-7注汽第3天溫度分布云圖圖4-8注汽第4天溫度分布云圖圖4-9注汽第5天溫度分布云圖圖4-10注汽第6天溫度分布云圖圖4-11注汽第7天溫度分布云圖圖4-12注汽第8天溫度分布云圖圖4-13注汽第9天溫度分布云圖圖4-14注汽第10天溫度分布云圖油層熱影響區(qū)域隨時間的變化如圖4-15。圖4-15油層熱影響區(qū)域根據(jù)以上各圖分析可知，油層高溫區(qū)從第3天開始逐漸向地層擴散，造成熱損失。油層熱影響區(qū)域先隨著注汽時間的增大而增大，大約以0.2m/d的速度增大，第5天熱影響區(qū)域達到最大，而后熱影響區(qū)域基本保持不變，維持在大約17.4米左右。5結(jié)論A.井筒熱損失（1）井筒單位長度熱損失整體上隨著隔熱油管視導熱系數(shù)的增大而增大。視導熱系數(shù)從0.07W/m?K提高至0.03W/m?K后，單位長度熱損失的變化速率隨著視導熱系數(shù)的減小而出現(xiàn)明顯的降低，建議選用隔熱油管視導熱系數(shù)在0.03~0.005W/m?K；（2）單位長度井筒的熱損失隨著注汽溫度的升高幾乎呈線性增大的趨勢。每提高50℃，井筒單位長度熱損失相差約2W/m（3）視導熱系數(shù)對井筒單位長度熱損失的影響大于注汽溫度的影響。B.油層熱影響區(qū)域（1）高溫區(qū)從第3天開始逐漸向地層擴散，造成熱損失；（2）油層熱影響區(qū)域先隨著注汽時間的增大而逐漸增大，在第5天達到最大，而后熱影響區(qū)域基本保持不變，大約在17.4米左右；（3）綜合1、2兩點，建議注汽時間在3~5天。參考文獻[1]RameyH.J.WellboreHeatTransmission[J].JPT,1962,14（4）:427-435[2]HasanAR.KabirCS.HeatTransferDuringTwo-PhaseFlowinWellboresPartIFormationTemperature[A].SPE22866,1991[3]HasanAR.KabirCS.HeatTransferDuringTwo-PhaseFlowinWellboresPareIIWellboreFluidTemperature[A].SPE22948,1991[4]HasanAR.KabirCS.AStudyofMultiphaseFlowBehaviorinVerticalWells[A].SPE264,1988[5]HasanAR.KabirCS.Two-PhaseFlowinVerticalandInclinedAnnuli[J].Intl.J.MultiphaseFlow,1992,18（2）:274[6]HasanAR.KabirCS.AspectsofWellboreHeatTransferDuringTwo-PhaseFlow[A].SPE22948,1994[7]FarouqAliSM.AComprehensiveWellboreSteam/WaterFlowModelforSteamInjectionandGeothermalApplications[A].SPE527,1981[8]SetterA.HeatlossesDuringFlowofSteamDownaWellbore[J].JPT,1965,17（7）:845-851[9]SquierDP,SmithDD,DoughertyEL.CalculatedTemperatureBehaviorofHot-WaterInjectionWells[J].JPT,1962,（4）:436-440[10]HuygenHH.HuittJL.WellboreHeatLossesandCasingTemperaturesDuringSteamInjection[J].Prod.Montlily,1966,30（8）:2[11]GristionS.WillhiteGP.