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

可修編可修編-石油化工大學(xué)畢業(yè)設(shè)計(jì)(論文)GraduationProject(Thesis)forUndergraduateofLSHU題目TITLE學(xué)院School專業(yè)班級(jí)一Major&Class姓名一Name指導(dǎo)教師一Supervisor20年月日論文獨(dú)創(chuàng)性聲明本人所呈交的論文，是在指導(dǎo)教師指導(dǎo)下，獨(dú)立進(jìn)行研究和開發(fā)工作所取得的成果。除文中已特別加以注明引用的容外，論文中不包含任何其他個(gè)人或集體已經(jīng)發(fā)表或撰寫過的研究成果。對(duì)本文的工作做出重要貢獻(xiàn)的個(gè)人和集體，均已在文中以明確方式標(biāo)明并致。本聲明的法律結(jié)果由本人承擔(dān)。特此聲明。論文作者（簽名）：年月日摘要針對(duì)稠油熱采過程中注入大量高溫蒸汽而造成井筒部分熱損失難以控制的問題，從熱工角度初步分析了高溫蒸汽在井筒部分的傳熱過程，建立了井筒及周圍地層部分的二維、瞬態(tài)傳熱數(shù)學(xué)模型，根據(jù)所求得的溫度場(chǎng)以及定量熱損失經(jīng)驗(yàn)公式，探討了隔熱油管視導(dǎo)熱系數(shù)和蒸汽溫度對(duì)單位長(zhǎng)度井筒熱損失的影響。研究表明：井筒單位長(zhǎng)度熱損失隨著隔熱油管視導(dǎo)熱系數(shù)和蒸汽溫度的增加而增加，且隔熱油管視導(dǎo)熱系數(shù)對(duì)井筒單位長(zhǎng)度熱損失的影響更大。高溫高壓蒸汽流經(jīng)井筒段后進(jìn)入油層段，蒸汽通過射孔進(jìn)入油層，并與油層進(jìn)行傳質(zhì)傳熱從而形成油層熱影響區(qū)域。本文利用有限元分析軟件ANSYS,對(duì)遼河某熱采井建立了井筒、地層和油層部分的二維、非穩(wěn)態(tài)數(shù)學(xué)模型，并模擬分析了油層和周圍地層溫度場(chǎng)，在此基礎(chǔ)之上探討了注汽時(shí)間對(duì)油層熱影響區(qū)域的影響。研究表明：油層熱影響區(qū)域先隨著注汽時(shí)間的增大而增大，而后基本不變，且注汽時(shí)間在3~5天最佳。關(guān)鍵詞：視導(dǎo)熱系數(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目錄緒論11.1目的意義11.2國(guó)外的研究現(xiàn)狀2井筒熱損失的研究2油層熱影響區(qū)域的研究31.3本文的研究容及創(chuàng)新點(diǎn)4研究容4創(chuàng)新點(diǎn)4研究方法62.1有限元分析法簡(jiǎn)介62.2ANSYS簡(jiǎn)介6ANSYS發(fā)展過程7ANSYS技術(shù)特點(diǎn)7ANSYS使用環(huán)境8ANSYS程序功能8ANSYS熱分析概述9ANSYS求解步驟9井筒熱損失113.1問題描述113.2數(shù)學(xué)模型的建立11能量平衡方程11連續(xù)介質(zhì)的守恒方程13傳熱過程分析14基本假設(shè)16控制方程16邊界條件163.3ANSYS模擬16建立有限元模型16施加載荷計(jì)算17求解18后處理183.4模擬結(jié)果與分析18隔熱油管視導(dǎo)熱系數(shù)對(duì)井筒熱損失的影響18注汽溫度對(duì)單位長(zhǎng)度井筒熱損失的影響204油層熱影響區(qū)域224.1問題描述224.2數(shù)學(xué)模型的建立22連續(xù)介質(zhì)的守恒方程23多孔介質(zhì)的守恒方程25傳熱過程分析26基本假設(shè)27邊界條件274.3ANSYS模擬27建立有限元模型27施加載荷計(jì)算28求解29后處理294.44.4模擬結(jié)果與分析29可修編可修編-5結(jié)論32參考文獻(xiàn)33辭351.11.1目的意義可修編可修編-緒論稠油在世界石油資源中儲(chǔ)量豐富，據(jù)統(tǒng)計(jì)，其總量約為1000x108噸，中國(guó)是世界上稠油資源豐富的國(guó)家之一，預(yù)測(cè)資源量約為198億噸。稠油由于粘度大，流動(dòng)性差等特點(diǎn)，常規(guī)開采方法并不適用，基于稠油對(duì)溫度的敏感性特點(diǎn)，熱力采油已成為國(guó)外普遍認(rèn)同的開采方法，即向地下油層注入高溫高壓蒸汽，一方面可以提高油層溫度、降低稠油粘度；另一方面可使油層壓力升高、增加驅(qū)油能力。常見的熱力采油方法有“蒸汽吞吐采油”和“蒸汽驅(qū)采油”兩種方式?！罢羝掏虏捎汀笔窍认蛴途⑷胍欢康恼羝?，燜井一段時(shí)間，待蒸汽的熱能向油層擴(kuò)散后，再開井生產(chǎn)的一種開采稠油的增產(chǎn)方法。蒸汽吞吐作業(yè)的過程分為三個(gè)階段，即注汽、燜井及回采?！罢羝?qū)采油”是稠油油藏經(jīng)過蒸汽吞吐采油之后，為進(jìn)一步提高采收率而采取的一項(xiàng)熱采方法，因?yàn)檎羝掏虏捎椭荒懿捎透鱾€(gè)油井附近油層中的原油，在油井與油井之間還留有大量的死油區(qū)?！罢羝?qū)采油”是由注入井連續(xù)不斷地往油層中注入高干度的蒸汽，蒸汽不斷地加熱油層，從而大大降低了地層原油的粘度。注入的蒸汽在地層中變?yōu)闊岬牧黧w，將原油驅(qū)趕到生產(chǎn)井的周圍，并被采到地面上來。在注蒸汽的過程中，由于蒸汽與地層之間存在溫差，蒸汽流動(dòng)過程中又有摩擦阻力，因此注汽過程中存在著能量損失，注汽過程中的能量損失，特別是井筒中的能量損失，直接影響熱采效果。為了有效地減少這部分能量損失，提高熱采效果，需要對(duì)注汽井筒進(jìn)行傳熱分析，進(jìn)而探討隔熱油管視導(dǎo)熱系數(shù)和蒸汽溫度對(duì)單位長(zhǎng)度井筒熱損失的影響，這為熱采過程中熱量的高效、合理利用提供了一定的理論依據(jù)。蒸汽流經(jīng)井筒后進(jìn)入目的層即油層，高溫高壓的蒸汽會(huì)通過射孔與油層進(jìn)行傳質(zhì)傳熱，將熱量釋放給原油，使其溫度升高，粘度降低，同時(shí)增加油層壓力，達(dá)到驅(qū)油效果。但隨著蒸汽與油層熱量交換的進(jìn)行，蒸汽自身溫度降低，逐漸冷凝成液體，阻止熱量繼續(xù)擴(kuò)散，會(huì)形成熱影響區(qū)域。研究油層熱影響區(qū)域?qū)侠泶_定注汽時(shí)間和鉆井位置具有重要指導(dǎo)意義。可修編可修編-國(guó)外的研究現(xiàn)狀井筒熱損失的研究自本世紀(jì)初以來，就有許多研究者致力于井筒熱損失的研究。井筒傳熱研究工作的先驅(qū)Ramey以能量平衡為基礎(chǔ)，忽略了流體流動(dòng)過程中動(dòng)能和摩擦的影響，最先提出了考慮時(shí)間和井深影響時(shí)的井筒液體溫度計(jì)算的理論模型。該模型作為井筒熱損研究的雛型，給出了井筒注入氣體和流體兩種工質(zhì)關(guān)于井深和時(shí)間的溫度解析表達(dá)式，為后來的井筒熱損研究奠定了基礎(chǔ)。但該模型中的時(shí)間函數(shù)要求注入流體的時(shí)間不得小于7天，所以此模型的精度在時(shí)間小于7天的時(shí)候受到限制。另外該模型僅適用于單相流，并在建立模型過程中假設(shè)井筒流體的物性參數(shù)不隨溫度變化，這與實(shí)際物理過程不相符［1］。國(guó)外許多研究者在Ramey模型的基礎(chǔ)上繼續(xù)做了多方面的研究，如Hasan和心回「建立起井筒兩相流模型，改善了Ramey模型僅適用于單相流的不足，此模型可用于計(jì)算井筒液體溫度的分布［2~6］。Pacheco和FarouqAli以流體流動(dòng)的能量平衡方程為基礎(chǔ)，建立了能同時(shí)計(jì)算蒸汽壓力和蒸汽干度的數(shù)學(xué)模型［7］。該模型彌補(bǔ)了Ramey模型中沒有考慮水泥殼的影響的缺點(diǎn)。另外此模型結(jié)構(gòu)形式較復(fù)雜。但它可以用來較精確地計(jì)算蒸汽物性和熱損失隨深度和時(shí)間的變化。不過此模型還用到了Ramey模型中不完善的理論，如地層時(shí)間函數(shù)￡(t)的計(jì)算，所以適用精度也受到一定限制。Satter在使用Ramey模型時(shí)，考慮了注汽過程中的蒸汽相變影響，進(jìn)而提出了一種預(yù)測(cè)冷凝流體干度隨深度及時(shí)間變化的計(jì)算方法，模型中提出的關(guān)系式［8］可用于計(jì)算飽和蒸汽注入時(shí)的井筒熱損失。Squier等人對(duì)熱水通過井筒注入時(shí)提出了一個(gè)完整的分析方法⑼；Huygen和Huit?］給出了在蒸汽注入期間的井筒熱損失理論和實(shí)驗(yàn)結(jié)果，同時(shí)指出輻射熱損失的重要性；Willhite提出了完整的計(jì)算總傳熱系數(shù)的方法⑺。Holst和Flock以及Earlougher總結(jié)前面的一些成果，并考慮了蒸汽壓力的計(jì)算［12,13］；Eickmeier等人采用有限差分模型計(jì)算注入初期過渡過程特性/。Hasan等人又對(duì)井筒壓力降和熱損失計(jì)算進(jìn)行了完善。但這些模型對(duì)于目前井筒注汽熱損失的計(jì)算尚有一些不足，首先它們?cè)谟?jì)算過程中，把某些參數(shù)看作是不隨溫度變化的常量，這與實(shí)際并不完全相符；另外由于隔熱管的發(fā)展，上述模型中的總傳熱系數(shù)不能較真實(shí)地反映井筒現(xiàn)在結(jié)構(gòu)，給計(jì)算精度的提高帶來困難。我國(guó)對(duì)稠油熱采的相關(guān)研究比較晚。胡智勉、王彌康對(duì)注汽井筒總傳熱系數(shù)的計(jì)算方法提出了一些見解［忸16］。王彌康在Ramey模型的基礎(chǔ)上，從傳熱機(jī)理出發(fā)，給出了蒸汽井井筒熱傳遞定量計(jì)算的數(shù)學(xué)模型，該模型避免了Ramey公式僅適用于長(zhǎng)時(shí)間注汽的不足，擴(kuò)大了計(jì)算時(shí)間的適用圍，從而改善了井筒傳熱計(jì)算的精度［17］?；莘粚⒌孛婀芫€與井筒進(jìn)行了綜合考慮，探討了計(jì)算蒸汽壓力降、干度降及熱量損失的計(jì)算模型［18］。景勤、艷華、王志國(guó)等人建立了計(jì)算井筒熱損數(shù)學(xué)模型，并提出了熱損合理分布的概念［19~21］。倪學(xué)鋒提出了井筒參數(shù)計(jì)算模型［22］。文章提出了用物理模擬方法確定井筒總傳熱系數(shù)的方法［23］。油層熱影響區(qū)域的研究1968年SpilletteNielsen通過建立一個(gè)二維模型，對(duì)油層注冷水和注熱水后溫度分布進(jìn)行了模擬計(jì)算，把流場(chǎng)分布和溫度分布藕合在了一起進(jìn)行迭代［24］。1970年埃克米爾等人在研究井筒換熱和溫度分布時(shí)，建立了一套關(guān)于注液和產(chǎn)液期間，液體與井筒周圍地層間熱交換的有限差分模型。計(jì)算中將油管、環(huán)空液、套管、水泥環(huán)及地層情況全部考慮在。該模型可以有效地模擬實(shí)際井筒條件，計(jì)算短時(shí)間和長(zhǎng)時(shí)間注液、產(chǎn)液時(shí)井筒溫度的變化，模型只是基于熱傳導(dǎo)的模式進(jìn)行分析的，沒有考慮注入層的對(duì)流傳熱和井筒中的對(duì)流傳熱。1970年，Smith，R.C.，Steffensen，口".等人提出了利用溫度測(cè)井的計(jì)算機(jī)分析估計(jì)流動(dòng)剖面附。他們通過改變影啊溫度特性的參數(shù)，對(duì)注水井的溫度測(cè)井進(jìn)行了詳細(xì)的研究，研究結(jié)果表明溫度測(cè)井對(duì)地層的熱學(xué)特性以及井眼的幾何形狀具有敏感性:另外，指出了注入水的溫度和小流量的漏失對(duì)溫度測(cè)井的影響，并對(duì)溫度測(cè)井解釋提出了指導(dǎo)性意見。1972年維脫霍特與梯克薩通過改變注入速率、注入深度和注入時(shí)間，來研究注入期間井下的溫度分布，并對(duì)注入層中的溫度分布進(jìn)行了研究。1973年Steffenson.RJ，SmlthR.C.等人提出在油層中，由于流動(dòng)的存在，勢(shì)必發(fā)生熱交換，要計(jì)算井眼的溫度必須確定井筒周圍油藏的大片區(qū)域的溫度分布，并提出了著名的Joul-Thomson效應(yīng)。1982年JohnFagley，H.SeottFogler，C.BrentDavenport，口6.乂3m^6等人對(duì)關(guān)井前后井筒的熱傳遞進(jìn)行了數(shù)值模擬閩，提出關(guān)井前一段時(shí)間的熱注是一種解釋老井注入剖面的潛在的重要工具。1988年袁益壯針對(duì)在油藏中水驅(qū)油的流動(dòng)和傳熱問題提出了多孔介質(zhì)中完全可壓縮、可混溶驅(qū)動(dòng)問題的差分方法。1994年FredericMaubeuge，MichelP.Didek，MartynB.Beardsen等人提出了一個(gè)多層注入模型［27］。1995年胡耀江等人建立了水飽和的一維多孔介質(zhì)在一端注入蒸汽后，熱液區(qū)溫度分布的數(shù)學(xué)模型，提出了運(yùn)用拉氏變換及數(shù)值反演求解此問題的計(jì)算方法［28］。1996年恩民等人在建立的溫度模型中，除了考慮熱傳導(dǎo)外，特別強(qiáng)調(diào)了水井中與油層傳質(zhì)傳熱，構(gòu)造了關(guān)于導(dǎo)熱系數(shù)、比熱等參變量的優(yōu)化識(shí)別模型［29］。還有部分學(xué)者利用解析法求解地層中的溫度分布。本文的研究容及創(chuàng)新點(diǎn)研究容本文研究容主要分為兩個(gè)部分，即井筒熱損失和油層熱影響區(qū)域。4井筒熱損失①利用傳熱學(xué)、熱力學(xué)、流體力學(xué)等知識(shí)分析隔熱油管管至地層的傳熱過程，建立連續(xù)介質(zhì)的守恒方程，并建立數(shù)學(xué)模型；②利用有限元軟件ANSYS建立隔熱油管管至地層的有限元模型，并利用其熱分析模塊對(duì)模型進(jìn)行熱分析，求得注汽過程中井筒的溫度場(chǎng)分布云圖；③探討注汽溫度和油管視導(dǎo)熱系數(shù)對(duì)井筒熱損失的影響。BUM層熱影響區(qū)域①利用傳熱學(xué)、熱力學(xué)、流體力學(xué)等知識(shí)分析套管至油層的傳熱過程，建立連續(xù)介質(zhì)的守恒方程和多孔介質(zhì)的守恒方恒，并建立數(shù)學(xué)模型；②利用有限元軟件ANSYS建立套管至油層的有限元模型，并利用其熱分析模塊對(duì)模型進(jìn)行熱分析，求得不同時(shí)刻油層的溫度場(chǎng)分布云圖；③探討注汽時(shí)間與油層熱影響區(qū)域之間的關(guān)系。創(chuàng)新點(diǎn)（1）在以往的井筒溫度場(chǎng)及熱損失研究中，以數(shù)值解法為主，本文利用有限元軟件ANSYS的熱分析模塊求得了井筒溫度場(chǎng)分布云圖，并結(jié)合熱損失計(jì)算經(jīng)驗(yàn)公式，探討了注汽溫度和油管視導(dǎo)熱系數(shù)對(duì)熱損失的影響。與數(shù)值法相比，本文研究方法可以清楚、直觀地觀察到溫度場(chǎng)的分布情況。(2)利用有限元軟件ANSYS的熱分析模塊模擬分析了油層熱影響區(qū)域，并探討了注汽時(shí)間與油層熱影響區(qū)域之間的關(guān)系。2研究方法有限元分析法簡(jiǎn)介有限元法（finiteelementmethod）是一種高效能、常用的計(jì)算方法。有限元法在早期是以變分原理為基礎(chǔ)發(fā)展起來的，所以它廣泛地應(yīng)用于以拉普拉斯方程和泊松方程所描述的各類物理場(chǎng)中（這類場(chǎng)與泛函的極值問題有著緊密的聯(lián)系）。自從1969年以來，某些學(xué)者在流體力學(xué)中應(yīng)用加權(quán)余數(shù)法中的迦遼金法（6她「麻）或最小二乘法等同樣獲得了有限元方程，因而有限元法可應(yīng)用于以任何微分方程所描述的各類物理場(chǎng)中，而不再要求這類物理場(chǎng)和泛函的極值問題有所聯(lián)系。其原理是將連續(xù)的求解域離散為一組單元的組合體，用在每個(gè)單元假設(shè)的近似函數(shù)來分片的表示求解域上待求的未知場(chǎng)函數(shù)，近似函數(shù)通常由未知場(chǎng)函數(shù)及其導(dǎo)數(shù)在單元各節(jié)點(diǎn)的數(shù)值插值函數(shù)來表達(dá)。從而使一個(gè)連續(xù)的無限自由度問題變成離散的有限自由度問題。有限元法作為目前工程中應(yīng)用較為廣泛的一種數(shù)值計(jì)算方法，以其獨(dú)有的計(jì)算優(yōu)勢(shì)得到了廣泛的發(fā)展和應(yīng)用，ANSYS軟件以其多物理場(chǎng)耦合分析功能而成為CAE（工程設(shè)計(jì)中的計(jì)算機(jī)輔助工程）軟件的應(yīng)用主流，在熱分析工程應(yīng)用中得到了較為廣泛的應(yīng)用。ANSYS簡(jiǎn)介ANSYS程序是融結(jié)構(gòu)、熱、流體、電磁、聲學(xué)于一體的大型通用有限元商用分析軟件，可廣泛應(yīng)用于核工業(yè)、鐵道、石油化工、航空航天、機(jī)械制造、能源、電子、造船、汽車交通、國(guó)防軍工、土木工程、生物醫(yī)學(xué)、輕工、地礦、水利、日用家電等各個(gè)領(lǐng)域。該軟件可在大多數(shù)計(jì)算機(jī)及操作系統(tǒng)中運(yùn)行，從PC到工作站直至巨型計(jì)算機(jī)，ANSYS文件在其所有的產(chǎn)品系列和工作平臺(tái)上均兼容；該軟件基于Motif的菜單系統(tǒng)使用戶能夠通過對(duì)話框、下拉式菜單和子菜單進(jìn)行數(shù)據(jù)輸入和功能選擇，此舉大大方便了用戶操作。它由世界上著名的有限元分析軟件公司美國(guó)ANSYS公司開發(fā)的，能與多數(shù)CAD軟件配合使用，實(shí)現(xiàn)數(shù)據(jù)的共享和交換，是現(xiàn)代產(chǎn)品設(shè)計(jì)中的高級(jí)CAD工具之一。在ANSYS公司相繼收購ICEM、CENTURYDYNAMICS、AAVIDTHERMAL、FLUENT等世界著名有限元分析程序制造公司并將產(chǎn)品整合之后，ANSYS實(shí)際上已成為世界上最通用和有效的商用有限元軟件。ANSYS發(fā)展過程ANSYS公司是由美國(guó)匹茲堡大學(xué)的JohnSwanson博士在1970年創(chuàng)建的，其總部位于美國(guó)賓夕法尼亞州的匹茲堡，目前是世界CAE行業(yè)最大的公司。自創(chuàng)建以來，ANSYS公司一直致力于設(shè)計(jì)分析軟件的開發(fā)，不斷吸取新的計(jì)算方法和計(jì)算技術(shù)，領(lǐng)導(dǎo)著世界有限元技術(shù)的發(fā)展。ANSYS程序的最初版本與本文用的版本ANSYS10.0相比有很大區(qū)別。最初版本僅僅提供了熱分析和線性分析功能，是一個(gè)批處理程序，而且只能在大型計(jì)算機(jī)上使用。20世紀(jì)70年代初，隨著非線性、子結(jié)構(gòu)以及更多的單元類型的加入。ANSYS程序發(fā)生了很大的變化，新技術(shù)的融入進(jìn)一步滿足了用戶的需求。20世紀(jì)70年代末。交互方式的加入是該軟件最為顯著的變化，此舉使得模型生成和結(jié)果評(píng)價(jià)大為簡(jiǎn)化。ANSYS技術(shù)特點(diǎn)與其他有限元計(jì)算軟件相比，ANSYS具有以下技術(shù)特征：能實(shí)現(xiàn)多場(chǎng)及多場(chǎng)耦合功能。集前后處理、分析求解及多場(chǎng)分析于一體。獨(dú)一無二的優(yōu)化功能，唯一具有流場(chǎng)優(yōu)化功能的CFD軟件。具有強(qiáng)大的非線性分析功能。具備快速求解器。最早采用并行計(jì)算技術(shù)的FEA軟件。支持從微機(jī)、工作站到巨型機(jī)的所有硬件平臺(tái)。從微機(jī)、工作站、大型機(jī)直至巨型機(jī)，與所有硬件平臺(tái)的全部數(shù)據(jù)文件兼容。從微機(jī)、工作站、大型機(jī)直至巨型機(jī)，所有硬件平臺(tái)具有統(tǒng)一用戶界面?？膳c大多數(shù)的CAD軟件集成并具備接口?？蛇M(jìn)行智能網(wǎng)格劃分。具有多層次多框架的產(chǎn)品系列。具備良好的用戶開發(fā)環(huán)境。ANSYS使用環(huán)境ANSYS程序是一個(gè)功能強(qiáng)大的有限元計(jì)算分析軟件包。它可運(yùn)行于PC、NT工作站、UNIX工作站以及巨型計(jì)算機(jī)等各類計(jì)算機(jī)及操作系統(tǒng)中，其數(shù)據(jù)文件在其所有的產(chǎn)品系列和工作平臺(tái)上均兼容。該軟件有多種不同版本，本文應(yīng)用版本為ANSYS10.0版，其微機(jī)版本要求的操作系統(tǒng)為Windows95/98/2000/XP,也可運(yùn)行與UNIX系統(tǒng)中；微機(jī)版的基本硬件要求為：顯示分辨率為1024*768像素，顯示存為2MB以上，硬盤大于350乂8,推薦使用17英寸顯示器。ANSYS多物理場(chǎng)耦合的功能，允許在同一模型上進(jìn)行各種耦合計(jì)算，例如：熱-結(jié)構(gòu)耦合、熱-電耦合、磁-結(jié)構(gòu)耦合以及熱-電-磁-流體耦合，同時(shí)在PC上生成的模型可運(yùn)行于工作站及巨型計(jì)算機(jī)上，所有這一切就保證了ANSYS用戶對(duì)多領(lǐng)域工程問題的求解。ANSYS可與多種先進(jìn)的CAD（如AutoCAD、Pro/Engineer、NASTRAN、Alogor、I-DEAS等）軟件共享數(shù)據(jù)，利用ANSYS的數(shù)據(jù)接口，可以精確地將在CAD系統(tǒng)下生成的幾何模型數(shù)據(jù)傳輸?shù)紸NSYS，并通過必要的修補(bǔ)可準(zhǔn)確地在模型上劃分網(wǎng)格并進(jìn)行求解，這樣就可以節(jié)省用戶在創(chuàng)建模型的過程中所花費(fèi)的大量時(shí)間，使用戶的工作效率大幅度提高。ANSYS程序功能ANSYS程序主要包括3個(gè)部分網(wǎng)］：前處理模塊、求解模塊和后處理模塊。前處理模塊提供了一個(gè)強(qiáng)大的實(shí)體建模及網(wǎng)格劃分工具，用戶可以方便地構(gòu)造有限元模型；求解模塊包括結(jié)構(gòu)分析（結(jié)構(gòu)線性分析、結(jié)構(gòu)非線性分析和結(jié)構(gòu)高度非線性分析）、熱分析、流體動(dòng)力學(xué)分析、電磁場(chǎng)分析、聲場(chǎng)分析、壓電分析以及多物理場(chǎng)的耦合分析，可模擬多種物理介質(zhì)的相互作用，具有靈敏度分析及優(yōu)化分析能力；后處理模塊可將計(jì)算結(jié)果以彩色等值線顯示、梯度顯示、矢量顯示、粒子流顯示、立體切片顯示、透明及半透明顯示等圖形方式顯示出來，也可將計(jì)算結(jié)果以圖表、曲線形式顯示或輸出。ANSYS程序提供了100種以上的單元類型，用來模擬實(shí)際工程中的各種結(jié)構(gòu)和材料。ANSYS熱分析概述在實(shí)際生產(chǎn)過程中，常常會(huì)遇到多種多樣的熱量傳遞問題：如計(jì)算某個(gè)系統(tǒng)或部件的溫度分布、熱量的獲取或損失、熱梯度、熱流密度、熱應(yīng)力、相變等。所涉及的部門包括：能源、化工、冶金、建筑、電子、航空航天、農(nóng)業(yè)、制冷、船舶等。以機(jī)械加工為例，往往需要估算和控制工件溫度場(chǎng)；分析不同工作條件下，不同材料及幾何形狀對(duì)溫度場(chǎng)變化的影響以及防止加工過程工件中缺陷的產(chǎn)生。ANSYS作為新穎的有限元分析軟件在熱分析問題方面具有強(qiáng)大的功能，而且界面友好，易于掌握。用戶可以隨心所欲地選擇圖形用戶界面方式(GUI)或命令流方式進(jìn)行計(jì)算。ANSYS的熱分析模塊一般包括ANSYS/Muliphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五種，其中ANSYS/FLOTRAN不含相變熱分析。ANSYS熱分析基于能量守恒原理的熱平衡方程，用有限元法計(jì)算物體部各節(jié)點(diǎn)的溫度，并導(dǎo)出其他熱物理參數(shù)。運(yùn)用ANSYS軟件可進(jìn)行熱傳導(dǎo)、熱對(duì)流、熱輻射、相變、熱應(yīng)力及接觸熱阻等問題的分析求解。ANSYS求解步驟ANSYS求解模型的過程中的主要步驟：①定義單元類型及選項(xiàng)，根據(jù)所分析的問題在單元庫中選擇分析問題的單元類型；②定義單元是常數(shù)，單元是常數(shù)是與特定單元有關(guān)的量；③定義材料屬性，在材料屬欄中選擇性的輸入熱導(dǎo)率、泊松比、密度等所需參數(shù)；④建立幾何模型；⑤定義網(wǎng)格并劃分參數(shù)；⑥網(wǎng)格劃分以創(chuàng)建的模型。ANSYS在后處理器中可以查看結(jié)果：①顯示等值線；②顯示分析結(jié)果數(shù)據(jù)；③對(duì)結(jié)果數(shù)據(jù)進(jìn)行計(jì)算和路徑操作；④繪制云圖。結(jié)合本文研究容及ANSYS10.0軟件的ANSYS/Thermal熱分析模塊功能，ANSYS10.0適用于本文研究容。3井筒熱損失問題描述以遼河油田某鉆井為研究對(duì)象，注汽時(shí)間T=5d,油套環(huán)空導(dǎo)熱系數(shù)2yt為0.023W/m-K,套管導(dǎo)熱系數(shù)2tg為43.2川/川小，水泥環(huán)導(dǎo)熱系數(shù)2cem為0.81W/m-兒地層導(dǎo)熱系數(shù)2cm為1.745川/m小，地層熱擴(kuò)散系數(shù)c為0.027川加，地層初始溫度1；40￡,蒸汽與隔熱油管壁對(duì)流換熱系數(shù)a為10kw/m2?oC。數(shù)學(xué)模型的建立隔熱油管管半徑r1為0.031m，外半徑「2為0.0365m；隔熱油管外管半徑r3為0.05015m,外半徑r4為0.05715m；套管半徑［為0.0809印，外半徑r6為0.0889印；水泥環(huán)外半徑r7為。1239江；地層半徑18為。2478m，詳見圖3-1井筒及周圍地層截面示意圖。隔熱八、、層水泥環(huán)圖3-1井筒及其周圍地層截面示意圖3.2.1能量平衡方程依據(jù)熱平衡原理，選取井筒的一微元流體單元體為研究對(duì)象，如圖3-2所示。圖3-2流體單元能量平衡模型當(dāng)蒸汽縱向流入或流出該單元體的控制面1、2時(shí)，攜帶的能量可表示為：E=m(h++gz).i=1,2(3-1)

