廖柯熹,馬曉磊,張中放,周衛(wèi)軍,張 瑤,任馳遠(yuǎn)(. 西南石油大學(xué) 石油與天然氣工程學(xué)院，成都 60500； . 中國(guó)石油天然氣集團(tuán)公司 塔里木油田分公司油氣運(yùn)銷(xiāo)部，庫(kù)爾勒 84000； . 西南石油大學(xué) 理學(xué)院，成都 60500)

基于邊界元法的站場(chǎng)陰極保護(hù)數(shù)值模擬軟件設(shè)計(jì)與開(kāi)發(fā)

廖柯熹1,馬曉磊1,張中放2,周衛(wèi)軍2,張 瑤2,任馳遠(yuǎn)3
(1. 西南石油大學(xué) 石油與天然氣工程學(xué)院，成都 610500； 2. 中國(guó)石油天然氣集團(tuán)公司 塔里木油田分公司油氣運(yùn)銷(xiāo)部，庫(kù)爾勒 841000； 3. 西南石油大學(xué) 理學(xué)院，成都 610500)

邊界元法在陰極保護(hù)電位分布計(jì)算方面有著明顯的優(yōu)勢(shì)，采用邊界元法對(duì)數(shù)值模型邊界進(jìn)行離散，計(jì)算節(jié)點(diǎn)少，計(jì)算速率快。站場(chǎng)陰極保護(hù)數(shù)值模擬軟件基于邊界元法對(duì)站場(chǎng)陰極保護(hù)電位分布進(jìn)行計(jì)算，得到可視化結(jié)果，通過(guò)該結(jié)果可以發(fā)現(xiàn)屏蔽和干擾位置，對(duì)輔助陽(yáng)極位置進(jìn)行優(yōu)化，從而解決站場(chǎng)陰極保護(hù)系統(tǒng)腐蝕區(qū)域不能預(yù)測(cè)的困難，使站場(chǎng)陰極保護(hù)系統(tǒng)工程更加合理有效。

邊界元法；站場(chǎng)；陰極保護(hù)；電位分布模擬

20世紀(jì)60年代國(guó)外開(kāi)始進(jìn)行區(qū)域性陰極保護(hù)的研究和應(yīng)用，我國(guó)從20世紀(jì)70年代末、80年代初開(kāi)始在油田和部分輸油站嘗試采用區(qū)域性陰極保護(hù)技術(shù)，到90年代中期該技術(shù)已相對(duì)成熟并開(kāi)始推廣應(yīng)用。2001年11月，鄯烏輸氣管道，鄯善首站成為第一個(gè)實(shí)施區(qū)域性陰極保護(hù)的壓氣站[1]。

目前，現(xiàn)場(chǎng)陰極保護(hù)施工過(guò)程中普遍采用試湊法，保護(hù)效果不太理想，往往需要進(jìn)行后期整改。隨著數(shù)值模擬技術(shù)在陰極保護(hù)領(lǐng)域的深入應(yīng)用，使得復(fù)雜環(huán)境中陰極保護(hù)效果的預(yù)測(cè)和低成本分析研究成為可能。與傳統(tǒng)方法相比，數(shù)值模擬技術(shù)在管道和站場(chǎng)陰極保護(hù)各種影響因素的干擾趨勢(shì)和規(guī)律研究方面具有無(wú)法比擬的優(yōu)勢(shì)，同時(shí)可用于陽(yáng)極位置優(yōu)化、陰極保護(hù)效果的預(yù)測(cè)和全面評(píng)價(jià)[2]。

目前，可以進(jìn)行區(qū)域陰極保護(hù)系統(tǒng)電位分布計(jì)算的軟件主要是基于有限元法和邊界元法進(jìn)行計(jì)算，基于有限元法進(jìn)行計(jì)算的軟件有FLUENT、COMSOL等，網(wǎng)格劃分復(fù)雜，計(jì)算節(jié)點(diǎn)多，不適用于現(xiàn)場(chǎng)情況。而B(niǎo)EASY是基于邊界元法進(jìn)行計(jì)算的軟件，在進(jìn)行網(wǎng)格劃分時(shí)，只需對(duì)邊界進(jìn)行劃分，計(jì)算量小，可以節(jié)省計(jì)算時(shí)間。

