綜放開(kāi)采頂煤采動(dòng)應(yīng)力場(chǎng)演化路徑

秦海初，邵 斌，程海星，霍昱名

（1.中煤華晉集團(tuán)有限責(zé)任公司，山西 運(yùn)城043300；2.中煤華晉集團(tuán)有限責(zé)任公司 王家?guī)X煤礦，山西 運(yùn)城043300；3.中煤能源研究院有限責(zé)任公司，陜西 西安710054；4.太原理工大學(xué) 原位改性采礦教育部重點(diǎn)實(shí)驗(yàn)室，山西 太原030024）

綜放開(kāi)采自20世紀(jì)80年代引入我國(guó)后，得到了較為充分的發(fā)展，已經(jīng)成為了我國(guó)煤炭開(kāi)采領(lǐng)域的重要標(biāo)志性成果[1,2]。綜放開(kāi)采是在綜合機(jī)械化采煤工藝的基礎(chǔ)上，歷經(jīng)高位放頂煤、中位放頂煤、低位放頂煤工藝多次變革，成為我國(guó)目前最為普遍的綜采主導(dǎo)技術(shù)。它具有高產(chǎn)高效、安全、降耗等諸多優(yōu)點(diǎn)。對(duì)于綜放開(kāi)采技術(shù)的基礎(chǔ)理論研究，主要集中在支架圍巖關(guān)系[3-4]、頂煤破壞機(jī)理[5-6]、煤巖運(yùn)移規(guī)律等方面[7-8]。頂煤的破壞由其所受的應(yīng)力狀態(tài)決定[4-9]，因此，綜放開(kāi)采工作面應(yīng)力場(chǎng)演化規(guī)律是關(guān)鍵基礎(chǔ)問(wèn)題之一[10-11]，國(guó)內(nèi)外學(xué)者對(duì)其展開(kāi)了大量的研究工作。陳忠輝[12]等通過(guò)損傷力學(xué)理論和實(shí)驗(yàn)驗(yàn)證相結(jié)合方法[13]，對(duì)綜放工作面前方支承壓力分布特性進(jìn)行了研究；吳建[14]等成莊礦采場(chǎng)圍巖應(yīng)力分布進(jìn)行了研究，得出煤體和上覆巖層采動(dòng)應(yīng)力的增加是同步的；劉長(zhǎng)友[15]等對(duì)孤島工作面支承壓力分布特征進(jìn)行了分析，得出超前支承壓力存在2個(gè)峰值區(qū)域，頂煤支承壓力峰值大于底煤；劉金海[16]等對(duì)深井特厚煤層綜放工作面支承壓力進(jìn)行了分析，將回采期間煤體垂直應(yīng)力的變化沿走向進(jìn)行了分區(qū)；趙鯤鵬[8]分析了綜放工作面不同傾角下采場(chǎng)圍巖的應(yīng)力狀態(tài)，與現(xiàn)場(chǎng)實(shí)測(cè)結(jié)果相吻合；王德超[17]等分析了回采過(guò)程中超前支承壓力變化規(guī)律，認(rèn)為工作面“見(jiàn)方”階段垂直應(yīng)力最大；劉欽節(jié)[18]等在朱仙礦地應(yīng)力實(shí)測(cè)的基礎(chǔ)上，對(duì)綜放工作面采動(dòng)應(yīng)力場(chǎng)演化規(guī)律進(jìn)行了研究，得出了地應(yīng)力與綜放采場(chǎng)采動(dòng)應(yīng)力場(chǎng)演化的關(guān)系；黃志增[19]等在實(shí)驗(yàn)室模擬了頂煤真實(shí)應(yīng)力環(huán)境，得出了巨厚煤層頂煤破壞形式及對(duì)不同層位頂煤的冒放性做出了評(píng)價(jià)；魏世明[20]等對(duì)支承壓力研究現(xiàn)狀進(jìn)行了綜述，得出了有關(guān)支承壓力的研究的特點(diǎn)及存在的問(wèn)題；王家臣[21]等對(duì)綜放開(kāi)采頂煤裂隙擴(kuò)展的應(yīng)力驅(qū)動(dòng)機(jī)制進(jìn)行了研究，得出綜放開(kāi)采過(guò)程中主應(yīng)力大小及方向的演化規(guī)律。大部分學(xué)者將其簡(jiǎn)化為最大主應(yīng)力為垂直應(yīng)力、中間主應(yīng)力與最小主應(yīng)力相等的情況，著重考慮了垂直應(yīng)力演化特性而忽略了原巖應(yīng)力狀態(tài)及中間主應(yīng)力。為此，研究以王家?guī)X煤礦12309綜放工作面為工程背景，在原始地應(yīng)力測(cè)試基礎(chǔ)上，針對(duì)工作面前方頂煤中主應(yīng)力場(chǎng)演化規(guī)律展開(kāi)。

