導(dǎo)語：如何才能寫好一篇流體的力學(xué)特征，這就需要搜集整理更多的資料和文獻(xiàn)，歡迎閱讀由公文云整理的十篇范文，供你借鑒。

篇1

【關(guān)鍵詞】  粒細(xì)胞; 流行病學(xué)研究; 感染

epidemiological characteristic of epidemic situation about the humaninfected granulocytic anaplasmosis in southern anhui cheng zhouxiang1, yang xiaoxiang1, li qun2, he jiangang1, dou zhengdong1, wu jiabing3, yu yanlin4, wang fei1. 1.wuhu center for disease control and prevention, wuhu 241000, china; 2.office for disease control and emergency response, chinese center for disease control and prevention, beijing 100050, china; 3.anhui center for disease control and prevention, hefei 230061, china; 4.yijishan hospital of wannan medical college, wuhu 241001, china

【abstract】 objective to describe the epidemiological and clinical characteristics of epidemic situation about the humaninfected granulocytic anaplasmosis. methods the whole cases were collected which had occurred in this epidemic. then the epidemiological and clinical characteristics were analyzed by descriptive epidemiology, compared with the clinical characteristic of previous tickinfected human granulocytic anaplasmosis. at the same time, the risk factor of infection was analyzed. results the cases had the typical clinical characteristic, such as hyperpyrexia, progressive falling of wbc and plt. the fever process was 1～6 days, the incubation period being 4～13 days and its median 7. the infection rate of the persons who had close contact with the first case during his or her illness was 14.3%. the attack cases among the persons who had close contact with the first case when he or she was in a critical condition were 23.1%. the blood sample from 8 secondary attack cases was tested positive about the 16s rrna gene of granulocytic anaplasma. the specific serum antibody was seroconversion in convalescence and its titer was 1∶64～1∶128. conclusions this is not only the first report about humaninfected granulocytic anaplasmosis but also the first one that granulocytic anaplasmosis can result in iatrogenic infection. there is even no infectiousness in the early light situation of disease, but the possibility of infection increases if people have close contact with patients when their illness is serious. therefore, isolation and protection have to be practiced.

【key words】 granulocytes; epidemiologic studies; infection

“人粒細(xì)胞無形體病”(human granulocytic anaplasmosis，hga)(以前稱“人粒細(xì)胞埃立克體病”[1])是一種人類新發(fā)傳染病。1994年在美國明尼蘇達(dá)州杜魯士首次發(fā)現(xiàn)一種不同于查菲埃立克體的感染，病原體主要侵犯中性粒細(xì)胞，被命名為人粒細(xì)胞埃立克體(human granulocytic ehrlichiosis,hge)[2]，此后美國[3]、斯洛文尼亞[4]、比利時(shí)[5]、澳大利亞[6]、意大利[7]、瑞典[8]、德國[9]等歐美國家都有病例報(bào)道或有流行病學(xué)資料證實(shí)該病原體感染。2001年高東旗等[10]在大興安嶺地區(qū)開展人群埃立克體感染的調(diào)查，結(jié)果首次在亞洲從人血中擴(kuò)增出人粒細(xì)胞埃立克體16s rrna基因片段。初步認(rèn)為我國大興安嶺地區(qū)存在人粒細(xì)胞埃立克體感染人群。但以上報(bào)道所有的病例或感染者都有蜱叮咬史，傳染源是野外動(dòng)物。2006年11月安徽省蕪湖市y醫(yī)院收治10例“人粒細(xì)胞無形體病”病例，由于其中后9例均是與首例病人有密切接觸的醫(yī)務(wù)人員或陪護(hù)親屬，且無蜱叮咬史等可能的傳播途徑，目前認(rèn)為這起疫情是全球首次發(fā)現(xiàn)的人傳“人粒細(xì)胞無形體病”疫情。使人類對(duì)人粒細(xì)胞無形體病的流行病學(xué)特征又有了進(jìn)一步了解和認(rèn)識(shí)?，F(xiàn)將這起疫情的流行病學(xué)和臨床特征做一描述。

1 對(duì)象與方法

1.1 病例臨床資料 來源于蕪湖市y醫(yī)院。

1.2 現(xiàn)場流行病學(xué)調(diào)查

1.2.1 病例定義 因急性發(fā)病，發(fā)熱(t≥38 ℃)，外周血白細(xì)胞(white blood cell,wbc)和血小板(platelet count,plt)均下降的病人，可合并多臟器功能損害，具有可疑動(dòng)物或首例病例接觸史并發(fā)病的病人，以及與這些病例在診療過程中接觸并具有以上類似病癥的病例，實(shí)驗(yàn)室檢驗(yàn)結(jié)果符合粒細(xì)胞埃立克體感染特征[11]。

1.2.2 個(gè)案調(diào)查 根據(jù)以上病例定義搜索病例，并對(duì)符合定義的10例病人或其家屬開展調(diào)查，調(diào)點(diǎn)內(nèi)容為發(fā)病前可能的暴露因素、發(fā)病經(jīng)過、診療過程、接觸者、接觸方式、接觸時(shí)間、個(gè)人防護(hù)情況、年齡、職業(yè)、生活環(huán)境和首例病人可能的動(dòng)物接觸史。

1.2.3 病例密切接觸者調(diào)查 內(nèi)容基本與病例相同。

1.2.4 病例資料描述性流行病學(xué)分析 發(fā)病時(shí)間分布，病例的性別、年齡、職業(yè)分布，地區(qū)分布。

1.3 樣品采集和實(shí)驗(yàn)室檢驗(yàn)

1.3.1 首例病人(死亡病例)血清 送安徽省疾病預(yù)防控制中心，免疫熒光法檢測流行性出血熱igg。

1.3.2 繼發(fā)病例和對(duì)照樣本 繼發(fā)病例病后每隔3 d采取抗凝血2.5 ml和不抗凝血5 ml，采取病例早期咽拭子樣3 ml，病程中后期尿液2 ml。密切接觸者血清標(biāo)本216人份(153人)，自然人群血清標(biāo)本40份(40人)，檢測特異抗體和核酸。送安徽省疾病預(yù)防控制中心和中國疾病預(yù)防控制中心進(jìn)行病原學(xué)檢測。

1.3.3 病例白細(xì)胞形態(tài)學(xué)檢查 吉姆薩染色觀察白細(xì)胞內(nèi)包涵體。

2 結(jié)果

2.1 臨床特征

2.1.1 臨床癥狀 所有病例均有畏寒和發(fā)熱，最高體溫都超過39 ℃，如果不進(jìn)行對(duì)癥處理，則高熱不退。大多體溫39.5 ℃，最高40 ℃。病例的熱程在1～6 d，中位數(shù)為4 d。大多數(shù)病例(7/10)有腹瀉。在發(fā)病后1～4 d出現(xiàn)，多為1～2 d,大多數(shù)每天腹瀉3～4次，多的有5～6次。見表1。1例進(jìn)展成急性呼吸窘迫綜合征(aucte respiratory distress syndrome,ards)，合并彌散性血管內(nèi)凝血(disseminated intravascular coagulation,dic)死亡。表1 皖南地區(qū)人粒細(xì)胞無形體病疫情臨床癥狀出現(xiàn)頻次

2.1.2 臨床檢驗(yàn)結(jié)果 所有病例均先出現(xiàn)血液系統(tǒng)的改變，wbc和plt均出現(xiàn)進(jìn)行性下降，wbc多為2.0×109/l～3.0×109/l、中性粒細(xì)胞比例減少多見 ，部分可見異形細(xì)胞，plt多為30×109/l左右;尿常規(guī)檢查4例尿蛋白陽性，可有血尿和顆粒管型;7例出現(xiàn)酶學(xué)改變，肌酶和肝酶升高，如乳酸脫氫酶、轉(zhuǎn)氨酶(alt和ast)異常，心肌酶譜升高(同工酶正常)等;腎功能異常出現(xiàn)較晚，且不是很嚴(yán)重(包括死亡病例)。

2.2 流行病學(xué)調(diào)查資料分析

2.2.1 發(fā)病時(shí)間分布 首例病人10月18日發(fā)現(xiàn)自己被野兔身上的草蜱咬傷右踝關(guān)節(jié)內(nèi)側(cè)，10月21日上午在家衣櫥中捕殺一窩老鼠，并清理和清洗部分被老鼠污染的衣物。10月30日晚20：00發(fā)病，10月31日在該村衛(wèi)生室就診，11月3日晚21:00到x縣人民醫(yī)院治療，11月4日中午轉(zhuǎn)往蕪湖市y醫(yī)院就診，初步診斷為“流行性出血熱(少尿期)”。11月5日凌晨患者存在多臟器衰竭，11月5日6：55死亡。

從11月9日起陸續(xù)出現(xiàn)新病例，在12日出現(xiàn)發(fā)病高峰(4例)，后逐漸下降，見圖1。

2.2.2 性別、年齡、職業(yè)分布 10例病例中，男性6例，女性4例。職業(yè)為：農(nóng)民5例，醫(yī)務(wù)人員4例，教師1例。年齡在25～67歲之間，其中20～年齡組4例，30～歲組3例，40～、50～、60～年齡組各1例。

2.2.3 地區(qū)分布 10例病例中，宣城市x縣新杭鎮(zhèn)6例，蕪湖市y醫(yī)院4例。

2.2.4 密切接觸者情況 經(jīng)宣城、蕪湖兩市疾病預(yù)防控制中心對(duì)所有密切接觸者進(jìn)行追蹤、隔離觀察及醫(yī)學(xué)觀察，結(jié)果共搜索與死亡病例(首例病人)病后有過密切接觸人員共63人，親屬21人，醫(yī)護(hù)人員42人(y醫(yī)院24人、x縣醫(yī)院16人、村醫(yī)2人)，其中9人陸續(xù)發(fā)病，4人為y醫(yī)護(hù)人員，5人為病人親屬，感染率14.3%。在y醫(yī)院有過密切接觸者39人中，發(fā)病9人，發(fā)病率為23.1%。與首例病人僅在病情尚不危重，未轉(zhuǎn)入y醫(yī)院之前密切接觸者24人沒有發(fā)現(xiàn)續(xù)發(fā)病例(χ2=4.7,p<0.05)。提示此階段病人不具備傳染性或傳染性不強(qiáng)，在y醫(yī)院與死亡病例危重時(shí)有過密切接觸是發(fā)病的危險(xiǎn)因素。

與9例續(xù)發(fā)病例的密切接觸人員有90人，其中親屬朋友48人，醫(yī)護(hù)人員42人(y醫(yī)院20人、x縣醫(yī)院22人)。所有密切接觸者經(jīng)醫(yī)學(xué)觀察14 d，未再發(fā)現(xiàn)類似病例。

2.2.5 潛伏期估計(jì) 首例病人10月18日被蜱叮咬,21日接觸老鼠,30日發(fā)病。因?yàn)槿肆＜?xì)胞無形體病通常是蜱傳疾病，所以推算潛伏期為12 d。

9例續(xù)發(fā)病例在y醫(yī)院與死亡病例危重時(shí)有過密切接觸時(shí)間在11月4日或5日，發(fā)病在此接觸后4～13 d之間，中位數(shù)為7 d，主要集中在7～9 d，占66.7%(6/9)，見圖2。

2.3 實(shí)驗(yàn)室結(jié)果

2.3.1 首例病人(死亡病例)血清流行性出血熱免疫熒光抗體igg(-)

2.3.2 繼發(fā)病例實(shí)驗(yàn)室結(jié)果

2.3.2.1 排除可能的病原體感染的有關(guān)檢驗(yàn) 11月15日～17日共5例病人血清標(biāo)本流行性出血熱igm抗體陰性(elisa法)。11月17日采集的7例患者血清除一人流感igg(+)外，檢測sars、流感、副流感、柯薩奇、呼吸道合胞病毒、腺病毒、支原體、衣原體igm、igg抗體，均為陰性。7例患者咽拭子樣本流感通用型、b型、h5核酸檢測均陰性。斑疹傷寒、斑點(diǎn)熱、恙蟲病核酸檢測均無特異條帶擴(kuò)增。 能引起出血熱的病毒如沙粒病毒屬、絲狀病毒(埃博拉病毒、馬爾堡病毒)、黃病毒、布尼亞病毒(漢坦病毒、cchf病毒、立夫特谷熱 )核酸和抗體均陰性。

2.3.2.2 病原體特異檢測 11月24日，9位患者早期(發(fā)病后3～5 d)血清抗粒細(xì)胞無形體igm(1∶20)、igg(1∶64)抗體檢測均為陰性，病后3 w對(duì)其中的5例病人恢復(fù)期血清進(jìn)行了特異抗體igg檢測，5例病人血清igg抗體均為陽性,其中3人陽性滴度1∶64，2人陽性滴度1∶128(呈4倍增高)。8位患者的血標(biāo)本的無形體和埃立克體屬16s rrna基因檢測陽性,pcr產(chǎn)物經(jīng)過測序分析驗(yàn)證與嗜粒細(xì)胞無形體16s rrna序列100%同源。白細(xì)胞形態(tài)學(xué)檢查，1例病人白細(xì)胞內(nèi)可見嗜粒細(xì)胞無形體桑椹包涵體。

3 討論

3.1 這起疫情是嗜粒細(xì)胞無形體感染引起的 從臨床特征上分析，10例病人均有急性起病、寒戰(zhàn)、高熱、體溫超過39 ℃、相對(duì)緩脈、乏力不適、外周血wbc計(jì)數(shù)下降、plt下降，部分患者出現(xiàn)腹瀉、肝功能損害、腎功能損害，少見皮疹。符合嗜粒細(xì)胞無形體病臨床特點(diǎn)[2，12，13]。實(shí)驗(yàn)室檢測排除斑疹傷寒、斑點(diǎn)熱、恙蟲病、能引起出血熱的一些病毒等病原體感染。繼發(fā)病例早期血清抗粒細(xì)胞無形體抗體陰性，恢復(fù)期轉(zhuǎn)為陽性,血標(biāo)本的嗜粒細(xì)胞無形體16s rrna基因檢測陽性,pcr產(chǎn)物經(jīng)過測序分析驗(yàn)證與嗜粒細(xì)胞無形體16s rrna序列100%同源。病人白細(xì)胞內(nèi)可見特異性桑椹包涵體。這些資料表明這些病人符合嗜粒細(xì)胞無形體病的診斷標(biāo)準(zhǔn)[11]。雖然首例病人在早期臨床上診斷為“流行性出血熱”，也沒有實(shí)驗(yàn)室結(jié)果支持為嗜粒細(xì)胞無形體病，但該患者有蜱叮咬史，潛伏期為12 d，符合該病潛伏期7～14 d(平均9 d)的特征[12～14]。其病程發(fā)展符合該病臨床特點(diǎn)。與其密切接觸者發(fā)生與其臨床特點(diǎn)相符的病癥?？梢酝茢嗍桌∪艘矐?yīng)該是嗜粒細(xì)胞無形體感染。