由于陰極保護(hù)計(jì)算通常只關(guān)心材料外表面和地表的參數(shù)分布情況，相較于有限元法，邊界元法具有只對(duì)邊界進(jìn)行離散化的特點(diǎn)，更適用于計(jì)算陰極保護(hù)等與腐蝕相關(guān)的問(wèn)題。為此，本工作設(shè)計(jì)開(kāi)發(fā)了基于邊界元法的陰極保護(hù)電位分布計(jì)算軟件——站場(chǎng)陰極保護(hù)數(shù)值模擬軟件(Numerical Simulation Software Of Cathodic Protection，簡(jiǎn)稱(chēng)NSSCP)，以期可靠分析與評(píng)價(jià)站場(chǎng)陰極保護(hù)系統(tǒng)的有效性。

1 陰極保護(hù)電位數(shù)值計(jì)算模型

埋地管道陰極保護(hù)電位控制方程和儲(chǔ)罐底板陰極保護(hù)電位控制方程的推導(dǎo)過(guò)程類(lèi)似，在此僅考慮儲(chǔ)罐底板陰極保護(hù)的電位控制方程，可以用Poisson方程表示為

2φ

式中：δ(X-Xe)表示以陽(yáng)極點(diǎn)X(x,y,z)為中心的狄拉克函數(shù)，其性質(zhì)見(jiàn)式(2)。

式中：φ是施加陰極保護(hù)后的極化電位值；x，y和z分別是計(jì)算點(diǎn)與儲(chǔ)罐底板中心的距離；Ie為陽(yáng)極點(diǎn)Xe(x,y,z)輸出的電流(即點(diǎn)源強(qiáng)度)。當(dāng)所研究的恒定電場(chǎng)中存在恒定的電源，電位滿(mǎn)足poisson方程，在給定邊界條件時(shí)，可通過(guò)求解該方程得到電位的分布。

公式(1)僅用來(lái)描述儲(chǔ)罐底板陰極保護(hù)的電位分布，但是如果想要求得具體解，就需要附加邊界條件。數(shù)值模型的邊界為Ω(研究區(qū)域)和總邊界Γ(Γ表面包圍3 Ω),總邊界Γ包括無(wú)限大地面表面邊界Γd、儲(chǔ)罐底板外表面邊界Γc、虛設(shè)半無(wú)限域土壤球冠邊界?！?。

邊界條件簡(jiǎn)化：地表面邊界Γd與大氣相接觸，外加陰極保護(hù)電流只能在土壤介質(zhì)中流動(dòng)，不能流入大氣，所以把地表面當(dāng)作絕緣面；假設(shè)無(wú)窮遠(yuǎn)處的電位為零，忽略陽(yáng)極對(duì)半無(wú)限大土壤表面的影響，所以在無(wú)窮遠(yuǎn)處表面的電流密度也為零；陰極的邊界條件采用描述電流密度J和電位φ之間函數(shù)關(guān)系的極化函數(shù)來(lái)定義

J=f(φ

式中：φeq是儲(chǔ)罐底板在土壤環(huán)境中的自腐蝕電位，f(φ-φeq)為陰極極化函數(shù)。

綜上所述得到數(shù)學(xué)模型及邊值問(wèn)題如下

2 邊界積分方程推導(dǎo)與處理

邊界元方法的求解是用格林公式對(duì)微分方程積分后再離散處理，最終通過(guò)求解關(guān)于電位和電流密度的線(xiàn)性方程組系統(tǒng)。可表示為

式中：H、G分別為通過(guò)邊界元方法定義的系數(shù)矩陣；E、I分別為節(jié)點(diǎn)電位矩陣和電流密度矢量，其矩陣中每一項(xiàng)分別為待求的電位和電流密度[3]。

2.1 邊界積分方程推導(dǎo)

以P′為圓心，以微小ε為半徑，在邊界上建立微小半球模型，設(shè)Q′點(diǎn)位于邊界上。通過(guò)積分方程式把域內(nèi)的點(diǎn)用P、Q表示，邊界上的點(diǎn)用P′、Q′的位勢(shì)值φ以及通量值q聯(lián)系起來(lái)，對(duì)區(qū)域Ω內(nèi)任意點(diǎn)均成立。在陰極保護(hù)電位分布模型中，邊界積分方程式為

C(P′)φ(P′)

所研究的模型屬于半無(wú)限域，很難在虛設(shè)半無(wú)限域土壤球冠邊界Γ∞和無(wú)限大地面Γd表面邊界上劃分網(wǎng)格。因此根據(jù)模型已知的邊界條件，對(duì)模型的邊值問(wèn)題進(jìn)行簡(jiǎn)化。在邊界?！尢?，滿(mǎn)足函數(shù)值和函數(shù)的法向?qū)?shù)都為零，即：