1 工程背景

1.1 工作面概況

王家?guī)X煤礦位于山西省河津市，12309工作面位于其西翼，煤層平均6.1 m，為近水平煤層，結(jié)構(gòu)簡(jiǎn)單，賦存狀況穩(wěn)定，平均埋深360 m，沿東西方向由西向東進(jìn)行回采。工作面位置如圖1。

工作面直接頂為粉砂巖，平均厚度3.54 m，基本頂為細(xì)粒砂巖，平均厚度4.2 m。直接底為細(xì)粒砂巖與泥巖互層，平均厚度2.77 m，基本底為中粒砂巖，平均厚度24.29 m。采用綜采放頂煤工藝進(jìn)行回采，采高3.1 m，放煤高度3.0 m，采用ZFY12000/23/34D型兩柱放頂煤液壓支架，額定工作阻力12 000 kN，初撐力9 420 kN。

圖1 工作面位置Fig.1 Location of panel 12309

1.2 地應(yīng)力測(cè)試

運(yùn)用空心包體類三維應(yīng)力測(cè)量法對(duì)王家?guī)X煤礦原始應(yīng)力進(jìn)行現(xiàn)場(chǎng)測(cè)量，于12309鄰近工作面選取3個(gè)測(cè)站，分別測(cè)量其地應(yīng)力大小及方向。首先在測(cè)站自巷幫用鉆機(jī)鉆進(jìn)應(yīng)力解除孔（直徑130 mm），至不受除鉆孔以外的其他工程影響的位置（1.25倍巷道寬度b）；磨平大孔孔底后，再同心鉆進(jìn)小直徑孔（直徑36 mm）穿過(guò)待測(cè)點(diǎn)；在小測(cè)量孔內(nèi)安裝環(huán)氧樹(shù)脂三軸應(yīng)變計(jì)，24 h后，待環(huán)氧樹(shù)脂完全固化，用直徑130 mm的取心鉆頭對(duì)含有應(yīng)變計(jì)的巖心進(jìn)行應(yīng)力解除，解除過(guò)程引起的巖石變形可由應(yīng)變計(jì)中不同方向的應(yīng)變片通過(guò)數(shù)據(jù)采集器檢測(cè)和記錄，地應(yīng)力測(cè)試過(guò)程如圖2。

圖2 地應(yīng)力測(cè)試過(guò)程Fig.2 Testing process of ground stress

由現(xiàn)場(chǎng)地應(yīng)力測(cè)試結(jié)果可知，礦區(qū)水平應(yīng)力大于垂直應(yīng)力，最大水平應(yīng)力σx為南北方向，平均側(cè)壓系數(shù)為λx=1.52，最小水平應(yīng)力σy為東西方向，平均側(cè)壓系數(shù)為λx=0.45，地應(yīng)力測(cè)試結(jié)果見(jiàn)表1。

表1 地應(yīng)力測(cè)試結(jié)果Table 1 Measurement results of ground stress

2 數(shù)值模擬

2.1 數(shù)值模型及參數(shù)

2.1.1 數(shù)值模型及測(cè)線

按照?qǐng)D2中各巖層厚度建立數(shù)值模型，以正南北向?yàn)閤軸，正東西向?yàn)閥軸，垂直方向?yàn)閦軸，全部由六面體單元組成。因旨在研究工作面中部推進(jìn)方向主應(yīng)力演化規(guī)律，固不考慮工作面長(zhǎng)度對(duì)頂煤體中主應(yīng)力空間演化規(guī)律的影響，因此模型尺寸為長(zhǎng)×寬×高=200 m×200 m×80 m，工作面模擬尺寸為走向150 m、傾向100 m，數(shù)值模型及測(cè)線位置如圖3。

