華北地區剪切波速與深度之間統計關系的通用模型特征研究1

閆振軍 呂悅軍 黃雅虹 方 怡

華北地區剪切波速與深度之間統計關系的通用模型特征研究1

閆振軍 呂悅軍 黃雅虹 方 怡

（中國地震局地殼應力研究所，北京 100085）

剪切波速是區別土動力學和靜力學的重要參數，其影響因素包括土層埋深、顆粒形狀、顆粒比重、壓縮模量、孔隙比、含水率和密度等，其中土層埋深對剪切波速的變化影響較大。本文搜集整理了華北地區10個城市的928個鉆孔共10703個測點的剪切波速與土層埋深之間的經驗統計關系，探討華北地區剪切波速隨深度變化的特征，并從巖性條件、沉積環境等方面分析其原因。通過對比分析，給出了華北地區黏性土和砂類土剪切波速隨深度變化的最佳擬合經驗統計關系，并進行實例驗證，所得結果可為缺乏數據的區域提供一定參考。

華北地區 剪切波速 統計關系 土層埋深 通用模型

引言

剪切波速是土動力學和靜力學的重要參數，在場地類別判定、砂土液化判別、場地卓越周期計算、土動剪切模量和阻尼比計算以及在地震安全性評價、地震小區劃、土木工程抗震、巖土工程勘察、地震地質災害評價等方面均有著重要應用（汪聞韶，1994；高玉峰等，2003；中華人民共和國建設部等，2004；王強等，2014）。在實際工程中，由于受塌孔、測試費用、測試時間、測試儀器以及場地條件等因素的影響，有時無法取得完整的波速資料或缺少波速資料，因此，關于剪切波速相關性的研究受到學者們的關注，也因此產生了很多經驗公式。其中，剪切波速與深度統計關系的研究最多（程祖鋒等，1997；陳國興等，1998；高飛等，2005；陶小三等，2009；溫一波等，2012；張小平等，2012；李鐵飛等，2016；王金艷等，2017；榮棉水等，2017），應用也最廣泛。但是，現有的剪切波速隨深度變化的經驗統計關系式，大多是基于1個城市、開發區、甚至1個工程場地的測試資料，利用數理統計原理得到的。這些統計關系式的場地條件針對性強，但缺乏較大區域的普適性應用價值。

雖然地殼在橫向和縱向的結構與物性方面均存在很強的非均勻性，但是影響中國大陸內部地貌格局和環境演變的構造運動具有分區性特征，成為構造分區的基礎。因此，土質狀態、地層結構、地形、地貌等地質條件也具有地區和區域性分布規律和特點（劉玉海等，1998；張培震等，2003）。因此，人們能夠從更大范圍的宏觀上分析剪切波速隨深度變化的統計關系，并從巖性條件、沉積環境等方面分析其變化規律特征，更好地服務工程建設。

根據資料的完整程度，并考慮地形地貌、地質構造條件，選擇數據量相對較多、涉及區域相對較廣的華北地區為研究對象，搜集整理了華北地區10個城市工程場地的928個鉆孔共10703個測點剪切波速與土層埋深之間的經驗統計模型關系。

通過對這些統計模型關系進行歸納梳理，探討大華北地區各區域剪切波速隨深度變化的特征，并從巖性條件、沉積環境等方面分析出現該特征的原因。通過對比分析，給出了大華北區主要巖性（黏性土和砂類土）剪切波速隨深度變化的最佳擬合模型類型，且得出了華北地區黏性土和砂類土剪切波速隨深度變化的最佳擬合經驗統計關系的具體模型公式。經實例驗證均取得較好的結果，所得結果能夠為缺乏數據的區域提供一定參考。

1 資料搜集與整理

圖1 資料來源城市（地區）分布

表1 原始資料概況

2 統計分析

地貌及地層巖性是決定剪切波速大小的主要因素。在不同的地貌單元下，即使相同巖性的地層，因其沉積環境、沉積厚度等條件影響，剪切波速也會有較大差異（董菲蕃等，2013）。因此，分析剪切波速變化特征時主要考慮地貌及地層巖性2種因素，表2為各統計區地貌和主要地層概況。

表2 城市（地區）地貌和主要地層巖性

華北地區第四系覆蓋層較厚（大于100m），巖性主要為黏性土和砂類土，局部夾粉土和砂礫石互層（李守禮等，2005）。由于各地鉆孔深度不一，為便于分析，鉆孔深度統一取50m；模型關系僅考慮主要巖性，即黏性土和砂類土。根據統計歸并得到了9個統計模型，見表3、表4以及圖2、圖3。

2.1 黏性土