代入后可以發(fā)現(xiàn)，不需要在邊界Γ∞上劃分單元和計(jì)算。在邊界Γd上：

代入后不能完全消除該邊界的積分項(xiàng)，因此需要對(duì)邊界Γd進(jìn)行單元?jiǎng)澐帧Ｓ捎讦是邊界無(wú)限大地面表面，不能確定應(yīng)該取多大的邊界來(lái)劃分單元，同時(shí)劃分單元后也會(huì)增大計(jì)算量。通過(guò)鏡像法可在Γd邊界，尋找一個(gè)特定的基本解φ*，使邊界Γd無(wú)需劃分。

根據(jù)Poisson方程的基本解得到的模型為

φ*(b)=φ*(A,B)+φ*(A′,B′)=

根據(jù)模型基本解，邊界積分方程可簡(jiǎn)化為

選用模型的基本解時(shí)，無(wú)需對(duì)邊界Γd劃分單元，地表面的邊界條件自動(dòng)滿(mǎn)足，所以只需對(duì)陰極的邊界進(jìn)行劃分。

2.2 邊界離散

邊界元法將邊界上的積分方程式離散化，并歸結(jié)為代數(shù)方程組求得近似解。通過(guò)插值函數(shù)得到各個(gè)邊界單元上的函數(shù)值和函數(shù)的法向?qū)?shù)值。

按單元求和的離散形式改寫(xiě)為按節(jié)點(diǎn)求和的邊界積分方程離散格式，即

由于j點(diǎn)也是邊界上的點(diǎn)，所以必然與某I節(jié)點(diǎn)為同一點(diǎn)，所以合并同類(lèi)項(xiàng)后，公式簡(jiǎn)化為

式(14)中:

對(duì)于m個(gè)節(jié)點(diǎn)，可以得到聯(lián)立的一次方程組為

式中:H與G表示系數(shù)矩陣，Φ表示邊界單元節(jié)點(diǎn)處的函數(shù)值的列向量，Q表示邊界單元節(jié)點(diǎn)處函數(shù)法向?qū)?shù)值的列向量，按照已給定的邊界條件，由式(16)可以計(jì)算出邊界上節(jié)點(diǎn)的未知函數(shù)值與其法向?qū)?shù)值。

2.3 輔助陽(yáng)極處理

根據(jù)δ函數(shù)的性質(zhì)可知，域內(nèi)存在一個(gè)集中源點(diǎn)X1時(shí)

式中:φ*(P′,A1)是對(duì)應(yīng)A1的基本解，得

當(dāng)域內(nèi)存在多個(gè)集中源點(diǎn)An(n=1LN)時(shí)，得

φ*(P′,An)b(Q)dΩ(Q)=

積分項(xiàng)Bj不會(huì)產(chǎn)生未知量。

2.4 非線(xiàn)性邊界條件處理

在數(shù)值模擬過(guò)程中，邊界條件的確定及其處理方式直接影響計(jì)算結(jié)果的準(zhǔn)確性，軟件總共有4種邊界條件，包括極化曲線(xiàn)邊界、電位邊界、電流邊界、阻抗邊界，其中最主要的邊界條件是非線(xiàn)性極化曲線(xiàn)邊界。

由于陰極和陽(yáng)極的電位和電流密度的關(guān)系并非簡(jiǎn)單的線(xiàn)性關(guān)系，而是非線(xiàn)性關(guān)系，在模型處理過(guò)程中必須對(duì)其進(jìn)行處理。儲(chǔ)罐底板表面受到多種不同因素的影響，極化函數(shù)不能只用單獨(dú)的環(huán)境或者材料的影響函數(shù)來(lái)表達(dá)，為得到極化函數(shù)定義的陰極邊界條件J=f(φ-φeq)，通過(guò)試驗(yàn)測(cè)得極化曲線(xiàn)作為邊界條件是非常有必要的。

通過(guò)對(duì)測(cè)得的極化曲線(xiàn)進(jìn)行擬合，可以確定電流密度與極化電位之間的定量關(guān)系。主流的擬合方法有Newton-Raphson迭代法和分段線(xiàn)性擬合法。