圖3 數(shù)值模型及測(cè)線位置Fig.3 Numerical model and line layout

根據(jù)λx=1.52、λy=0.45，設(shè)置模型側(cè)壓系數(shù)，進(jìn)行初始應(yīng)力反演，再對(duì)工作面一次開(kāi)挖，計(jì)算至平衡（以模型最大不平衡力比率小于10-5為平衡條件），導(dǎo)出工作面中部測(cè)線位置的3個(gè)主應(yīng)力大小及方向。測(cè)線位于工作面中部頂煤內(nèi)，測(cè)線起始坐標(biāo)為（100，100，21.5），終止坐標(biāo)為（100，200，12.5），測(cè)線中每隔1 m選取1個(gè)測(cè)點(diǎn)，共設(shè)置測(cè)點(diǎn)100個(gè)，待模型開(kāi)挖平衡后，通過(guò)FISH語(yǔ)言中內(nèi)置“z_pstress”函數(shù)，提取每個(gè)測(cè)點(diǎn)的3個(gè)主應(yīng)力大小、方向，測(cè)線位置如圖3（b）。

2.1.2 邊界條件及本構(gòu)模型

1）應(yīng)力邊界。數(shù)值模型總高80 m，煤層底板在數(shù)值模型中位于27.06 m處，即數(shù)值模型中煤層埋深為52.94 m，由于實(shí)際煤層平均埋深為360 m，因此需要在模型上邊界施加相當(dāng)于307.06 m厚度的巖層重力，按照巖層平均密度為2.2 t/m3計(jì)算，得模型上邊界應(yīng)施加應(yīng)力σz=6.76 MPa。

2）位移邊界。對(duì)數(shù)值模型前、后、左、右4個(gè)邊界面施加法向位移約束，對(duì)模型底面施加全位移約束。

2.1.3 本構(gòu)模型及參數(shù)

數(shù)值模型均為煤、巖石等材料，因此選用FLAC3D內(nèi)置Mohr-Coulomb本構(gòu)模型進(jìn)行模擬，數(shù)值計(jì)算中各巖層物理力學(xué)參數(shù)匯總見(jiàn)表2。

表2 巖層物理力學(xué)參數(shù)Table 2 Lithologic parameters in numerical model

3 結(jié)果分析

3.1 主應(yīng)力大小演化規(guī)律

工作面開(kāi)挖后，由于開(kāi)采擾動(dòng)，頂煤體內(nèi)應(yīng)力場(chǎng)重新分布，但采動(dòng)影響并不會(huì)波及到整個(gè)測(cè)線。測(cè)線30 m范圍內(nèi)的3個(gè)主應(yīng)力值大小分布曲線如圖4。

圖4 主應(yīng)力大小分布曲線Fig.4 Principle stresses distribution

圖4中（0，0）點(diǎn)位置處3個(gè)主應(yīng)力大小均為0；0～4 m范圍為模擬工作面液壓支架控頂區(qū)，該區(qū)域3個(gè)主應(yīng)力大小基本呈線性逐漸升高，最大主應(yīng)力σ1從0升高到3.20 MPa，中間主應(yīng)力σ2從0升高到2.05 MPa，最小主應(yīng)力σ3從0升高到0.69 MPa；最大主應(yīng)力σ1、中間主應(yīng)力σ2在11 m左右達(dá)到峰值，而最小主應(yīng)力σ3則在13 m左右達(dá)到峰值，其中最大主應(yīng)力σ1峰值為25.09 MPa，應(yīng)力集中系數(shù)為1.90，中間主應(yīng)力σ2峰值為18.17 MPa，應(yīng)力集中系數(shù)為2.13，最小主應(yīng)力σ3峰值為8.13 MPa，應(yīng)力集中系數(shù)為1.87。

中間主應(yīng)力集中系數(shù)最大，最大主應(yīng)力與最小主應(yīng)力集中系數(shù)相近；在16 m左右，最大主應(yīng)力與中間主應(yīng)力大小相等，即σ1=σ2，該現(xiàn)象的產(chǎn)生是由于王家?guī)X煤礦原巖應(yīng)力場(chǎng)最大水平應(yīng)力大于垂直應(yīng)力，受到超前采動(dòng)應(yīng)力的影響，垂直應(yīng)力不斷升高，逐漸大于了水平應(yīng)力，在這個(gè)過(guò)程中最大主應(yīng)力向垂直方向旋轉(zhuǎn)，中間主應(yīng)力向水平方向旋轉(zhuǎn)。

