橫坡敷設管道力學響應規律研究

摘要：

山體滑坡是威脅油氣管道安全運行的常見地質災害之一。橫坡敷設的管道最易受到滑坡影響發生破壞，研究滑坡作用下埋地管道的力學響應極為重要。基于彈性地基梁理論，考慮管道變形非線性和土體抗力物理非線性，建立橫坡敷設管道力學模型，推導管道撓度、軸向應力及Mises應力解析表達式，并對工程算例進行分析，得到管道的力學響應規律。通過ABAQUS軟件建立山體和管道的三維有限元模型，對比分析仿真結果和理論結果，驗證了理論模型的可靠性。在此基礎上，考慮并分析了坡面凹凸程度對管道力學響應規律的影響。研究結果表明：滑坡區域內橫坡敷設管道在主滑面處撓度最大，處在滑坡邊界附近和主滑面的管道易發生破壞；坡面凹陷程度的增加會使管道撓度和Mises應力的峰值先減小再增大；凸起坡面內橫坡敷設管道產生的最大撓度和受到的峰值Mises應力較平面坡體更小，敷設埋地管道時應盡可能避開凹陷區域。根據坡面凹凸程度對管道力學響應規律的影響，對所建立的理論模型引入形狀修正公式，建立了考慮坡面形狀的Mises應力計算表達式。研究結果可為管道安全運行和維護提供理論依據。

關鍵詞：

山地管道；橫坡敷設；力學響應規律；凹凸程度；撓度；山體坡面

中圖分類號：TE973

文獻標識碼：A

202405082

Mechanical Responses of Pipeline Laid on Cross Slope

Han Xuelong1 Wei Kai1 Li Yukun2 Xu Guoqing3 Xi Youqiang4 Guo Yong5 Zhou Peng6

（1.School of Petroleum Engineering，Yangtze University；2.College of Pipeline and Civil Engineering，China University of Petroleum （East China）；3.Dehong Branch，PipeChina Yunnan Company；4.PipeChina Beijing Pipeline Co.，Ltd.；5.Tarim Oil and Gas Transportation Branch，PipeChina West Pipeline Company Limited；6.Shandong Airlines Co.，Ltd.）

Landslide is one of the common geological hazards that threaten the operation of oil and gas pipelines.Pipelines laid on cross slopes are most susceptible to landslide and damage.It is extremely important to understand the mechanical responses of buried pipelines under the action of landslides.Based on the theory of beam on elastic foundation，and considering the nonlinear deformation of pipeline and the physical nonlinearity of soil resistance，a mechanical model of pipeline laid on cross slope was built to derive the analytical expressions of pipeline deflection，axial stress and Mises stress，analyze the engineering examples，and obtain the mechanical responses of the pipeline.Then，the ABAQUS software was used to build a 3D finite element model of mountain mas and pipeline.Next，the simulation and theoretical results were compared and analyzed to verify the reliability of the theoretical model.Finally，the influence of slope roughness on the mechanical response of pipeline was considered and analyzed.The research results show that the pipeline laid on the cross slope in the landslide area exhibits the greatest deflection at the main slide face，and is prone to damage near the landslide boundary and the main slide face.The increase in the degree of slope depression first reduces and then increases the peak values of pipeline deflection and Mises stress.The maximum deflection generated and peak Mises stress received by the pipeline laid on the convex cross slope are smaller those on the flat cross slope，and every effort should be made to avoid the concave area while laying the buried pipelines.Moreover，based on the influence of the slope roughness on the mechanical response of the pipeline，a shape correction formula was introduced into the built theoretical model to establish a Mises stress calculation expression considering the slope shape.The research results provide a theoretical basis for ensuring the safe operation and maintenance of the pipeline.

mountain pipeline；cross slope laying；mechanical response；roughness；deflection；mountain slope

0 引 言

埋地管道不可避免地穿越山區［1］，一旦發生山體滑坡，管道將處于易損或受損狀態，可能導致管道斷裂、泄漏甚至爆炸，造成嚴重人員傷亡，危害巨大［2-6］。垂直于土體下滑方向的橫坡敷設管道最易受到滑坡破壞［7-8］。我國油氣管網里程將在2025年達到24萬km［9-10］，隨著大量埋地管道的建設完成，更多的山區管道將面臨這種潛在危險。因此，掌握滑坡作用下橫坡敷設管道的力學響應規律十分必要。