i oi2i式中：E—蒸汽攜帶的能量，W;im—蒸汽流量，kg/s；oh—蒸汽焓值，J/kg;iv一蒸汽流速，m/s；iz—高度，m;ig—重力加速度，9.81m/S2。該流體單元的能量平衡表達(dá)式為：v2 v2Q+m?(h+j+g?z)=m?(h+—+g-z)(3-2)022 2 012 1式中：Q—單位時(shí)間的徑向散熱量，W。蒸汽的物性參數(shù)在井筒截面上是不變的，但不同截面處的物性參數(shù)是隨溫度變化的。由式(3-2)看出蒸汽沿井筒縱向的焓值在減少，蒸汽的能量是按遞減規(guī)律傳至地層的。3.2.2連續(xù)介質(zhì)的守恒方程一井筒流體的能量平衡方程井筒油管高溫蒸汽在注入過程中遵循熱力學(xué)第一規(guī)律，其能量平衡方程為：v2 v2H+gz+p—z—H+g(z-dz)+p zdz,+Q(3-3)z z2 z-dz z-dz2dz式中：H、H—z及z-dz處的焓，kJ；z z-dzvz、v-dz—z及z-dz處的速度，m/h；p、p—2及2^2處的流體密度，kg/m3；zz-dzQ一蒸汽散熱量，kJ。dz公式(3-3)表示包括蒸汽焓、勢(shì)能和動(dòng)能在2處的蒸汽能量與2+dz處的蒸汽能量及蒸汽散熱量的平衡方程。8.井筒流體的連續(xù)性方程高溫蒸汽在油管流動(dòng)時(shí)，質(zhì)量流量保持不變，因此滿足流體的連續(xù)性方程，即pAv—pAv—m(3-4)zzz z+dz.z+dz.z+dz. 0式中：A、A—z及z+dz處的橫截面積，m2；z z+dz.vz、vz+dz—z及z+dz處的速度，m/h；p、p—z及z+dz處的流體密度，kg/m3。z z+dz。井筒流體的動(dòng)量平衡方程把蒸汽在油管的流動(dòng)看作是兩相流的流動(dòng)，為描述其流動(dòng)特性，通常從兩相流壓力降的研究出發(fā)。管穩(wěn)態(tài)均質(zhì)兩相流滿足動(dòng)量平衡方程，即dP=pgdz-fpmvdz-pvdv(3-5)