3.2 這起疫情是一次罕見的、首次發(fā)現(xiàn)的人傳嗜粒細(xì)胞無形體病疫情和醫(yī)源性感染暴發(fā) 在此以前，嗜粒細(xì)胞無形體病已確認(rèn)是蜱叮咬傳播的人畜共患的新發(fā)傳染病[1，8，13]。傳播媒介是硬蜱屬的一些蜱種[12，15]。在我國已經(jīng)報(bào)道過的嗜粒細(xì)胞無形體病(埃立克體病)病例或流行病學(xué)調(diào)查發(fā)現(xiàn)的感染者中，也都有蜱叮咬史[10，14]。皖南這起疫情的9例續(xù)發(fā)病例都沒有蜱叮咬史，也沒有其他野生動(dòng)物接觸史，只是與首例病人在醫(yī)院中有過密切接觸，且發(fā)生醫(yī)務(wù)人員感染。由此來看，嗜粒細(xì)胞無形體病的傳染源可以是病人，傳播途徑可以是與危重病人近距離接觸而非一定是蜱傳。

3.3 人傳嗜粒細(xì)胞無形體病臨床特征和潛伏期與蜱傳嗜粒細(xì)胞無形體病基本相同 這起疫情病人臨床表現(xiàn)和體征以及實(shí)驗(yàn)室檢驗(yàn)結(jié)果與以往報(bào)道蜱傳嗜粒細(xì)胞無形體病沒有特異性。病人潛伏期多為7～9 d，以往報(bào)道蜱傳嗜粒細(xì)胞無形體病潛伏期為7～14 d(平均9 d)。病原體不同侵入部位和方式所引起的病理過程是否有差異有待研究。與危重病人近距離接觸傳播最可能的侵入途徑是吸入或沾染粘膜和破損皮膚，而蜱叮咬多在四肢暴露部位。

3.4 病人作為嗜粒細(xì)胞無形體病的傳染源的傳染力大小因其處在不同病程階段而不同 這起疫情與死亡病例(首例病人)病后有過密切接觸人員感染率14.3%。但其中與首例病人僅在病情尚不危重，未轉(zhuǎn)入y醫(yī)院之前密切接觸者24人沒有發(fā)現(xiàn)續(xù)發(fā)病例。與病人在危重時(shí)有過密切接觸者發(fā)病率為23.1%。提示這類病人在危重時(shí)應(yīng)嚴(yán)格隔離，禁止病房內(nèi)陪護(hù)，醫(yī)務(wù)人員應(yīng)加強(qiáng)衛(wèi)生防護(hù)。

【參考文獻(xiàn)】

 

[1] oteo ja, brouqui p. ehrlichiosis and human anaplasmosis [j]. enferm infecc microbiol clin, 2005,23(6):375380.

[2] 曹務(wù)春. 人埃立克體病的發(fā)現(xiàn)與研究進(jìn)展 [j]. 中國人獸共患病雜志, 1997,13(4):5760.

[3] gardner sl, holman rc, krebs jw, et al. national surveillance for the human ehrlichioses in the united states, 19972001, and proposed methods for evaluation of data quality [j]. ann n y acad sci, 2003,990:8089.

[4] lotricfurlan s, petrovec m, avsiczupanc t, et al. human granulocytic ehrlichiosis in slovenia [j]. ann n y acad sci, 2003,990:279284.

[5] heyman p, cochez c, bigaignon g, et al. human granulocytic ehrlichiosis in belgium: an underestimated cause of disease [j]. j infect, 2003,47(2):129132.

[6] walder g, falkensammer b, aigner j, et al. first documented case of human granulocytic ehrlichiosis in austria [j]. wien klin wochenschr, 2003 ,115(78):263266.

[7] ruscio m, cinco m. human granulocytic ehrlichiosis in italy: first report on two confirmed cases [j]. ann n y acad sci， 2003，990:350352.

[8] strle f. human granulocytic ehrlichiosis in europe [j]. int j med microbiol, 2004，293:2735.

篇2

關(guān)鍵詞：油氣成藏動(dòng)力，學(xué)油氣運(yùn)移油，油氣成藏機(jī)理

 

1.油氣成藏動(dòng)力學(xué)研究方法

成藏動(dòng)力學(xué)研究是在綜合分析區(qū)域鉆探、地球物理、分析測試和地質(zhì)地化等資料的基礎(chǔ)上, 采用靜態(tài)描述和動(dòng)態(tài)模擬相結(jié)合的方法, 其中計(jì)算機(jī)模擬方法可以定量地、動(dòng)態(tài)地刻劃各種因素相互作用的歷史過程, 從而更深刻地揭示其內(nèi)在規(guī)律性, 因此是成藏動(dòng)力學(xué)過程研究的一項(xiàng)關(guān)鍵技術(shù)。成藏動(dòng)力學(xué)模擬實(shí)質(zhì)上是成藏動(dòng)力學(xué)過程模擬, 是一項(xiàng)高度復(fù)雜的系統(tǒng)工程, 它需要以當(dāng)代最先進(jìn)的地質(zhì)學(xué)和石油地質(zhì)學(xué)理論為基礎(chǔ), 全面利用各種地質(zhì)、物探資料, 采用最先進(jìn)的盆地描述和盆地模擬技術(shù)方可進(jìn)行[1]。，油氣成藏機(jī)理。盆地描述部分用于刻劃盆地現(xiàn)今的構(gòu)造、沉積巖性和各種地質(zhì)參數(shù)的空間展布特征, 為盆地模擬奠定基礎(chǔ)。盆地模擬方面包括構(gòu)造、沉積、儲(chǔ)層、古水動(dòng)力場、古地溫、生烴、排烴、圈閉演化和油氣運(yùn)移聚集等各個(gè)部分。其中, 從生烴到運(yùn)移的模擬構(gòu)成成藏動(dòng)力學(xué)過程模擬的主體, 而其他的描述和模擬則是成藏動(dòng)力學(xué)過程模擬必不可少的重要基礎(chǔ)。成藏動(dòng)力學(xué)過程模擬的最終結(jié)果體現(xiàn)在油氣資源量計(jì)算部分上, 包括計(jì)算出盆地的生烴量、排烴量、烴碳轉(zhuǎn)換量、油氣損失量, 最后要計(jì)算出盆地中聚集的油氣資源量[2]。，油氣成藏機(jī)理。

2.油氣成藏動(dòng)力學(xué)系統(tǒng)的劃分及類型

田世澄(1996) 提出將受地球深部動(dòng)力學(xué)控制的盆地構(gòu)造2沉積旋回作為一個(gè)成藏動(dòng)力學(xué)系統(tǒng), 把改變地下成藏動(dòng)力學(xué)條件, 影響成藏動(dòng)力學(xué)過程的區(qū)域不整合和區(qū)域分布的異?？紫读黧w壓力界面作為不同成藏動(dòng)力學(xué)系統(tǒng)的界面。并據(jù)動(dòng)力學(xué)特征將成藏動(dòng)力學(xué)系統(tǒng)分為開放型、封閉型、半封閉型3 種類型, 據(jù)油源特征又區(qū)分為自源成藏動(dòng)力學(xué)系統(tǒng)和他源成藏動(dòng)力學(xué)系統(tǒng)。因此共可劃分出6 種油氣成藏動(dòng)力學(xué)系統(tǒng)[3-6]。康永尚(1999) 根據(jù)系統(tǒng)動(dòng)力的來源、去向和系統(tǒng)的演化方式將油氣成藏流體動(dòng)力系統(tǒng)分為重力驅(qū)動(dòng)型、壓實(shí)驅(qū)動(dòng)型、封存型和滯留4 種。，油氣成藏機(jī)理。實(shí)際上重力驅(qū)動(dòng)型對(duì)應(yīng)開放型, 壓實(shí)驅(qū)動(dòng)型對(duì)應(yīng)半開放型, 封存型和滯留型則對(duì)應(yīng)封閉型。，油氣成藏機(jī)理。，油氣成藏機(jī)理。因此二者是一致的。這種以油氣成藏的動(dòng)力因素來劃分油氣系統(tǒng)的方法比經(jīng)典的含油氣系統(tǒng)的一套源巖對(duì)應(yīng)一個(gè)油氣系統(tǒng)的粗略劃分方法更深入, 更能體現(xiàn)油氣作為一種流體的運(yùn)動(dòng)分布規(guī)律, 從而有效指導(dǎo)我國陸相含油氣盆地的勘探[7]。

3.油氣成藏主要?jiǎng)恿σ蛩氐难芯?/p>

沉積盆地實(shí)際上是一個(gè)低溫?zé)峄瘜W(xué)反應(yīng)器, 油氣的富集是由溫度、力和有效受熱時(shí)間控制的化學(xué)動(dòng)力學(xué)過程, 及由壓力、地應(yīng)力、浮力和流體勢控制的流體動(dòng)力學(xué)過程的綜合結(jié)果, 也是盆地中各個(gè)成藏動(dòng)力學(xué)系統(tǒng)中的油、氣、水三相滲流過程的結(jié)果。張厚福(1998) 認(rèn)為: 地溫場、地壓場、地應(yīng)力場等“三場”系受地球內(nèi)能控制, 是地球內(nèi)部能量在地殼上的不同表現(xiàn)表現(xiàn)形式。“三場”相互之間彼此影響與聯(lián)系。“三場”的作用使地殼上形成海盆、湖盆等各種水域, 才衍生出水動(dòng)力場, 有了水體才能出現(xiàn)化學(xué)場與生物場, 后二者也相互聯(lián)系與相互制約。綜合這些場的作用, 在含油氣盆地內(nèi)才出現(xiàn)油氣成藏動(dòng)力系統(tǒng)與流體壓力封存箱等地質(zhì)實(shí)體, 后二者之間互有聯(lián)系和影響。油氣從烴源巖生成并排出到相鄰的輸導(dǎo)層經(jīng)運(yùn)移聚集而成藏及成藏后發(fā)生的物理化學(xué)變化這一系列過程都始終貫穿“三場”的作用[8-10]。

4.含油氣系統(tǒng)和油氣成藏動(dòng)力學(xué)的關(guān)系探討

目前對(duì)含油氣系統(tǒng)和油氣成藏動(dòng)力系統(tǒng)之間的關(guān)系眾說紛紜。主要有3 種說法。(1) 含油氣系統(tǒng)研究是油氣成藏動(dòng)力學(xué)研究的起點(diǎn)。(2) 油氣成藏動(dòng)力學(xué)研究是含油氣系統(tǒng)研究的基礎(chǔ)。王英民(1998) 認(rèn)為含油氣系統(tǒng)劃分是成藏動(dòng)力學(xué)研究的結(jié)果。，油氣成藏機(jī)理。(3) 含油氣系統(tǒng)和油氣成藏動(dòng)力學(xué)系統(tǒng)是交叉關(guān)系。筆者認(rèn)為由油氣運(yùn)聚的物質(zhì)空間和動(dòng)力因素控制的流體輸導(dǎo)系統(tǒng)的研究是油氣成藏動(dòng)力學(xué)研究的核心內(nèi)容, 油氣成藏動(dòng)力學(xué)研究應(yīng)按照從源巖到圈閉這一歷史主線, 側(cè)重于油氣成藏的動(dòng)力學(xué)與運(yùn)動(dòng)學(xué)機(jī)制的研究。但油氣成藏動(dòng)力系統(tǒng)對(duì)應(yīng)的狀態(tài)空間是油氣藏。而含油氣系統(tǒng)是從油氣顯示開始, 而不考慮其是否具有工業(yè)價(jià)值。因此油氣成藏動(dòng)力系統(tǒng)是在大的合油氣系統(tǒng)研究基礎(chǔ)上進(jìn)一步按油氣運(yùn)聚動(dòng)力學(xué)條件而追蹤油氣分布規(guī)律。因此筆者傾向于第一種說法, 認(rèn)為在含油氣系統(tǒng)宏觀研究思路基礎(chǔ)上進(jìn)行油氣成藏動(dòng)力學(xué)過程的系統(tǒng)研究, 并根據(jù)成藏動(dòng)力源泉進(jìn)一步劃分油氣成藏動(dòng)力系統(tǒng), 才能弄清我國陸相盆地的成藏機(jī)理和油氣分布規(guī)律并建立當(dāng)代高等石油地質(zhì)理論, 從而更好地指導(dǎo)21 世紀(jì)的油氣勘探[11]。

參考文獻(xiàn)

[]佩羅東1石油地質(zhì)動(dòng)力學(xué)[M]1北京:石油工業(yè)出版社,19931

[2]孫永傳1石油地質(zhì)動(dòng)力學(xué)的理論與實(shí)踐[J]1地學(xué)前緣,1995,2(224):92141

[3]康永尚,等1油氣成藏流體動(dòng)力系統(tǒng)分析原理及應(yīng)用[J]1沉積學(xué)報(bào),1998,16(3):802841

[4]康永尚,等1油氣成藏流體動(dòng)力學(xué)[M]1北京:地質(zhì)出版社,19991

[5]王英民1殘余盆地成藏動(dòng)力學(xué)過程研究方法[J]1成都理工學(xué)院學(xué)報(bào),1998,25(3):38523921

[6]胡朝元,等1成油系統(tǒng)概念在中國的提出及應(yīng)用[J]1石油學(xué)報(bào),1996,17(1)1

[7]龔再升,等1南海北部大陸邊緣盆地分析與油氣聚集[M]1北京:科學(xué)出版社,19971

[8]張厚福,等1石油地質(zhì)學(xué)[M]1北京:石油工業(yè)出版社,19931

[9]田世澄,等1論成藏動(dòng)力學(xué)系統(tǒng)[J]1勘探家,1996,1(2):202241