3.2 主應(yīng)力方向演化規(guī)律

3.2.1 最大主應(yīng)力方向演化規(guī)律

綜放開(kāi)采過(guò)程中，最大主應(yīng)力與3個(gè)坐標(biāo)平面的夾角演化曲線如圖5（截取測(cè)線前60 m區(qū)間）。

圖5 最大主應(yīng)力與各坐標(biāo)平面夾角Fig.5 Angels between maximum principle stress and planes

在原巖應(yīng)力狀態(tài)下，最大主應(yīng)力與y-z平面夾角為90°，與x-z、x-y 2個(gè)平面的夾角均為0°，與實(shí)際情況相符。最大主應(yīng)力與x-z、x-y、y-z 3個(gè)平面的夾角變化開(kāi)始于17 m左右，因此將曲線0～17 m范圍定義為最大主應(yīng)力偏轉(zhuǎn)區(qū)，進(jìn)一步地，以煤壁為分界線，將偏轉(zhuǎn)區(qū)劃分為壁前偏轉(zhuǎn)區(qū)和壁后偏轉(zhuǎn)區(qū)。

對(duì)于x-y平面，最大主應(yīng)力偏轉(zhuǎn)始于17 m位置處，17～15 m范圍，最大主應(yīng)力與x-y平面夾角從1.75°迅速偏轉(zhuǎn)至82.10°；15 m至最大主應(yīng)力峰值位置（11 m）區(qū)間內(nèi)，最大主應(yīng)力方向與x-y平面夾角基本保持不變；11 m（最大主應(yīng)力峰值位置）到4 m（煤壁位置）區(qū)間內(nèi)，最大主應(yīng)力方向與x-y平面夾角先減小后增大，在5 m位置降低到50.27°，5～4 m區(qū)間內(nèi)，最大主應(yīng)力方向與x-y平面夾角增加10°，達(dá)到60.04°；4～0 m區(qū)間為壁后旋轉(zhuǎn)區(qū)，最大主應(yīng)力方向與x-y平面夾角先增大至86.01°然后保持基本穩(wěn)定。

對(duì)于x-z平面，最大主應(yīng)力偏轉(zhuǎn)也開(kāi)始于17 m位置處，在17～16 m范圍內(nèi)，最大主應(yīng)力與x-z平面夾角從0.27°偏轉(zhuǎn)至7.62°；16 m至最大主應(yīng)力峰值位置（11 m）區(qū)間內(nèi)，最大主應(yīng)力方向與x-z平面夾角基本保持不變；在11 m（最大主應(yīng)力峰值位置）到4 m（煤壁位置）區(qū)間內(nèi)，最大主應(yīng)力方向與x-z平面夾角先增大后減小，在5 m位置增大到39.37°，5～4 m區(qū)間內(nèi)，最大主應(yīng)力方向與x-z平面夾角減小10°，至29.63°左右；4～0 m區(qū)間為壁后旋轉(zhuǎn)區(qū)，最大主應(yīng)力方向與x-z平面夾角先減小至3.99°，然后保持基本穩(wěn)定。

對(duì)于y-z平面，最大主應(yīng)力方向偏轉(zhuǎn)從距離17 m左右處開(kāi)始，在15 m左右，最大主應(yīng)力方向與yz平面夾角為1.51°，近似0°，即在17～15 m內(nèi)的2 m左右區(qū)間內(nèi)，最大主應(yīng)力與y-z平面旋轉(zhuǎn)了約90°；從15 m處至0 m位置，最大主應(yīng)力方向與y-z平面夾角保持0°左右，基本穩(wěn)定。

3.2.2 中間主應(yīng)力方向演化規(guī)律

綜放開(kāi)采過(guò)程中，中間主應(yīng)力與3個(gè)坐標(biāo)平面的夾角演化曲線如圖6（截取測(cè)線前60 m區(qū)間）。