勃洛達夫金［11］研究了山體橫向滑坡、縱向滑坡時管道的受力情況。B.B.RAJANI等［12］采用線彈性體模型分析了橫向滑坡的管道力學狀態，并指出管道滑動方向的土體對管道產生阻礙作用。M.J.O’ROURKE等［13］在線彈性體模型的基礎上分析了縱向滑坡下管道的力學狀態，指出管道最大主應力與管道埋深、土壤黏聚力等因素的影響關系。張東臣等［14］建立了斜向滑坡埋地管道受力分析模型，研究發現管道最大應力隨滑坡推力與管軸線夾角的增加而增大。G.COCCHETTI等［15］將管道考慮為埋地梁模型，得到三維管土作用下的橫向滑坡管道力學響應規律。梁政等［16］基于勃洛達夫金提出的模型建立了橫向滑坡模型，獲得了滑坡區域管道的變形及應力解。鄧道明等［17］、吳銳等［18］、尚玉杰等［19］、王榮有等［20］、謝強等［21］、張杰等［22］分別將管道簡化為大變形梁模型，推導計算了橫向滑坡載荷作用下埋地管道變形及內力的計算解析解，開展了管道力學響應的分析。

實際坡體存在凸起或凹陷的情況，目前關于坡體凹凸形狀這一因素對管道力學響應規律的研究較少。為此，基于Winkler彈性地基梁理論，構建橫坡敷設管道力學模型，提出管道控制方程、協調條件，推導用于計算管道變形和內力的解析解。通過解析解和有限元仿真對算例進行分析，總結出管道力學響應規律。在此基礎上，分析了不同坡面形狀（凸起、凹陷和平整）對管道變形、受力的影響，并基于數值模擬結果對力學模型解析解引入形狀修正公式，彌補了理論模型無法考慮坡面凹凸程度的缺陷。研究結果可為管道的安全運行和維護提供參考。

1 橫坡敷設管道力學模型與分析

1.1 橫坡敷設管道力學模型

橫坡敷設管道所處的山體發生滑坡時，由于土體剩余下滑力的推移作用，管道后端受到動載，土體變形分布直接影響管道受到的載荷分布。統計表明，滑坡平面形態圓弧形滑坡較多，實際滑坡橫斷面上剩余下滑力隨著靠近滑坡邊界逐漸減小至0［23］。為方便計算，本文采用拋物線分布的動載q，建立橫坡敷設埋地管道在滑坡體推移作用下的管道受力模型，如圖1所示，其中管道用加粗虛線表示。

對該力學模型做出如下說明及假設：①土體物性均勻，土體變形分布影響管道受到的載荷分布，管道自重、前端土壓力對管道變形的影響遠小于滑坡下滑力，忽略管道自身、內部輸送介質的重力及管道前端土體的阻擋作用；②管道材料沿軸向、剛度均勻，將管道視作一根無限長梁；③橫向土抗力符合Winkler彈性地基梁假設，縱向土抗力符合雙線性假設；④管道變形、受力關于滑坡區域管道中心截面對稱。不同管道地質災害（滑坡、塌陷以及沉降）作用下管道力學模型相似，此模型只針對滑坡作用下橫坡敷設管道開展討論。

圖1中：l0為滑坡區域寬度，m；M0為管道滑坡區域兩側交界處截面的彎矩，N·m；Q0為管道滑坡區域兩側交界處截面的剪力，N；ν0為管道滑坡區域兩側交界處截面的撓度，m；θ0為管道滑坡區域兩側交界處截面的轉角，rad；f為滑坡區域管道最大撓度，m；q0為管道受到的滑坡體主滑面處土體的載荷（即滑坡區域中間最大的土體載荷），N/m；N0為滑坡區域管道兩端受到的軸力，N。

1.2 橫坡敷設管道靜力學分析

將滑坡區域外部管道視為放在Winkler地基上的半無限長彈性地基梁，其受力和變形關于主滑面對稱，取B側滑坡區域外管道進行分析，管道只受地基反力的作用。滑坡區域內部管道受到拋物線分布的載荷，邊界處管道受到剪力、軸力、彎矩。

1.2.1 管道控制方程

基于圖1b坐標系，可以得到滑坡體外管道滿足如下控制方程：

EId4ydx14=-Ky=-ky0Dy（1）

式中：E為管道的彈性模量，Pa；I為管道的截面慣性矩，m4；K為地基彈性系數，N/m2；kyo為土體橫向阻力綜合系數，N/m3；D為管道外徑，m；y為管道撓度，m。

針對滑坡區域內管道，考慮到其端部當量軸力存在呈現拉力和壓力的可能性，基于力學模型，建立滑坡區域內管道兩端當量軸力分別為拉力、壓力和當量軸力為0時的控制方程：

EId2ydx2=M0-Q0x+S0y-v0+2q0x33l0-q0x43l02（2）

EId2ydx2=M0-Q0x-S0y-v0+2q0x33l0-q0x43l02（3）

EId2ydx2=M0-Q0x+2q0x33l0-q0x43l02（4）

S0=N0-pπ4d2 （S0為拉力時）