表3 黏性土剪切波速最優擬合經驗關系統計

注：SD（Standard Deviation）為標準差。

從表3可以看出，華北區10個城市（地區）共分為5種擬合模型：線性函數（衡水）、冥函數（北京、唐山、煙臺、焦作）、二次函數（天津、邯鄲）、三次函數（阜陽）和分段函數（東營、下遼河）；擬合優度方面均較高（大于0.8）。而從圖2的曲線形狀看出，雖然剪切波速的總體變化隨深度增加而增大，但不同地區模擬曲線的變化趨勢存在一定差異，除煙臺、阜陽擬合模型曲線形狀和變化趨勢差異較大外，其它8條曲線的變化趨勢以邯鄲為界分為2部分：天津、衡水、東營和北京、唐山、焦作。

（1）煙臺、阜陽擬合模型曲線的形狀接近，其變化趨勢均與其它8條曲線明顯不同，表現為隨著深度增大，剪切波速起初增大較快，但在深度達10m和20m后增加較緩慢。說明該地區黏性土隨深度增加，軟硬程度變化較小。由于煙臺地區地貌類型為沖積—海積平原，相應的沉積層類型為第四紀殘坡積層（Q4eld）和沖海積層（Q4fm），地層較軟，因此剪切波速隨深度變化較小。阜陽地區在地貌單元上屬于淮北沖積平原，地層有第四系沖積層（Q4f）和晚更新世沖湖積層（Q3fl），受古河道影響，地層較軟，故剪切波速增大也較緩。

（2）冥函數模型中，除煙臺地區的模型曲線外，其它3條（北京、唐山、焦作）擬合模型非常接近，符合隨深度增加、土層硬度逐漸增加的趨勢。這與3個地方的沉積層成因均為第四系沖洪積層（Q4fp）有很大關系。

圖2 黏性土剪切波速與深度變化關系的擬合模型曲線類型

圖3 砂類土剪切波速與深度變化關系的擬合模型曲線類型

（3）雖然衡水模型曲線為線性，但僅從曲線形狀上看與天津（二次曲線）、東營（分段式曲線）非常接近，剪切波速隨深度變化的趨勢基本符合線性增大趨勢。

由上述分析可知，沉積環境對剪切波速的變化影響較大，相同或相似的沉積環境，其模型曲線的形狀和變化趨勢均相近，如沖積—海積平原的煙臺與沖湖積平原的阜陽，另外還有同為沖洪積的北京、唐山和焦作。

2.2 砂類土

與黏性土相比，華北地區砂類土擬合模型具有以下特點：

（1）9條擬合曲線中只涉及冥函數模型和一元二次模型，其中一元二次曲線2條（天津、邯鄲），其余7條為冥函數，如圖3所示。

（2）從圖3曲線分布可以看出，2種函數的變化趨勢較一致，該區砂類土的擬合模型具有較好的規律性；邯鄲砂類土統計擬合模型曲線的規律性與其它曲線差異較大，其原因與黏性土相同。

（3）從擬合優度看，冪函數的擬合優度明顯高于一元二次曲線，且模型曲線之間的差異值基本在50—150m/s，因此在該區域內砂類土的剪切波速隨深度變化的擬合模型中，冥函數較為理想。

表4 砂類土剪切波速最優擬合模型統計

由上述分析可知，即使相同模型、相同巖性的土類，由于沉積環境相差大，其模型曲線的差異仍較大，如邯鄲（位于滑積物組成的沖洪積平原）與天津（位于沖積物組成的沖積平原）的模型曲線。沉積環境越相近，模型曲線之間的差異性越小。

3 最佳經驗關系擬合

3.1 最佳經驗關系

由上述統計對比分析可知，無論是黏性土還是砂類土，其剪切波速與深度之間的規律性較好且較一致的均為冪函數。黏性土與砂類土不同的是，黏性土冪函數模型的代表性區域較小（北京、唐山、焦作），而砂類土的代表性區域相對較大（除東營外的其它地區）。

為了找出便于工程應用的通用冥函數模型，對3條黏性土和7條砂類土的冥函數曲線進一步優化擬合，分別得到2種類型土的區域通用擬合曲線（圖4）。其中黏性土的擬合關系式為=122.60.3595，擬合優度2=0.9765；砂類土擬合關系式為=92.330.3525，擬合優度2=0.9436。

根據得到的擬合優度，認為上述擬合結果作為華北平原區沖洪積成因的黏性土和砂類土的剪切波速隨深度變化通用計算公式是可行的。

3.2 實例驗證