在原巖應(yīng)力狀態(tài)下，中間主應(yīng)力與x-y平面夾角為90°，與x-z、y-z 2個(gè)平面的夾角均為0°，與實(shí)際情況相符。中間主應(yīng)力與x-z、x-y、y-z 3個(gè)平面的夾角變化開(kāi)始于53 m左右（以角度變化大于5°為旋轉(zhuǎn)分界線），因此將曲線0～53 m范圍定義為中間主應(yīng)力偏轉(zhuǎn)區(qū)，進(jìn)一步地，以煤壁為分界線，將偏轉(zhuǎn)區(qū)劃分為壁前偏轉(zhuǎn)區(qū)和壁后偏轉(zhuǎn)區(qū)。

圖6 中間主應(yīng)力與各坐標(biāo)平面夾角Fig.6 Angels between intermediate principle stress and planes

對(duì)于x-y平面，中間主應(yīng)力偏轉(zhuǎn)開(kāi)始于53 m位置處，在53～17 m（σ1=σ2）區(qū)間中，中間主應(yīng)力與x-y平面夾角緩慢減小，從原始90°減小到80.17°；中間主應(yīng)力與x-y平面夾角從17 m處開(kāi)始劇烈變化，至15 m處時(shí)，中間主應(yīng)力與x-y平面夾角由17 m處的80.17°降低至1.50°，接近0°；15～0 m區(qū)間內(nèi)，中間主應(yīng)力與x-y平面夾角保持0°基本穩(wěn)定。

對(duì)于x-z平面，中間主應(yīng)力偏轉(zhuǎn)開(kāi)始于53 m位置處，在53～17 m（σ1=σ2）區(qū)間中，中間主應(yīng)力與x-z平面夾角緩慢增大，從原始0°增大到8.65°；中間主應(yīng)力與x-z平面夾角從17 m處開(kāi)始減小，至15 m處時(shí)，中間主應(yīng)力與x-z平面夾角由17 m處的8.65°降低至0.20°；15～0 m區(qū)間內(nèi)，中間主應(yīng)力與x-z平面夾角保持0°基本穩(wěn)定。

對(duì)于y-z平面，中間主應(yīng)力方向偏轉(zhuǎn)從距離17 m左右處開(kāi)始，偏轉(zhuǎn)過(guò)程劇烈，在15 m左右，最大主應(yīng)力方向與y-z平面夾角為近似88.49°，即在17～15 m內(nèi)的2 m左右區(qū)間內(nèi)，中間主應(yīng)力與y-z平面旋轉(zhuǎn)了近90°；從15 m處至0 m位置，最大主應(yīng)力方向與y-z平面夾角保持0°左右，基本穩(wěn)定。

3.2.3 最小主應(yīng)力方向演化規(guī)律

綜放開(kāi)采過(guò)程中，最小主應(yīng)力與3個(gè)坐標(biāo)平面的夾角演化曲線如圖7（截取測(cè)線前60 m區(qū)間）。

在原巖應(yīng)力狀態(tài)下，最小主應(yīng)力與x-z平面夾角為90°，與x-y、y-z 2個(gè)平面的夾角均為0°，與實(shí)際情況相符。最小主應(yīng)力與x-z、x-y、y-z 3個(gè)平面的夾角變化開(kāi)始于53 m左右（以角度變化大于5°為旋轉(zhuǎn)分界線），因此將曲線0～53 m范圍定義為最小主應(yīng)力偏轉(zhuǎn)區(qū)，進(jìn)一步地，以煤壁為分界線，將偏轉(zhuǎn)區(qū)劃分為壁前偏轉(zhuǎn)區(qū)和壁后偏轉(zhuǎn)區(qū)。

圖7 最小主應(yīng)力與各坐標(biāo)平面夾角演化規(guī)律Fig.7 Angels between minimum principle stress and planes

對(duì)于x-y平面，最小主應(yīng)力偏轉(zhuǎn)開(kāi)始于53 m位置處，在53～11 m（最大主應(yīng)力峰值位置）區(qū)間中，最小主應(yīng)力與x-y平面夾角緩慢增大又緩慢減小，該區(qū)間內(nèi)，最小主應(yīng)力方向與x-y平面夾角變化浮動(dòng)較?。ㄎ闯^(guò)10°）；最小主應(yīng)力與x-y平面夾角從11 m處開(kāi)始，發(fā)生較顯著變化，在11～4 m區(qū)間（煤壁位置）內(nèi)，最小主應(yīng)力與x-y平面夾角先增大后減小，先由11 m處的7.31°增加至5 m處的39.73°，在減小到4 m處的29.96°；壁后旋轉(zhuǎn)區(qū)間內(nèi)，最小主應(yīng)力與x-y平面夾角由4 m處的29.96°減小至3 m處的3.99°，后保持相對(duì)穩(wěn)定。