S0=pπ4d2-N0 （S0為壓力時）（5）

式中：S0為管道在x=0、x=l0截面的當量軸力，N；p為管道受到的內壓，Pa；d為管道的內徑，m。

滑坡區域內外交界截面的剪力Q0=q0l03。

1.2.2 管道協調條件

滑坡區域外管道考慮為半無限長彈性地基梁，可采用龍馭球［24］提到的邊界條件。根據管道變形對稱性和連續性，滑坡區域內外邊界處管道橫向變形邊界條件為：

y0=v0

y′0=-θ0

y′l0/2=0（6）

滑坡區域外管道在模型中視作半無限長梁。管道截面切應力的分布密度小于一定的數值時，管道與土壤處于彈性相互作用階段。達到一定數值時后，管道與土壤處于塑性相互作用階段，管道位移存在極限平衡段。需確定管道與土壤是否處于塑性相互作用階段，因此應判定是否存在極限平衡段長度。判定依據為：

N0-Nrγtfgt;1，管道存在極限平衡段

N0-Nrγtf≤1，管道不存在極限平衡段（7）

式中：Nr為管道嵌固端受到的軸向作用力，N，Nr=-αΔTE+μσhA；α為管道線膨脹系數，m/（m·℃）；ΔT為管道運行溫度與敷設溫度差，℃；μ為管道泊松比；σh為管道環向應力，Pa，σh=pd2t；t為管道壁厚，m;γ為中間計算系數，m-1，γ=πDkx0EA；kx0為土體縱向阻力綜合系數，N/m3;A為管道金屬橫截面積，m2;tf為滑坡區域外土壤的極限抗剪強度，Pa。

滑坡區域外管道存在極限平衡段時，縱向位移量u01為：

u01=12EAN0-Nr2tf+tfγ2（8）

滑坡區域外管道無極限平衡段時，縱向位移量u01為：

u01=N0-NrEAγ（9）

滑坡區域內管道由于受軸力、溫度、壓力、橫向變形影響，產生縱向位移u02，其量化公式為：

u02=12-N0LEA-αΔTL+μσhLE+π2f24L"" （10）

滑坡邊界處管道縱向變形連續，滑坡區域內外變形一致，即u01=u02。

1.2.3 管道當量軸力拉壓性質的判別

結合管道變形協調條件，求解當量軸力為0時滑坡區域內管道控制方程，以管道當量軸力為0時滑坡區域外縱向位移作為的判別依據。為方便計算，引入中間計算系數η=4K4EI，m-1。滑坡區域內管道撓曲線方程為：

yx=1EI12M00x2-118l0q0x3+q030l0x5-

q090l02x6+x30q0l03-x2M00l0+2l0q0η3K+M002EIη2（11）

M00=-20EIl0q0η3+Kηq0l0315Kηl0+30K（12）

式中：M00為滑坡區域內外交界處管道兩端當量軸力為0時的彎矩，N·m。

滑坡區域內管道最大撓度f為：

f=1EI-18M00l02+615 760q0l04+2q0ηL3K+M002EIη2（13）

當量軸力拉壓性質的判別依據為：

（1）當π2f24l0+μσhl0E-αΔTl0-πpd2l04EAgt;Uw時，S0為拉力；

（2）當π2f24l0+μσhl0E-αΔTl0-πpd2l04EA=Uw時，S0為0；

（3）當π2f24l0+μσhl0E-αΔTl0-πpd2l04EAlt;Uw時，S0為壓力。

pπ4d2-Nrγtf≤1時，Uw=2pπ4d2-NrEAγ

pπ4d2-Nrγtfgt;1時，

Uw=1EApπ4d2-Nr2q0+q0γ2（14）

式中：Uw為滑坡區域外管道相對縱向位移，m。

1.2.4 管道變形及內力計算

由控制方程和協調條件得到滑坡區域外管道變形方程為：

yx1=e-ηx12Q0ηK+M02EIη2cosηx1-

2M0η2Ksinηx1（15）

滑坡區域內外交界處截面管道撓度和轉角為：

v0=2Q0ηK+M002EIη2

θ0=-2Q0η2K-M00EIη（16）

為方便計算，設中間計算系數k=S0EI，m-1。滑坡區域內管道受到的當量軸力為拉力時，其撓曲線方程為：

y=1S0M0-8q0l02k4coshkx-

tanhkl02sinhkx+q03S0l02x4-2q03l0S0x3+

4EIq0S20l02x2+q0l03S0-4EIq0l0S20x-1S0M0-8q0l02k4+v0（17）

M0=8q0l02S0k3tanhkl02+q0l03S0-4q0l0S0k2-q0l06η2EIkS0tanhkl02+1EIη（18）

