劉開顏,付 湘,龔來紅,謝亨旺,劉博懿
(1. 武漢大學(xué) 水資源工程與調(diào)度全國重點(diǎn)實(shí)驗(yàn)室,湖北 武漢 430072; 2. 江西省灌溉試驗(yàn)中心站,江西 南昌 330201)
河流涉水工程可以滿足防洪、灌溉、發(fā)電、供水等方面效益,在地區(qū)經(jīng)濟(jì)社會(huì)中有重要地位。但是涉水工程的建設(shè)會(huì)占用河道原有的過水?dāng)嗝?,壅高水位,?duì)河道行洪能力造成不利影響,可能導(dǎo)致洪水漫溢。而洪水沖擊下對(duì)工程本身造成破壞,可能導(dǎo)致工程潰決。兩種災(zāi)害形式都將對(duì)當(dāng)?shù)氐纳钌a(chǎn)和周邊生態(tài)環(huán)境造成不可估計(jì)的損失。因此,我們不僅需要對(duì)涉水工程對(duì)河道行洪能力的影響進(jìn)行深入研究,還應(yīng)該關(guān)注洪水作用下的工程安全風(fēng)險(xiǎn)。
在涉水工程影響行洪能力的研究中,張大茹等[1]利用MIKE21 建立平面二維水流數(shù)學(xué)模型,考慮建筑物設(shè)計(jì)標(biāo)準(zhǔn)、橋梁尺寸、結(jié)構(gòu)等影響因素,分析跨河工程對(duì)河道行洪能力的影響。何貞俊等[2]采用物理模型試驗(yàn)不同水流條件、不同橋梁間距及碼頭群布置形式下,涉水工程疊加所引起的水位壅高值變化規(guī)律。郭立兵等[3]對(duì)涉水建筑物概化處理,依據(jù)圣維南方程基本原理,采用有限差分法離散求解方程組,建立了河道一維數(shù)值仿真模型,得出涉水建筑物在不同運(yùn)行工況下對(duì)河道行洪產(chǎn)生不利影響。劉麗紅等[4]構(gòu)建非結(jié)構(gòu)化水動(dòng)力學(xué)模型,將進(jìn)洪退洪過程作為整體對(duì)濛洼蓄洪區(qū)洪水演進(jìn)過程進(jìn)行實(shí)時(shí)模擬,分析橋梁工程建設(shè)對(duì)蓄洪區(qū)行洪能力的影響。WANG 等[5]采用MIKE21 建立二維數(shù)值模型,從最大水位和最大流速兩方面定量分析橋墩對(duì)中國嘉陵江行洪的影響。HAN 等[6]針對(duì)瓶頸河段河道寬度小、斷面窄深的問題,采用水平二維數(shù)值模型計(jì)算河流平面形狀和斷面,模擬不同洪水條件下洪水水力因子和施工后階段的變化橋梁建設(shè)對(duì)河流防洪的影響。TOMASZ等[7]使用程序模擬洪水水面線,進(jìn)而用淹沒范圍和淹沒深度評(píng)估橋梁對(duì)波蘭瓦爾塔河行洪的影響。PERUZZI 等[8]建立了包含復(fù)雜形狀水工結(jié)構(gòu)的三維流體動(dòng)力學(xué)模型,從水位、流量系數(shù)、流線、能量損失等方面探究水工結(jié)構(gòu)對(duì)水流的影響。
在洪水作用下的工程安全風(fēng)險(xiǎn)評(píng)估中,王仲玨等[9]基于水文、氣象、地質(zhì)、工程結(jié)構(gòu)的隨機(jī)特征,采用蒙特卡羅法及一階二次矩法分析洪水作用下南水北調(diào)中線工程渡槽滑移失效的風(fēng)險(xiǎn)。殷丹等[10]基于糙率的不確定性,采用情景模擬法評(píng)估河道防洪能力以及堤段設(shè)計(jì)洪水標(biāo)準(zhǔn)下的河道行洪風(fēng)險(xiǎn)。琚烈紅等[11]針對(duì)海堤防洪問題,采用故障樹法進(jìn)行風(fēng)險(xiǎn)識(shí)別,并運(yùn)用直接積分法和蒙特卡洛方法評(píng)估海堤洪水災(zāi)害風(fēng)險(xiǎn)。
鑒于各類涉水工程具有不同的結(jié)構(gòu)特點(diǎn)和功能特性,對(duì)于河道行洪能力的影響需要做針對(duì)性研究,并且綜合考慮河道行洪能力和工程安全,可以更加全面地評(píng)估涉水建筑物修建的合理性以及管理的必要性。本文提出針對(duì)修建于河道內(nèi)的跨水渡槽,考慮工程安全及阻水效果,評(píng)估其對(duì)河道行洪風(fēng)險(xiǎn)的影響。以江西省贛撫平原崗前渡槽為研究對(duì)象,從渡槽本身的結(jié)構(gòu)穩(wěn)定性以及河道水面高度是否超越堤防高度兩方面,探究其對(duì)清豐山溪行洪風(fēng)險(xiǎn)的影響。以清豐山溪不同設(shè)計(jì)頻率的洪水為輸入,通過概化渡槽結(jié)構(gòu),依據(jù)洪水特性分析并計(jì)算渡槽所受到的洪水水平?jīng)_擊力、洪水豎向上托力、浮力、渡槽自身重力、渡槽內(nèi)水體重力以及風(fēng)力,進(jìn)而求解維護(hù)渡槽抗滑抗傾穩(wěn)定性所允許清豐山溪行洪的流量閾值。其次,依托一維水動(dòng)力模型和基于拉丁超立方抽樣改進(jìn)的蒙特卡羅法隨機(jī)模擬法,計(jì)算在有無渡槽的兩種情況下清豐山溪行洪風(fēng)險(xiǎn),以期為崗前渡槽和清豐山溪的后續(xù)管理提供科學(xué)參考。
崗前渡槽位于江西省贛撫平原灌區(qū),是江西省最大的渡槽,槽身及槽身支撐段均為鋼筋混凝土結(jié)構(gòu),其運(yùn)行六十余年后,逐漸出現(xiàn)混凝土老化剝蝕、槽身開裂漏水等安全問題,被定性為四類渡槽。
清豐山溪位于江西省中部,發(fā)源于升華山,區(qū)域內(nèi)水系復(fù)雜,流經(jīng)崗前渡槽下接撫河故道,最終流入鄱陽湖。
崗前渡槽跨越清豐山溪排洪道,二者具體的地理位置見圖1。清豐山溪存在河道淤堵嚴(yán)重,行洪不暢的問題,所以研究崗前渡槽對(duì)清豐山溪行洪能力的影響,進(jìn)而權(quán)衡河道治理與工程安全問題具有重大意義,可以為灌區(qū)管理局的決策提供參考。
圖1 清豐山溪與崗前渡槽地理位置示意圖Fig.1 Sketch of Gangqian aqueduct
工程及地質(zhì)數(shù)據(jù)來自《江西省贛撫平原灌區(qū)“十四五”續(xù)建配套與現(xiàn)代化改造工程二期可研報(bào)告》、《江西省贛撫平原灌區(qū)崗前渡槽安全鑒定報(bào)告》和《贛撫平原初設(shè)報(bào)告》(CS3-2002.06),工程規(guī)范采用《灌溉與排水渠系建筑物設(shè)計(jì)規(guī)范》(SL482-2011)、《水利水電工程等級(jí)劃分及洪水標(biāo)準(zhǔn)(SL252-2017)》。水文數(shù)據(jù)來自清豐山溪的兩個(gè)控制站——吳石水文站(115°58'E, 28°21'N)和崗前渡槽下水位站(115°57'E,28°21'N)。
河道行洪風(fēng)險(xiǎn)與河道水面線緊密聯(lián)系,而水面線的推求過程中存在諸多不確定性因子。將主要的不確定性因子納入考慮,構(gòu)建河道行洪風(fēng)險(xiǎn)模型,通過蒙特卡羅方法隨機(jī)模擬水面線,抽取大量樣本后得到較為準(zhǔn)確的河道行洪風(fēng)險(xiǎn)。
2.1.1 河流洪峰流量的不確定性
超標(biāo)準(zhǔn)洪水是影響河道行洪風(fēng)險(xiǎn)的主要因素。對(duì)于某種工況,不同標(biāo)準(zhǔn)的洪水對(duì)河道行洪風(fēng)險(xiǎn)有不同的影響。根據(jù)《報(bào)告》建議,采用P-Ⅲ分布模擬清豐山溪的洪峰流量。
P-Ⅲ分布表達(dá)式為:
式中:Q為洪峰流量;m3/s,α、β、a0為皮爾遜Ⅲ型分布的形狀、尺度和位置參數(shù);Γ(α)為α的伽馬函數(shù)。
3個(gè)參數(shù)α、β、a0的計(jì)算方法如下:
式中:EX為總體均值;Cs為偏態(tài)系數(shù);Cv為變差系數(shù)。
2.1.2 河道初值水位的不確定性
根據(jù)上述河流洪峰流量的不確定性,運(yùn)用清豐山溪的水位—流量關(guān)系曲線,可得到河流初值水位的不確定性。
2.1.3 河道糙率的不確定性
河道糙率作為反映河流阻力的綜合性系數(shù),是河道水力計(jì)算中重要而敏感的參數(shù)[12]。若缺乏河道糙率具體數(shù)值,可以根據(jù)天然河道糙率表,假設(shè)其服從正態(tài)分布[13]。
2.1.4 河道兩岸堤防高程的不確定性
堤防高程由于設(shè)計(jì)誤差及人類活動(dòng)存在高低不平的情況,造成其不確定性。一般可以假設(shè)其符合正態(tài)分布,根據(jù)實(shí)測情況,計(jì)算均值和方差。
根據(jù)《灌溉與排水渠系建筑物設(shè)計(jì)規(guī)范》(SL482-2011),渡槽抗滑穩(wěn)定性安全系數(shù)K1需滿足下式要求:
式中:N表示作用于槽身的豎向力,N;P表示作用于槽身的水平力,N;fb表示支座的摩擦系數(shù)。
渡槽抗傾覆穩(wěn)定性安全系數(shù)公式K2應(yīng)滿足:
式中:K2表示抗傾覆穩(wěn)定安全系數(shù);∑MV表示所有垂直力繞背水面支點(diǎn)轉(zhuǎn)動(dòng)的傾覆力矩總和,kN·m;∑MP表示所有水平力繞背水面支點(diǎn)轉(zhuǎn)動(dòng)的力矩總和,kN·m。
渡槽受力圖如圖2所示。
圖2 渡槽受力示意圖Fig.2 Sketch of the force on aqueduct
其中橫向風(fēng)壓力、渡槽自重、渡槽內(nèi)水重、靜水壓力、靜水浮托力、洪水橫向沖擊力、洪水豎向上托的計(jì)算公式如下:
(1)作用于槽身的橫向風(fēng)壓力。
式中:WZ表示風(fēng)荷載標(biāo)準(zhǔn)值,KN/m2;W0表示基本風(fēng)壓,KN/m2;μS表示風(fēng)載體型系數(shù);μZ表示風(fēng)壓高度變化系數(shù);βZ表示風(fēng)振系數(shù);S表示迎風(fēng)面面積,m2。
(2)渡槽自重。
式中:m表示渡槽自身重量,kg。
(3)渡槽內(nèi)水重。
式中:Gw表示水重,N;B表示渡槽寬度,m;b表示渡槽槽身厚度,m;h表示渡槽內(nèi)水深,m。
(4)靜水壓力。
式中:Fs表示靜水壓力,N;;H表示清風(fēng)山溪淹沒崗前渡槽水深,m;L表示渡槽長度,m。
(5)靜水浮托力。
式中:Ff表示靜水浮托力,N;VP表示渡槽淹沒部分的體積,m3。
(6)洪水橫向沖擊力。MALAVASI 和GUADAGNIN[14]根據(jù)流水具有的動(dòng)能平衡原理提出水流對(duì)橋梁的水平作用力計(jì)算公式,以此為參考計(jì)算渡槽承受的洪水沖擊力:
式中:δ為渡槽阻水系數(shù);A為渡槽在垂直于水流方向的投影面積,m2;v為斷面平均流速,m/s。
(7)洪水豎向托力。洪水不是穩(wěn)定的層流,存在一定的波動(dòng)。洪水中任一點(diǎn)的瞬時(shí)流速均可視為由時(shí)均流速-v和脈動(dòng)流速u'組成,時(shí)均流速為穩(wěn)定的層流,其方向平行于渡槽底面,不對(duì)渡槽產(chǎn)生豎向上托力。而由于波動(dòng)引起的脈動(dòng)流速u'方向是不確定的,將其分解為水平分量ux'和豎向分量uy',渡槽底部受到由uy'產(chǎn)生的上托力[15]。
研究表明,梁體底面豎向流速水壓的分布呈直線衰減,在迎水面底部水壓為ρgΔh,在背水面底部水壓為0,具體見圖2。故渡槽底部中間的豎向上托力為:
式中:Δh為渡槽迎水面的壅水高度,m。
D’ Aubuisson 公式是目前工程中常用的壅水高度經(jīng)驗(yàn)計(jì)算公式[16]:
式中:μ為測收縮系數(shù);hb表示渡槽后水深,m;Δβ表示渡槽排架總寬度,m。
本文研究洪水使渡槽失穩(wěn)以及洪水超過堤防高程的河道行洪風(fēng)險(xiǎn),風(fēng)險(xiǎn)率可用下式表示:
式中:Z(x)表示行洪水位高程,m;M(x)表示堤防高程,m;K1由式(3)求得;K2由式(4)求得。
Z(x)由圣維南方程組求解,不考慮恒定性的影響,忽略時(shí)間因素的作用,進(jìn)而得到如下形式的恒定非均勻流水面線基本方程[17]。
式中:α為流速水頭系數(shù);ξ為局部損失系數(shù);ω為流量模數(shù);n為糙率;R為水力半徑。
河道行洪風(fēng)險(xiǎn)采用常規(guī)蒙特卡羅方法與基于拉丁超立方體的蒙特卡羅方法進(jìn)行計(jì)算。常規(guī)蒙特卡羅方法通過計(jì)算機(jī)模擬產(chǎn)生處于(0,1)區(qū)間的隨機(jī)數(shù),在已知概率分布的情況下,通過反變換法得到各變量的隨機(jī)數(shù),代入到相應(yīng)的風(fēng)險(xiǎn)率計(jì)算函數(shù)中,從而求得風(fēng)險(xiǎn)率的無偏估計(jì)量。該方法不受極限狀態(tài)函數(shù)復(fù)雜程度的影響,也不需要將狀態(tài)函數(shù)線性化和隨機(jī)變量等效正態(tài)化,當(dāng)抽樣數(shù)量足夠多時(shí),可以達(dá)到非常高的計(jì)算精度[18],是預(yù)測和估算可靠度常用的方法之一,但也存在計(jì)算次數(shù)多耗時(shí)長的缺點(diǎn)[19]。
基于拉丁超立方體的蒙特卡羅方法對(duì)常規(guī)蒙特卡羅法的隨機(jī)抽樣技術(shù)進(jìn)行改,具有更高的運(yùn)算效率和計(jì)算精度。拉丁超立方抽樣是一種多維分層抽樣法,由MCKAY 等人[20]于1979年提出,其本質(zhì)是通過產(chǎn)生分布更加均勻的樣本,在抽取更少樣本的同時(shí)獲得更高的精度[21],該方法下每個(gè)子區(qū)間僅產(chǎn)生一個(gè)隨機(jī)數(shù),避免重復(fù)的隨機(jī)數(shù)降低計(jì)算效率。公式如下:
式中:i= 1,2,…,N;U為[0,1]區(qū)間內(nèi)均勻分布的隨機(jī)數(shù);Ui是第i個(gè)子區(qū)間中的隨機(jī)數(shù)。
根據(jù)吳石站1961-2018 年共58 年流量資料系列,采用公式(1)~(2)計(jì)算清豐山溪洪峰流量的P-Ⅲ型分布,統(tǒng)計(jì)參數(shù)為:EX=490 m3/s,CV=0.63,CS=3.0CV,累積密度函數(shù)(CDF)見圖3。
圖3 清豐山溪洪峰流量概率分布圖Fig.3 PDF of peak discharge of Qingfeng Stream
圖3表明清豐山溪洪峰流量理論分布概率值和經(jīng)驗(yàn)分布概率值相近,進(jìn)一步通過SPSS軟件分析得到二者R2為0.994,可知P-Ⅲ分布對(duì)清豐山溪洪峰流量的擬合較好。
根據(jù)《報(bào)告》中“現(xiàn)狀崗前渡槽為四類渡槽,當(dāng)其運(yùn)行水深不超過2.0 m時(shí),渡槽結(jié)構(gòu)滿足安全運(yùn)行要求”,設(shè)定渡槽的4種運(yùn)行工況,即控制渡槽內(nèi)水位h分別為0.5、1.0、1.5、2.0 m。運(yùn)用Matlab R2020b 編程得到偽隨機(jī)數(shù)后采用反變換法得到符合P-Ⅲ分布的清豐山溪洪峰流量,根據(jù)公式(3)~(12)計(jì)算渡槽的抗滑穩(wěn)定系數(shù)K1和抗傾穩(wěn)定系數(shù)K2,計(jì)算結(jié)果見圖4。
圖4 不同槽內(nèi)水位下渡槽抗滑抗傾穩(wěn)定系數(shù)Fig.4 Coefficient of anti-sliding and anti-overturn stability of aqueduct under different h
由圖4(a)可知,當(dāng)h分別為0.5、1.0、1.5、2.0 m 時(shí),符合渡槽的抗滑穩(wěn)定性系數(shù)規(guī)范值[K1]=1.05 的流量閾值依次為759.9、865.6、972.1以及1 079.4 m3/s。
由圖4(b)可知,當(dāng)h分別為0.5、1.0、1.5、2.0 m 時(shí),符合渡槽的抗傾穩(wěn)定性系規(guī)范值[K2]=1.1 的流量閾值依次為868.5、1 024.3、1 192.8以及1 361.7 m3/s。
相同清豐山溪洪水流量的情況下,隨著渡槽內(nèi)水位升高,抗滑穩(wěn)定系數(shù)和抗傾穩(wěn)定系數(shù)均變大,渡槽具有更好的穩(wěn)定性,這說明在洪水來臨時(shí),應(yīng)適當(dāng)抬高槽內(nèi)水位保證渡槽結(jié)構(gòu)的穩(wěn)定性。
綜合抗滑抗傾穩(wěn)定性系數(shù)對(duì)應(yīng)的流量閾值,若上游發(fā)生了20 年一遇的洪水,洪峰流量達(dá)到1 200 m3/s,在渡槽內(nèi)水位分別為0.5,1.0,1.5,2.0 m 的4 種工況下,清豐山溪(崗前渡槽段)的行洪能力分別下降36.68%、27.87%、19.00%、10.05%;
根據(jù)清豐山溪兩岸堤防高程實(shí)測值得到均值和方差,將其看作服從N(26.1,12)。河床地質(zhì)成分主要有粉質(zhì)黏土、細(xì)砂及淤泥質(zhì)粉質(zhì)黏土,因?yàn)槿狈暮拥啦诼实膶?shí)測值,故依據(jù)《天然河道糙率表》,以及附近水域相關(guān)文獻(xiàn)所使用的糙率值[22,23],假定其服從N(0.022,0.0012)。運(yùn)用Matlab R2020b 編程得到堤防高程和糙率的隨機(jī)模擬值。
根據(jù)公式(13)~(16),基于水面線計(jì)算結(jié)果,采用常規(guī)MC法和LHS-MC 法,計(jì)算在沒有渡槽和有渡槽且槽內(nèi)水位控制在0.5 m 兩種情況下,不同模擬次數(shù)所對(duì)應(yīng)的河道行洪風(fēng)險(xiǎn),具體見圖5。
圖5 常規(guī)MC法和LHS-MC法下不同模擬次數(shù)對(duì)應(yīng)的河道行洪風(fēng)險(xiǎn)Fig.5 Flood risk corresponding to different simulation times under MC and LHS-MC
結(jié)合圖5和模擬數(shù)值,MC法在50 000次左右時(shí)可以達(dá)到誤差小于0.001,LHS-MC 法在38 000 次左右可以達(dá)到誤差小于0.001,因此,LHS-MC 法比常規(guī)MC 法所需的模擬次數(shù)更少,風(fēng)險(xiǎn)收斂更快。
無渡槽時(shí),河道行洪風(fēng)險(xiǎn)收斂后的穩(wěn)定值約為0.007;有渡槽時(shí),河道行洪風(fēng)險(xiǎn)收斂后的穩(wěn)定值約為0.34。說明渡槽對(duì)河道行洪能力有負(fù)面影響,致使河道防洪標(biāo)準(zhǔn)低于設(shè)計(jì)值。
運(yùn)用相同方法計(jì)算有渡槽情況下,渡槽內(nèi)水位分別為1.0,1.5及2.0 m 時(shí)的河道行洪風(fēng)險(xiǎn),結(jié)果依次為0.045、0.016、0.012,這說明提高渡槽內(nèi)水位可以大幅降低河道行洪風(fēng)險(xiǎn)。
控制渡槽內(nèi)水位為2.0 m,計(jì)算不同清豐山溪流量下的河道行洪風(fēng)險(xiǎn),具體見圖6。
圖6 當(dāng)h=2.0 m時(shí),清豐山溪流量和行洪風(fēng)險(xiǎn)關(guān)系Fig.6 h=2.0 m, relationship between discharge and flood risk of Qingfeng Stream
由圖6 可知,盡管控制渡槽內(nèi)水位為其所能承擔(dān)的最高水位2.0 m,當(dāng)清豐山溪流量大于1 600 m3/s時(shí),河道行洪風(fēng)險(xiǎn)仍會(huì)顯著上升,此時(shí)可以考慮報(bào)廢渡槽以保證清豐山溪的行洪能力,盡量避免河道兩旁的農(nóng)田和房屋承受洪澇災(zāi)害損失。
綜合研究所得結(jié)果,渡槽本身結(jié)構(gòu)穩(wěn)定性的要求會(huì)降低清豐山溪的可行洪流量,適當(dāng)抬升渡槽內(nèi)水位可以減少這種不利影響。此外,渡槽使得河道過水?dāng)嗝鏈p小,水面線升高,造成河道行洪風(fēng)險(xiǎn)上升。所以當(dāng)清豐山溪上游發(fā)生超設(shè)計(jì)流量1 200 m3/s 的洪水時(shí),應(yīng)盡可能抬高渡槽內(nèi)水位,降低河道行洪風(fēng)險(xiǎn);當(dāng)洪水流量超過1 600 m3/s時(shí)應(yīng)考慮報(bào)廢渡槽。
本文提出了輸水渡槽對(duì)河道行洪能力的影響分析方法,以江西省贛撫平原灌區(qū)崗前渡槽和清豐山溪為例,考慮洪峰流量、初值水位、河道糙率以及堤防高程的不確定性,計(jì)算河道行洪風(fēng)險(xiǎn),結(jié)果如下。
(1)基于洪水力學(xué)特性,從維護(hù)渡槽抗滑抗傾穩(wěn)定性角度出發(fā),反推清豐山溪行洪流量閾值。當(dāng)渡槽內(nèi)水位分別為0.5、1.0、1.5、2.0 m 時(shí),行洪流量閾值依次為759.9、865.6、972.1 和1 079.4 m3/s。相比于清豐山溪設(shè)計(jì)流量1 200 m3/s,河道行洪能力分別下降36.68%、27.87%、19.00%、10.05%。
(2)沒有渡槽的情況下,河道行洪風(fēng)險(xiǎn)約0.007;有渡槽的情況下,考慮渡槽壅水高度影響,在(1)所述4種工況下,河道行洪風(fēng)險(xiǎn)依次為0.34、0.045、0.016、0.012,說明渡槽對(duì)河道行洪能力有負(fù)面影響,致使河道防洪標(biāo)準(zhǔn)低于設(shè)計(jì)值,但提高渡槽內(nèi)水位可以有效削減負(fù)面影響。
(3)當(dāng)清豐山溪流量超過1 200 m3/s 時(shí),應(yīng)適當(dāng)提升渡槽內(nèi)水位降低行洪風(fēng)險(xiǎn);當(dāng)清豐山溪流量超過1 600 m3/s 時(shí),即使渡槽內(nèi)水位處于安全范圍內(nèi)可承受的最高值,行洪風(fēng)險(xiǎn)仍隨著河道流量增大顯著增加,此時(shí)應(yīng)考慮報(bào)廢渡槽。
(4)基于拉丁超立方體的蒙特卡羅隨機(jī)模擬法比常規(guī)蒙特卡羅隨機(jī)模擬法具有更好的收斂性,但兩種方法達(dá)到同一精度的模擬次數(shù)差別不大,對(duì)于渡槽對(duì)河道行洪風(fēng)險(xiǎn)影響評(píng)估,兩種隨機(jī)模擬方法均適用。