對(duì)于x-y平面，最小主應(yīng)力偏轉(zhuǎn)開(kāi)始于53 m位置處，在53～11 m（最大主應(yīng)力峰值位置）區(qū)間中，最小主應(yīng)力與x-y平面夾角緩慢減小又緩慢增大，該區(qū)間內(nèi)，最小主應(yīng)力方向與x-y平面夾角變化浮動(dòng)較?。ㄎ闯^(guò)10°）；最小主應(yīng)力與x-y平面夾角從11 m處開(kāi)始，發(fā)生較顯著變化，在11～4 m區(qū)間（煤壁位置）內(nèi)，最小主應(yīng)力與x-y平面夾角先減小后增大，先由11 m處的82.69°減小至5 m處的50.27°，在增大到4 m處的60.04°；壁后旋轉(zhuǎn)區(qū)間內(nèi)，最小主應(yīng)力與x-y平面夾角由4 m處的60.04°增大至3 m處的86.01°，后保持相對(duì)穩(wěn)定。

對(duì)于y-z平面，最小主應(yīng)力方向與y-z平面夾角在整個(gè)綜放過(guò)程中均為0°，保持穩(wěn)定。

3.3 主應(yīng)力分區(qū)及演化路徑

將主應(yīng)力場(chǎng)分布曲線進(jìn)行分區(qū)，將主應(yīng)力場(chǎng)由遠(yuǎn)及近地分為4個(gè)區(qū)域：原巖應(yīng)力區(qū)、應(yīng)力升高區(qū)、應(yīng)力峰后降低區(qū)、支架控頂區(qū)。

在上述各區(qū)域內(nèi)取特征點(diǎn)，提取其應(yīng)力大小作為應(yīng)力路徑，其中，原巖應(yīng)力區(qū)特征點(diǎn)1個(gè)、應(yīng)力升高區(qū)特征點(diǎn)3個(gè)、應(yīng)力峰后降低區(qū)3個(gè)、支架控頂區(qū)2個(gè)，各分區(qū)特征點(diǎn)位置如圖8。

圖8 各分區(qū)特征點(diǎn)位置Fig.8 Location of characteristic points

4結(jié)論

1）王家?guī)X礦原始地應(yīng)力場(chǎng)中，最大水平應(yīng)力>垂直應(yīng)力>最小水平應(yīng)力，其中，最大水平應(yīng)力方向?yàn)槟媳毕颍骄鶄?cè)壓系數(shù)為1.52，最小水平應(yīng)力方向?yàn)闁|西向，平均側(cè)壓系數(shù)為0.45。

2）頂煤各主應(yīng)力在煤壁前方一定距離均存在應(yīng)力集中現(xiàn)象，其中，最大主應(yīng)力、中間主應(yīng)力峰值均出現(xiàn)在測(cè)線11 m位置處，應(yīng)力集中系數(shù)為1.90、2.13，最小主應(yīng)力應(yīng)力峰值出現(xiàn)在測(cè)線13 m位置處，應(yīng)力集中系數(shù)為1.87。

3）最大主應(yīng)力方向旋轉(zhuǎn)的特征點(diǎn)位置有3個(gè)，分別出現(xiàn)在位置、最大主應(yīng)力峰值位置、煤壁位置，中間主應(yīng)力方向旋轉(zhuǎn)的特征點(diǎn)位置有1個(gè)，出現(xiàn)在位置；最小主應(yīng)力方向旋轉(zhuǎn)的特征點(diǎn)位置有2個(gè)，依次出現(xiàn)在最大主應(yīng)力峰值位置、煤壁位置。

4）根據(jù)主應(yīng)力場(chǎng)演化特征，將頂煤分為原巖應(yīng)力區(qū)、應(yīng)力升高區(qū)、應(yīng)力峰后降低區(qū)及支架控頂區(qū)，并在4個(gè)區(qū)域內(nèi)共選取了9個(gè)特征點(diǎn)，組成了頂煤主應(yīng)力演化路徑，為后續(xù)頂煤破壞規(guī)律的研究提供應(yīng)力基礎(chǔ)。