坡滑區域內管道受到的當量軸力為壓力時，其撓曲線方程為：

y=1S08q0l02k4-M0coskx+tankl02sinkx-

q03S0l02x4+2q03l0S0x3+4EIq0S20l02x2-q0l03S0+4EIq0l0S20x+

v0-1S08q0l02k4-M0（19）

M0=8q0l02S0k3tankl02-q0l03S0-4q0l0S0k2-q0l06η2EIkS0tankl02+1EIη

（20）

滑坡區域內管道的轉角、彎矩和剪力方程可分別通過θ=dydx、M=-EId2ydx2和Q=-EId3ydx3得到。管道前端和后端的軸向應力為：

σf=N0A+MW=N0A+MD2I

σa=N0A-MW=N0A-MD2I（21）

式中：σf為沿滑坡方向管道前端的軸向應力，Pa；σa為沿滑坡方向管道后端的軸向應力，Pa;W為管道截面抗彎模量，m3。

根據Mises屈服條件，當量應力σMises的表達式為：

σMises=σ1-σ22+σ1-σ32+σ2-σ322（22）

式中：σ1、σ2、σ3分別為第一、第二、第三主應力，Pa。

2 橫坡敷設管道力學響應規律分析

2.1 工程算例

以忠-武線沿程某一坡體作為算例，分析埋地管道在滑坡條件下的力學響應。該坡體表面較為平緩，斜坡地層構成主要為頁巖，地層傾角為20°～55°，坡體角度約為22°。滑坡體表面總面積約240 m2，滑坡體厚度為2～6 m、寬度約為40 m。管道采用規格為φ711 mm×7.9 mm的X65型鋼，線膨脹系數為1.2×10-5 m/（m·℃），泊松比0.3，運行壓力4 MPa，土體極限抗剪強度24.6 kPa，縱向綜合系數9.45 N/cm3，橫向綜合系數3.43 N/cm3，主滑面土體滑坡推力3.316 kN/m。

2.2 理論模型計算分析

基于文中推導得到的解析解，對算例管道、山體等參數進行編程計算，得到滑坡區域內橫坡敷設管道撓度、管道前后位置軸向應力與Mises應力沿軸向分布的曲線圖（見圖2）。

由圖2a可知，滑坡區域內管道變形整體呈常態分布趨勢，最大撓度為0.20 m，滑坡區域兩端邊界管段受到滑坡區域外土體嵌固作用，僅有少量變形。管道前端軸向應力沿管道呈常態分布，管道后端的分布則幾乎與之完全相反。管道前端在靠近主滑面處（滑坡區域中心）呈軸向受拉的狀態，而管道后端在靠近主滑面區域軸向受壓。在靠近滑坡區域邊界處，管道前端軸向受壓，后端軸向受拉。管道前端和后端分別在距離滑坡區域邊界兩端約l0/10和l0/4處發生軸向應力拉壓性質的轉換。從圖2b可知，管道前端在滑坡區域邊界附近受到較大的Mises應力，Mises應力沿管道軸向整體呈現靠近邊界處大，靠近主滑面附近小的分布規律，但在主滑面的管道前端存在極大值，Mises應力值達到171 MPa。管道后端Mises應力沿軸向分布呈現靠近主滑面大，其余區域整體偏小的分布，最大Mises應力出現在主滑面。

2.3 有限元數值模擬分析

基于ABAQUS軟件建立滑坡作用下管道與土體相互作用的三維有限元模型（見圖3a），對滑坡作用下管道的力學影響進行數值模擬分析。首先，采用FENG W.K.等［25］的管道橫穿滑坡區試驗結果，對數值模型進行驗證分析。試驗［21］采用φ325 mm×8 mm的L245NB鋼管，埋設于距地表1.5 m滑動尺寸為500 m3的土體中。土體參數：土壤容重為20.98 kN/m3，泊松比為0.33，內摩擦角為3.5°，彈性模量2.88 MPa，黏聚力為5.03 kPa。管-土相互作用模型采用非線性接觸模型，該模型能充分考慮管道與周圍土體的實際接觸狀態，可以模擬兩者之間的非線性摩擦現象［26］。切向作用通過設置“罰函數”來模擬管-土的摩擦與相對運動，法向作用則采用“硬接觸”模擬管-土間的接觸與分離。模型底面完全約束，其周邊設置法向約束；模型整體受到重力作用，管道受內壓作用，采用等效位移法模擬滑坡對管道的影響。管道本構模型選用Ramberg-Osgood模型［27］，其應力應變關系式為：

ε=σE+ασEσσsr（24）

式中：ε為應變；σ為應力，MPa；α和r均為R-O模型參數；σs為屈服強度MPa。

