• 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

      基于MIKE 21的二維水沙模型研究與應(yīng)用

      2021-03-31 01:40:54張宇倩
      四川建筑 2021年1期
      關(guān)鍵詞:河床泥沙沖刷

      張宇倩

      (西南交通大學(xué)土木工程學(xué)院,四川成都 610031)

      隨著管道建設(shè)規(guī)模的擴大,越來越多的輸油、氣管道將穿越江河等流域。水下穿越管道由于具有施工周期短、造價低廉等優(yōu)勢而被廣泛應(yīng)用于實際工程中[1]。穿越管道在長期服役過程中,由于自然、人為等因素容易在沖刷過程中發(fā)生破壞,造成巨大的財產(chǎn)安全損失。河床的極限沖刷深度將直接影響管道的設(shè)計埋深,因此研究管道附近河床的極限沖刷深度將對管道的埋深設(shè)計和安全運營提供重要參考價值。國內(nèi)外許多學(xué)者對河床沖刷深度均有一定研究。Sumer[2]等就海底管線在沙床中的沖刷開展了物理模型試驗,得出了管道局部沖刷的經(jīng)驗公式。趙彥波[3]等通過縮尺物理試驗歸納出了溝床平均沖刷深度的經(jīng)驗公式。常懷民[4]等發(fā)現(xiàn)用不同經(jīng)驗公式計算出的最大沖刷深度結(jié)果有明顯差異,選用具體公式時應(yīng)綜合考慮河床水文條件等因素。夏云峰[5]建立了非交錯網(wǎng)格的二、三維水沙數(shù)學(xué)模型,為預(yù)測河流的沖刷深度及沖淤變化提供了重要依據(jù)。

      目前,在對河床沖刷深度的研究中,物理試驗和經(jīng)驗公式法應(yīng)用較多,而數(shù)學(xué)模型計算法應(yīng)用較少。因此,本文基于MIKE 21軟件建立了與物理試驗1∶1的平面二維水沙數(shù)學(xué)模型,將數(shù)學(xué)模型計算結(jié)果與試驗結(jié)果進行對比分析,并將該模型運用到實際工程中進行計算,運用經(jīng)驗公式進行驗證。為預(yù)測天然河床的沖刷深度及管道埋深設(shè)計提供了重要參考價值。

      1 物理模型建立

      進行試驗前,根據(jù)弗勞德相似準(zhǔn)則確定試驗幾何比尺λL=80,其余參數(shù)比尺見表1。試驗在長9.1 m、寬0.4 m、高0.7 m的單向循環(huán)水槽內(nèi)進行,如圖1所示,模型試驗段布置在水槽的中下游,全長2 m,高0.8 m,試驗段為中值粒徑d50=0.55mm的模型沙,試驗段前后均鋪設(shè)有0.3 m的粗粒卵石護坡。

      試驗研究了流量分別為4.05L/s、5.40L/s和6.75L/s下的三組河床沖刷深度。試驗流速v=0.25m/s大于試驗?zāi)P蜕称饎恿魉賉v0]=0.2m/s,保證沖刷可進行,三組流量對應(yīng)的水深h分別為4.1 cm、5.4 cm、6.8 cm。

      表1 水槽模型比尺關(guān)系表

      圖1 試驗段模型布置

      2 數(shù)學(xué)模型建立

      2.1 水動力模型控制方程

      MIKE21水動力模塊是解決水動力及泥沙運輸問題的基礎(chǔ)模塊,其主要控制方程如下[6]:

      連續(xù)性方程:

      (1)

      X方向動量方程:

      (2)

      Y方向動量方程:

      (3)

      2.2 泥沙運輸模型求解公式

      FM模型中泥沙平衡方程如下[7]:

      (4)

      式中:z為河底高程;Sx和Sy是平均泥沙輸運矢量因子;tmor是時間。

      泥沙輸運算法采用如下公式:

      qT=qb+qs

      (5)

      式中:qb是推移質(zhì)泥沙運輸方程,qs是懸移質(zhì)泥沙運輸方程。推移質(zhì)泥沙公式為:

      (6)

      式中:P為所有推移質(zhì)泥沙都啟動的概率;θ′表示底表面摩阻有關(guān)的無量綱剪切應(yīng)力;θc為泥沙啟動臨界剪切力;γs表示泥沙的相對密度d50表示中值粒徑。

      (7)

      啟動概率定義為:

      (8)

      式中:β為動摩阻參數(shù)。

      懸移質(zhì)泥沙公式為:

      (9)

      (10)

      式中:c為瞬時懸移質(zhì)泥沙濃度;t為時間;d50表示中值粒徑;z為瞬時水深;D為積分范圍。

      2.3 模型構(gòu)建與網(wǎng)格劃分

      根據(jù)水槽試驗?zāi)P徒?∶1的平面二維河道模型進行數(shù)值模擬,矩形河道長2 m,寬0.4 m,模型建立完畢后利用SMS軟件進行結(jié)構(gòu)化網(wǎng)格劃分,并定義上下游邊界,最后導(dǎo)入MIKE21FM模型建立模塊,生成的網(wǎng)格如圖2所示。

      圖2 數(shù)學(xué)模型網(wǎng)格(單位:m)

      2.4 模型參數(shù)設(shè)置

      模型建立完畢后,進行水動力模塊和泥沙運輸模塊的參數(shù)設(shè)置,具體參數(shù)設(shè)置見表2,數(shù)值計算時入口邊界條件取三組試驗設(shè)計流量,出口邊界條件根據(jù)試驗水深換算為相應(yīng)水位進行設(shè)置。參數(shù)設(shè)置完畢后,進行數(shù)學(xué)模型的模擬計算。

      表2 MIKE 21 FM參數(shù)設(shè)置

      3 驗證結(jié)果與分析

      取河床下游中部某點A的沖刷深度進行分析,河床斷面及取點位置見圖3。

      圖3 河床取樣位置(單位:m)

      各工況下A點沖刷深度的數(shù)學(xué)模型計算結(jié)果與試驗對比結(jié)果見圖4,由圖可知,計算結(jié)果基本一致,由此可驗證該數(shù)學(xué)模型計算結(jié)果的正確性。通過觀察沖刷過程中河床附近的演變規(guī)律,對比分析各工況下的試驗數(shù)據(jù),結(jié)合數(shù)學(xué)模型計算出的河床沖刷深度曲線,可知河床沖刷大致可分為3個階段,分別是:緩慢沖刷階段、極速沖刷階段和沖刷平衡階段。在初始緩慢沖刷階段,河床沖刷深度變化不大,沖刷速率較小,沖刷深度曲線的斜率較為平緩。各工況下緩慢沖刷階段持續(xù)時間分別為:380 min、630 min、950 min,沖刷深度分別為:5 mm、11 mm、16 mm,由此可知該階段持續(xù)時間隨著流量的增大而增大,沖刷深度也逐漸增加,這是因為在初始河面流速一定時,隨著流量的增大,河道水深加深,河床附近的流速減小,初始沖刷速率減緩,緩慢沖刷階段持續(xù)的時間延長,隨之對應(yīng)的沖刷深度逐漸加深。隨著河床沖刷深度逐漸加深,沖刷速率越來越快,進入極速沖刷階段,該階段持續(xù)時間隨著流量的增大而減小,沖刷速度和沖刷深度隨著流量的增大而增大,河床沖刷深度變化顯著。這是因為隨著河床地形的變化,沖刷深度逐漸加深,對應(yīng)的沖刷速率逐漸增大,如此往復(fù)循環(huán)影響,河床沖刷深度發(fā)生極速變化。當(dāng)沖刷深度增大到一定值時,河床沖刷速率開始逐漸減緩,進入沖刷平衡階段,此時沖刷速率逐漸減小至0,沖刷深度達到極限值,河床地形將穩(wěn)定不變。各工況下的極限平衡沖刷深度分別為:3.084 cm、3.363 cm、4.174 cm,到達沖刷平衡歷時分別為:1 830 min、2 290 min、2 970 min,由此可得在流速一定時,隨著流量的增大,河床的沖刷深度逐漸加深,到達沖刷平衡的歷時越長,河床變化時間越長。

      (a) Q=4.05L/s, v=0.25m/s

      4 實際工程應(yīng)用

      為了驗證該數(shù)學(xué)模型在實際工程中的可行性,根據(jù)實際工程資料,選取北干線中河河流進行河床沖刷數(shù)值模擬,并將計算結(jié)果與64-1修正公式計算出的計算結(jié)果進行比對分析。由于缺乏實際河道地形圖,建立中河穿越段簡化矩形河道數(shù)學(xué)模型,河道寬105 m,長800 m,按照實際工程地質(zhì)條件設(shè)置水沙模型參數(shù),計算河段河床中值粒徑d50=44.25 mm,起動流速[v0]=1.9m/s。根據(jù)規(guī)范,北干線中河穿越管道工程等級為中型,設(shè)計洪水頻率P=2%。根據(jù)水文資料得到該洪水頻率下的水文參數(shù),見表3。

      表3 中河河道斷面P=2%水文參數(shù)

      4.1 沖刷深度經(jīng)驗公式

      根據(jù)北干線中河河道的地形特點,選用64-1修正式進行沖刷深度的計算,計算公式如下:

      H=hp-hcq

      (11)

      式中:hp為沖刷后的最大水深(m);Bz為該洪水頻率下的河槽寬度(m);Hz為該洪水頻率下的河槽平均水深(m);Ad為單寬流量的集中系數(shù);Q2為通過河槽的設(shè)計流量(m3/s);μ為水流側(cè)向壓縮系數(shù),取μ=0.98;Bcj為河槽部分的過水凈寬(m);hcm為河槽最大水深(m);hcq為河槽平均水深(m);E為與汛期含沙量有關(guān)的系數(shù),取0.46;d為河槽泥沙中值粒徑(mm),H為一般沖刷深度。

      4.2 計算結(jié)果分析對比

      經(jīng)數(shù)學(xué)模型計算,中河河道在50年一遇的洪峰流量下沖刷6 d左右達到?jīng)_刷平衡,整個河道的平衡沖刷深度變化情況見圖5,結(jié)果表明,河床沖刷深度沿河床順?biāo)鞣较虺侍荻仍龃蟮淖兓?guī)律,在河床最下游處達到最深,最大沖刷深度為1.5 m。

      圖5 河床沖刷深度變化(單位:m)

      通過64-1經(jīng)驗公式計算出的中河最大沖刷深度為1.41 m,略小于數(shù)學(xué)模型計算結(jié)果,考慮到數(shù)學(xué)模型是在絕對理想條件下進行數(shù)值計算,與經(jīng)驗公式計算條件存在一定偏差,且二者誤差較小,因此認為該數(shù)學(xué)模型的計算結(jié)果可靠,可應(yīng)用于實際工程河流最大沖刷深度的預(yù)測。

      5 結(jié)論

      (1)通過建立與物理模型試驗段1∶1的MIKE21平面二維水沙簡化數(shù)學(xué)模型計算河床沖刷深度,計算結(jié)果與試驗結(jié)果吻合較好,驗證了該數(shù)學(xué)模型的正確性。

      (2)天然河床的沖刷一般可分為3個階段:緩慢沖刷階段、極速沖刷階段和沖刷平衡階段。三個階段的持續(xù)時間與沖刷深度值均隨著流量、水深的不同而有所差異,當(dāng)流速一定時,流量越大、水深越深,三個階段的持續(xù)時間越長,河床沖刷深度越深。

      (3)對于實際工程河流,該數(shù)學(xué)模型計算結(jié)果與經(jīng)驗公式計算結(jié)果誤差較小,可應(yīng)用于實際工程河流計算。

      猜你喜歡
      河床泥沙沖刷
      泥沙做的父親
      新疆多泥沙河流水庫泥沙處理措施
      土壤團聚體對泥沙沉降速度的影響
      透射槽波探測技術(shù)對煤層沖刷帶的研究與應(yīng)用
      走過冬日河床
      水庫壩區(qū)沖刷漏斗的形成機理
      泥沙滅火
      兒童繪本(2015年2期)2015-05-25 18:10:15
      在沙漠中淹死的人比渴死的多
      基于CFD的液固兩相流沖刷腐蝕預(yù)測研究
      ArcGIS在河床沖淤量分析中的應(yīng)用
      河南科技(2014年23期)2014-02-27 14:18:45
      广德县| 扶余县| 和政县| 深水埗区| 兴安盟| 安徽省| 五指山市| 汝州市| 衡阳市| 沐川县| 怀化市| 永寿县| 花莲市| 红桥区| 漾濞| 门头沟区| 晋中市| 石家庄市| 沽源县| 东城区| 昌江| 凤阳县| 呼和浩特市| 桦甸市| 安平县| 新野县| 荣成市| 鸡东县| 青冈县| 乌拉特中旗| 三河市| 绥江县| 玉田县| 拉萨市| 安塞县| 赤水市| 永城市| 景宁| 诸城市| 青海省| 松溪县|