將得到的擬合函數應用于實際工程，并將鉆孔的擬合剪切波速與實測剪切波速進行比較，以驗證利用通用公式進行剪切波速與埋深關系回歸分析的可靠性。

黏性土實例數據選自潘國營等（2000）給出的43個鉆孔在50m范圍內的不同深度（2m、4m、6m、8m、10m、12m、28m、32m、34m、36m和38m）的實測剪切波速，測試方法為單孔檢層法；砂類土實例數據則選自何仲太等（2014）給出的82個鉆孔在100m范圍內的不同深度（4m、8m、12m、16m、20m、24m、28m、56m、60m、64m、68m、72m、76m、80m、84m、88m、92m、96m和100m）的實測剪切波速，該波速同樣采用單孔檢層法測試得到，采樣間隔為1m。

圖4 剪切波速隨深度變化冥函數統計模型及區域最優擬合曲線

黏性土實例采用焦作城區分布的主要巖土類型，包括人工填土、黏性土（粉土、粉質黏土和黏土）、砂類土和碎石土，選取黏土的實測剪切波速，分別用該統計區的擬合公式及區域通用擬合公式進行剪切波速擬合，并與實測剪切波速比較進行誤差分析，結果見表5。

砂類土實例為衡水市某工程項目鉆孔實測剪切波速，土層巖性主要為填土、黏性土（粉土、粉質黏土和黏土）、粉砂和細砂，選取粉砂進行驗算，結果見表6。

表5 黏性土剪切波速統計模型計算值、通用模型計算值與實測值對比

續表

表6 砂類土剪切波速統計模型計算值、通用模型計算值與實測值對比

由表5、表6可以看出：

（1）砂類土通用模型的剪切波速擬合誤差與統計模型差別不大，這與砂類土各區總體變化規律一致性較好有很大關系；也與實例數據與通用模型接近有關。

（2）黏性土的通用模型計算結果與統計模型差別較大，由圖4（a）也可明顯看出此差別，大部分測點的誤差率小于20%，個別點的差值較大，但由此計算出的剪切波速不會使場地的軟硬程度相差1級，故誤差在可接受范圍內。

4 討論與結論

雖然地貌條件與地層巖性是影響土體剪切波速大小的決定因素，但土的結構、狀態、土體中存在的裂隙、孔隙度及孔隙中的充填物等也對剪切波速有一定影響，因此，在同一地貌單元內，相同土層的剪切波速也存在差別。由于地層巖性生成的地質環境不同，導致其土層密度、埋深、顆粒成分含量等物理性質各不相同，因此土層剪切波速具有明顯的差異性。

通過上述統計分析可知，雖然華北地區黏性土和砂類土剪切波速隨深度變化的擬合模型類型較多，但最佳模型類型均為冪函數模型，且砂類土相比于黏性土更具規律性。沉積環境相似的區域模型曲線形狀和變化趨勢更接近（如煙臺和阜陽），沉積環境相差較大的區域其曲線形狀和變化趨勢差異均較大（如邯鄲和天津），且沉積環境對黏性土的影響相比于砂類土更大。經實例驗證，優化擬合后得到的華北區黏性土和砂類土的通用函數模型的誤差均在合理范圍之內，可作為該大區缺乏剪切波速資料時的估算公式。然而，通用模型也有一定的適用條件，黏性土的剪切波速受沉積環境影響較大，其通用模型的代表性有限，對于沉積條件接近的平原區適用性較強；砂類土的通用模型對于全新世堆積物、沖積層和沖洪積層的砂類土適用性較好。

由于統計數據量的不足，導致區域覆蓋面不夠廣，因此今后還需要更多的數據統計分析，進一步充實和完善相關的研究工作。

陳國興，徐建龍，袁燦勤，1998. 南京城區巖土體剪切波速與土層深度的關系. 南京建筑工程學院學報，（2）：32—37.

程祖鋒，李萍，李燕等，1997. 深圳地區部分巖土類型剪切波速與深度的關系分析. 工程地質學報，5（2）：163—168.

董菲蕃，陳國興，金丹丹，2013. 福建沿海3個盆地的土層剪切波速與深度的統計關系. 巖土工程學報，35（S2）：145—151.

高飛，孫小剛，2005. 上海地區場地剪切波速的特征分析. 上海地質，（2）：27—29，36.

高玉峰，劉漢龍，2003. 合肥膨脹土剪切波速的特征分析. 巖土工程學報，25（3）：371—373.

何仲太，馬保起，盧海峰，2014. 衡水市巖土體剪切波速與土層深度的關系. 科學技術與工程，14（35）：90—96.

侯頡，塔拉，2017. 利用深度回歸分析天津市區土體剪切波速. 華北地震科學，35（2）：48—53.