通過仿真得到管道變形分布曲線（見圖3b），可以看出，位移曲線呈現相同的類“正態分布”趨勢，最大位移誤差僅為2.2%，表明建立的有限元模型可靠。

在此基礎上，基于驗證的有限元模擬方法，結合忠-武線沿程某一坡體的基本特征，對滑坡作用下管道的力學響應進行分析。模型網格劃分、邊界條件、載荷施加方式等均與上述驗證模型一致。忠-武線沿程坡體土體材料參數：滑動面以上土體容重為19.11 kN/m2，彈性模量0.035 GPa，泊松比0.35，摩擦角14.5°，黏聚力18 kPa；滑床彈性模量12 GPa，泊松比0.35，容重25.7 kN/m2。圖4展示了模擬計算得到的管道應力、變形的情況。由圖4可以看出，管道撓度沿軸向呈現常態分布。管道軸向應力在滑坡區域內沿軸向關于主滑面呈對稱分布。管道前端在滑坡區域中間部分約l0/2長度受到軸向拉應力，在靠近滑坡區域邊界附近受軸向壓應力，主滑面上管道前端受到最大的軸向拉應力。管道后端在主滑面附近軸向受到壓應力，在滑坡邊界附近軸向受到拉應力。

管道Mises應力在前端沿軸向呈現類似“W”形分布，主滑面管道前端出現Mises應力的極大值，后端出現Mises應力最大值。滑坡區域內，管道前端Mises應力整體呈現邊界附近較大，中間區域小的趨勢；管道后端則呈現Mises應力主滑面大，邊界附近小的分布特征。

2.4 管道力學行為的對比分析

對比圖2和圖4知，管道撓度沿軸向分布趨勢一致，且最大撓度數值均為0.2 m。理論計算模型和有限元模型得到的管道軸向應力沿管道軸向分布的規律大致相同，管道前、后端均出現2個拉壓性質轉換點，在滑坡區域邊界附近和主滑面存在較大的軸向應力。管道最大Mises應力值數值相近，且均出現在主滑面的管道后端處。基于理論模型和有限元模型計算的結果，根據第四強度理論，可發現在滑坡區域邊界附近的管道前端和主滑面的管道后端為應力集中區，易發生破壞，在管道的運行維護中應著重監測其應變數值和變化速率，加強管道附近土體的抗滑治理措施。對比上述有限元模型與力學模型的案例結果，二者得到的管道力學響應規律基本一致。

3 坡面凹凸程度對管道力學響應規律的影響

3.1 坡面凹凸形狀參數

在有限元模型的計算中，為方便計算，往往將山體坡面近似為平面斜坡進行處理，但實際山體可能由于自然或人為因素，出現坡面凹陷或凸起的情況。為了探明相同滑坡區域下坡面凹凸形狀對管道力學響應規律的作用機理，首先需確立描述斜坡凹凸形狀的參數。現有的研究并沒有明確的描述參數，因此針對同一滑坡區域設定形狀參數N值（見圖5）來描述斜坡的凹凸程度。綜合考慮凹凸高度、坡面曲率，將坡面的凹凸程度N值定義為：

N=100a-bl0（25）

式中：a為坡面圓弧中心距地表最高處的垂直距離，m；b為坡面圓弧中心距地表最低處的垂直距離，m。

坡面凹陷時，N lt; 0；坡面凸起時，N gt; 0；坡面平整，則N =0。

3.2 管道力學響應規律分析

針對相同工況下不同坡面形狀特征山體及管道，采用和2.3節相同的建模思路，分別建立N值為-2.50、-1.25、1.25、2.50的三維山體-管道模型。管道敷設于斜坡中部（見圖6），對斜坡施加同樣的位移邊界條件，固定相同的區域進行仿真計算。

提取仿真結果并繪制軸向應力隨管道位置變化的曲線圖，如圖7所示。

由圖7可以看出，坡體凹凸形狀的變化并未改變管道軸向應力沿管道軸向不同位置分布的整體特征，但是隨著凹凸程度的變化，軸向應力的大小會發生明顯的改變。坡面凸起越明顯，受到的軸向應力整體越小。當N值由-2.50增加到2.50時，管道前端軸向拉應力峰值由148 MPa減小到119 MPa，后端軸向拉應力峰值由167 MPa減小到122 MPa。拉壓性質轉換點為平直虛線與各曲線交點，坡面越凸起，管道前端拉壓性質轉換點越遠離主滑面，管道后端拉壓性質轉換點更加靠近主滑面。這意味著，N值越大，坡面突起情況越嚴重，滑坡區域內管道前端受軸向拉應力的區域更大且有向滑坡區域邊界發展的趨勢，相應的受軸向壓應力的區域越小。滑坡區域內管道后端受到軸向拉應力的區域也隨著N值的變大而變大，有向主滑面處擴大的趨勢，軸向壓應力區域隨著N值的增大在主滑面附近減小。需注意的是，拉壓轉換點位置隨N值變化的規律并不一直如此，例如，N值從-2.50增大到0的過程，管道軸向受拉區域存在先增大再減小的情況。