[10]張厚福，石油地質(zhì)學(xué)新進(jìn)展[M]1北京:石油工業(yè)出版社,19981

[11]費(fèi)琪,等1成油體系與成藏動(dòng)力學(xué)論文集[C]1北京:地震出版社,19991

篇3

完成單位：東華大學(xué)

“纖維/高速氣流兩相流體動(dòng)力學(xué)及“應(yīng)用基礎(chǔ)研究”項(xiàng)目解決了纖維在高速氣流場中的耦合相互作用問題，這在國內(nèi)外同類研究中未見報(bào)道。作為基礎(chǔ)科學(xué)研究項(xiàng)目，東華大學(xué)紡織學(xué)院獲得紡織之光2012年度科技獎(jiǎng)一等獎(jiǎng)的項(xiàng)目，為紡織新技術(shù)、新工藝的產(chǎn)生與發(fā)展奠定了基礎(chǔ)。

紡織氣流問題是紡織科學(xué)中的重要基礎(chǔ)性課題。自上世紀(jì)80年代高速氣流在噴氣織機(jī)、噴氣紡紗機(jī)得到商業(yè)化應(yīng)用以來，在陸續(xù)產(chǎn)生了變形絲、網(wǎng)絡(luò)絲、空氣捻接、熔噴非織造和噴氣渦流紡紗等現(xiàn)代紡織新技術(shù)中，高速氣流技術(shù)已成為現(xiàn)代紡織加工的主流技術(shù)之一。

高速氣流與其加工對(duì)象——纖維或紗線之間的耦合作用特性是這些技術(shù)共有的基本特征。但是，國內(nèi)外的相關(guān)研究還主要局限在低速氣流，纖維/高速氣流兩相流體動(dòng)力學(xué)的基礎(chǔ)研究還相對(duì)薄弱，難以為高速氣流技術(shù)的紡織應(yīng)用提供有力的支撐。

項(xiàng)目研究團(tuán)隊(duì)于1995年起，在國家自然科學(xué)基金委、教育部和上海市等科研計(jì)劃項(xiàng)目的支持下，圍繞纖維/高速氣流兩相流動(dòng)力學(xué)開展了系統(tǒng)的研究工作，實(shí)現(xiàn)了五個(gè)方面的技術(shù)創(chuàng)新。

作為項(xiàng)目組成員之一，東華大學(xué)教授曾泳春為我們做了更加詳細(xì)的介紹。首先，構(gòu)建了基于柔彈性特征的纖維模型?！爸椤獥U”鏈?zhǔn)胶突谟邢迒卧ǖ睦w維模型，將纖維的柔彈性物理特征納入其中，實(shí)現(xiàn)了纖維在高速氣流中位置、取向及變形的合理描述。二是揭示了纖維在高速氣流場中的耦合作用特性與運(yùn)動(dòng)變形規(guī)律。在國際上首次實(shí)現(xiàn)了纖維在噴氣紡、噴氣渦流紡、氣流減羽等噴嘴中運(yùn)動(dòng)的數(shù)值模擬，獲得了纖維運(yùn)動(dòng)、變形特征及其與高速氣流場的相互作用規(guī)律。三是實(shí)現(xiàn)了紡紗噴嘴內(nèi)高速氣流場流動(dòng)特性的數(shù)值模擬與實(shí)驗(yàn)測試。四是對(duì)纖維，高速氣流兩相流體動(dòng)力學(xué)理論研究成果在高速氣流紡紗中進(jìn)行應(yīng)用。揭示了噴氣紡與噴氣渦流紡加捻、氣流噴嘴減少紗線毛羽的機(jī)理，設(shè)計(jì)了具有自主知識(shí)產(chǎn)權(quán)的紡紗噴嘴，實(shí)現(xiàn)了工藝的系統(tǒng)優(yōu)化與成紗質(zhì)量的精確預(yù)測。五是將纖維，高速氣流兩相流體動(dòng)力學(xué)模型拓展應(yīng)用于超細(xì)纖維紡絲拉伸技術(shù)中。

基礎(chǔ)科學(xué)研究不僅僅是實(shí)驗(yàn)室中的數(shù)據(jù)和論文上的文字，它的投產(chǎn)使用才是研發(fā)的初衷。紡織科學(xué)國家重點(diǎn)學(xué)科現(xiàn)代紡織技術(shù)研究方向?qū)W科帶頭人之一，本項(xiàng)目的第一完成人，東華大學(xué)教授郁崇文讓我們看到的是這項(xiàng)研究的前景。

篇4

關(guān)鍵詞：液壓支架；流體；動(dòng)力響應(yīng)；非穩(wěn)態(tài)

1液壓支架前連桿三維結(jié)構(gòu)模型

利用大型有限元分析軟件Solidworks對(duì)液壓支架前連桿進(jìn)行建模，并且進(jìn)行簡化處理，以便進(jìn)行流體動(dòng)力學(xué)分析。簡化的具體情況：（1）忽略倒角圓角等細(xì)小特征，該特征將直接影響收斂速度，從而影響結(jié)果的精度；（2）忽略建立銷軸零件模型，在載荷設(shè)置中可以設(shè)置相關(guān)的約束，否則會(huì)影響前連桿動(dòng)力學(xué)參數(shù)曲線的輸出精度；（3）忽略焊縫的影響。材料屬性為：彈性模量210GPa，泊松比μ=0.3，密度7850kg/m3，前連桿三維結(jié)構(gòu)模型如圖1所示。圖1前連桿三維結(jié)構(gòu)模型在網(wǎng)格劃分過程中，單元選擇SHELL63，該單元是一種應(yīng)用于大變形和應(yīng)力剛化的四節(jié)點(diǎn)彈性殼單元，每個(gè)節(jié)點(diǎn)自由度具有繞X、Y、Z軸的轉(zhuǎn)動(dòng)和沿X、Y、Z方向的平動(dòng)6個(gè)自由度。通過賦予4個(gè)節(jié)點(diǎn)來建立不同厚度分布的模型，對(duì)于不規(guī)則形狀的結(jié)構(gòu)具有較高精度。網(wǎng)格化的前連桿如圖2所示，節(jié)點(diǎn)數(shù)14537，單元總數(shù)7089。

2基本控制方程

液壓支架前連桿附近流體為黏性不可壓縮流體流動(dòng)，運(yùn)動(dòng)類型為湍流運(yùn)動(dòng)，雷諾數(shù)為4200，并且滿足以下控制方程：連續(xù)性方程。

3前連桿流體動(dòng)力學(xué)分析

邊界條件及數(shù)值計(jì)算方法為前連桿承受的載荷主要為液壓支架的擺動(dòng)及轉(zhuǎn)動(dòng)。在工程應(yīng)用中，一般轉(zhuǎn)速為0.8m/s，本文模型的計(jì)算域大小為長5m、寬3m、高4m，前連桿位于計(jì)算域的中心位置。數(shù)值計(jì)算方法為Simple算法，周圍流體滿足的連續(xù)性方程、動(dòng)量方程、能量方程等的殘差均小于10-6，當(dāng)參數(shù)趨于漸近值時(shí)，達(dá)到收斂，之后采用Quick格式求解，再次進(jìn)行迭代直至收斂。圖3、圖4為前連桿在非穩(wěn)態(tài)情況下切面云圖變化。由圖3可知，前連桿切面周圍溫度由同一方向293.17~293.25K過渡，最終趨于293.19K，呈穩(wěn)態(tài)變化，這是由于該系統(tǒng)符合能量守恒定律，表面有一定的溫度，流體接觸前連桿形成散射波，并且繞流在前連桿周圍，導(dǎo)致溫度升高。流體距離前連桿越遠(yuǎn)，繞流逐漸散去，溫度降低。圖4給出了塔架切面靜壓等高線，經(jīng)過流體的前連桿靜壓變化由1.01278×105Pa向1.01348×105Pa過渡，呈非線性增長，最終趨于1.01320×105Pa，壓力變化平緩，避免了由于壓力變化幅度過大引起結(jié)構(gòu)失效。圖3、圖4說明了前連桿結(jié)構(gòu)能夠適應(yīng)在非穩(wěn)態(tài)流體載荷的溫度及壓力變化，數(shù)值模擬結(jié)果中的溫度及壓力變化沒有出現(xiàn)極值幅度，周圍不易形成漩渦，保證了前連桿在液壓支架中的正常運(yùn)行。圖5、圖6、圖7分別為湍流動(dòng)能、平均動(dòng)能、相對(duì)壓力變化曲線。從圖5可知，前連桿的湍流動(dòng)能最大值為4.8J/kg，隨著長度的逐漸增大先減小后增大，最后達(dá)到最大值，這是由于前連桿機(jī)械能增大，隨后逐漸減小，趨于一個(gè)定值。從圖6可知，平均運(yùn)動(dòng)能在0~12s出現(xiàn)振幅較小的振蕩，12~24s出現(xiàn)幅度較大的振蕩，最后趨于定值，完成了一個(gè)脈動(dòng)周期，說明能量耗散變化趨勢和機(jī)械能與內(nèi)能之間轉(zhuǎn)換的關(guān)系，12~24s能量耗散是最大的，內(nèi)能轉(zhuǎn)換為機(jī)械能也是最大的，符合能量守恒定律，輻射阻尼使系統(tǒng)的部分能量向四周輻射出去，使前連桿平均運(yùn)動(dòng)能呈現(xiàn)振蕩。從圖7可知，隨著長度的增加呈現(xiàn)減小與增加的交替現(xiàn)象，這是由于前連桿在剛啟動(dòng)時(shí)，相對(duì)壓力與平均力矩即扭矩呈正比關(guān)系，隨著前連桿的轉(zhuǎn)動(dòng)，扭矩的勁逐漸增大，然后減小，呈現(xiàn)交替現(xiàn)象，符合扭矩變化規(guī)律。圖8、圖9分別為前連桿應(yīng)力和位移變化云圖，應(yīng)力變化在范圍允許之內(nèi)，整個(gè)前連桿應(yīng)力幅度變化不大，承受主要載荷的區(qū)域顯示淺綠色，其余部分由圖中應(yīng)力的顏色可知在安全范圍內(nèi)，符合強(qiáng)度要求。從位移變化云圖可知，位移在0.02~0.15mm浮動(dòng)，符合穩(wěn)定性要求，穩(wěn)定性尤為重要，直接關(guān)系前連桿與其他配合零件相對(duì)位置，進(jìn)而影響整個(gè)系統(tǒng)的穩(wěn)定。

4結(jié)語

經(jīng)過液壓支架前連桿的流體動(dòng)力學(xué)分析，可以得出，液壓支架前連桿承受非穩(wěn)態(tài)流體載荷的靜壓與溫度變化，符合流體動(dòng)力學(xué)規(guī)律，壓力變化平緩，呈現(xiàn)非線性增長。平均動(dòng)能出現(xiàn)振幅較小的振蕩，最后趨于定值，符合能量守恒定律，應(yīng)力與位移變化幅度較小，位移變化量有10%以上的裕度，穩(wěn)定性較好，符合強(qiáng)度要求，從而為前連桿在運(yùn)行過程中動(dòng)力學(xué)參數(shù)的監(jiān)測提供一定的參考。

參考文獻(xiàn)：

［1］陸丹丹,陸金桂.基于灰熵關(guān)聯(lián)的液壓支架前連桿疲勞壽命影響因素分析［J］.機(jī)械制造與自動(dòng)化,2016,45(5):37-39.

［2］廖華林,李根生,李敬彬,等.徑向水平鉆孔直旋混合射流噴嘴流場特性分析［J］.煤炭學(xué)報(bào),2012,37(11):1895-1900.

［3］王永龍,王振鋒,孫玉寧,等.水力割縫閥體實(shí)現(xiàn)切換的流體動(dòng)力學(xué)分析［J］.煤礦機(jī)械,2015,36(3):123-125.

篇5

關(guān)鍵詞：工業(yè)廠房；煤氣泄漏；火災(zāi)模擬；煙氣流動(dòng)

中圖分類號(hào)：X92 文獻(xiàn)標(biāo)識(shí)碼:A

1 概述

隨著工業(yè)化發(fā)展的不斷提高，工業(yè)用火、用電、用氣和化學(xué)物品的應(yīng)用也日益廣泛，工業(yè)建筑物火災(zāi)的危險(xiǎn)性和危害性大大增加，伴隨而來的火災(zāi)發(fā)生的頻率也越來越高，尤其在工業(yè)廠房內(nèi)，一旦有火災(zāi)發(fā)生就會(huì)對(duì)人員的生命財(cái)產(chǎn)造成直接的危害，嚴(yán)重的情況下會(huì)導(dǎo)致巨大的人員傷亡和財(cái)產(chǎn)損失。

例如以簡化的廠房空間為研究對(duì)象，借助大型專業(yè)化計(jì)算流體力學(xué)軟件FLUENT, 對(duì)單室內(nèi)由于煤氣泄漏引起的火災(zāi)及煙氣蔓延過程進(jìn)行數(shù)值模擬，得到了火場中溫度場和燃燒產(chǎn)物二氧化碳濃度場的直觀顯示，為人們了解火災(zāi)發(fā)生和發(fā)展的過程提供了新的方法和手段，也為工業(yè)建筑防火設(shè)計(jì)和消防安全評(píng)估提供了新的科學(xué)工具，是消防安全工程學(xué)和性能化設(shè)計(jì)的重要基礎(chǔ)。

2本文設(shè)定的場景與幾何建模

實(shí)驗(yàn)房間幾何尺寸為長6.0m，寬3.3m，高4.5m，房間墻壁有一高1.2m，寬2.0m的窗戶，窗戶對(duì)面是一高2.0m，寬0.8m的房門，在一側(cè)壁中央位置距地面高1.0m有一處煤氣泄漏點(diǎn)，本文考慮的泄漏危險(xiǎn)源是沿墻鋪設(shè)的煤氣匯流排管道有一處泄漏，形成直徑為1cm的圓形泄漏噴口，泄漏流量為0.1413m3/h（即泄漏速度為0.5m/s），泄漏煤氣的主要成分為甲烷。如圖2.1所示為在GAMBIT幾何建模軟件中建立的實(shí)驗(yàn)房間幾何模型。

2.1 實(shí)驗(yàn)房間幾何模型