李守禮，李迎春，2005. 東北與華北鐵路工程地質特征與成就. 鐵道工程學報，（S1）：140—150.

李鐵飛，陳學良，高孟潭，2016. 玉溪盆地淺層土層顆粒特征與剪切波速的關系研究. 地震學報，38（6）：934—941.

劉玉海，陳志新，倪萬魁，1998. 地震工程地質學. 北京：地震出版社.

劉之雨，徐振霖，2016. 阜陽市潁東區剪切波速的特征分析. 赤峰學院學報（自然科學版），32（9）：57—59.

潘國營，許昭平，趙元松，2000. 焦作市區巖土剪切波速和場地卓越周期研究. 焦作工學院學報（自然科學版），19（1）：1—4.

齊靜靜，2011. 山東東營城區剪切波速回歸分析及應用. 勘察科學技術，（2）：23—26.

齊鑫，丁浩，2012. 下遼河平原區剪切波速與土層埋深關系分析. 世界地震工程，28（3）：151—156.

榮棉水，呂悅軍，蔣其峰等，2017. 渤海常見土類剪切波速與埋深關系分析. 震災防御技術，12（2）：288—297.

陶小三，楊偉林，高志兵等，2009. 南京河西地區巖土體剪切波速與土層深度的關系. 防災減災工程學報，29（3）：320—324.

王幫圓，李世峰，王幫團等，2013. 邯鄲市區場地土剪切波速與埋深相關性分析. 河北工程大學學報（自然科學版），30（3）：92—95.

王廣軍，蘇經宇，1986. 連云港堿廠層狀土剪切波速沿深度變化的推測. 勘察科學技術，（3）：35—38.

王金艷，侯瑩，2017. 常州市城區土層剪切波速與土層深度關系分析. 地震工程學報，39（S1）：236—240.

王強，王蘭民，吳志堅等，2014. 天水市巖土體剪切波速與埋深的變化關系. 地震工程與工程振動，34（S1）：247—252.

汪聞韶，1994. 土工地震減災工程中的一個重要參量——剪切波. 水利學報，（3）：80—84.

溫一波，趙洞明，陳堯，2012. 呈貢新區剪切波特征與分布. 云南大學學報（自然科學版），34（S2）：221—226.

張培震，鄧啟東，張國民等，2003. 中國大陸的強震活動與活動地塊. 中國科學（D輯），33（S1）：12—20.

張小平，馬順，金源等，2012. 大連地區工程場地各類巖土剪切波速的變化特征分析. 防災減災學報，28（4）：7—11.

中華人民共和國建設部，國家質量監督檢驗檢疫總局，2004. GB 50011—2010 建筑抗震設計規范（2008年版）. 北京：中國建筑工業出版社.

Research on the Characteristics of Statistical Relationship between Shear Wave Velocity and Depth in North China

Yan Zhenjun, Lü Yuejun, Huang Yahong and Fang Yi

（Institute of Crustal Dynamics, China Earthquake Administration, Beijing 100085, China）

Shear wave velocity is an important parameter to distinguish soil dynamics and statics. The influencing factors include soil depth, particle shape, particle specific gravity, compressive modulus, void ratio, water content and density, etc. Among them, the soil depth has a great influence on the change of shear wave velocity. Many scholars have given statistical relationships between soil depth and the shear wave velocity of different lithological soils such as clay, silt, sandy soil and gavel soil in different regions of the country. This paper collects and compiles the empirical statistical relationship between the shear wave velocity and soil depth of 16741 points in 753 holes in 11 cities of North China. The statistical relationship is summarized and the characteristics of shear wave velocity variation in North China are discussed from aspects of lithological conditions, sedimentary environment and so on. Through comparison and analysis, the best fitting empirical relationship between the shear wave velocity and the depth of clay soil and sand soil in North China is given and verified. The obtained results can provide certain reference for the areas lacking data, which has certain engineering practical significance.

North China region; Shear wave velocity; Statistical relationship; Soil depth; General model

10.11899/zzfy20190205

中央級公益性科研院所基本科研業務專項（ZDJ2017-28），國家自然科學基金（51708516）

2018-08-01

閆振軍，男，生于1992年。碩士研究生。主要從事三維地質建模與可視化研究。E-mail：1321140479@qq.com

呂悅軍，男，生于1966年。研究員。主要從事工程地震研究。E-mail：luyj1@263.net

閆振軍，呂悅軍，黃雅虹，方怡，2019．華北地區剪切波速與深度之間統計關系的通用模型特征研究．震災防御技術，14（2）：304—313．