霍魁 曾金艷 丁學文

(山西省地震局， 太原 030006)

運城市城區(qū)土層剪切波速與埋深關系分析

霍魁 曾金艷 丁學文

(山西省地震局， 太原 030006)

根據運城市城區(qū)地震小區(qū)化建設項目中72個土層的剪切波速測試資料，采用冪函數模型對粉土、粉砂、粉黏土、細砂的土層剪切波速與埋深進行加權回歸分析，得到不同土類的模型參數及擬合優(yōu)度指標。將土層剪切波速預測值與實測值進行對比，驗證了回歸關系式的合理性和適用性。土層剪切波速與埋深呈較好的冪函數關系，擬合優(yōu)度均在0.9以上。

運城市城區(qū)； 土層埋深； 回歸分析； 剪切波速

土層剪切波速是巖土地震工程中的一個重要物理量，在場地土類劃分、地基飽和土液化判別、土層地震反應分析等方面得到了廣泛應用[1-2]。通常采用現場原位測試方法確定土層剪切波速，現場原位測試方法分為單孔檢層法、跨孔法和表面波法。其中，單孔檢層法能直接給出土層的準確位置、厚度及波速，很少受測試場地環(huán)境限制，是目前應用最廣泛的方法[3]。受土體本身固有的變異性及測試技術水平等影響，測試結果具有較大的變異性，誤差甚至可能超過30%[4-5]。本次研究中，將利用實測土層剪切波速資料建立土層剪切波速的統計關系，并對統計估計值、現場實測值和土層特征進行對比分析，綜合確定土層剪切波速值。

1 城區(qū)工程地質條件

根據土層剪切波速與埋深及其標準貫入值的相關研究[6-11]，可知：

(1) 除土層埋深外，影響土層剪切波速的因素還包括土質類別、形成條件、天然重度、物理性質等。

(2) 土層剪切波速與埋深的關系可通過線性模型、一元二次多項式模型、冪函數模型進行擬合。總體上，一元二次多項式或冪函數模型的擬合精度較高。

(3) 不同地區(qū)土層剪切波速與埋深的經驗關系具有較大差異，表現出很明顯的地域性。

(4) 以土層動力性質基本一致的工程地質單元為統計單元，可得到較好的土層剪切波速與埋深經驗關系。

運城市城區(qū)主要建設規(guī)劃用地位于具巨厚松散沉積物的斷陷區(qū)，覆蓋層厚度多為74～76 m?？傮w上來看，場地內第四紀松散層厚度較穩(wěn)定，土層厚度變化明顯受斷裂控制：層序也比較穩(wěn)定，地層比較連續(xù)，呈近水平狀，不同時代地層為整合接觸。覆蓋層厚度涉及地層自上而下有第四系全新統、上更新統、中更新統地層[12]。

全新統沉降中心主要在運城市區(qū)以南的鹽池一帶，厚度為10～30 m，中條山北麓斷裂和鹽湖北岸斷裂全新世時期活動強烈。巖性下部為深灰、灰綠色粉砂質黏土，上部為黑色粉砂質黏土，多呈淤泥狀，并出現多層石鹽晶體，湖積成因。鹽池外市區(qū)全新統厚度約0～10 m，地層為褐黃色 — 淺黃色粉土，稍密，含白色菌絲、鈣核，湖積成因；上更新統在鹽池一帶厚度約30～60 m，為一套灰綠色粉砂質黏土，具有水平層理，含有豐富的鹽類礦物，呈沉積旋回，湖積成因。鹽池外市區(qū)上更新統湖相沉積較薄，厚度約5～30 m，地層為褐黃色-淺黃色粉土、粉質黏土，稍密-中密，可塑，土質不均勻，含少量細砂、云母碎片；中更新統以河湖相地層為主，為湖相的粉砂質黏土、粉土層，淺灰、褐灰色，中密-密實，可塑，含少量細砂、云母碎片、銹斑點，偶見蝸牛殼。運城市城區(qū)土層水平范圍內變化不大，垂直范圍內結構簡單，可作為一個工程地質單元，用于確定土層剪切波速與土層埋深的統計關系。

2 分析資料收集

本次工作中，收集了運城市地震小區(qū)化建設項目中72個鉆孔，圖1為工程場地位置分布圖。從圖中可以看出，場地均勻分布于運城市區(qū)主要建設用地規(guī)劃范圍，保證了本次工作統計結果的適用性。每個鉆孔波速測試深度按《工程場地地震安全性評價》(GB 17741 — 2005)、《巖土工程勘察規(guī)范》(GB 50021 — 2001)等相關規(guī)范要求，達到假想基巖面。結合土層結構測點垂直間取2 m， 分析各數據測試點。選取72個鉆孔和5個收集鉆孔共2 160個測點的剪切波數據，各土類測試點深度間隔為10 m，其分布情況見表1。

圖1 運城市小區(qū)化工程場地位置分布示意圖

土類剪切波速測試點分布數0-10m10-20m20-30m30-40m40-50m50-60m60-70m70-80m80-85m粉土23510617515713615718916514粉砂151285934624725367粉黏土9213252374642481細砂22821163091441

3 回歸分析