3建立數(shù)學(xué)物理模型

火災(zāi)的場模擬理論基礎(chǔ)是基于質(zhì)量守恒、動(dòng)量守恒、能量守恒、化學(xué)組分平衡的微分方程，而場模擬的過程就是對(duì)這些方程組進(jìn)行數(shù)值求解。然而這些方程組是不封閉的，為了使方程組封閉，就必須引入湍流的附加方程，然后借助計(jì)算機(jī)對(duì)封閉方程組進(jìn)行求解，就可以得出火災(zāi)過程中各種參數(shù)的詳細(xì)空間分布及其隨時(shí)間的變化。

3.1湍流平均量方程

火災(zāi)中燃燒的湍流過程影響著整個(gè)流體的流動(dòng)，要想求解方程組就必須正確處理湍流過程。計(jì)算湍流粘性系數(shù)Ut的方法就是所謂的湍流模型，湍流模型種類很多，本文引用應(yīng)用范圍較廣的k-epslion雙方程模型：

3.2邊界條件

FLUENT軟件包含了更為廣泛的邊界條件設(shè)置類型，本文僅結(jié)合所研究的實(shí)際問題進(jìn)行討論。

（1）流體進(jìn)口邊界條件：本文對(duì)于著火廠房設(shè)置了兩個(gè)流體進(jìn)口，一個(gè)是燃料進(jìn)口，另一個(gè)是空氣進(jìn)口。其中甲烷的泄漏速度為0.5m/s，溫度為300K；窗戶進(jìn)風(fēng)速度始終為0.5m/s，溫度300K。

（2）流體出口界面：本文場景中廠房的房門為出流口。

（3）壁面：本文的研究中采用絕熱、無滲透、無滑移壁面。

4模擬結(jié)果及分析

數(shù)值模擬計(jì)算采用FLUENT軟件進(jìn)行，使用分離解方法。求解過程中，燃燒控制方程具有較好的收斂記錄。

圖4.1 60秒時(shí)溫度等值線圖

圖4.2 60秒時(shí)CO2濃度等值線圖

結(jié)果分析

高溫?zé)煔鈱?duì)人體的危害如果煙氣層面高于人眼的特征高度時(shí)（人眼的特征高度通常為1.2～1.8m，一般取1.5m），其煙氣層的熱輻射強(qiáng)度就會(huì)傷害人，而當(dāng)煙氣層面低于人眼的特征高度時(shí)，對(duì)人的危害將是直接燒傷或由于吸入熱氣體而對(duì)呼吸系統(tǒng)的破壞。由60s時(shí)圖4.1可以看出，房間窗戶下沿1.8m處煙氣溫度已經(jīng)達(dá)到180℃左右，此時(shí)可對(duì)人體造成輻射傷害，在房間中后部煙氣高度已接近1.5m，將對(duì)人體構(gòu)成直接傷害。所以，人置身于煤氣泄漏的火場時(shí)，要在短時(shí)間內(nèi)逃離，以防熱煙氣層的沉積造成高溫?zé)煔獾妮椛?、灼傷而阻礙逃生。

CO2濃度在60秒時(shí)有明顯的增長趨勢。主要是廠房空間內(nèi)原有的氧氣和窗戶的進(jìn)風(fēng)使煤氣燃燒充分，生成較多的CO2，幾乎在短時(shí)間內(nèi)就彌漫了廠房的上半層空間，圖4.2可以看出，廠房的大部分空間已經(jīng)彌漫了濃度較高的CO2，并且在生成CO2的同時(shí)消耗了大量的氧氣，CO2對(duì)人體的危害主要是窒息作用。

結(jié)論

采用專業(yè)化計(jì)算流體力學(xué)軟件FLUENT對(duì)工業(yè)廠房的煤氣泄漏火災(zāi)進(jìn)行了模擬，把火災(zāi)數(shù)值模擬技術(shù)應(yīng)用于實(shí)際，實(shí)現(xiàn)了火場中的物理量，如溫度、生成物組分二氧化碳濃度隨泄漏時(shí)間變化的直觀顯示，為廠房內(nèi)危險(xiǎn)區(qū)域的劃分及災(zāi)害控制提供了參考依據(jù)，也為廠房結(jié)構(gòu)空間的性能化設(shè)計(jì)和煤氣泄漏火災(zāi)的防治提供了理論依據(jù)。

參考文獻(xiàn)

[1]王福軍.計(jì)算流體動(dòng)力學(xué)分析[M].清華大學(xué)出版社，2004.

[2]韓占先，徐寶林，霍然等.火災(zāi)科學(xué)與消防工程[M].山東科學(xué)技術(shù)出版社，2001.

[3]霍然等.建筑火災(zāi)中煙氣蔓延的動(dòng)態(tài)顯示[J].火災(zāi)科學(xué)，1995，Vol.4，No.3.

[4]范維澄，萬躍鵬.流體及燃燒的模型與計(jì)算[M].北京：中國科學(xué)技術(shù)大學(xué)出版社，1992.

篇6

所謂的流體力學(xué)，是一門專門研究像液體、氣體等流體運(yùn)動(dòng)形態(tài)以及其運(yùn)動(dòng)規(guī)律和應(yīng)用的學(xué)科。流體在力的作用下，其本身的物理性質(zhì)會(huì)發(fā)生一定的變化，比如流體的固體壁面，流體和固體之間的相互作用，流體的壓強(qiáng)、溫度，還有流體的密度等等變量，都會(huì)發(fā)生變化。在選礦的過程中，有許多地方都會(huì)用到流體力學(xué)的知識(shí)，比如一些設(shè)備的設(shè)計(jì)，就是按照流體力學(xué)的原理來設(shè)計(jì)的，所以流體力學(xué)在選礦中起著極其重要的作用。本文就流體力學(xué)的特點(diǎn)做了詳細(xì)的闡述，以及對(duì)流體力學(xué)在選礦中所存在的問題進(jìn)行了深入的探討。

前言

流體力學(xué)的發(fā)展起源于阿基米德，他是一位古希臘人，建立了液體平衡理論，在這個(gè)基礎(chǔ)上，流體力學(xué)的運(yùn)用才得以發(fā)展。而如今流體力學(xué)的應(yīng)用范圍越來越廣泛，在空氣動(dòng)力學(xué)、電磁計(jì)算流體學(xué)、氣液相流體學(xué)等方面都有十分重大的突破。

冶金和煤礦的生產(chǎn)都需要進(jìn)行選礦，在工業(yè)生產(chǎn)上，選礦起著非常重要的作用。所謂選礦，就是通過化學(xué)或者是物理等方法把有用的礦物和無用的礦物分開，最終獲取能夠被生產(chǎn)中得以利用的礦物的過程。選擇有用物料的方式分為兩種情況。

一種是針對(duì)單體解離狀態(tài)的，可以采用不同的選礦方法，把我們所需要的有用的礦物篩選出來；

另一種是針對(duì)不能達(dá)到單體解離狀態(tài)的，這就需要將物料用機(jī)械破碎和磨碎的方法來取得當(dāng)中的有用礦物。回看整個(gè)世界的選礦發(fā)展史，不難看出，流體力學(xué)在整個(gè)選礦發(fā)展過程中起到了不可替代的作用。眾所周知，在工業(yè)生產(chǎn)上選礦的方法有很多，比如利用重力來進(jìn)行篩除選礦，還有利用浮力來選礦等等。隨著科技的發(fā)展，應(yīng)用流體力學(xué)的原理來進(jìn)行選礦，已經(jīng)越來越普及，現(xiàn)在許多大型選礦機(jī)器都有流體力學(xué)的應(yīng)用。也就是說，要進(jìn)行高質(zhì)量、高效率的選礦，流體力學(xué)的運(yùn)用是必不可少的環(huán)節(jié)。

1.選礦中所用的設(shè)備

1.1旋流器

所謂的旋流器，在液體的分級(jí)分離中經(jīng)常用到它。它的工作原理是利用離心的方法，使得不同質(zhì)量不同密度的物質(zhì)分離沉降。旋流器整體是一個(gè)圓錐形，采出來的礦漿通過管道進(jìn)入旋流器內(nèi)部，進(jìn)入的礦漿在旋流器強(qiáng)大的離心力作用下，產(chǎn)生旋轉(zhuǎn)。礦漿中含有不同質(zhì)量和大小的礦物質(zhì)顆粒，在旋流器離心力的作用下，大的礦物質(zhì)顆粒和密度高的顆粒就被拋向旋流器內(nèi)壁，在重力的作用下，這些顆?？恐鴥?nèi)壁逐漸的下沉，最后在旋流器的底部流出。而那些比較細(xì)小的顆粒和液體則經(jīng)過溢流管口流出來。所以，礦漿經(jīng)過旋流器，就把大的礦物顆粒和細(xì)小的礦物顆粒分離開來，離心力的篩選是很有效的，相比于重力場，它具有快速、高效、穩(wěn)定分離的特點(diǎn)。目前的選礦旋流器主要有兩種，

一是用水作為介質(zhì)的水力旋流器，第二種就是專門配制的重介質(zhì)液體，密度比水更高，以此作為介質(zhì)的重介質(zhì)旋流器。這兩種旋流器的作用各不相同，

第一種是專門用于脫泥以及脫浮選藥劑的水力旋流器，第二種是專門用于最后分選冶煉。旋流器的原理簡單，但是液體在其內(nèi)部的流動(dòng)卻是十分復(fù)雜的，這要涉及到流體力學(xué)中的湍流運(yùn)動(dòng)。所謂的湍流運(yùn)動(dòng)就是將流體的流動(dòng)看作是一個(gè)三維空間內(nèi)的運(yùn)動(dòng)，在這個(gè)三維空間內(nèi)，任何一點(diǎn)的速度都可以分解成橫向速度、縱向速度和切向速度。就是這樣的運(yùn)動(dòng)特征和旋流器內(nèi)復(fù)雜的邊界，使得內(nèi)部的流體力學(xué)并沒有一定的規(guī)律可循，自然就沒有具體的方法來研究，只能通過實(shí)際的測量來估算和分析。通過實(shí)際測量，旋流器內(nèi)的運(yùn)動(dòng)實(shí)際是一個(gè)組合運(yùn)動(dòng)，它是通過自由渦和強(qiáng)制渦組合在一起共同形成的。在旋流器工作的時(shí)候，由于切向速度最大，導(dǎo)致了旋流器內(nèi)壁的壓力小于其中問的壓力而形成一個(gè)負(fù)壓區(qū)，空氣從底部進(jìn)入，旋流器中間區(qū)域就成了真空的空氣柱。那些比較重的礦物質(zhì)直接由旋流器內(nèi)壁被分離了下來，從底部排出，那些輕的礦物質(zhì)就由上面的口排出。

1.2浮選機(jī)

浮選礦物質(zhì)實(shí)際就是利用浮力，將礦物質(zhì)浮到水面來進(jìn)行篩選。這就要求礦漿中必須要充滿大量的氣泡，必要時(shí)，還要在礦漿中添加一種捕收劑，這種藥劑會(huì)作用在礦物質(zhì)顆粒表面，使得礦粒的疏水性增加，利于捕收。有些礦粒是疏水性，而有些是親水性的，疏水性的礦粒和產(chǎn)生的氣泡相結(jié)合，由于浮力原因而浮出水面，那些清水性的礦粒則不會(huì)與氣泡結(jié)合而下沉。產(chǎn)生氣泡的方法有很多，一般采用機(jī)械攪拌和電解水來產(chǎn)生氣泡。

1.3浮選柱

浮選柱也是一種專門用來選擇礦物質(zhì)的設(shè)備。浮選柱早在上世紀(jì)六十年代就開始應(yīng)用，只不過由于氣泡發(fā)生器設(shè)計(jì)的并不好，始終無法很好的產(chǎn)生氣泡，就使得浮選柱并沒有廣泛的運(yùn)用。浮選柱的氣泡產(chǎn)生器位于柱體下端，礦物由上部倒入。其原理和上面介紹的浮選機(jī)差不多，也是通過重力和浮力的相互作用來實(shí)現(xiàn)礦粒的分離。所以氣泡發(fā)生器才是浮選柱工作的關(guān)鍵所在，如果氣泡發(fā)生器無法正常工作產(chǎn)生氣泡，浮選柱也就無法分離礦物質(zhì)。

2.流體力學(xué)在選礦設(shè)備中的應(yīng)用

浮選機(jī)的工作原理是利用氣泡產(chǎn)生的浮力和礦物質(zhì)本身的重力實(shí)現(xiàn)礦粒分離，分離的過程中就需要大量的氣泡，氣泡發(fā)生器正是大量產(chǎn)生氣泡的設(shè)備。通常我們都把氣泡發(fā)生器分為兩類，即外部發(fā)生器和內(nèi)部發(fā)生器這兩類。

外部氣泡發(fā)生器也有兩種，通常用的是氣水噴射型。它是把壓縮空氣和水一起混合，通過噴頭直接噴入浮選機(jī)內(nèi)，再通過分流孔將水流分開，從而產(chǎn)生氣泡。

顧名思義，內(nèi)部氣泡發(fā)生器是將有孔介質(zhì)固定在浮選機(jī)內(nèi)部，然后沖入空氣，通過介質(zhì)的孔向外噴出氣泡。內(nèi)部氣泡發(fā)生器常常出現(xiàn)問題，原因在于有孔介質(zhì)放在浮選機(jī)內(nèi)部，很容易被礦粒堵塞和結(jié)垢，使浮選機(jī)無法正常工作，現(xiàn)在已經(jīng)很少采用這種方法了。

上述兩種氣泡發(fā)生器都是以往一直被使用的。但畢竟它的結(jié)構(gòu)太過復(fù)雜容易出現(xiàn)問題，而且能量消耗也不小。目前被廣泛采用的是另一種更加先進(jìn)的設(shè)備――微泡發(fā)生器。這種氣泡發(fā)生器利用了流體力學(xué)的原理，首先流體在其內(nèi)部分散開來，在后半段又重新成為連續(xù)相。微泡發(fā)生器動(dòng)力消耗小，產(chǎn)生氣泡效率高，礦化充分，而且使用壽命長，大大節(jié)省了維護(hù)費(fèi)用。