由于被保護(hù)金屬表面受到多種不同因素的影響，所以極化函數(shù)不能只用單獨(dú)的環(huán)境或者材料的影響函數(shù)來(lái)表達(dá)，為得到極化函數(shù)定義的陰極邊界條件，綜合考慮現(xiàn)場(chǎng)實(shí)際狀況，通過(guò)試驗(yàn)測(cè)得被保護(hù)金屬在土壤中的電位分布，為了取得良好的收斂性和計(jì)算精度，該軟件采用分段擬線(xiàn)性化法處理非線(xiàn)性邊界條件。擬合得到電流密度J與電位φ的之間的線(xiàn)性關(guān)系

J=aφ

推導(dǎo)得

將式(22)代入下列矩陣方程

移項(xiàng)整理得

矩陣最終化為AX=F的形式，輔以迭代的方法，可以求解得到電位值大小。

3 軟件開(kāi)發(fā)

基于邊界元算法開(kāi)發(fā)了“站場(chǎng)陰極保護(hù)數(shù)值模擬軟件(NSSCP)”，該軟件可以實(shí)現(xiàn)三維站場(chǎng)陰極保護(hù)系統(tǒng)的模擬，獲得可視化電位分布規(guī)律，從而解決現(xiàn)有站場(chǎng)和新改擴(kuò)建站場(chǎng)內(nèi)埋地金屬構(gòu)筑物密集所造成的屏蔽和干擾問(wèn)題，并防止站場(chǎng)內(nèi)陰極保護(hù)電流對(duì)未納入陰極保護(hù)管道的干擾影響。相較于有限元法計(jì)算,該軟件的優(yōu)點(diǎn)就是簡(jiǎn)單快捷。

3.1 軟件開(kāi)發(fā)流程和計(jì)算流程

在軟件的設(shè)計(jì)開(kāi)發(fā)過(guò)程中，使用Access 2007數(shù)據(jù)庫(kù)作為后臺(tái)數(shù)據(jù)庫(kù)，并基于Visual Studio C# 2010設(shè)計(jì)了軟件用戶(hù)界面。軟件共劃分出陰極保護(hù)系統(tǒng)參數(shù)存儲(chǔ)、邊界條件處理、電位計(jì)算及數(shù)據(jù)導(dǎo)出四個(gè)主模塊。圖1和圖2分別是軟件開(kāi)發(fā)的流程圖和軟件計(jì)算流程圖。

圖1 軟件開(kāi)發(fā)流程Fig. 1 Software development process

圖2 軟件計(jì)算流程Fig. 2 Software calculation process

3.2 陰極保護(hù)系統(tǒng)參數(shù)存儲(chǔ)和建模

用戶(hù)可以在“工程參數(shù)”界面輸入工程名稱(chēng)、站場(chǎng)名稱(chēng)、土壤電阻、站場(chǎng)簡(jiǎn)介、模擬目的、精度選擇(高精度、一般精度、較低精度)、精度參數(shù)(最大迭代次數(shù)和收斂誤差閾值)等基礎(chǔ)數(shù)據(jù)。土壤電阻率會(huì)隨地域、季節(jié)和天氣發(fā)生變化，在設(shè)置土壤電阻率時(shí)，應(yīng)考慮這些因素的影響。在“實(shí)體對(duì)象”界面中可以進(jìn)行腐蝕實(shí)體對(duì)象建模，儲(chǔ)罐、陽(yáng)極、埋地管道需要設(shè)置保護(hù)層系數(shù)、直徑、埋深、位置坐標(biāo)、邊界條件等參數(shù)。

3.3 電位分布模擬計(jì)算

采用對(duì)于求解區(qū)域進(jìn)行自動(dòng)剖分網(wǎng)格的方法來(lái)劃分單元，將形成的邊界單元的坐標(biāo)保存在相應(yīng)的輸入文件中。計(jì)算系數(shù)矩陣，調(diào)用邊界條件，求解線(xiàn)性方程組，完成迭代。