m2dm1式中：dP—dz段的壓力降，MPa；p—蒸汽密度kg/m3；mf—摩阻系數(shù)；d—井筒徑，m。1加V2pgdz表示dz段的勢(shì)能變化，"d示dz段的摩擦損失，Pvdv表示dz段的動(dòng)能變m 2d m1化。也就是說管兩相流的壓力降是摩擦損失、勢(shì)能變化和動(dòng)能變化的綜合結(jié)果。D.油管濕蒸汽的焓由于濕蒸汽為汽水兩相混合物，所以它的焓可用相同溫度或壓力下的飽和水、飽和蒸汽及濕蒸汽的干度來表示，即h=hx+h(1-x)(3-6)gl式中：x—蒸汽干度；h—濕蒸汽的焓，kJ/kg；h—飽和蒸汽的焓，kJ/kg；gh—飽和水的焓，kJ/kg。l3.2.3傳熱過程分析根據(jù)熱力學(xué)［31］和傳熱學(xué)［32］知識(shí)對(duì)井筒傳熱進(jìn)行分析，由于水泥環(huán)外側(cè)的地層體積和熱容量比水泥環(huán)側(cè)的井筒大得多，水泥環(huán)與地層交界處的溫度t隨注入時(shí)間的推移緩慢上升。k表征熱量在地層中擴(kuò)散速率的參數(shù)一熱擴(kuò)散系數(shù)a的數(shù)量級(jí)僅為10-72m/s，因而t與未k受熱影響的遠(yuǎn)處地層溫度t的差值不大；相反，隔熱管蒸汽溫度t與t的差別較大，有一個(gè)ci fk數(shù)量級(jí)的差別。因而，可以假設(shè)水泥環(huán)側(cè)(簡(jiǎn)稱井筒)熱量傳遞是穩(wěn)態(tài)的，而水泥層外側(cè)(簡(jiǎn)稱地層)熱量傳遞則按非穩(wěn)態(tài)處理。一井筒穩(wěn)態(tài)傳熱井筒傳熱過程主要包括：蒸汽與隔熱油管壁之間的對(duì)流換熱、隔熱油管多孔隔熱材料中伴有熱輻射的導(dǎo)熱、油套管環(huán)空中對(duì)流、輻射同時(shí)作用的換熱以及套管、水泥環(huán)和地層的導(dǎo)熱等環(huán)節(jié)。對(duì)于注汽井一微元段dz的徑向熱損失，可表示為：