根據收集到的剪切波速測試資料，分別對粉土、粉砂、粉黏土、細砂進行了回歸分析。首先繪制了各類散點圖(圖2)，從圖中可以看出呈非線性規(guī)律，參照以往研究成果，采用冪函數模型Vs=a*Hb+c(Vs為剪切波速，H為土層埋深，a、b、c為回歸系數)表示土層剪切波速與埋深的關系。

首先，以埋深分檔樣本數量的倒數作為加權平均數相對值wr；然后，為了估計剪切波速的變化，對加權系數進行標準化處理，即wn=wrmean(wr)；式中wn為向量平均值；最后，利用Matlab軟件提供的nlinfit和nlpredci函數進行加權回歸分析，分別計算出粉土、粉砂、粉黏土、細砂回歸參數，判定系數R2、誤差標準差及95%預測區(qū)間。圖2中虛線表示95%預測區(qū)間，實線為預測函數Vs=a×Hb+c，式中各參數見表2。

圖2 土層剪切波速與埋深關系擬合曲線

土類abcR2σ粉土31.40.576114.10.96320.1粉砂36.80.54794.80.92624.4粉黏土22.80.640136.90.92421.9細砂24.20.611143.40.90022.4

4 工程檢驗

選取運城市城區(qū)幾個工程場地的鉆孔，利用上述回歸關系式進行剪切波速預測，計算其相對誤差，驗證所得關系式的預測精確度。具體對比結果見表3?？梢钥闯觯翆蛹羟胁ㄋ偻茰y值相對誤差均小于7.8，15個測試點中只有4個測試點的相對誤差超過了5%。這表明，此回歸關系式可以較準確地預測該地區(qū)的土層剪切波速。

表3 土層剪切波速預測值與實際值對比

5 結 語

(1) 根據城區(qū)場地土層結構分析，認為城區(qū)土層結構相對簡單，水平變化不大，可將其作為一個工程地質單元，并確定土層剪切波速與埋深具有關聯性。

(2) 通過對回歸方法和回歸系數的顯著性檢驗，得到冪函數模型Vs=a*Hb+c，可以較好地擬合出城區(qū)粉土、粉砂、粉黏土、細砂土層剪切波速與埋深關系。

(3) 通過某工程場地剪切波速實測值對比分析，認為本次所給出的回歸公式可較準確地預測該地區(qū)土層剪切波速。

[1] 廖振鵬.地震小區(qū)劃：理論與實踐[M].北京：地震出版社，1989：5-12.

[2] 周錫元，王廣軍，蘇經宇.場地·地基·設計地震[M].北京：地震出版社，1990:57-61.

[3] 汪聞韶.土工地震減災工程中的一個重要參量：剪切波速[J].水利學報，1994，25(3)：80-84.

[4] 孫銳，劉紅帥，劉德東.土體動力參數不確定性及影響研究[G].哈爾濱： 中國地震局工程力學研究所，2009：22-30.[5] 劉紅帥，鄭桐，薄景山，等.黏性土剪切波速不確定性的統計分析[J].世界地震工程，2010，26(增刊1)：99 -103.

[6] 李凡，李大華.天津軟土剪切波速影響因素的灰色關聯分析[J].勘察科學技術，2004(3)：9-11.

[7] 程祖鋒，李萍，李燕，等.深圳地區(qū)部分巖土類型剪切波速與深度的關系分析[J].工程地質學報，1997，5(2)：163-168.

[8] 陳國興，徐建龍，袁燦勤.南京城區(qū)巖土體剪切波速與土層深度的關系[J].南京建筑工程學院學報，1998，45 (2)：32-37.

[9] 李存志，李向新，姚明波.昆明盆地剪切波速與地基特性相關分析研究[J].昆明冶金高等?？茖W校學報，2006，22(3)：1-5.

[10] 齊文浩，劉德東，蘭景巖，等.西安閻良區(qū)土層剪切波速統計分析[J].防災科技學院學報，2008，10(4)：10-12.

[11] 張忠利，袁媛 ，羅云. 北海市土層剪切波速與埋深的關系[J].華南地震，2015,9(3)：40-45.

[12] 王乃樑. 山西地塹系新生代沉積與構造地貌[M]. 北京：科學出版社, 1996：206-217.

Relationship Between Soil Shear Wave Velocity and Soil Depth in Urban Areas of Yuncheng City

HUOKuiZENGJinyanDINGXuewen

(Earthquake Administration of Shanxi Province, Taiyuan 030006, China)

Based on data of shear wave velocities measured from 72 Yuncheng districts, the paper does the weighted regression analysis for shear wave velocity and depth of soils with various of soil type, including silt, silty sand, silty clay, fine sand, by using the power function model, and gets the model parameters of different soil types and goodness of fit index. Finally, the paper verifies the rationality and applicability of the regression equations via comparing measured values with the predicted values. The results show that there is power function relation between the sheaf wave velocities and the depths of soils in Yuncheng, and the goodness of fit is over 0.9.

Yuncheng City; soil depth; regression analysis; shear wave velocity

2017-05-02

山西省地震局科研項目“運城市土層剪切波速與埋深的關系”(SBK-1633)

霍魁(1987 — )，男，碩士，助理工程師，研究方向為巖土地震、強震臺網。

TU411

A

1673-1980(2017)04-0112-04