軟件在使用時(shí)最重要的步驟就是進(jìn)行參數(shù)收集與導(dǎo)入，在軟件的主界面應(yīng)先進(jìn)行陰極極化曲線(xiàn)邊界參數(shù)的導(dǎo)入；其次是土壤電阻率的導(dǎo)入，需要考慮不同季節(jié)導(dǎo)致的土壤電阻率的變化，根據(jù)實(shí)測(cè)結(jié)果進(jìn)行導(dǎo)入；最后導(dǎo)入的是實(shí)體腐蝕對(duì)象的參數(shù)(包括地理位置、幾何尺寸等)，通過(guò)一系列的計(jì)算便可以得到電位分布。

干擾分析菜單欄可以實(shí)現(xiàn)對(duì)計(jì)算結(jié)果查看和分析，得到整個(gè)站場(chǎng)節(jié)點(diǎn)電位分布數(shù)據(jù)列表,見(jiàn)圖3。

干擾問(wèn)題是站場(chǎng)陰極保護(hù)系統(tǒng)實(shí)施的一個(gè)難題。

圖3 數(shù)值計(jì)算結(jié)果查看界面Fig. 3 View interface of numerical results

對(duì)油氣輸送站場(chǎng)進(jìn)行陰極保護(hù)，由于是將整個(gè)站區(qū)內(nèi)的所有地下金屬結(jié)構(gòu)全部納入保護(hù)系統(tǒng)，因此一般不會(huì)產(chǎn)生內(nèi)部干擾，但是由于站場(chǎng)陰極保護(hù)電流通常遠(yuǎn)大于干線(xiàn)保護(hù)電流，故常常對(duì)管道干線(xiàn)等外部結(jié)構(gòu)及其陰極保護(hù)系統(tǒng)造成干擾。一般系統(tǒng)電流輸出越大，干擾也會(huì)越嚴(yán)重。通過(guò)該軟件干擾分析產(chǎn)生的干擾數(shù)據(jù)可以有效判斷出是否發(fā)生發(fā)生干擾，提高了站場(chǎng)陰極保護(hù)系統(tǒng)的可靠性。

4 實(shí)例驗(yàn)證

某站場(chǎng)使用一套陰極保護(hù)系統(tǒng)對(duì)四個(gè)儲(chǔ)罐進(jìn)行保護(hù)，6組深井陽(yáng)極埋設(shè)45 m分別位于儲(chǔ)罐的兩側(cè)每組陽(yáng)極由8支研究塊組成,共8 m,直徑為25 cm，儲(chǔ)罐直徑為60 m，恒電位儀的輸出電位為8 V。對(duì)儲(chǔ)罐周邊測(cè)試樁進(jìn)行檢測(cè)，測(cè)得電位均處于-0.85～-1.2 V(相對(duì)銅/硫酸銅參比電極)，達(dá)到保護(hù)狀態(tài)，保護(hù)效果良好。

通過(guò)采用COMSOL軟件模擬結(jié)果進(jìn)行比對(duì)分析，模擬時(shí)，所取土壤電阻率為50 Ω·m，砂墊層200 Ω·m，瀝青砂層1 000 Ω·m。采用COMSOL軟件模擬得到的儲(chǔ)罐保護(hù)電位分布云圖，通過(guò)模擬計(jì)算，并將COMSOL軟件模擬結(jié)果和實(shí)測(cè)數(shù)據(jù)進(jìn)行統(tǒng)計(jì),見(jiàn)表1。

表1 計(jì)算結(jié)果比較Tab. 1 Comparison of calculated results

COMSOL基于有限元方法進(jìn)行計(jì)算，計(jì)算時(shí)間較長(zhǎng)，一般認(rèn)為其結(jié)果較為精確，但是存在外邊界截?cái)嗾`差；而NSSCP是基于邊界元法進(jìn)行開(kāi)發(fā)的，不存在邊界截?cái)嗾`差。由表1可見(jiàn)，NSSCP的計(jì)算精度接近COMSOL的，證明該軟件計(jì)算結(jié)果可靠，可以作為站場(chǎng)陰極保護(hù)系統(tǒng)的檢測(cè)維護(hù)軟件。

5 結(jié)束語(yǔ)

(1) NSSCP軟件在計(jì)算結(jié)果精度方面接近COMSOL軟件，并且計(jì)算過(guò)程中收斂速率快，計(jì)算時(shí)間短，提高了使用效率。

(2) NSSCP軟件具有較高的實(shí)用性，實(shí)際應(yīng)用效果較好，能夠用于指導(dǎo)站場(chǎng)陰極保護(hù)系統(tǒng)的設(shè)計(jì)計(jì)算和科學(xué)維護(hù)。