dQ-2nrdzk(t-1) (3-7)2fk式中：dQ—dz段井筒的徑向熱損失，W;r一隔熱油管外半徑,m；2k—以隔熱油管管外表面為基準(zhǔn)面積的總傳熱系數(shù)，W/(m2-k)。當(dāng)井筒隔熱管抽真空時(shí)，k可表示為:k-[k-[R+R+R123+ 1R4+1R5+R+R11(3-8)67式中：R1—蒸汽凝結(jié)換熱熱阻;R—隔熱油管熱阻；2R—隔熱管隔熱層換熱熱阻；3R—環(huán)空輻射換熱熱阻；4R—環(huán)空對(duì)流熱熱阻；5R—套管壁導(dǎo)熱熱阻；6R—水泥環(huán)導(dǎo)熱熱阻。78.地層非穩(wěn)態(tài)傳熱在井筒穩(wěn)態(tài)傳熱的同一微元段dz，地層非穩(wěn)態(tài)導(dǎo)熱方程如下：(3-9)d21 1dt 1dt(3-9) e-+ e- e

dr2rdradT初始條件T-0時(shí):t-t-t+aZ(3-10)

ecibT邊界條件工―0時(shí):匕-0(3-11)

drr-r(3-12)7式中：t—地層某處溫度，t-f(r,T),oC；ee九一地層導(dǎo)熱系數(shù)，W/(m^K);ea—地層擴(kuò)散系數(shù)，a6/Vc，m2/s。epe基本假設(shè)(1)從油管的蒸汽到地層的熱量傳遞為一維穩(wěn)態(tài)傳熱，地層為非穩(wěn)態(tài)導(dǎo)熱；(2)視油管的蒸汽流為一維均質(zhì)兩相流，且無相變及泄露現(xiàn)象；(3)忽略地層溫度梯度的影響?？刂品匠谈鶕?jù)傳熱基本理論，二維瞬態(tài)傳熱問題的控制方程為：巴+更」區(qū)(3-13)。x2dy2a沆式中：x,y分別指的是x和y方向；t是時(shí)間，s；a是熱擴(kuò)散率，m，s。邊界條件蒸汽與油管壁之間屬于對(duì)流換熱，所以油管壁為第三類邊界，應(yīng)滿足：at 、，一一—k—=(T—T)(3-14)

an frr在不考慮地溫梯度的情況下，地層溫度為定值，因此地層邊界處為第一類邊界條件，其他邊界均為絕熱邊界。3.3ANSYS模擬利用ANSYS10.0熱分析軟件，做出在不同隔熱油管視導(dǎo)熱系數(shù)和注汽溫度時(shí)的井筒溫度場(chǎng)分布情況。建立有限元模型(1)打開ANSYS10.0有限元分析軟件，過濾圖形用戶界面；(2)確定工作文件名76^lossinwell”,單位制默認(rèn)為國(guó)際單位制；(3)進(jìn)入PREP7前處理器，選擇plane55熱分析單元，并設(shè)置成軸對(duì)稱單元；(4)將已知物性參數(shù)定義為材料熱物性參數(shù)；(5)根據(jù)圖3-1建立有限元模型，如圖3-3；(6)采用自由網(wǎng)格劃分方法對(duì)模型進(jìn)行網(wǎng)格劃分，結(jié)果如圖3-4。施加載荷計(jì)算人定義分析類型定義新的分析類型為瞬態(tài)分析“1口門5?6鵬”。8.設(shè)置瞬態(tài)熱分析的初始條件由于初始溫度場(chǎng)是不均勻的且又是未知的，所以首先要進(jìn)行穩(wěn)態(tài)熱分析，為穩(wěn)態(tài)熱分析建立初始條件：在井筒壁處施加對(duì)流載荷，在地層處施加原始地層溫度場(chǎng)載荷，其他邊界為絕熱邊界。圖3-3井筒有限元模型E!JJj￡TSANw1mt]-9:33:93：圖3-4網(wǎng)格劃分結(jié)果求解采用GUI命令：Main|Solution|LoadStepOpts|TimeFrequenc|Time-Time51?口，設(shè)置計(jì)算終止時(shí)間為0.01s,時(shí)間步長(zhǎng)為0.0公，并將時(shí)間積分效果設(shè)定為?！?對(duì)模型進(jìn)行穩(wěn)態(tài)熱分析。穩(wěn)態(tài)熱分析結(jié)束后，將時(shí)間積分效果設(shè)定為On,并刪除所有節(jié)點(diǎn)的溫度載荷，計(jì)算終止時(shí)間設(shè)定為5天，時(shí)間步長(zhǎng)設(shè)定為1小時(shí)，對(duì)模型進(jìn)行瞬態(tài)熱分析。后處理對(duì)于瞬態(tài)熱分析，ANSYS提供了兩種后處理方式，即POST1和POST26。本文采用POST1進(jìn)行后處理，利用GUI命令：MainMenu|GeneralPostproc|PlotResults顯示溫度場(chǎng)彩色云圖,并拾取套管徑邊緣處溫度。3.4模擬結(jié)果與分析取水泥環(huán)與地層之間截面的平均溫度作為溫度參數(shù),并結(jié)合文獻(xiàn)［33］中給出的熱損失經(jīng)驗(yàn)公式,分別算出在不同的隔熱油管視導(dǎo)熱系數(shù)及注汽溫度下井筒熱損失情況。隔熱油管視導(dǎo)熱系數(shù)對(duì)井筒熱損失的影響現(xiàn)在使用的隔熱油管多為抽真空隔熱油管，結(jié)合預(yù)應(yīng)力隔熱油管標(biāo)準(zhǔn)（5丫/15324-94）中提出的視導(dǎo)熱系數(shù)概念,把隔熱管的隔熱層以導(dǎo)熱、對(duì)流和輻射三種方式傳遞的熱量視為以純導(dǎo)熱方式傳遞的熱量。其中,視導(dǎo)熱系數(shù)數(shù)值劃分為五個(gè)圍等級(jí),如表3-1所示,因此,分別選取視導(dǎo)熱系數(shù)為0.07、0.05、0.03、0.01和0.005W/m-K。表3-1隔熱油管等級(jí)分類單位：W/m^K隔熱性能等級(jí)ABCDE視導(dǎo)熱系數(shù)0.06<k<0.080.04<k<0.060.02<人<0.040.006<卜<0.020.002<卜<0.006在實(shí)際的稠油熱采過程中,注汽溫度一般為300~350℃之間,基于此,本文以蒸汽溫度為312℃情況下進(jìn)行模擬,圖3-5為在視導(dǎo)熱系數(shù)為0.07時(shí)的溫度分布云圖,結(jié)合熱損失經(jīng)驗(yàn)公式,得到單位長(zhǎng)度井筒熱損失與視導(dǎo)熱系數(shù)之間的關(guān)系,如表3-2和圖3-6所示。