圖8為不同斜坡N值對應的管道前、后端Mises應力沿管道軸向分布的曲線。圖9展示了不同N值下管道最大變形量和最大Mises應力。結果表明：隨著N值的變化，管道最大撓度和最大Mises應力的變化趨勢趨于一致。N值由-2.50到0的過程，管道最大變形量、峰值Mises應力最大變化幅度分別為6%和4%，相較于N≤0的部分，N值從0增加2.50時，最大變形量和峰值Mises應力變化更大，分別為36%和11%。說明凸起坡面比凹陷坡面對管道變形受力的影響大。坡面形狀的變化并不會引起管道Mises應力沿管道軸向整體分布規律的改變，在主滑面的管道后端和滑坡區域邊界的管道前端均出現Mises應力較大的區域。當凹陷坡面N越小即坡面凹陷深度越大時，橫坡敷設管道產生的最大撓度和受到的最大Mises應力表現出先減小再增大的趨勢。實際斜坡由于風化、雨水沖刷等因素，表層土出現缺失，造成坡面凹陷，由于主滑面管道上方土層厚度減少，管道后方受到的土體推力也相對減弱，管道產生的變形和受到的Mises應力也會整體減小。隨著凹陷情況加劇，管道更進一步接近地表，斜坡滑動會更直接影響到管道本體。使得管道更容易產生變形，受到的應力增大；隨著凸起坡體N值越大即坡體凸起高度越高時，橫坡敷設管道的最大撓度和最大Mises應力都越來越小，這意味著坡面發生滑動時，管道橫坡敷設在凸起坡體比在平坡中發生破壞風險的可能性更低。

3.3 管道應力計算的修正

前述對比可知，力學理論模型與有限元模型在管道變形一致的情況下，得到的管道力學響應規律一致，最大Mises應力接近。根據第四強度理論，基于Mises應力可以對管道強度進行校核，是評價管道安全情況的一個標準。為研究管道承受的最大Mises應力，基于有限元結果提取不同坡體N值對應的滑坡區域管道前、后端最大Mises應力值，擬合可得到形狀參數N值與管道前后端最大Mises應力的量化關系。

管道前端最大Mises應力與N值的量化關系為：

σMises-fmax=446 571l0N1002-200 000l0N+σfmm（26）

在坡面凹陷時，管道后端最大Mises應力與N值的量化關系為：

σMises-amax=1 000（l0N）2-90 000l0N+σamm（28）

在坡面凸起時，管道后端最大Mises應力與N值的量化關系為：

σMises-amax=-200 000l0N+σamm（29）

式中：σfmm為坡面是平坡時解析公式得到的管道前端最大Mises應力，Pa；σamm為坡面是平坡時解析公式得到的管道后端最大Mises應力，Pa；σMises-fmax為考慮形狀參數修正后的管道前端最大Mises應力，Pa；σMises-amax為考慮形狀參數修正后的管道后端最大Mises應力，Pa。

將量化關系引入橫坡敷設管道力學模型的計算，用以修正坡面不同凹凸形狀對管道最大Mises應力的影響，有利于形成基于應力的反映坡面形狀對管道安全情況影響的評價方法，進一步為管道的運維提供理論支撐。

4 結 論

（1）基于彈性地基梁理論，提出考慮管道變形非線性和土體抗力物理非線性的橫坡敷設管道力學模型，根據變形協調與當量應力的判別條件求解控制方程，推導得到管道撓度、彎矩等解析表達式。并結合工程算例，將本文得出的解析解編程，計算得到管道變形、軸向應力、當量應力沿管道軸向變化規律。

（2）基于ABAQUS軟件中建立了于坡體中部的管道和山體的三維有限元模型。仿真得到滑坡作用下管道力學響應規律與理論模型結果一致，即管道在滑坡區域主滑面和邊界附近受到較大的Mises應力，主滑面處管道變形最大，滑坡區域內管道前后端均存在拉壓性質轉化點，且管道主滑面和滑坡邊界附近存在較大的軸向應力。

（3）設立坡面形狀參數N值并分析了坡面凹凸情況對管道力學響應的影響規律，發現相同邊界條件下，隨著凹陷程度的增加，管道的最大變形和峰值Mises應力會先減小再增大，峰值Mises應力的最大變化幅度達到4%。坡面凸起情況的加劇會使得管道變形和受到的Mises應力逐漸變小，峰值Mises應力最大可降低11%。在此基礎上，根據N值對管道最大Mises應力的影響規律，修正了理論模型計算當量應力。