[1] 陳洪源,范志剛,劉玲莉,等. 區(qū)域性陰極保護(hù)技術(shù)在輸氣站場(chǎng)中的應(yīng)用[J]. 油氣儲(chǔ)運(yùn),2005,24(5):41-44.

[2] 張豐,陳洪源,李國(guó)棟,等. 數(shù)值模擬在管道和站場(chǎng)陰極保護(hù)中的應(yīng)用[J]. 油氣儲(chǔ)運(yùn),2011,30(3):208-212.

[3] 陳石,秦朝葵. 邊界元法介紹及其在陰極保護(hù)領(lǐng)域應(yīng)用[J]. 上海煤氣,2015(3):7-10.

“2017材料物理測(cè)試新技術(shù)研討會(huì)”征 文 通 知

由上海材料研究所等單位聯(lián)合主辦，《理化檢驗(yàn)-物理分冊(cè)》編輯部、《機(jī)械工程材料》編輯部共同承辦的“物理測(cè)試新技術(shù)研討會(huì)”已連續(xù)成功舉辦兩屆，受到了國(guó)內(nèi)廣大材料物理測(cè)試專(zhuān)家、學(xué)者及從業(yè)人員的熱烈歡迎和積極參與。在此基礎(chǔ)上，“2017材料物理測(cè)試新技術(shù)研討會(huì)”將繼續(xù)圍繞“新的測(cè)試技術(shù)與方法”、“新的測(cè)試儀器和設(shè)備”、“新的試驗(yàn)標(biāo)準(zhǔn)與規(guī)范”等方面展開(kāi)討論，旨在為廣大物理測(cè)試工作者以及設(shè)備廠(chǎng)商提供一個(gè)溝通、交流和展示的平臺(tái)，以總結(jié)經(jīng)驗(yàn)，促進(jìn)創(chuàng)新，進(jìn)而推動(dòng)我國(guó)材料物理測(cè)試技術(shù)的發(fā)展與進(jìn)步。本次會(huì)議將于2017年10月底在上海舉行，由上海材料研究所、中國(guó)機(jī)械工程學(xué)會(huì)理化檢驗(yàn)分會(huì)、中國(guó)機(jī)械工程學(xué)會(huì)材料分會(huì)聯(lián)合主辦，《理化檢驗(yàn)-物理分冊(cè)》編輯部、《機(jī)械工程材料》編輯部和上海市工程材料應(yīng)用與評(píng)價(jià)重點(diǎn)實(shí)驗(yàn)室共同承辦。會(huì)議將邀請(qǐng)業(yè)內(nèi)權(quán)威專(zhuān)家作專(zhuān)題報(bào)告，同時(shí)熱忱歡迎行業(yè)專(zhuān)家、學(xué)者及從業(yè)人員踴躍投稿，積極參會(huì)。所征論文將以《理化檢驗(yàn)-物理分冊(cè)》正刊專(zhuān)欄或增刊形式發(fā)表，評(píng)選出的優(yōu)秀論文將在會(huì)議期間進(jìn)行口頭交流。

1 征文內(nèi)容

(1) 材料物理測(cè)試的新技術(shù)、新方法、新標(biāo)準(zhǔn)和新儀器設(shè)備的介紹及應(yīng)用；

(2) 物理測(cè)試技術(shù)在新材料開(kāi)發(fā)、研究中的應(yīng)用；

(3) 物理測(cè)試技術(shù)在產(chǎn)品質(zhì)量控制過(guò)程中的應(yīng)用；

(4) 物理測(cè)試技術(shù)在材料表面處理工藝、材料冷熱加工工藝、材料改性技術(shù)以及先進(jìn)制造技術(shù)方面的應(yīng)用；

(5) 物理測(cè)試技術(shù)在產(chǎn)品在役安全監(jiān)控、壽命估算、可靠性評(píng)估等方面的應(yīng)用；

(6) 機(jī)械產(chǎn)品及其構(gòu)件失效分析的最新思路、技術(shù)、方法及典型的失效分析案例介紹；

(7) 試驗(yàn)設(shè)備的升級(jí)改造技術(shù)；

(8) 試驗(yàn)室的管理及在全面質(zhì)量管理中的成功經(jīng)驗(yàn)；

(9) 其他與材料物理測(cè)試技術(shù)相關(guān)的論文。