■m10.0EHaMU時(shí)CfiLMLCZICtfrrEP=iJOB=11?57D€E-15ZW1ZZHTI1VG-.IELfXS=DEH^ezarjfiLcaEBUTEI=1MTRE5KU七SHI-ID5KC=a&t.ti44E-B凡9%M值“BLJlit.tifiB15&.5！?MlEt.tl?LJ=TLON圖3-6隔熱油管視導(dǎo)熱系數(shù)對(duì)井筒熱損失的影響圖3-5視導(dǎo)熱系數(shù)為0.07W/m圖3-6隔熱油管視導(dǎo)熱系數(shù)對(duì)井筒熱損失的影響表3-2不同視導(dǎo)熱系數(shù)下的計(jì)算結(jié)果序號(hào)視導(dǎo)熱系數(shù)（川向水）套管溫度（￡）熱損失（W/m）10.00541.5882.71520.0142,8484.87030.0346,04610,33840.0547,79713,33250.0748,90215,222根據(jù)圖3-6分析可知，井筒單位長(zhǎng)度熱損失整體上隨著隔熱油管視導(dǎo)熱系數(shù)的增大而增大。其中，對(duì)于4級(jí)的隔熱油管，單位長(zhǎng)度井筒熱損失約達(dá)15W/m左右，而對(duì)于￡級(jí)的隔熱油管，單位長(zhǎng)度井筒熱損失僅約為3.5W/m,兩者單位長(zhǎng)度井筒熱損失相差約11.5W/m,因此，選用視導(dǎo)熱系數(shù)高級(jí)別的隔熱油管對(duì)于降低井筒單位長(zhǎng)度的熱損失具有明顯的影響。此外，當(dāng)視導(dǎo)熱系數(shù)從0.07川佃水提高至0.03川向小后，單位長(zhǎng)度熱損失的變化速率隨著視導(dǎo)熱系數(shù)的減小而出現(xiàn)明顯的降低。綜上所述，在實(shí)際的應(yīng)用過程中，建議選用隔熱油管視導(dǎo)熱系數(shù)在0.03~0.005W/m?Ko注汽溫度對(duì)單位長(zhǎng)度井筒熱損失的影響上述視導(dǎo)熱系數(shù)的建議圍，在視導(dǎo)熱系數(shù)為0.03川佃水情況下，分別研究以注汽溫度分別為312、324、333、346和350oC五種情況下，注汽溫度對(duì)單位長(zhǎng)度井筒熱損失的影響，詳細(xì)情況如表3-3和圖3-7所示，圖3-8為當(dāng)注汽溫度為324oC時(shí)的溫度分布云圖。表3-3不同注汽溫度下的計(jì)算結(jié)果序號(hào)注汽溫度（。0套管溫度（。0熱損失（w/m）131248.90215.222232449.29715.897333349.59216.401434649.98717.077535050.08217.239根據(jù)上圖分析所知，單位長(zhǎng)度井筒的熱損失隨著注汽溫度的升高幾乎呈線性增大的趨勢(shì)。在注汽溫度為312oC情況下，單位長(zhǎng)度井筒熱損失約為15.2川佃，在注汽溫度在350oC情況下，單位長(zhǎng)度井筒熱損失約為17.2W/m,由此可見，注汽溫度的提高對(duì)于井筒單位長(zhǎng)度熱損失的影響較小，每提高50。*井筒單位長(zhǎng)度熱損失相差約2W/m。注氣溫度代圖3-7注汽溫度對(duì)井筒熱損失的影響Id.?jujt5注氣溫度代圖3-7注汽溫度對(duì)井筒熱損失的影響Id.?jujt52d12■35:54;39NOIWjgOLrrriOKI37E.P-=1SOB-12DT：3G-3223JCi■呂B呂二■an?1.3231?3,6461E5.Z911?6.227.33^259.2€29C'.5G2321.905圖3-8注汽溫度為3240c時(shí)的溫度分布云圖4油層熱影響區(qū)域問題描述高溫高壓的蒸汽流經(jīng)井筒后進(jìn)入封隔器以下的油層區(qū)，通過射孔與油層進(jìn)行傳質(zhì)傳熱，傳熱過程如圖4-1所示。一一環(huán)形空間-『陶器#出魏婷 油窟嵬，孔段,熱水帶nniinntoi冷油帶下部地層一一環(huán)形空間-『陶器#出魏婷 油窟嵬，孔段,熱水帶nniinntoi冷油帶下部地層圖4-1油層傳熱過程示意圖以遼河油田某鉆井為研究對(duì)象，注汽時(shí)間T=10d,套管導(dǎo)熱系數(shù)atg為43.2川/m小，水泥環(huán)導(dǎo)熱系數(shù)acem為1.2W/m小，油層綜合導(dǎo)熱系數(shù)acm為2.53W/m-七地層導(dǎo)熱系數(shù)為1.968川/四小，地層熱擴(kuò)散系數(shù)c為0.027川辦，地層初始溫度T前=40￡,蒸汽與隔熱油管壁對(duì)流換熱系數(shù)a為10kw/m2?oC。數(shù)學(xué)模型的建立套管半徑r為0.0809印，外半徑r為0.0889??；水泥環(huán)外半徑r為0.1239印；地層和油層半徑「4為20m,詳見圖4-2油層截面示意圖。圖4-2油層截面示意圖連續(xù)介質(zhì)的守恒方程求解與流體中動(dòng)量、熱量和質(zhì)量傳遞有關(guān)的問題時(shí),首先應(yīng)根據(jù)基本的守恒定律對(duì)該物理系統(tǒng)建立一組方程。這些方程包括:①能量平衡方程；②連續(xù)性方程；③焓方程。A.連續(xù)介質(zhì)的能量平衡方程井筒油管高溫蒸汽在注入過程中遵循熱力學(xué)第一規(guī)律,其能量平衡方程為：v2 v2H+gz+p—=H+g(z-dz)+p —z-dz+Q(4—1)z z2 z-dz z-dz2dz式中：Hz、Hz-d—z及z—dz處的焓，kJ；vz、vzdz—z及z-dz處的速度，m/h；p、p—2及2^2處的流體密度，kg/m3；zz-dzQ一蒸汽散熱量，kJ。dz上式表示包括蒸汽焓、勢(shì)能和動(dòng)能在2處的蒸汽能量與2+dz處的蒸汽能量及蒸汽散熱量的平衡方程。B.連續(xù)介質(zhì)的連續(xù)性方程高溫蒸汽在油管流動(dòng)時(shí)，質(zhì)量流量保持不變，因此滿足流體的連續(xù)性方程，即pAv=pAv=m(4-2)zzz z+dz.z+dz.z+dz. 0式中：A、A—z及z+dz處的橫截面積，m2；z z+dz.vz、vz+dz—z及z+dz處的速度，m/h；p、p—z及z+dz處的流體密度，kg/m3。z z+dzC.連續(xù)介質(zhì)的動(dòng)量平衡方程把蒸汽在油管的流動(dòng)看作是兩相流的流動(dòng)，為描述其流動(dòng)特性，通常從兩相流壓力降的研究出發(fā)。管穩(wěn)態(tài)均質(zhì)兩相流滿足動(dòng)量平衡方程，即dP=pgdz-fpm^dz-pvdv(4-3)

