本發(fā)明海洋工程地質(zhì)，具體涉及到一種波浪對直墻式巖質(zhì)岸坡沖擊力的評估方法。

背景技術(shù)：

1、隨著國家經(jīng)濟社會的高速發(fā)展，海洋資源的探測和建設(shè)規(guī)模在不斷擴大，同時逐漸從淺海區(qū)向深海區(qū)發(fā)展。在工程建設(shè)運營的過程中對周圍海洋環(huán)境可能發(fā)生的地質(zhì)災(zāi)害的評估必不可缺，其中由波浪沖擊導(dǎo)致海洋工程周邊巖質(zhì)邊坡失穩(wěn)的地質(zhì)災(zāi)害需要引起高度重視。

2、海洋地質(zhì)災(zāi)害的水動力觸發(fā)機制的成因之一，波浪是由風(fēng)暴等海洋氣候自然因素產(chǎn)生的，在其行進過程中會不斷對海洋岸坡施加正應(yīng)力和剪應(yīng)力，影響巖質(zhì)岸坡的穩(wěn)定性，最終導(dǎo)致岸坡內(nèi)部結(jié)構(gòu)產(chǎn)生破壞并產(chǎn)生滑坡。波浪沖擊作用對巖質(zhì)岸坡的穩(wěn)定性影響主要由兩點：一是波浪對巖質(zhì)岸坡坡面的循環(huán)沖擊作用，在該作用下，巖體內(nèi)部產(chǎn)生瞬態(tài)的孔隙水壓力以及殘余孔隙水壓力，內(nèi)部超孔隙水壓力積累使得內(nèi)部應(yīng)力發(fā)生顯著變化；二是波浪對對結(jié)構(gòu)面的弱化作用，體現(xiàn)在波浪對巖質(zhì)岸坡的低頻率、長時間的反復(fù)作用，在長期對巖體與其結(jié)構(gòu)面的循環(huán)應(yīng)力影響下，其物理、化學(xué)性能發(fā)生了改變，降低了巖石的強度與變形模量。波浪主要影響巖質(zhì)岸坡的內(nèi)部結(jié)構(gòu)，不易觀測且一旦產(chǎn)生破壞危害性極大，嚴(yán)重威脅與岸坡相關(guān)的海洋工程和工作人員的生命財產(chǎn)安全。

3、但是，目前對于在波浪影響下的岸坡穩(wěn)定性分析的研究對象主要集中于土質(zhì)岸坡上，雖然在該方面取得了一系列規(guī)范計算分析與模型試驗的研究成果，但是波浪對巖質(zhì)岸坡的沖擊力的評估內(nèi)容相對較少，也缺少對規(guī)范計算結(jié)果和數(shù)值模擬結(jié)果的比對研究。波浪對直墻式巖質(zhì)岸坡沖擊力的評估需要一種更為規(guī)范、有效的評估方法來指導(dǎo)實際的海洋工程建設(shè)工作。

技術(shù)實現(xiàn)思路

1、針對現(xiàn)有技術(shù)存在的不足，本發(fā)明提出一種波浪對直墻式巖質(zhì)岸坡沖擊力的評估方法，采用規(guī)范計算法和數(shù)值模擬法對巖質(zhì)岸坡受到的波浪作用進行計算分析，最后開展比對分析選取計算結(jié)果，以解決相關(guān)評估內(nèi)容較少且缺少比對參照的問題，從而有效、快速地確定波浪對直墻式巖質(zhì)邊坡的作用力。

2、本發(fā)明采用的技術(shù)方案是，一種波浪對直墻式巖質(zhì)岸坡沖擊力的評估方法，該方法通過構(gòu)建三維數(shù)值模型和規(guī)范計算相結(jié)合，精確評估波浪對巖質(zhì)岸坡的沖擊力，包括以下步驟：

3、步驟s1：根據(jù)直墻式巖質(zhì)岸坡的幾何形態(tài)、結(jié)構(gòu)特征及其地質(zhì)組成，構(gòu)建三維數(shù)值模型，模擬岸坡結(jié)構(gòu)的幾何形態(tài)和物理特性；

4、步驟s2：根據(jù)實際海洋環(huán)境條件和波浪特性，設(shè)定數(shù)值模型的邊界條件，包括波浪參數(shù)、岸坡材料屬性以及周圍環(huán)境因素；

5、步驟s3：選擇適合波浪與巖質(zhì)岸坡相互作用的數(shù)值計算方法，設(shè)定數(shù)值模擬程序的計算方法和參數(shù)，包括使用的數(shù)值方法、時間步長；