目前較為先進(jìn)的浮選機(jī)是德國制造的切向噴射浮選機(jī)，其也是利用流體力學(xué)的原理進(jìn)行選料。它有兩個(gè)槽室，一個(gè)收集精礦，一個(gè)用作分離室。礦漿沿著切線方向射入容器，從而產(chǎn)生旋轉(zhuǎn)，和氣泡結(jié)合的礦物質(zhì)浮力大，向上運(yùn)動(dòng)。沒有結(jié)合水的礦粒則向下沉降，最后通過底部排出。

篇7

關(guān)鍵詞： 巖石變形巖土模型試驗(yàn)透明模型

中圖分類號(hào)：F470.1 文獻(xiàn)標(biāo)識(shí)碼：A

0 引言

地下巖體在施工過程中存在很多不確定性因素，由于這些風(fēng)險(xiǎn)無法用傳統(tǒng)的經(jīng)驗(yàn)公式法、理論計(jì)算法和數(shù)值模擬分析法進(jìn)行求解，因此研究人員只有通過一些模型試驗(yàn)來進(jìn)行類比推理。可以說，運(yùn)用模型試驗(yàn)法來探索巖體內(nèi)部的變形相對(duì)于其他方法來說比較直接有效。

1 傳統(tǒng)的物理模型試驗(yàn)法

物理模型試驗(yàn)以實(shí)驗(yàn)力學(xué)為基礎(chǔ)，作為一種研究手段，發(fā)揮了非常大的作用[1]，這種研究方法在模擬巖石的應(yīng)力變化和分布情況、巖石結(jié)構(gòu)的破壞情況以及發(fā)展變化等方面都非常有效[2]。

相似模型試驗(yàn)是在實(shí)驗(yàn)室內(nèi)按相似原理制作與原型相似的模型，借助測試儀表觀測模型內(nèi)力學(xué)參數(shù)及其分布規(guī)律，利用在模型上研究的結(jié)果，借以推斷原型中可能發(fā)生的力學(xué)現(xiàn)象以及巖體壓力分布的規(guī)律，以便為巖土工程設(shè)計(jì)和施工方案的選擇提供依據(jù)。這種研究方法直觀且實(shí)驗(yàn)周期短，而且能夠根據(jù)需要，通過固定某些參數(shù)，改變另一些參數(shù)來研究隧道圍巖應(yīng)力在空間和時(shí)間上的分布規(guī)律和變化情況，這在現(xiàn)場條件下是難以實(shí)現(xiàn)的。而且根據(jù)模型實(shí)驗(yàn)的結(jié)果，可以推斷原型可能出現(xiàn)的力學(xué)現(xiàn)象與變形狀態(tài)以及工程的穩(wěn)定性與安全程度。

在巖土模型試驗(yàn)中，研究人員用來觀察量測巖體內(nèi)部變形的方法主要有在土體中預(yù)埋傳感器或電阻應(yīng)變片法、X-射線、γ-射線、計(jì)算機(jī)軸向斷層掃描（CAT）和核磁共振成像（MRI）法、數(shù)字圖像技術(shù)法及光學(xué)測量技術(shù)法等，但是這些方法都有一定的局限性，不能解決隧道模型內(nèi)部變形過程可觀測的問題。

2 數(shù)字照相變形量測技術(shù)

數(shù)字照相量測就是利用數(shù)碼相機(jī)、數(shù)碼攝像機(jī)、CCD攝像機(jī)、CT掃描、激光照相等圖像采集手段，獲得觀測目標(biāo)的數(shù)字圖像后，再利用計(jì)算機(jī)數(shù)字圖像處理與分析技術(shù)，對(duì)觀測目標(biāo)進(jìn)行變形或特征識(shí)別分析的一種現(xiàn)代量測技術(shù)。

位移測量和模型內(nèi)部變形情況是相似材料模型試驗(yàn)的主要觀測內(nèi)容，也是進(jìn)一步研究模型內(nèi)部力學(xué)形態(tài)的基礎(chǔ)。土的變形測量通常有兩種方式：（1）在模型邊界處設(shè)置位移探頭以測定某點(diǎn)的位移；（2）利用光學(xué)技術(shù)對(duì)整個(gè)研究區(qū)域的變形場進(jìn)行非插入式測量。傳統(tǒng)的測量方法有物理測量和機(jī)械測量，其缺點(diǎn)就是觀測裝置或者傳感器安裝比較麻煩、工作量大，而且采樣點(diǎn)有限，只能測量土體塊邊界上的變形，隨著數(shù)碼產(chǎn)品的高性能化和計(jì)算機(jī)圖像處理技術(shù)的不斷發(fā)展，現(xiàn)在已經(jīng)能夠測量并展現(xiàn)土體內(nèi)部的變形過程，數(shù)字照相量測技術(shù)也成為地下工程相似模型試驗(yàn)中越來越重要的測試手段。李元海[5]開發(fā)出了一套數(shù)字圖像變形量測無標(biāo)點(diǎn)分析軟件PhotoInfor程序和后處理程序PostViewer，該項(xiàng)技術(shù)有助于研究巖土模型全程漸進(jìn)性變形破壞過程。

數(shù)字照相量測技術(shù)可以幫助研究人員進(jìn)一步認(rèn)識(shí)土體的內(nèi)部變形情況。在對(duì)巖體模型進(jìn)行數(shù)字照相變形量測時(shí)，利用數(shù)碼相機(jī)和激光切片試樣可以獲得土體內(nèi)部清晰的變形圖像。無損測量試樣的變形是依靠由計(jì)算機(jī)、激光、CCD相機(jī)和幀存儲(chǔ)器構(gòu)成的光學(xué)系統(tǒng)獲得的?；蛘撸覀円部梢岳眉す獬上?，由土壤和激光共同作用來生成特別的散斑圖，由兩幅散斑圖可以獲得二維變形場，綜合分析多個(gè)二維變形場即可獲得巖土體的三維變形場。

3 透明材料模型試驗(yàn)法

透明巖體模型試驗(yàn)法克服了巖土體內(nèi)部變形不可視的不足，實(shí)現(xiàn)了內(nèi)部變形的可視化觀測，這種方法不需要接觸模型，因此排除了觀測方法對(duì)實(shí)驗(yàn)結(jié)果的干擾，同時(shí)該試驗(yàn)方法較以往的光學(xué)觀測實(shí)驗(yàn)而言更為便捷，測量準(zhǔn)確度也得到了很大提高[6]。

自然土體由兩部分組成：土顆粒和孔隙流體。在用透明材料進(jìn)行模擬時(shí)需要分別選定性質(zhì)相近的材料將兩者進(jìn)行合成，選擇時(shí)為了保證合成材料的透明性，務(wù)必使兩種透明材料的折射率近似相等。用來模擬土顆粒的材料必須滿足以下兩點(diǎn)要求：物理性質(zhì)與自然土相近、具有一定的透明度。

模擬材料主要有無定型硅石粉末、無定型硅石凝膠和熔融石英砂。用來模擬孔隙流體的材料，必須滿足以下幾點(diǎn)要求：物理性質(zhì)與孔隙流體相近、化學(xué)性質(zhì)穩(wěn)定且透明性良好，更為重要的一點(diǎn)是，這種模擬孔隙流體的材料務(wù)必具有與模擬土顆粒的材料相似的折射系數(shù)，而且不能與模擬土顆粒的材料發(fā)生反應(yīng)。目前為止，孔隙流體的模擬材料主要有兩種：（1）石蠟和白色礦物油混合而成，二者的重量比為1:1；（2）溴化鈣和水混合而成，即所謂的鹵水[7]。

吳明喜[9]用折射率為1.458、粒徑為0.1-1mm的熔融石英砂模擬砂土顆粒，用濃度為75.5%的溴化鈣水溶液模擬孔隙流體，根據(jù)石英砂的折射率采用折射率-濃度法配制與石英砂相同折射率的溴化鈣溶液，然后將石英砂溶解在溴化鈣溶液中并進(jìn)行真空排氣得到透明土試樣，對(duì)試樣進(jìn)行三軸壓縮試驗(yàn)。研究結(jié)果表明，該法得到的砂土試樣具有良好的透明度，當(dāng)不同應(yīng)力作用其上時(shí)，該合成砂土與自然砂土的應(yīng)力應(yīng)變特征相近，可在某種程度上模擬砂土的巖土特征。

許國安[10]選用液體石蠟與正十三烷按質(zhì)量比0.85混合而成的溶液模擬孔隙流體，用細(xì)度為300目的硅石凝膠模擬土顆粒，將硅膠粉溶解在混合溶液中并進(jìn)行真空排氣得到透明土試樣。對(duì)試樣進(jìn)行三軸試驗(yàn)，測定了固結(jié)壓力不同時(shí)透明材料的相關(guān)力學(xué)參數(shù)。

用透明材料來模擬巖體，同時(shí)結(jié)合數(shù)字照相量測系統(tǒng)對(duì)巖體內(nèi)部的變形進(jìn)行觀測分析，這樣可以更加方便和形象地得到外界荷載作用下巖體內(nèi)部漸進(jìn)變化情況，從而彌補(bǔ)常規(guī)巖體相似材料的不透明性缺點(diǎn)，可為巖土工程施工的可視化試驗(yàn)研究提供一定的幫助。因此，借助于透明相似材料和數(shù)字照相量測方法以及相似模型試驗(yàn)技術(shù)對(duì)隧道圍巖變形破裂規(guī)律進(jìn)行試驗(yàn)研究具有一定的學(xué)術(shù)價(jià)值，也會(huì)對(duì)實(shí)際工程有一定的指導(dǎo)意義。

4 研究不足

由于透明材料的特殊性，其強(qiáng)度有一定上限，合成的透明材料只能模擬軟巖的性質(zhì)，暫時(shí)無法模擬硬巖。

合成的透明材料顆粒之間的粘結(jié)力不夠，與實(shí)際情況有一定差距。

透明巖體的透明度有限，實(shí)驗(yàn)過程中對(duì)外界環(huán)境要求較高。

參考文獻(xiàn)：

徐挺.相似方法及其引用[M].北京：機(jī)械工業(yè)出版社，1995.

趙文，謝強(qiáng)，詹志鋒.北盤江大橋岸坡位移特征模型試驗(yàn)研究[J].四川大學(xué)學(xué)報(bào)，2002，34（4）：60-63.

吳明喜.人工合成透明砂土及其三軸試驗(yàn)研究[D].大連理工大學(xué)，2006.

許國安.深部巷道圍巖變形損傷機(jī)理及破裂演化規(guī)律研究[D].中國礦業(yè)大學(xué)，2011.

李元海，靖洪文，曾慶有，等.巖土工程數(shù)字照相量測軟件系統(tǒng)研發(fā)與應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào)，2006，25（增刊2）：3859-3866.

隋旺華，高岳，LIU Jin-yuan.透明土實(shí)驗(yàn)技術(shù)現(xiàn)狀與展望[J].煤炭學(xué)報(bào)，2011，36（4）：577-582.

王秀華，宰金珉，蔣剛，陳亞東.用于透視土體內(nèi)部變形的透明材料的研究[J].常州工學(xué)院學(xué)報(bào)，2008，21：270-273.

佘躍心.基于透明介質(zhì)和顆粒圖像技術(shù)的土體變形測量研究進(jìn)展[J].勘察科學(xué)技術(shù)，2005，（6）：7-10.

篇8

關(guān)鍵詞：儲(chǔ)層砂巖；出砂；三維顆粒流程序；力學(xué)響應(yīng)；不同物性

中圖分類號(hào)：TU45

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：16744764（2014）02002807

Abstract：The reservoir sandstone is the sedimentary rock with many cemented sand particles. The different physical property of the reservoir even in the same area results in different mechanical response and sand production of the sandstone. Taking two types of the reservoir sandstones as research subjects， a numerical model based on 3Dimensional Particlebased Distinct Element （PFC3D） under cylindrical coordinate system was used to simulate the micromicro response of the sandstone considering the given confining pressure and the oil flow rate. Meanwhile， the sanding initiation and the process of the development were analyzed. The macro stress indicated that the reservoir sandstones with weaker cemented sand particles and less percent of the cemented materials would yield and fail more easily， the sand production was initiated more easily as well. Meanwhile， the stress of the parallel bonds indicated that the reservoir sandstones with granule and less cemented materials dislodged from the sandstone more easily， and the force on the particle contact was larger and the sandstone failure was more serious.Hence， the probability of the dislodged particles flowing into the wellbore was also much more. In addition， the particle displacement and rotation indicated that the physical property of the sandstone played a significant influence on the mechanical response and the sand production， the results also agreed with the above results. Therefore， the sand mechanical response and sand production of reservoir sandstone are both different due to different physical property of the reservoir. As a result， the reliable measures of the sand prediction and sand control will be adopted based on the reservoir physical property and conditions.

Key words：reservoir sandstone； sand production； 3Dimensional particlebased distinct element； mechanical response； different physical property

儲(chǔ)層開采過程中，油井出砂與產(chǎn)能提高的平衡問題一直受到不少學(xué)者和工程人員的關(guān)注，解決該課題對(duì)石油的可持續(xù)發(fā)展具有重要的實(shí)踐意義。目前，隨著油藏勘探難度的增大，優(yōu)質(zhì)儲(chǔ)量的減少，不同物性特征的砂巖儲(chǔ)層陸續(xù)投入開發(fā)。不同物性的儲(chǔ)層由于沉積時(shí)間、埋深、巖性、砂巖顆粒的粒度分布、砂巖顆粒的形狀和排列以及砂巖顆粒間的黏結(jié)物質(zhì)等存在差異，直接影響著砂巖的受力性能，間接干擾了油藏開采中的油井出砂對(duì)生產(chǎn)產(chǎn)能的作用效應(yīng)。由此可見，出砂機(jī)理的研究[1]能提供有效的出砂預(yù)測和防砂措施，具有重要的實(shí)踐價(jià)值。但不同物性的儲(chǔ)層即使在相同的賦存環(huán)境中，也會(huì)出現(xiàn)不同的出砂率和出砂量，文獻(xiàn)[25]基于室內(nèi)試驗(yàn)和理論分析了儲(chǔ)層的物性特征，其結(jié)果表明儲(chǔ)層物性的差異形成不同的巖石內(nèi)部結(jié)構(gòu)，砂巖顆粒在油藏流動(dòng)過程中的運(yùn)動(dòng)軌跡就會(huì)有差異。與此同時(shí)，有些學(xué)者[67]分析了不同的顆粒形狀和黏結(jié)物質(zhì)對(duì)儲(chǔ)層砂巖的孔隙率和滲透性的影響較大，為研究細(xì)觀結(jié)構(gòu)對(duì)儲(chǔ)層巖石力學(xué)的影響提供了重要的技術(shù)手段。上述研究成果說明不同物性儲(chǔ)層的組成結(jié)構(gòu)不同，外力的擾動(dòng)導(dǎo)致儲(chǔ)層結(jié)構(gòu)發(fā)生變化，巖層的力學(xué)響應(yīng)不同，出砂機(jī)理也有較大的差異。因此，需要對(duì)不同物性儲(chǔ)層的出砂機(jī)理進(jìn)行研究，以解決出砂與產(chǎn)能提高之間的矛盾。