m 2d m1式中：dP—dz段的壓力降，MPa；p—蒸汽密度kg/m3；mf—摩阻系數(shù)；d—井筒徑，m。1fpv2pgdz表示dz段的勢(shì)能變化，匕一超小dz段的摩擦損失，pvdv表示dz段的動(dòng)能變m 2d m1化。D.油管濕蒸汽的焓由于濕蒸汽為汽水兩相混合物，所以它的焓可用相同溫度或壓力下的飽和水、飽和蒸汽及濕蒸汽的干度來表示，即h=hx+h(1-x)(4-4)

gl式中：X—蒸汽干度；h—濕蒸汽的焓，kJ/kg；h—飽和蒸汽的焓，kJ/kg；gh—飽和水的焓，kJ/kg。l4.2.2多孔介質(zhì)的守恒方程一多孔介質(zhì)的質(zhì)量守恒方程多孔介質(zhì)的質(zhì)量守恒方程為：史+1vpv=0(4-5)dt①式中：p—密度，kg/m3；①一介質(zhì)的孔隙度;V—滲流速度，m/d。8.多孔介質(zhì)的動(dòng)力學(xué)方程@.單相流在只有飽和牛頓流體的空隙介質(zhì)中，達(dá)西定律如下：V=-—(Vp-pgvD)(4-6)