6、步驟s4：運行數(shù)值模擬程序，模擬波浪作用下巖質(zhì)岸坡的應(yīng)力應(yīng)變反應(yīng)，獲取孔隙水壓力、應(yīng)力分布數(shù)據(jù)；

7、步驟s5：結(jié)合規(guī)范計算結(jié)果，分析數(shù)值模擬結(jié)果，對比和驗證數(shù)值模擬結(jié)果的準(zhǔn)確性，調(diào)整模型參數(shù)以優(yōu)化結(jié)果；

8、步驟s6：根據(jù)模擬結(jié)果和實際工況，評估波浪對直墻式巖質(zhì)岸坡的沖擊力，并提供針對性建議用于工程設(shè)計和安全防護。

9、優(yōu)選的，步驟s1中，所述構(gòu)建三維數(shù)值模型的方法，包括根據(jù)直墻式巖質(zhì)岸坡的主要形態(tài)，考慮實際條件使用數(shù)值模擬軟件(例如flow-3d程序、mike?21程序及swan程序等)建立波浪傳播和直墻式巖質(zhì)岸坡的1：1還原的數(shù)值模型，同時輸入波浪的特性、直墻式巖質(zhì)岸坡的幾何形狀和材料屬性參數(shù)。

10、優(yōu)選的，步驟s2中，所述數(shù)值模型的邊界條件，包括波浪源的輸入、水深分布、直墻式巖質(zhì)岸坡的幾何形狀、底床摩擦系數(shù)的重要數(shù)值模擬設(shè)置選項。

11、優(yōu)選的，步驟s3中，所述計算方法和參數(shù)的設(shè)定為選定的規(guī)范和公式，結(jié)合波浪參數(shù)和模型的幾何形狀，計算波浪對模型的沖擊力。

12、優(yōu)選的，步驟s3中，所述計算方法選擇的公式為基于規(guī)范《海港水文規(guī)范》(jts145-2-2013)的沖擊力計算部分按以下步驟進行計算：

13、判定作用于直墻式巖質(zhì)岸坡的波浪的波態(tài)；

14、計算波峰與波谷作用下直墻式巖質(zhì)岸坡受到的沖擊作用力。

15、優(yōu)選的，步驟s3中，所述作用于直墻式巖質(zhì)岸坡的波浪的波態(tài)為根據(jù)岸坡和波浪條件判定，立波的產(chǎn)生條件應(yīng)滿足波峰線與岸坡大致平行，且岸坡長度大于一個波長，同時還需要滿足下列公式：

16、

17、

18、式中：——波浪平均周期(s)；d——岸坡前水深(m)；h——岸坡所在處進行波的波高(m)；g——重力加速度(m/s2)。

19、優(yōu)選的，所述波峰與波谷作用下直墻式巖質(zhì)岸坡受到的沖擊力，包括根據(jù)岸坡與波浪的特征，選擇對應(yīng)的公式分別計算：

20、d≥1.8h且d/l＝0.05～0.12時(l為波長)，波峰作用下的沖擊力可根據(jù)規(guī)范《海港水文規(guī)范》(jts145-2-2013)的波峰作用時沖擊波波壓力分布圖，按下列規(guī)定計算：

21、①波面高程根據(jù)下列公式進行計算：

22、

23、

24、

25、

26、式中：ηc——波面高程(m)；d——岸坡前水深(m)；bη、m——修正系數(shù)；t*——無因次周期(s)；——平均周期(s)；g——重力加速度(m/s2)；

27、②在靜水面以上hc處的坡面波壓力強度按式(7)～(9)計算：

28、

29、n＝[0.636618+4.23264(h/d)1.67,1.0]??????(9)

30、式中：hc——波壓力強度pac在靜水面以上的作用點位置(m)；pac——與hc對應(yīng)的坡面波壓力強度(kpa)；n——靜水面以上波浪壓力強度分布曲線的指數(shù)，其值取式中兩數(shù)的大值；poc——靜水面上的波壓力強度(kpa)；γ——水的重度(kn/m3)；

31、③靜水面上的波壓力強度poc及坡面上其他各特征點的波壓力強度按式(10)計算，當(dāng)pbc>poc，取pbc＝poc，其中pbc為處于d/2水深對應(yīng)的坡面波壓力強度(kpa)：

32、

33、式中：p——坡面上各特征點的波壓力強度(kpa)；ap、bp、q——波峰作用系數(shù)，按相應(yīng)參數(shù)表確定；