目前，關(guān)于出砂機(jī)理的研究主要基于室內(nèi)試驗(yàn)、理論計(jì)算和數(shù)值分析。由于室內(nèi)試驗(yàn)的局限性和理論分析的簡化，目前的數(shù)值分析大多基于宏觀力學(xué)理論，如汪緒剛[8]建立連續(xù)介質(zhì)模型研究油井出砂，未能反映出砂過程中砂巖顆粒的運(yùn)動(dòng)特征。隨著數(shù)值計(jì)算理論和計(jì)算機(jī)性能的提高，一些學(xué)者從細(xì)觀角度研究砂巖顆粒的運(yùn)動(dòng)，如Papamichos等[9]、Jensen等[10]、Li等[1112]建立了笛卡爾坐標(biāo)系下的出砂模型，描述二維達(dá)西滲流作用下的砂巖顆粒運(yùn)動(dòng)，未考慮顆粒間的膠結(jié)對(duì)砂巖顆粒運(yùn)動(dòng)的影響，但為細(xì)觀角度研究砂巖出砂提供了重要的技術(shù)手段。

對(duì)于實(shí)際儲(chǔ)層的射孔試驗(yàn)，基于笛卡爾坐標(biāo)系模擬射孔的力學(xué)邊界和滲流邊界較復(fù)雜，為了提高計(jì)算結(jié)果的可靠性，基于劉先珊[1314]提出的三維顆粒流數(shù)值模型分析油藏開采過程中的砂巖力學(xué)響應(yīng)。為了更準(zhǔn)確地反映不同物性的儲(chǔ)層對(duì)出砂的影響，主要考慮砂巖顆粒的粒徑、粒度分布、顆粒間的膠結(jié)物質(zhì)等不同的兩種砂巖儲(chǔ)層，在圍壓和油藏流速相同的情況下，模擬油藏流動(dòng)過程中的流固耦合效應(yīng)，分析砂巖的宏細(xì)觀力學(xué)響應(yīng)以及出砂的發(fā)生發(fā)展過程，其成果可為實(shí)際儲(chǔ)層開采過程中的出砂預(yù)測和砂控優(yōu)化提供新的研究思路。

1數(shù)值建模

模擬兩種儲(chǔ)層砂巖，基于射孔的特征建立射孔的物理模型，如圖1所示，兩個(gè)模型的幾何尺寸為：外徑75 mm，內(nèi)徑6 mm和高度12 mm。根據(jù)砂巖顆粒的分布曲線（如圖2所示）生成砂巖顆粒，建立柱坐標(biāo)系下的三維顆粒流數(shù)值模型[1314]，如圖3所示的俯視圖，兩種儲(chǔ)層的結(jié)構(gòu)特性如表1所示。圖2、圖3的顆分曲線、顆粒分布圖及表1的顆粒特性比較顯示：第1種儲(chǔ)層砂巖的粒徑較小，均勻的粒徑使得這些小顆粒的接觸顆粒過少，外力作用下易成為“漂浮”顆粒，而第2種儲(chǔ)層砂巖的粒徑分布范圍較大，初始孔隙率相對(duì)較大，外荷載作用下形成較為密實(shí)的狀態(tài)；另外，第1種儲(chǔ)層砂巖的膠結(jié)程度和膠結(jié)百分含量（采用膠結(jié)半徑比α=RA表示，為膠結(jié)平均半徑，RA為接觸處顆粒中的較小半徑）、膠結(jié)百分含量（采用顆粒膠結(jié)數(shù)/模型總接觸數(shù)表示）均較低，說明砂巖顆粒間的黏結(jié)物質(zhì)較少，且黏結(jié)性較差，同樣的外力作用下膠結(jié)物質(zhì)承擔(dān)的應(yīng)力更大，顆粒間的膠結(jié)物質(zhì)更容易破壞，砂巖顆粒成為離散的顆粒，這些離散的顆粒成為出砂的幾率更大。

由于儲(chǔ)層砂巖是由砂巖顆粒膠結(jié)而成的沉積巖，顆粒間的膠結(jié)物質(zhì)使得砂巖的受力不同于無黏結(jié)砂巖，數(shù)值計(jì)算中采用具有一定膠結(jié)程度的平行黏結(jié)模型[1314]來模擬儲(chǔ)層砂巖的膠結(jié)特性。另外，油藏開采過程中，油藏的流動(dòng)與砂巖顆粒之間相互作用，即流固耦合作用，該作用力對(duì)整個(gè)砂巖模型的受力產(chǎn)生重要的影響，數(shù)值計(jì)算中要重點(diǎn)模擬。數(shù)值計(jì)算中，基于圖1設(shè)置上、下兩個(gè)水平邊界以及兩個(gè)內(nèi)外半徑組合而成的環(huán)形邊界，將整個(gè)砂巖模型劃分為288個(gè)單元[1314]，計(jì)算參數(shù)如表1所示，其參數(shù)是對(duì)同尺寸的圓柱模型進(jìn)行數(shù)值模擬，不斷調(diào)整參數(shù)使之與試驗(yàn)曲線一致，如圖4所示，兩種儲(chǔ)層砂巖顯示數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果相吻合，其參數(shù)可作為后續(xù)分析。

2計(jì)算結(jié)果分析

對(duì)于砂巖儲(chǔ)層，即使在相同的區(qū)域，由于地質(zhì)構(gòu)造的不同，其物性特征也會(huì)有差異，油藏開采過程中砂巖的力學(xué)響應(yīng)不同，出砂的發(fā)生和發(fā)展過程也有差異。本文模擬了兩種儲(chǔ)層砂巖在圍壓10 MPa、油藏流速4 m/s時(shí)的宏細(xì)觀力學(xué)響應(yīng)，比較分析了兩種儲(chǔ)層砂巖的出砂特性。

2.1砂巖宏觀力學(xué)特性的研究

射孔圍壓相同時(shí)，油藏流動(dòng)對(duì)顆粒的作用力導(dǎo)致砂巖應(yīng)力狀態(tài)的改變，顆粒間的膠結(jié)越弱，越容易失去膠結(jié)而成為離散顆粒，并被流體攜帶形成出砂。圖5比較了兩種不同儲(chǔ)層砂巖在相同圍壓和流速時(shí)的切向應(yīng)力σθσro隨r/Ri（r代表任意位置到射孔中心點(diǎn)的距離）的變化曲線。其結(jié)果顯示兩者的應(yīng)力曲線差別較大，第1種儲(chǔ)層砂巖相較于第2種砂巖的應(yīng)力峰值較大，且峰值應(yīng)力點(diǎn)對(duì)應(yīng)的位置離射孔中心點(diǎn)較遠(yuǎn)，說明在相同的賦存環(huán)境中，油藏運(yùn)動(dòng)對(duì)第1種儲(chǔ)層砂巖力學(xué)特性的影響較大。如表1所示第1種儲(chǔ)層砂巖的粒徑較小，與小粒徑顆粒相接觸的顆粒數(shù)較少，這些小粒徑顆粒在外力作用下更容易成為離散顆粒；另外，第1種儲(chǔ)層砂巖的膠結(jié)百分含量和膠結(jié)程度均較低，同樣的外載環(huán)境中顆粒膠結(jié)的破壞更容易一些。對(duì)于本文砂巖賦存環(huán)境，圍壓相同時(shí)，初始時(shí)刻相同的油藏流速對(duì)兩種砂巖顆粒的作用力是相同的，隨著流固耦合的進(jìn)行，由于第1種儲(chǔ)層砂巖膠結(jié)物質(zhì)的黏結(jié)性要差一些，膠結(jié)破壞更嚴(yán)重，顆粒的自由度較大，粒徑較小的顆粒隨流體進(jìn)入射孔中，使得承擔(dān)外力作用的顆粒數(shù)減少，砂巖的應(yīng)力增大更多，整個(gè)砂巖模型對(duì)應(yīng)的塑性區(qū)也更大一些。圖5的結(jié)果表明，在相同的外界條件下，由于儲(chǔ)層物性特征的不同，特別是砂巖顆粒間的膠結(jié)物質(zhì)對(duì)砂巖力學(xué)響應(yīng)的影響很大，由于顆粒間的弱膠結(jié)性，油藏的運(yùn)動(dòng)對(duì)顆粒的作用力更容易削弱顆粒的強(qiáng)度特性，顆粒間膠結(jié)破壞越嚴(yán)重，越容易在流體的作用下攜帶進(jìn)入油井而形成出砂。因此，不同物性儲(chǔ)層的出砂發(fā)生和發(fā)展過程有差異，實(shí)際的油藏開采中，即使是在同一區(qū)域，也會(huì)由于儲(chǔ)層砂巖細(xì)觀結(jié)構(gòu)的不同導(dǎo)致出砂不同，因此不能一概而論，需要針對(duì)具體的儲(chǔ)層出砂給出具體的出砂預(yù)測和防砂措施。

圖6和圖7的結(jié)果還可以通過分析整個(gè)砂巖的分區(qū)應(yīng)力來解釋。分區(qū)應(yīng)力以已有計(jì)算結(jié)果為基礎(chǔ)，區(qū)域內(nèi)的顆粒應(yīng)力加權(quán)平均后作為該區(qū)域的平均應(yīng)力。圖5和圖6顯示兩種砂巖的塑性區(qū)距離射孔中心點(diǎn)距離大約18 mm和12 mm，可選擇每個(gè)區(qū)域的厚度為3 mm，則將模型分成23個(gè)區(qū)域，如圖8所示。分區(qū)數(shù)量對(duì)塑性區(qū)范圍和砂巖受力特性無影響，但影響圖10曲線的比較：分區(qū)太多，每個(gè)區(qū)域的應(yīng)力增量相對(duì)較小，顯示在圖10中的應(yīng)力曲線較接近，影響了不同材料響應(yīng)曲線的比較，如第一種材料的塑性區(qū)距離中心點(diǎn)大約18 mm，每個(gè)區(qū)域的厚度若為1 mm，分區(qū)數(shù)量至少大于12個(gè)，才能獲得如圖10所示的塑性區(qū)附近的應(yīng)力曲線；而分區(qū)太少，塑性區(qū)太小的儲(chǔ)層力學(xué)響應(yīng)曲線在圖10中無法顯示。

根據(jù)數(shù)值計(jì)算結(jié)果，建立圖9～10來分析偏應(yīng)力σθ-σr與平均應(yīng)力p之間的關(guān)系，并與試驗(yàn)得到的屈服包絡(luò)線和峰值包絡(luò)線進(jìn)行比較研究。圖9顯示相同的油藏流速條件下，第1種儲(chǔ)層砂巖的偏應(yīng)力曲線與屈服曲線相交更為容易，且當(dāng)流速為5 m/s時(shí)，第1種儲(chǔ)層砂巖已與峰值曲線相交，而第2種儲(chǔ)層砂巖的偏應(yīng)力曲線與屈服曲線和峰值曲線還比較遠(yuǎn)，砂巖破壞的幾率較小。另外，圖10描述了砂巖模型不同分區(qū)的偏應(yīng)力分布，結(jié)果表明第1種儲(chǔ)層砂巖與屈服曲線的相交更為容易，距離射孔中心18 mm的區(qū)域均為塑性區(qū)，與圖5的結(jié)果一致，而第2種儲(chǔ)層砂巖與屈服曲線相差較遠(yuǎn)，砂巖屈服破壞的可能性較小，塑性區(qū)也只位于油井附近，主要是流體攜帶離散顆粒到射孔附近，使得射孔堵塞而導(dǎo)致砂巖應(yīng)力的增大，砂巖局部破壞。根據(jù)上述的偏應(yīng)力變化規(guī)律可知，由于儲(chǔ)層的物性特征不同，整個(gè)砂巖的屈服包絡(luò)線和峰值包絡(luò)線也有差異，對(duì)于砂巖顆粒膠結(jié)較弱和和膠結(jié)含量較少的儲(chǔ)層，相同外力作用下的顆粒膠結(jié)破壞越嚴(yán)重，砂巖骨架破壞越嚴(yán)重，形成的出砂通道越多，離散顆粒被流體攜帶的可能性越大，出砂越多。

2.2砂巖細(xì)觀力學(xué)特性的研究

上述砂巖的宏觀應(yīng)力圖形說明了在圍壓和油藏運(yùn)動(dòng)的共同作用下，砂巖更容易發(fā)生剪切破壞，且膠結(jié)較弱的砂巖模型更容易破壞，砂巖顆粒從砂巖母體上剝離更容易，出砂量更多。為了更清晰地解釋儲(chǔ)層的出砂機(jī)理，本文從細(xì)觀角度研究砂巖顆粒的力學(xué)特性。

圖11和圖12描述了砂巖顆粒粘結(jié)的平均張拉應(yīng)力和剪應(yīng)力的變化。兩個(gè)圖形均表明第1種儲(chǔ)層砂巖的應(yīng)力較大，說明油藏的運(yùn)動(dòng)對(duì)該儲(chǔ)層砂巖應(yīng)力的影響較大，主要在于該儲(chǔ)層顆粒較小，與之連接的顆粒數(shù)較少，外力作用下成為自由顆粒的可能性要大；另外，該儲(chǔ)層的膠結(jié)含量和膠結(jié)程度要小，其膠結(jié)物質(zhì)的黏結(jié)性相比第2種儲(chǔ)層要小，相同的外力作用下，第1種儲(chǔ)層砂巖的膠結(jié)物質(zhì)更容易破壞，并被流體攜帶至射孔內(nèi)，能承擔(dān)外力作用的砂巖顆粒逐漸減少，顆粒上的接觸應(yīng)力增強(qiáng)，則平行粘結(jié)上的應(yīng)力增大較多。圖12與圖11比較可能，在砂巖模型相同的位置，黏結(jié)上的剪應(yīng)力要比張拉應(yīng)力大，在離射孔中心點(diǎn)較遠(yuǎn)的位置，剪應(yīng)力仍然較大，說明第1種儲(chǔ)層砂巖的屈服破壞主要取決于剪應(yīng)力的狀態(tài)，與前述的圖9、圖10相吻合。