參考文獻［1］ 吳東容，李姝璇，王彬彬，等．油氣管道高后果區識別標準在中緬山區管道的應用［J］．油氣儲運，2023，42（12）：1376-1382．

WU D R，LI S X，WANG B B，et al.Application of identification standards for high consequence areas of oil and gas pipeline in Myanmar-China mountainous pipeline［J］.Oil amp; Gas Storage and Transportation，2023，42（12）：1376-1382.

［2］ 李華章，何文剛，李生紅．中緬油氣管道貴州段沿線滑坡變形特征及成因［J］．油氣儲運，2023，42（2）：178-187．

LI H Z，HE W G，LI S H.Characteristics and causes of landslide deformation along Guizhou section of Myanmar-China oil and gas pipelines［J］.Oil amp; Gas Storage and Transportation，2023，42（2）：178-187.

［3］ 唐俊杰．滑坡作用下管道變形特征與應變計算模型研究［D］．綿陽：西南科技大學，2022．

TANG J J.Study on deformation characteristics and strain calculation model of pipeline under landslide［D］.Mianyang：Southwest University of Science and Technology，2022.

［4］ 江珂，王東源，余志峰，等．管道受橫向滑坡作用破壞分析：以中緬管道貴州晴隆段兩次爆炸事故為例［J］．科學技術與工程，2023，23（21）：8988-8995．

JIANG K，WANG D Y，YU Z F，et al.Analysis of pipeline damage caused by lateral landslides：taking two explosion accidents in the Qinglong section of China-Myanmar pipeline in Guizhou as an example［J］.Science Technology and Engineering，2023，23（21）：8988-8995.

［5］ 唐柳怡，魏然然，王峰，等.高鋼級管道安全評價指標體系構建［J］.油氣儲運，2023，42（10）：1137-1148.

TANG L Y，WEI R R，WANG F，et al.Construction of safety evaluation index system for high-grade steel pipeline［J］.Oil amp; Gas Storage and Transportation，2023，42（10）：1137-1148.

［6］ 姜昌亮.油氣管道全生命周期質量管控與安全管理探討［J］.油氣儲運，2023，42（10）：1081-1091.

JIANG C L.Quality control and safety management of oil and gas pipelines during its whole life cycle［J］.Oil amp; Gas Storage and Transportation，2023，42（10）：1081-1091.

［7］ 王磊，鄧清祿．滑坡作用對輸氣管道危害的靜力學分析［C］∥2010年全國工程地質學術年會暨“工程地質與海西建設”學術大會．福州：中國地質學會工程地質專業委員會，2010：347-345．

WANG L，DENG Q L.Static analysis of the hazards of landslides on gas pipelines［C］∥2010 National Engineering Geology Academic Annual Conference and “Engineering Geology and Haixi Construction” Academic Conference.Fuzhou：Engineering Geology Professional Committee of the Geological Society of China，2010：347-345.

［8］ 劉嘯奔，胡匯霖，費凡，等.基于SBAS-InSAR與有限元的滑坡段管道應力計算方法［J］.油氣儲運，2024，43（2）：163-170.

LIU X B，HU H L，FEI F，et al.A computational method for assessing pipeline stress in landslide areas based on SBAS-InSAR and FEM［J］.Oil amp; Gas Storage and Transportation，2024，43（2）：163-170.

［9］ 李秋揚，趙明華，張斌，等．2020年全球油氣管道建設現狀及發展趨勢［J］．油氣儲運，2021，40（12）：1330-1337，1348．

LI Q Y，ZHAO M H，ZHANG B，et al.Current construction status and development trend of global oil and gas pipelines in 2020［J］.Oil amp; Gas Storage and Transportation，2021，40（12）：1330-1337，1348.

［10］ 黃維和，宮敬.天然氣管道與管網多能融合技術展望［J］.油氣儲運，2023，42（12）：1321-1328.

HUANG W H，GONG J.Prospect for the development of natural gas network and the multi-energy integration technology in pipeline networks［J］.Oil amp; Gas Storage and Transportation，2023，42（12）：1321-1328.

［11］ 勃洛達夫金．埋設管線［M］．北京：石油工業出版社，1980：149-162．

БOPOДABKNH П П.Buried pipeline［M］.Beijing：Petroleum industry press，1980：149-162.

［12］ RAJANI B B，ROBERTSON P K，MORGENSTERN N R.Simplified design methods for pipelines subject to transverse and longitudinal soil movements［J］.Canadian Geotechnical Journal，1995，32（2）：309-323.

［13］ O’ROURKE M J，LIU X J，FLORES-BERRONES R.Steel pipe wrinkling due to longitudinal permanent ground deformation［J］.Journal of Transportation Engineering，1995，121（5）：443-451.