34、④單位長度岸坡上的水平總波浪力pc(kn/m)按式(11)計算：

35、

36、式中：pc——單位長度岸坡上的水平總波浪力(kn/m)；pdc——坡底處波壓力強度(kpa)；pbc——與d/2水深對應(yīng)的坡面波壓力強度(kpa)；

37、⑤單位長度岸坡上的水平總波浪力矩mc(kn·m/m)按式(12)計算：

38、

39、式中：mc——單位長度岸坡上的水平總波浪力矩(kn·m/m)；

40、⑥單位長度坡底面上的波浪浮托力puc(kn/m)按式(13)計算：

41、

42、式中：puc——單位長度坡底面上的波浪浮托力(kn/m)；b——邊坡底寬(m)；

43、d≥1.8h且d/l＝0.05～0.12時，波谷作用下的沖擊作用力可根據(jù)規(guī)范《海港水文規(guī)范》(jts145-2-2013)的波谷作用時沖擊波波壓力分布圖按下列規(guī)定計算：

44、①波谷波面高程ηt按式(14)計算：

45、

46、式中：ηt——波谷波面高程(m)；ap、bp、q——波谷作用系數(shù)，可根據(jù)規(guī)范《海港水文規(guī)范》(jts145-2-2013)的參數(shù)表確定；

47、②所述坡面上各特征點的波浪壓力強度p均按(15)計算，系數(shù)ap、bp、q按相應(yīng)參數(shù)表確定，當(dāng)pdt>pot時，取pdt＝pot，其中pdt為坡底處波壓力強度(kpa)；pot為與波谷波面對應(yīng)的波壓力強度(kpa)：

48、p＝ap+bp(h/d)q????????(15)

49、式中：p——坡面上各特征點的波壓力強度(kpa)；

50、③所述單位長度坡身上的水平總波浪力pt按式(16)計算：

51、

52、式中：pt——單位長度坡身上的水平總波浪力(kn/m)；pdt——坡底處波壓力強度(kpa)；pot——與波谷波面對應(yīng)的波壓力強度(kpa)；

53、④所述單位長度坡底面上方向向下的波浪力put按式(17)計算：

54、

55、式中：put——單位長度坡底面上方向向下的波浪力(kn/m)；pdt——邊坡底處波壓力強度(kpa)；b——邊坡底寬(m)；

56、d≥1.8h、0.12≤d/l<0.139且8<t*≤9時，所述波浪壓強和波面高程等各量值可按下列公式(18)和(19)計算：

57、

58、

59、式中：t*——實際波況時的無因次周期；xt*——代表波浪壓強和波面高程等各量值；——取t*＝8和實際波況的h/d值；——取t*＝9和實際波況的h/d值。

60、優(yōu)選的，步驟s4中，所述數(shù)值模擬程序的運算，包括：

61、使用數(shù)值模擬程序模擬波浪在防波堤周圍的傳播、反射、折射和干擾；

62、使用數(shù)值模擬程序計算波浪高度、波浪壓力、波浪速度參數(shù)及波浪與直墻式巖質(zhì)邊坡之間的相互作用，最后得出具體的沖擊力結(jié)果。

63、優(yōu)選的，步驟s5中，所述分析數(shù)值模擬結(jié)果，包括比較規(guī)范計算結(jié)果和數(shù)值模擬計算結(jié)果之間差異，具體通過計算不同工況下的沖擊力結(jié)果，分析兩種方法計算的沖擊力結(jié)果之間的差異，以便進一步開展分析。

64、優(yōu)選的，步驟s6中，所述評估波浪對直墻式巖質(zhì)岸坡的沖擊力，包括結(jié)合實際工況和海洋工程建設(shè)需要評估分析波浪沖擊力的計算結(jié)果，并以此作為實際工程建設(shè)的參考依據(jù)。

65、由上述技術(shù)方案可知，本發(fā)明的有益技術(shù)效果如下：

66、本方法中通過建立高精度的數(shù)值模型對具體實例進行開展數(shù)值模擬計算，計算不同工況的結(jié)果與對應(yīng)規(guī)范計算的結(jié)果進行比對，以確保評估結(jié)果的準(zhǔn)確性和合理性，能更為規(guī)范、有效地評估波浪對直墻式巖質(zhì)岸坡沖擊力，數(shù)值模擬程序?qū)С龅乃畨毫υ茍D能更直觀分析沖擊力的影響，同時為當(dāng)?shù)氐暮Ｑ蠊こ探ㄔO(shè)提供參考依據(jù)。