圖13和圖14從細(xì)觀角度描述了顆粒的位移。圖13的結(jié)果可知射孔附近的砂巖顆粒位移較大，且第1種砂巖模型對(duì)應(yīng)的位移要大一些。其結(jié)果同樣說明第1種儲(chǔ)層砂巖由于顆粒間的膠結(jié)性較差，相同的賦存環(huán)境中，砂巖受剪破壞更容易，離散的顆粒越多，則出砂幾率越大。同時(shí)，圖14描述了不同粒徑的顆粒對(duì)應(yīng)的平均位移，其結(jié)果說明直徑越小的顆粒其位移值越大，且對(duì)于相同的砂巖粒徑，第1種儲(chǔ)存砂巖顆粒的位移要大，其結(jié)果同樣說明粒徑較小的顆粒周圍接觸的顆粒較少，顆粒與顆粒間的膠結(jié)性較差，砂巖更容易破壞，形成的顆粒運(yùn)移通道更多，小尺寸顆粒更容易在這些通道中移動(dòng)，出砂更容易。

另外，還分析了兩種儲(chǔ)層砂巖的顆粒轉(zhuǎn)動(dòng)，如圖15所示。其結(jié)果顯示第1種儲(chǔ)層砂巖顆粒的轉(zhuǎn)動(dòng)較大，且轉(zhuǎn)動(dòng)較大值的范圍大，而第2中砂巖顆粒轉(zhuǎn)動(dòng)較大值的范圍小，只出現(xiàn)在油井附近，兩者的比較說明第1種儲(chǔ)層砂巖更容易出砂。其結(jié)果同樣說明在外界條件相同時(shí)，儲(chǔ)層物性對(duì)砂巖出砂的影響較大，由于第1種儲(chǔ)層砂巖顆粒較小，且顆粒間的膠結(jié)性較弱，油藏流動(dòng)對(duì)顆粒的拖曳力較大，應(yīng)力增大使得砂巖屈服破壞范圍加大，則顆粒間的平行粘結(jié)破壞就越多，失去膠結(jié)的砂巖顆粒也越多，被流體攜帶至油井的砂巖顆粒越多，與圖13和圖14描述的變化規(guī)律一致。

上述圖形顯示外界條件一定時(shí)，儲(chǔ)層砂巖的物性特征不同表征的出砂力學(xué)響應(yīng)不同。特別是對(duì)小粒徑且膠結(jié)性較弱的儲(chǔ)層砂巖，顆粒與顆粒較少的連接使得顆粒的自由度更大，且較差的膠結(jié)性使得砂巖更容易破壞，從砂巖母體上剝離的顆粒越多，砂巖骨架破壞越嚴(yán)重，則油藏?cái)y帶砂巖顆粒至油井的幾率越大。上述數(shù)值模擬反映了油藏開采過程中的砂巖力學(xué)響應(yīng)和出砂特征，是可行的。在實(shí)際的油藏開采中，即使在同一區(qū)域，由于位置的不同也會(huì)導(dǎo)致儲(chǔ)層物性特征的不同，出砂的發(fā)生和發(fā)展過程也會(huì)有差異，因此需要針對(duì)具體的儲(chǔ)層開采進(jìn)行出砂預(yù)測和防砂優(yōu)化。

3結(jié)論

基于柱坐標(biāo)系下的三維顆粒流數(shù)值模型，分析了不同物性的儲(chǔ)層砂巖力學(xué)響應(yīng)，據(jù)此分析出砂的發(fā)生和發(fā)展過程。

1）砂巖的宏觀應(yīng)力結(jié)果顯示油藏開采過程中，顆粒較小且膠結(jié)性較差的第1種儲(chǔ)層砂巖應(yīng)力變化更明顯，且峰值應(yīng)力和塑性區(qū)也較大。其結(jié)果表明油藏的運(yùn)動(dòng)對(duì)該儲(chǔ)層的干擾較大，對(duì)顆粒產(chǎn)生的拖曳力較大，顆粒膠結(jié)破壞多，砂巖骨架破壞嚴(yán)重，形成的通道越多，離散顆粒在這些通道中的移動(dòng)更容易，出砂幾率較大。

2）砂巖模型的細(xì)觀力學(xué)圖形描述了砂巖黏結(jié)應(yīng)力的變化以及顆粒的位移和轉(zhuǎn)動(dòng)。其結(jié)果說明第1種儲(chǔ)層砂巖的應(yīng)力較大，且運(yùn)動(dòng)更劇烈。主要在于第1種儲(chǔ)層的顆粒較小，膠結(jié)性較差，油藏運(yùn)動(dòng)中顆粒膠結(jié)物質(zhì)破壞嚴(yán)重，從砂巖母體上剝離的顆粒增多，增大了被油藏?cái)y帶進(jìn)入油井的幾率。

上述研究結(jié)果表明，對(duì)于復(fù)雜賦存環(huán)境下的油藏開采，需要針對(duì)儲(chǔ)層物性特征的差異進(jìn)行出砂機(jī)理分析，才能提出有效的出砂預(yù)測和防砂優(yōu)化方法，以實(shí)時(shí)改善油藏的開發(fā)效果。

參考文獻(xiàn)：

[1]范小軍， 黃勇， 厚東琳， 等. 儲(chǔ)層物性與產(chǎn)能的關(guān)系以元壩長興組礁灘相儲(chǔ)層為例[J]. 天然氣技術(shù)與經(jīng)濟(jì)， 2011，5（4）： 2932.

[2]吳斌， 王興志， 張帆，等. 川北地區(qū)飛仙關(guān)組儲(chǔ)層特征及其主控因素[J]. 現(xiàn)代地質(zhì)， 2012，26（1）： 168174.

[3]梁旭， 鄧宏文， 秦雁群，等. 大慶長垣泉三、四段扶余油層儲(chǔ)層特征與主控因素分析[J]. 特種油氣藏， 2012， 19（1）：5861.

[4]李彬， 賀凱， 何治亮，等. 準(zhǔn)東地區(qū)石炭系風(fēng)化殼儲(chǔ)層特征及影響因素[J]. 石油天然氣學(xué)報(bào)，2012， 34（1）： 2024.

[5]李嬌娜， 盧雙舫， 郝愛萍. 勝利油區(qū)車古20潛山儲(chǔ)層特征及儲(chǔ)層評(píng)價(jià)[J]. 復(fù)雜油氣藏， 2012， 5（1）： 15：18.

[6]周樺， 黃思靜， 蘭葉芳. 華慶地區(qū)長油層組黏土礦物類型及其對(duì)儲(chǔ)層物性的影響[J]. 巖性油氣藏， 2012， 24（3）： 6673.

[7]李龍， 鞠斌山， 江懷友，等. 油層微粒運(yùn)移及其對(duì)儲(chǔ)層物性的影響[J]. 中外能源， 2011， 16（12）：5054.

[8]汪緒剛， 鄒洪嵐， 李國誠， 等. 蘇丹6 區(qū)稠油有限攜砂冷采數(shù)值模擬及采油工藝技術(shù)[J]. 遼寧工程技術(shù)大學(xué)學(xué)報(bào)： 自然科學(xué)版， 2009， 28（Sup）： 131134.

[9]Papamichos E， Vardoulakis I， Tronvoll J， et al. Volumetric sand production model and experiment [J]. International Journal for Numerical and Analytical Methods in Geomechanics， 2001， 25（8）： 789808.

[10]Jensen R P， Preece D S， Cook B K， et al. Modeling sand production with darcyflow coupled with discrete elements [C]//Proceedings of the 10th International Conference on Computer Methods and Advances in Oeomechanics， Desai， EeL Tucson， Arizona， 2001： 819822.

[11]Li L， Holt R. Particle scale reservoir mechanics [J]. Oil and Gas Science and Technology， 2002，57（5）：525538.

[12]Li L， Papamichos E， Cerasi P. Investigation of sand production mechanics using DEM with fluid flow [C]//Eurock 2006： Multiphysics Coupling and Long Term Behavior in Rock Mechanics：241247.

篇9

流固耦合力學(xué)是一門比較新的力學(xué)邊緣分支，是流體力學(xué)與固體力學(xué)二者相互交叉而生成的。它的研究對(duì)象是固體在流場作用下的各種行為以及固體變形或運(yùn)動(dòng)對(duì)流場的影響。流固耦合力學(xué)的重要特征是兩相介質(zhì)之間的相互作用：固體在流體動(dòng)載荷作用下產(chǎn)生變形或運(yùn)動(dòng)，而固體的變形或運(yùn)動(dòng)又反過來影響到流場，從而改變流體載荷的分布和大小。流固耦合問題涉及流體和固體分析理論計(jì)算和其之間的耦合關(guān)系，所以流體力學(xué)構(gòu)成流固耦合理論的基礎(chǔ)，計(jì)算流體動(dòng)力學(xué)則是關(guān)于流體的數(shù)值計(jì)算，而建立流體力學(xué)的有限元方程求解卻極為困難，因此過去三十年人們一直在研究的流體數(shù)值分析的方法。流固耦合問題是強(qiáng)非線性的問題，求解方程的規(guī)模加大，還必須研究、借鑒有關(guān)計(jì)算機(jī)處理技術(shù)，如有限元并行算法的成果[2]。本文以攪拌機(jī)葉片系統(tǒng)作為研究對(duì)象，在SolidWorks中建立三維模型。在ANSYSWorkbench中進(jìn)行單向流固耦合，其基本思路是：先計(jì)算流場和固體結(jié)構(gòu)，然后通過中間平臺(tái)交換耦合量。每次大迭代中，進(jìn)行一次流體計(jì)算，并交換數(shù)據(jù)到固體計(jì)算，直到最終收斂。分析基本流程框架結(jié)構(gòu)如下。

2攪拌機(jī)葉片系統(tǒng)前處理

Pro/E、I-Deas、UG和SolidWorks是目前比較流行的CAD軟件，這些三維建模軟件其實(shí)原理都是相同的，各有利弊。本文采用SolidWorks軟件對(duì)攪拌機(jī)葉片系統(tǒng)進(jìn)行三維實(shí)體模型繪制，SolidWorks適合于通用機(jī)械的三維設(shè)計(jì),非常適合攪拌機(jī)這樣的通用機(jī)械產(chǎn)品設(shè)計(jì),具有與其他CAE、CAM軟件的良好接口，避免了轉(zhuǎn)換格式所造成的信息丟失和模型缺損。在三維建模軟件中,所建的三維模型完全尊重?cái)嚢铏C(jī)葉片系統(tǒng)的原型，這樣的優(yōu)點(diǎn)是，模擬更加接近真實(shí)情況。完成建模后，因?yàn)镾olidWorks軟件和ANSYSWorkbench共享數(shù)據(jù)，可以直接相連接，即在SolidWorks軟件主界面有ANSYSWorkbench菜單項(xiàng)，可直接啟動(dòng)ANSYSWorkbench12.0，模型導(dǎo)入后，一般不會(huì)發(fā)生曲面丟失、破面、形狀變形等問題。從SolidWorks導(dǎo)入到ANSYSWorkbench的攪拌機(jī)葉片系統(tǒng)流固耦合三維模型如圖2所示,內(nèi)部為攪拌機(jī)葉片系統(tǒng)，外部為管狀流體場，包括靜流場和旋轉(zhuǎn)流場，包裹著攪拌機(jī)葉片系統(tǒng)。ANSYSWorkbench工作允許添加自定義的模板和提供預(yù)先定義的自定義模板，如單向耦合分析和形變分析等。建立流固耦合分析，操作將Workbench12左側(cè)的Toolbox內(nèi)FSI:FluidFlow(CFX)>StaticStructural項(xiàng)直接拖到右邊的A2欄內(nèi)即可。ANSYSWorkbench里面劃分網(wǎng)格可選取系統(tǒng)默認(rèn)值，不需要操作者去選擇，只要控制網(wǎng)格劃分方法和有限元單元尺寸，模型局部也可以細(xì)化，這比ANSYS的經(jīng)典界面方便很多，而且復(fù)雜模型網(wǎng)格質(zhì)量也比較好。對(duì)于旋轉(zhuǎn)葉片的網(wǎng)格，需要對(duì)其局部的網(wǎng)格進(jìn)行細(xì)化，以使計(jì)算結(jié)果更加精準(zhǔn)。一般來說，葉片和流場的網(wǎng)格是一起被劃分的，但可以分開劃分進(jìn)行尺寸設(shè)定，本文所舉的例子是流固耦合分析，流體場尺寸遠(yuǎn)遠(yuǎn)大于葉片的尺寸，這樣設(shè)置的原因也是為了能夠更加真實(shí)地模擬流體場和加快分析計(jì)算速度。為便于后期的CFX操作和流固藕合分析，故對(duì)相關(guān)而先對(duì)有限元模型各表面進(jìn)行命名。網(wǎng)格劃分完畢，最后得有限元模型有190952個(gè)單元和35113個(gè)節(jié)點(diǎn)，有限元模型。

3CFX流場分析