［14］ 張東臣，БЫKOBЛ И．滑坡條件下埋地管道受力分析［J］．石油規劃設計，2001，12（6）：1-3，6．

ZHANG D C，БЫKOBЛ И.The force-summing analysis of buried pipeline under landslide condition［J］.Petroleum Planning amp; Engineering，2001，12（6）：1-3，6.

［15］ COCCHETTI G，DI PRISCO C，GALLI A，et al.Soil-pipeline interaction along unstable slopes：a coupled three-dimensional approach.Part 1：theoretical formulation［J］.Canadian Geotechnical Journal，2009，46（11）：1289-1304.

［16］ 梁政．滑坡地區管線應力和位移的分析［J］．天然氣工業，1991（3）：55-59，9．

LIANG Z.Analyses of pipeline stress and disolacement at landslide-prone areas［J］.Natural Gas Industry，1991（3）：55-59，9.

［17］ 鄧道明，周新海，申玉平．橫向滑坡過程中管道的內力和變形計算［J］．油氣儲運，1998，17（7）：18-22．

DENG D M，ZHOU X H，SHEN Y P.Calculation of pipeline inner force and distortion during transverse landslide body［J］.Oil amp; Gas Storage and Transportation，1998，17（7）：18-22.

［18］ 吳銳，梅永貴，鄧清祿，等．滑坡作用下輸氣管道受力分析［J］．建筑科學與工程學報，2014（3）：105-111．

WU R，MEI Y G，DENG Q L，et al.Stress analysis of buried gas pipeline under landslide load［J］.Journal of Architecture and Civil Engineering，2014（3）：105-111.

［19］ 尚玉杰，王殿龍，閆生棟，等．橫向滑坡作用下埋地管道力學響應分析［J］．安全與環境工程，2019，26（1）：155-161．

SHANG Y J，WANG D L，YAN S D，et al.Mechanical response analysis of buried pipeline under the action of transverse landslide［J］.Safety and Environmental Engineering，2019，26（1）：155-161.

［20］ 王榮有，張家銘，王聰．橫向滑坡作用下埋地管道靜力學與數值模擬對比分析［J］．科學技術與工程，2021，21（1）：326-333．

WANG R Y，ZHANG J M，WANG C.Comparative analysis of statics and numerical simulation of buried pipeline under the action of lateral landslide［J］.Science Technology and Engineering，2021，21（1）：326-333.

［21］ 謝強，王雄，張建華，等．不同滑坡形式下埋地管的縱向受力分析［J］．地下空間與工程學報，2012，8（3）：505-510．

XIE Q，WANG X，ZHANG J H，et al.Longitudinal stress analysis of pipeline buried in two different forms of landslide［J］.Chinese Journal of Underground Space and Engineering，2012，8（3）：505-510.

［22］ 張杰，梁博豐，潘波，等．牽引式滑坡作用下天然氣管道力學響應研究［J］．壓力容器，2020，37（6）：19-25，38，58．

ZHANG J，LIANG B F，PAN B，et al.Mechanical response analysis of gas pipeline under traction landslide［J］.Pressure Vessel Technology，2020，37（6）：19-25，38，58.

［23］ 劉金濤．管道橫穿滑坡相互作用大尺度模型試驗研究［D］．成都：成都理工大學，2012．

LIU J T.Study on the large scale model test of the interaction of pipeline across landslide［D］.Chengdu：Chengdu University of Technology，2012.

［24］ 龍馭球．彈性地基梁的計算［M］．北京：高等教育出版社，1981．

LONG Y Q.Calculation of elastic foundation beam［M］.Beijing：Higher Education Press，1981.

［25］ FENG W K，HUANG R Q，LIU J T，et al.Large-scale field trial to explore landslide and pipeline interaction［J］.Soils and Foundations，2015，55（6）：1466-1473.

［26］ ZHANG S Z，LI S Y，CHEN S N，et al.Stress analysis on large-diameter buried gas pipelines under catastrophic landslides［J］.Petroleum Science，2017，14（3）：579-585.

［27］ 曹平，顧曉婷，臧雪瑞，等．凍土區橫向滑坡對埋地管道應變的影響研究［J］．石油機械，2020，48（2）：141-146．

CAO P，GU X T，ZANG X R，et al.ZANG X R et al effect of lateral landslide in frozen Soil area on strain of the buried pipeline［J］.China Petroleum Machinery，2020，48（2）：141-146.

第一韓雪龍，在讀碩士研究生，生于2000年，研究方向為油氣管道安全。地址：（430100）湖北省武漢市。email：ak6123456@163.com。

通信作者：李玉坤，教授，博導。email：mliyk@qq.com。2024-05-252024-08-25劉 鋒