2 征文要求

(1) 未公開(kāi)發(fā)表過(guò)的原創(chuàng)性論文。

(2) 來(lái)稿要求論點(diǎn)明確、數(shù)據(jù)可靠、邏輯嚴(yán)密、文字精煉。每篇論文必須包括題目、作者姓名、作者單位、單位所在地及郵政編碼、摘要和關(guān)鍵詞、正文、參考文獻(xiàn)和第一作者及通信作者簡(jiǎn)介(包括姓名、性別、職稱(chēng)、出生年月、所獲學(xué)位、目前主要從事的工作和研究方向)，在文稿的首頁(yè)地腳處注明論文屬何項(xiàng)目、何基金(編號(hào))資助，沒(méi)有的不注明。論文題目、作者姓名、作者單位、單位所在地及郵政編碼、摘要和關(guān)鍵詞應(yīng)中英文對(duì)照。

(3) 論文摘要盡量寫(xiě)成報(bào)道性文摘，包括目的、方法、結(jié)果、結(jié)論4方面內(nèi)容(200字左右)，應(yīng)具有獨(dú)立性與自含性，關(guān)鍵詞選擇貼近文義的規(guī)范性單詞或組合詞(3～8個(gè))。

(4) 文稿篇幅(含圖表)一般不超過(guò)8 000字。文中量和單位的使用請(qǐng)參照中華人民共和國(guó)法定計(jì)量單位最新標(biāo)準(zhǔn)。外文字符必須分清大、小寫(xiě)，正、斜體，黑、白體，上下角標(biāo)應(yīng)區(qū)別明顯。

(5) 文中的圖、表應(yīng)有自明性，圖像要清晰，層次要分明；所有圖和表均應(yīng)標(biāo)明中英文圖題和表題。

(6) 參考文獻(xiàn)的著錄格式采用順序編碼制，請(qǐng)按文中出現(xiàn)的先后順序編號(hào)。所引文獻(xiàn)必須是作者直接閱讀參考過(guò)的、最主要的、公開(kāi)出版文獻(xiàn)。未公開(kāi)發(fā)表的、且很有必要引用的，請(qǐng)采用腳注方式標(biāo)明，參考文獻(xiàn)不少于5條。

3 征文截止時(shí)間

2017年8月31日。

4 論文提交

論文投稿至E-mail: pt@mat-test.com，請(qǐng)注明稿件為“物理測(cè)試會(huì)議論文”。

5 聯(lián)系方式

地址： 上海市邯鄲路99號(hào)，《理化檢驗(yàn)-物理分冊(cè)》編輯部、《機(jī)械工程材料》編輯部。

聯(lián)系人：李玲；金靜靜。

電話(huà)：021-65559079；021-65556775-361。E-mail：pt@mat-test.com。

(會(huì)議的具體時(shí)間與地點(diǎn)，及會(huì)務(wù)費(fèi)用等詳情將在第二輪通知中發(fā)出。)

Design and Development of Numerical Simulation Software for Cathodic Protection of Stations Based on Boundary Element Method

LIAO Kexi1, MA Xiaolei1, ZHANG Zhongfang2, ZHOU Weijun2, ZHANG Yao2, REN Chiyuan3
(1. College of Petroleum & Natural Gas Engineering, Southwest Petroleum University, Chengdu 610500, China;2. Tarim Oilfield Company Oil and Gas Distribution Unit, China National Petroleum Corporation, Korla 841000, China;3. College of Science, Southwest Petroleum University, Chengdu 610500, China)

Boundary element method in distributed computing of cathodic protection potential has a distinct advantage. When the boundary element method to scatter a numerical model was used, the calculation nodes were less and calculation speed was fast. The software based on the boundary element method calculated the cathodic protection potential distribution in cathodic station, and the results were visual, shielding and interference position could be found through the results, auxiliary anode position could be optimized in order to resolve the difficulty that the station cathodic corrosion protection system area could not be predicted, so the station cathodic protection system engineering could be more rational and effective.

boundary element method; station; cathodic protection; potential distribution simulation

2016-01-15

廖柯熹(1970-)，教授，博士，從事油氣儲(chǔ)運(yùn)系統(tǒng)完整性管理與安全評(píng)價(jià)工作,liaokxswpi@163.com

10.11973/fsyfh-201707015

TE88

A

1005-748X(2017)07-0551-06