N式中：p一流體的滲流速度，m/d;k—滲透率， 2;N一粘度，Pa^s；P一流體的密度，kg/m3；p一壓力，MPa；D一垂直高度，m；G—重力加速度，m/S2。團(tuán)多相流多相流的滲流微分方程為：kk av-[—4(Vp-pgVD)]+q=—(ps6 (4-7)從ii iatiil式中：5—1相流體的飽和度；lql一流體流入或流出的流體質(zhì)量流量，kg/s。4.2.3傳熱過程分析一井筒穩(wěn)態(tài)傳熱井筒傳熱過程主要包括：蒸汽與套管壁之間的對(duì)流換熱、套管與水泥環(huán)和水泥環(huán)與油層的導(dǎo)熱等環(huán)節(jié)。對(duì)于注汽井一微元段dz的徑向熱流量，可表示為：dQ=2兀rdzk(t-1) (4-8)2 fk式中：dQ—dz段井筒的徑向熱流量，W;r—套管外半徑,m；2k—以套管外表面為基準(zhǔn)面積的總傳熱系數(shù)，W/(m2-k)。k的表達(dá)式為：k=[R+R+R卜i(4-9)123式中：R—蒸汽與套管壁對(duì)流換熱熱阻；1R—套管壁導(dǎo)熱熱阻；2R—水泥環(huán)導(dǎo)熱熱阻。3B.油層傳熱分析井筒蒸汽與套管通過對(duì)流換熱向井筒周圍散熱，同時(shí)產(chǎn)生從射孔到油層的傳質(zhì)熱傳導(dǎo)。由于在有限元中無法解決傳質(zhì)問題，因此一般簡(jiǎn)化為具有一定半徑的加熱區(qū)模型。非穩(wěn)態(tài)無熱源軸對(duì)稱溫度場(chǎng)的微分方程為：D[T(x,r,t)]=kJr2+r巴+二]-pc三=0(4-10)ax2 ar2ar pat應(yīng)用6^「師法可得到軸對(duì)稱溫度場(chǎng)整體區(qū)域“變分”計(jì)算的基本方程為:kr+”竺kr+”竺+”三Sxdxx drdr+pcWr空plStdxdr』kWr支dS=0iSnT(4-11)利用溫度場(chǎng)變分方程對(duì)模型進(jìn)行剖分離散、單元分析及最終總體合成，即可完成模型溫度場(chǎng)的求解。4.2.4基本假設(shè)①沿井深方向上的熱能不發(fā)生傳遞；②視蒸汽流為一維均質(zhì)兩相流，且無相變；③忽略油層溫度梯度的影響；④以井眼中心為中心的同心圓上的溫度是相同的。4.2.5邊界條件蒸汽與套管壁之間屬于對(duì)流換熱，所以套管壁為第三類邊界，應(yīng)滿足：ST 、，一—k一=(T—T)(4-12)Sn frr加熱區(qū)為第一類邊界條件，由于不受油層溫度梯度的影響，所以油層為第一類邊界條件，其他為絕熱邊界。4.3ANSYS模擬建立有限元模型①打開ANSYS10.0有限元分析軟件，過濾圖形用戶界面；②確定工作文件名“youcengreyingxiangquyu”,單位制默認(rèn)為國(guó)際單位制；③進(jìn)入PREP7前處理器，選擇plane55熱分析單元，并設(shè)置成軸對(duì)稱單元；④將已知物性參數(shù)定義為材料熱物性參數(shù)；⑤根據(jù)問題描述建立有限元模型，結(jié)果如圖4-3所示；

圖4-3油層有限元模型（局部放大圖）⑥采用自由網(wǎng)格劃分方法對(duì)模型進(jìn)行網(wǎng)格劃分，結(jié)果如圖4-4所示。圖4-4網(wǎng)格劃分結(jié)果（局部放大圖）施加載荷計(jì)算本定義分析類型定義新的分析類型為瞬態(tài)分析機(jī)丁15?6鵬”。8.設(shè)置瞬態(tài)熱分析的初始條件由于初始溫度場(chǎng)是不均勻的且又是未知的，所以首先要進(jìn)行穩(wěn)態(tài)熱分析，為穩(wěn)態(tài)熱分析

建立初始條件：在套管壁處施加對(duì)流載荷，在油層和地層處施加原始溫度場(chǎng)載荷，其他邊界為絕熱邊界。求解采用GUI命令：Main|Solution|LoadStepOpts|TimeFrequenc|Time-Time51?口，設(shè)置計(jì)算終止時(shí)間為0.01s,時(shí)間步長(zhǎng)為0.0公，并將時(shí)間積分效果設(shè)定為。￡3對(duì)模型進(jìn)行穩(wěn)態(tài)熱分析。穩(wěn)態(tài)熱分析結(jié)束后，將時(shí)間積分效果設(shè)定為On,并刪除所有節(jié)點(diǎn)的溫度載荷，計(jì)算終止時(shí)間設(shè)定為10天，時(shí)間步長(zhǎng)設(shè)定為1小時(shí)，對(duì)模型進(jìn)行瞬態(tài)熱分析。后處理采用POST1進(jìn)行后處理，利用GUI命令：MainMenu|GeneralPostproc|PlotResults顯示溫度場(chǎng)彩色云圖。4.4模擬結(jié)果與分析模擬油層溫度場(chǎng)分布云圖如圖4-5至圖4-14所示?！鯩ik寫ID..'E1￡?￡291191D9E?KkiL30147710+1"：H==B?46dTEMP|liV5-.iRna-oEFJUCEI=i『■Mik寫ID..'E1￡?￡291191D9E?KkiL30147710+1"：H==B?46dTEMP|liV5-.iRna-oEFJUCEI=i『黑?呂口呂呂■-34B.F1S削?tM74,&111DB.?E1143.^3177.T29212.1197日口■聒敦。/*白NC-ZXIkL?LOTICO7IHZ=li3a'MZEMF|AVG|HSTS=G籃S呂呂口二Euve=芯邛SEJCEr-lfiVKS=Hai：=4籃S呂呂口二74..LC9.246Ida.5071.17.M?EL2.2S7H72加，的日315..MB手，口.￡DB圖4-5注汽第1天溫度分布云圖圖4-6注汽第2天溫度分布云圖aj?wfllo.4JVQ1NKN20=21:46HODALl]hG-：l5E^aIEM? |MFGIFdiU4313rB^lEUEF>￡ET-=1ma-1L.?51-3iqJ43U.JM4T"卻白140.GQ2iia.i7T170.2712L3.9fiG214.942HJGB3⑶傳319.F3Mns!i0.41ML/JOjJDj43KCZiiEimum工支Mmc=E:55M0?￡iifXEW=^EHpji-=ABFMblXTOSUbt.Sa-4d.:-95UE=S4fl_72fTT9.3EH109.lEtlie.l￡s￡：3--14；144.71S514H.7IT圖4-7注汽第3天溫度分布云圖圖4-8注汽第4天溫度分布云圖建;AMSXS1U.UJVH12C1220:23:10HMftLSOLUIICtfiike=4320￡iDTEMP5聞附￥3?QF=Ha=C=aphiEaEEMEI-lA'inU：5=HJitSOI-41.39SSOC=^43.75541.2u=73.5S-3=1M白恥|1-3-9.112■1陽型3售212.LJ3**皆相LJ2ai.2O7I^4fi.755M3T51U.Z?EHjQ址I9CUL5CIZTTIM=ee=EMMZ*WIA.7CIPSHI-E;■?1C31IZRIX：=1"=41.4U?3jg_7C341.4=-315.TM]n?,M6：-9-9.223*日■小212-TJG2-96.BS3洲1,亞3115?WE三■■□呂呂.圖4-9注汽第5天溫度分布云圖圖4-10注汽第6天溫度分布云圖AJ用與ICiQJWJ2U1Z3Di2^li39NEOKLSCUITjCG門您二助時(shí)時(shí)IEMT fin畫，PCMe-rtZBphlOSEF*￡E>1=41.501=543.7￡t41.5S113.B32ho.ptj144.317170.55B:iia.m247.00H1..3IS：13.32739369MTTJmsE12D=2GiDOHDIXU-3C1￡C[C?IIME-691KOTEMPIA7C)unrz?QEaw?EGs 2.CsE?*CEI-1-AJFEirHm-=4i.tlSTM.TTS41.1IS：73.M2IIO.ITL176.￡^a212.￡Ba317.Oft?ZEI1.31G315.M3aM.TTS圖4-11注汽第7天溫度分布云圖圖4-12注汽第8天溫度分布云圖■M1ZOLZ￡6=2￡=41TEHf lira卜R5T5=0Iwtrirafiii05EFJUCEr=lA7523-HW==41_114TE.-M1im，1TS.E922HT-LZ6ati.SiSnx*5Mt.TM；AH3坨10?0JUN1制;I戈2Ql37l^ZHEDAZ￡O1ZTTjOQZIME-E.F4OKIEHFU3es=0F口黃虺h恥工。、EHU：EI=LAVREfl-BiT.ELKH-11.5255MX=246-Ul7E.LS1119-S3T1H-M1114-TS21T.1E1S15-&1634P-TS1圖4-13注汽第9天溫度分布云圖圖4-14注汽第10天溫度分布云圖油層熱影響區(qū)域隨時(shí)間的變化如圖4-15。17617.417.217.016.816.616.416.21 2 5 4 5 6 7 6 9 10 11 12注汽時(shí)間，d圖4-15油層熱影響區(qū)域根據(jù)以上各圖分析可知，油層高溫區(qū)從第3天開始逐漸向地層擴(kuò)散，造成熱損失。油層熱影響區(qū)域先隨著注汽時(shí)間的增大而增大，大約以0.2m/d的速度增大，第5天熱影響區(qū)域達(dá)到最大，而后熱影響區(qū)域基本保持不變，維持在大約17.4米左右。5結(jié)論4井筒熱損失（1）井筒單位長(zhǎng)度熱損失整體上隨著隔熱油管視導(dǎo)熱系數(shù)的增大而增大。視導(dǎo)熱系數(shù)從0.07W/m-K提高至0.03W佃水后，單位長(zhǎng)度熱損失的變化速率隨著視導(dǎo)熱系數(shù)的減小而出現(xiàn)明顯的降低，建議選用隔熱油管視導(dǎo)熱系數(shù)在0.03~0.005W/m^K;（2）單位長(zhǎng)度井筒的熱損失隨著注汽溫度的升高幾乎呈線性增大的趨勢(shì)。每提高50℃，井筒單位長(zhǎng)度熱損失相差約2W/m；（3）視導(dǎo)熱系數(shù)對(duì)井筒單位長(zhǎng)度熱損失的影響大于注汽溫度的影響。BUM層熱影響區(qū)域（1）高溫區(qū)從第3天開始逐漸向地層擴(kuò)散，造成熱損失；（2）油層熱影響區(qū)域先隨著注汽時(shí)間的增大而逐漸增大，在第5天達(dá)到最大，而后熱影響區(qū)域基本保持不變，大約在17.4米左右；（3）綜合1、2兩點(diǎn)，建議注汽時(shí)間在3~5天。參考文獻(xiàn)[1]RameyH.J.WellboreHeatTransmission[J].JPT,1962,14（4）:427-435[2]HasanAR.KabirCS.HeatTransferDuringTwo-PhaseFlowinWellboresPartI---FormationTemperature[A].SPE22866,1991HasanAR.KabirCS.HeatTransferDuringTwo-PhaseFlowinWellboresPareII---WellboreFluidTemperature[A].SPE22948,1991HasanAR.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,1994FarouqAliSM.AprehensiveWellboreSteam/WaterFlowModelforSteamInjectionandGeothermalApplications[A].SPE527,1981SetterA.HeatlossesDuringFlowofSteamDownaWellbore[J].JPT,1965,17（7）:845-851SquierDP,SmithDD,DoughertyEL.CalculatedTemperatureBehaviorofHot-WaterInjectionWells[J].JPT,1962,（4）:436-440HuygenHH.HuittJL.WellboreHeatLossesandCasingTemperaturesDuringSteamInjection[J].Prod.Montlily,1966,30（8）:2GristionS.WillhiteGP.NumericalModelforConcentricSteamInjectionWells[