CFX軟件主要包括三個(gè)部分：前處理模塊（CFX-Pre）、求解模塊（CFX-Solver）和后處理模塊（CFX-Post）。CFX-Pre主要用于定義物理模型、材料屬性、邊界條件、初始條件和求解參數(shù)等，樹形結(jié)構(gòu)方便用戶修改參數(shù)；CFX-Solver為求解問題所有變量的整個(gè)過程，輸入文件有CFX-Pre生成，特別地，該solver采用耦合求解器，比傳統(tǒng)分離求解器需要更少的求解步數(shù)就能達(dá)到收斂，另外，CFX-solver還有一個(gè)控制管理計(jì)算任務(wù)的manager，方便用戶管理作業(yè)。CFX-Post可以將計(jì)算結(jié)果以彩色等值線、梯度、矢量等方式顯示，也可以將計(jì)算結(jié)果以圖表、曲線的形式顯示。用CFX對(duì)流場進(jìn)行數(shù)值模擬，流場的設(shè)置主要是對(duì)水的性質(zhì)進(jìn)行設(shè)置，例如溫度、密度、傳熱系數(shù)等等。本例所選的是k-Epsilon流體方程模型，室溫條件和isothermal熱傳遞模式。然后給定流場邊界條件：入口速度為0；出口因?yàn)槭桥c空氣相同，通常這里的壓強(qiáng)都設(shè)置成相對(duì)大氣的壓強(qiáng)；其它為固定的無滑移壁面(wall)[3]。葉片部分設(shè)置的邊界條件：本例中，葉片在風(fēng)場中模擬的轉(zhuǎn)速為80r/min；葉片的固體屬性。在CFX里，葉片表面設(shè)置成流固耦合接觸面（interface），這樣就更能符合實(shí)際，模擬出其固有屬性。在獲得CFX流場數(shù)值分析的結(jié)果后，如圖4所示，可以發(fā)現(xiàn)攪拌機(jī)旋轉(zhuǎn)葉片附近的速度變化，最高變?yōu)?.7m/s，葉片正面的壓強(qiáng)值比背面的要大很多，其最大壓強(qiáng)為1234pa。

4結(jié)構(gòu)靜力分析

依據(jù)在ANSYSWorkbench中建立的CAD模型（如圖3所示），在StaticStructural中定義材料性質(zhì)，密度7800kg/m^3、泊松比0.3和彈性模量為207GPa；更新網(wǎng)格。添加邊界條件，添加攪拌機(jī)的支架頂部為固定約束，添加CFX得出的瞬態(tài)壓強(qiáng)載荷[4]。計(jì)算結(jié)果如圖5、圖6所示，由于壓應(yīng)力場相對(duì)較小，攪拌機(jī)葉片最大應(yīng)力發(fā)生在葉片前端內(nèi)側(cè)邊緣處（中部偏下），其最大值為40.7MPa,遠(yuǎn)低于鋼材的屈服強(qiáng)度。根據(jù)以上計(jì)算分析，機(jī)身強(qiáng)度符合要求。這樣的分析結(jié)果，可以讓我們對(duì)葉片結(jié)構(gòu)進(jìn)行優(yōu)化設(shè)計(jì)，在后面的設(shè)計(jì)中對(duì)葉片加筋板，起到穩(wěn)定加固葉片結(jié)構(gòu)的作用。此外，在StaticStructural的分析基礎(chǔ)上，ANSYSWorkbench還可以進(jìn)行模態(tài)分析(即振動(dòng)分析)，對(duì)葉片的振行進(jìn)行分析，得出其固有頻率和相應(yīng)振型，對(duì)葉片結(jié)構(gòu)進(jìn)行更深入的分析[5]。

篇10

數(shù)學(xué)模型：1）控制方程：根據(jù)不可壓縮黏性流體非定常流動(dòng)的Navier-Stokes方程，選用kε雙方程湍流模型對(duì)環(huán)冷機(jī)內(nèi)流動(dòng)換熱規(guī)律進(jìn)行研究?？梢詫h(huán)冷機(jī)問題整體求解方程描述為：連續(xù)性方程：()=0+jjuxρτρ(1)動(dòng)量傳輸方程：ijiijijjigfxuuxu+=+(ρ)(ρ)pτ(2)式中：ρ為流體密度；ui為流體在i方向的速度；τ為冷卻時(shí)間；pij為表面壓力矢量，包括靜壓力和流體黏性壓力；gi為作用于單位體積流體在i方向的體積力；fi為作用于單位體積流體的反方向的阻力；u為床層顆粒間隙內(nèi)的氣體流速，由表觀流速ub與空隙率ε決定：u=ub/ε。采用壓力沿床層線性分布的假設(shè)，利用Darcy定律計(jì)算氣體的表觀流速：()bL0u=Kp/z=Kpp其中，pL和p0分別為臺(tái)車進(jìn)出口壓力；滲透系數(shù)fK=k/μ，滲透率k用Ergun關(guān)系式[6]計(jì)算：k=/[150/(1)]322εεpd。能量方程利用局部非熱力學(xué)平衡換熱理論，建立氣固兩相換熱雙方程，使用編寫的用戶自定義函數(shù)(UDF)進(jìn)行數(shù)值計(jì)算。2）局部非熱平衡能量雙方程：Coberly等[7]采用局部熱力學(xué)平衡方程對(duì)二維偽均質(zhì)模型進(jìn)行研究，忽略了氣固兩相之間的溫差；DeWasch等[89]研究表明，只有當(dāng)氣固兩相溫差很小且畢渥數(shù)小于0.05時(shí)，局部熱力學(xué)平衡方程可以用于簡化的一維和二維模型，但不能滿足環(huán)冷機(jī)中的氣固換熱問題。Wakao等[10]研究表明：氣固兩相熱容和熱導(dǎo)率相差較大時(shí)，各相局部溫度變化率會(huì)明顯不同。本研究將氣相溫度Tf和固相溫度Ts作為2個(gè)獨(dú)立的變量，分別表征同一特征單元每相的熱狀態(tài)，把多孔結(jié)構(gòu)內(nèi)的傳熱視為兩相之間的傳熱，得到通用方程組[1112]如式(4)和(5)所示：固相：=τερss(1)()Tc(1)()(1)()sssvsfελT+εqhTT(4)氣相：==ffff()(c)uTTcppρτερ()()fffvsfελT+εq+hTT(5)式中，qs和qf分別為固相和氣相發(fā)熱源項(xiàng)；Tf為氣相溫度；Ts為固相溫度；hv和h分別為固相骨架與流動(dòng)介質(zhì)之間的單位體積與單位表面積的對(duì)流傳熱系數(shù)。hv可由Achenbach準(zhǔn)則關(guān)系式確定：ph6h(1ε)/dv=(6)h由下式確定[13]：1/31/2fNuhd/2.00.6PrRep=λ=+ffPrcv/λp=，ffReεdu/λp=(7)其中：Nu，Pr和Re分別為始塞爾數(shù)，氣體普朗特?cái)?shù)和雷諾數(shù)；vf為流體的運(yùn)動(dòng)黏性系數(shù)；uf為流體速度；λf為流體熱導(dǎo)率；cp為流體的比熱容。4）邊界條件與初始條件：邊界條件：Logtenberg等[4,14]認(rèn)為應(yīng)將環(huán)冷機(jī)篦板壁面邊界條件設(shè)為流體溫度。流體出口溫度與壓力均滿足第二類邊界條件：0f=zT，=0zp。初始條件：當(dāng)環(huán)冷機(jī)運(yùn)行在余熱循環(huán)利用區(qū)時(shí)(即τ＜τ循環(huán))，氣相溫度Tf為循環(huán)風(fēng)溫，固相溫度Ts為常數(shù)；當(dāng)環(huán)冷機(jī)運(yùn)行在非循環(huán)區(qū)時(shí)(即τ＞τ循環(huán))，Tf為自然風(fēng)溫。

模型結(jié)果驗(yàn)證

考慮到現(xiàn)場測試條件較艱苦，且固相與氣相之間較強(qiáng)的對(duì)流換熱會(huì)對(duì)環(huán)冷機(jī)臺(tái)車內(nèi)物理場測量產(chǎn)生很大的影響，故文獻(xiàn)[15]選用環(huán)冷機(jī)處于不同時(shí)刻時(shí)，出口空氣平均溫度的現(xiàn)場測試值與仿真結(jié)果數(shù)值對(duì)本研究所采用模型的正確性進(jìn)行驗(yàn)證。從表1可以看出，在數(shù)值仿真結(jié)果和測試結(jié)果之間存在不同程度的誤差。該誤差主要來源于：(1)測試期間環(huán)冷機(jī)操作參數(shù)的波動(dòng)；(2)測試時(shí)在煙罩上進(jìn)行了開孔，對(duì)環(huán)冷機(jī)內(nèi)的溫度場、速度場和壓力場產(chǎn)生了干擾破壞作用；(3)環(huán)冷機(jī)存在漏風(fēng)。但是，環(huán)冷機(jī)出口空氣溫度的數(shù)值仿真結(jié)果與實(shí)驗(yàn)測試結(jié)果的最大誤差小于10%，環(huán)冷機(jī)內(nèi)燒結(jié)礦的溫度分布與實(shí)際趨勢也基本一致，因此，可以認(rèn)為本文所建立的模型及計(jì)算結(jié)果是可靠的。

計(jì)算結(jié)果與分析

由環(huán)冷機(jī)對(duì)流換熱控制方程可以看出：物料粒徑、空隙率、進(jìn)風(fēng)溫度、進(jìn)風(fēng)速度、料層高度等都會(huì)對(duì)環(huán)冷機(jī)溫度場、流場分布產(chǎn)生影響[15]。本文主要研究不同固相顆粒粒徑對(duì)余熱利用量的影響，3種粒徑的物料沿臺(tái)車高度方向按粒徑從小到大的順序布置于上、中、下3層，試驗(yàn)工況見表2（略）。

1）溫度場分布：冷卻時(shí)間為581s時(shí)，環(huán)冷機(jī)內(nèi)物料溫度如圖2所示。由圖2可見：經(jīng)過分層布料工藝后，環(huán)冷機(jī)內(nèi)出現(xiàn)高溫區(qū)與低溫區(qū)，除工況Ⅵ(體積換熱系數(shù)按料層高度由大向小分布)外，其余工況均有明顯的高、低溫區(qū)交錯(cuò)分布現(xiàn)象。環(huán)冷機(jī)下層物料均能得到很好的冷卻，但在工況Ⅰ，Ⅱ，Ⅴ中，環(huán)冷機(jī)中層或上層靠近壁面的區(qū)域出現(xiàn)部分高溫區(qū)域，這3種工況粒徑配置的共同點(diǎn)為：中層向上層過渡時(shí)，物料粒徑均減小，即流體自中層向上層流動(dòng)時(shí)，所受到的阻力增加，于是流體更多從中間區(qū)域流出，壁面區(qū)域的物料由于冷卻不充分而出現(xiàn)高溫區(qū)。工況Ⅲ和Ⅵ上層物料粒徑最大，故換熱效果較差，出現(xiàn)較明顯的高溫區(qū)；工況Ⅲ和Ⅳ中層物料粒徑最小，換熱效果較好，故臺(tái)車中層物料冷卻效果最好；工況Ⅰ和Ⅳ下層物料粒徑最大，但由于臺(tái)車結(jié)構(gòu)影響，下層流體的物理速度最大，氣固兩相溫差最大，故換熱效果較好。為反應(yīng)臺(tái)車內(nèi)溫度分布的均勻性，表3列出了環(huán)冷機(jī)不同截面處溫度的標(biāo)準(zhǔn)差。由表3可以看出：工況Ⅰ和Ⅳ中的物料溫度分布較均勻，有利于提高燒結(jié)礦冷卻質(zhì)量。環(huán)冷機(jī)不同工況下出口截面物料溫度分布如圖3所示。由圖3可以看出：不同工況下的燒結(jié)物料在1500s之前冷卻速度較快，整個(gè)循環(huán)過程中，物料溫度隨冷卻時(shí)間呈指數(shù)形式減小。對(duì)于余熱循環(huán)利用區(qū)出口空氣的溫度T，若選取無量綱溫度()/()fsfT=TTTT作為空氣的特征溫度，定義τ=τ/(H/u)為特征冷卻時(shí)間，其中，H為物料高度。圖4所示為出口截面空氣無量綱溫度隨時(shí)間的變化曲線。通過線性回歸分析，特征溫度隨特征時(shí)間滿足指數(shù)函數(shù)關(guān)系：τeBT=A(8)式中，A反映初始階段特征溫度隨特征時(shí)間的變化速率；B反映整個(gè)冷卻過程征溫度隨特征時(shí)間的變化速率。A與B隨H/d及22Nu(1ε)H/d變化關(guān)系分別如圖5和圖6所示。由圖5和6可以看出：擬合函數(shù)變量A隨H/d線性變化，隨22Nu(1ε)H/d對(duì)數(shù)變化；擬合函數(shù)變量B不隨H/d變化。#p#分頁標(biāo)題#e#

2）流場分布：環(huán)冷機(jī)流場分布如圖7所示。由圖7可以看出：由于篦板自身的影響，流場在環(huán)冷機(jī)中的分布并不是完全對(duì)稱的，流體更多從與篦板傾斜方向一致的區(qū)域流出；由于環(huán)冷機(jī)自身結(jié)構(gòu)影響，流體會(huì)經(jīng)歷先收縮再向四周擴(kuò)張的“8”字形趨勢，這樣易使中層與上層出現(xiàn)部分高溫區(qū)；分層工藝之后與標(biāo)準(zhǔn)未分層工況的流場分布基本一致。綜合以上特點(diǎn)，縱向布料時(shí)可以考慮將小粒徑物料分布在環(huán)冷機(jī)四周。Leong等[3]研究結(jié)果表明：環(huán)冷機(jī)內(nèi)流場分布與空隙率有關(guān)，空隙率對(duì)其影響主要表現(xiàn)在環(huán)冷機(jī)底部位置，與本試驗(yàn)結(jié)果一致。3）余熱回收量：根據(jù)表2所示實(shí)驗(yàn)工況，考慮到溫度對(duì)空氣體積的影響以及進(jìn)入臺(tái)車空氣質(zhì)量守恒，以下基于定壓質(zhì)量熱容進(jìn)行余熱回收量的計(jì)算。進(jìn)風(fēng)溫度404K，風(fēng)量569.4t/h。實(shí)驗(yàn)工況Ⅰ中余熱利用量為1.519779×108kJ/h。同理可以算出其他工況時(shí)的余熱量。標(biāo)準(zhǔn)工況中沒有采用分層布料工藝。Jang等[16]研究表明，物料粒徑與空隙率越小，換熱系數(shù)越大，換熱量越多。通過對(duì)余熱回收量的分析可以得出：大粒徑物料分布在下層，小粒徑物料分布在中層更加有利于余熱的回收利用；小粒徑物料分布在下層，大粒徑物料分布在中層不利于余熱回收利用；采用分層布料工藝可以提高14%的余熱利用量。

結(jié)論