模拟深海环境砂土地层的材料配比试验研究

沈翔; 袁大军; 曹宇陶; 高振峰

doi:10.3969/j.issn.0258-2724.20180285

模拟深海环境砂土地层的材料配比试验研究

doi: 10.3969/j.issn.0258-2724.20180285

沈翔^1,2,3,,
袁大军^1,2,3, ,,
曹宇陶^1,2,3,
高振峰^1,2,3

基金项目: 中央高校基本科研业务费专项资金资助（2017YJS151）；国家重点基础研究发展计划（973计划）（2015CB057802）

详细信息

作者简介:
沈翔（1991—），男，博士研究生，研究方向为隧道与地下工程，E-mail：xshenbjtu@126.com

通讯作者:
袁大军（1961—），男，教授，博士，研究方向为隧道与地下工程，E-mail：djyuan603@163.com

中图分类号: V221.3
计量
- 文章访问数: 463
- HTML全文浏览量: 338
- PDF下载量: 15
- 被引次数: 0
出版历程
- 收稿日期: 2018-09-05
- 修回日期: 2018-10-24
- 网络出版日期: 2018-10-31
- 刊出日期: 2020-06-01

Experiments on Material Proportions for Simulating Sandy Layer in Deep Sea

SHEN Xiang^{1,2,3
,},
YUAN Dajun^{1,2,3
, ,},
CAO Yutao^1,2,3,
GAO Zhenfeng^1,2,3

摘要

摘要: 为了研究可近似模拟琼州海峡深海环境下砂土地层的相似材料，基于正交试验法，以重晶石粉/标准砂、黏结剂浓度、凡硅比和石膏含量作为4个控制因素，对不同配比相似材料的物理力学参数的变化规律进行了分析. 首先确定标准砂、重晶石粉作为骨料，凡士林和硅油作为黏结剂，石膏粉作为调节剂，其次基于正交试验原理设计25组相似材料配比，分别进行密度试验、直剪试验、压缩试验以及渗流试验，得到了重度、内摩擦角、黏聚力、压缩模量以及渗透系数等物理力学指标，通过极差和显著性分析，研究了各个控制因素对材料参数的影响规律，并对25组试验结果进行了去异常点的多元线性回归分析. 研究结果表明：相似材料重度受重晶石粉/标准砂比值的影响最大，石膏的比例是控制压缩模量的主要因素，内摩擦角和黏聚力无主要控制因素，重晶石粉的比重主要控制着渗透系数的大小；适用于深海环境砂土地层条件下的相似材料配比为重晶石粉/标准砂0.4，黏结剂浓度4.5%，凡硅比3∶1，石膏含量3%.
- 材料配比 /
- 正交试验设计 /
- 盾构隧道 /
- 力学参数
Abstract: In order to obtain the material proportions that can simulate the sandy layer in the deep sea of Qiongzhou Strait, using on the orthogonal experiment method, the physical and mechanical parameters of different material proportions were analyzed with four control factors of barite powder/standard sand, binder concentration, silicon ratio and gypsum content. Firstly, standard sand or barite powder is used as aggregate, petroleum jelly and silicone oil as binders and gypsum powder as modifier. Secondly, 25 groups of material proportions were designed on the basis of orthogonal experiment principle, and types of tests of density, direct shear, compression and seepage were performed to obtain physical and mechanical indices such as weight density, internal friction angle, cohesion force, compression modulus, and permeability coefficient. Through the range analysis and significance analysis, the influence of control factors on the material parameters was explored, and a multivariate linear regression analysis of test results was carried to eliminate outliers. The results show that the weight density is most affected by the ratio of barite powder to standard sand; gypsum content is the leading factor to control the compression modulus; the internal friction angle and cohesion force have no major controlling factors; and the proportion of barite powder mainly affects the permeability coefficient. The appropriate material proportions for the sandy layer in deep sea environment are barite powder /standard sand of 0.4), binder concentration of 4.5%, silicon ratio of 3∶1, and gypsum content of 3%.
- material proportions /
- orthogonal experiment design /
- shield tunnel /
- mechanical parameter

HTML全文

场地效应被广泛定义为特定的近地表地质特征对地震波的影响，可以显著改变地震运动的峰值、持续时间和频谱幅度^[1]. 由于软土沉积物的各向异性或不均匀性、局部地形、沉积盆地深度和近地表衰减等多维效应的影响，量化场地效应具有挑战性.

地震危险性评估旨在预测特定区域或场地的地面运动强度，其核心手段是采用经验记录数据后进行回归分析以构建地震动预测方程（GMPEs）. 在GMPEs中，场地效应通常由场地卓越频率和地下30 m的平均剪切波速（V_s30）表征. V_s30参数可以反映一些场地特性，但不足以描述考虑场地下方土壤特性引起地面运动频率相关的场地放大^[2]. 随着概率地震危险性分析的不断深入，对场地效应进行更为精细的定量评估以改进地震危险性分析正逐渐成为研究焦点^[3-4].

近年来，基于地表观测记录与土层底部输入运动之比为场地放大的认识，提出了用不同的经验方法来获取频率相关的场地放大系数，如井上井下谱比法、传统谱比法^[5]、水平垂直谱比法（HVSR）^[6]等. 由于钻孔技术在世界各地并未完全得到发展，井上井下谱比法被限制. 传统谱比法假设地表基岩场地的记录相当于土层底部输入值，但在土层场地紧靠地表基岩场地时才有效^[5]，地震动广义反演技术（GIT）的提出克服了这一局限性^[7]. 广义反演技术又称线性反演法，通过对强震记录的联合反演，可以高效地分离强震记录中震源、路径、场地的3种贡献^[8-9]. HVSR假设地下基岩的各向是同性的^[6]，使用单台记录H/V即可得到场地效应^[10]. GIT需要大量的强震数据作为支持，而HVSR的适用性至今备受争议^[11].

该研究选取由 17 个强震动台站获取的21个汶川地震余震中的 233 组强震动记录，采用广义线性反演法、HVSR法对比研究该区域内的局部场地效应，讨论了地震波传播的方向性对HVSR的影响，同时使用线性反演法对HVSR进行经验修正. 本文结果可对该区域的地震危险性分析提供重要的参考.

1. 强震数据

1.1 数据选取

基于汶川余震强震动波形记录展开研究，HVSR方法只需要单个台站强震波形即可获取场地效应值，因此，采用线性反演法的数据选取原则：

1）选取地表岩石台站作为参考台站，同时不考虑未被参考台站记录到的地震事件.

2）选取10 km≤震源距≤150 km，考虑对区域内的衰减路径的影响^[12].

3）选择2 cm/s²≤PGA<150 cm/s²的记录，PGA为地面运动的峰值加速度. 考虑到一组强震记录包括3个不同方向的PGA值，采取各组中PGA的最大值作为评估标准，目的是一定程度地消除背景噪声与场地土层非线性效应的影响.

4）保证选取的每个台站记录到5次以上的事件，且每个事件被5个及以上台站记录到.

基于以上原则，该研究共选取了17个台站在21次地震中记录到的233组三分量强震记录，台站51MXT是线性反演的参考站，该台站布设于稳定的露头岩石^[13]. 强震动波形记录的PGA与震源距的关系如图1所示. 可以发现：大部分记录都小于50 cm/s²，大于100 cm/s²的记录只有7条，尽可能在保留地震动数据的同时，排除非线性反应的影响. 震中与台站的方位地形如图2.

图 1 PGA与震源距的关系

Figure 1. Relationships between PGA and hypocentral distance

下载: 全尺寸图片幻灯片

图 2 震中与台站的位置

Figure 2. Location of epicenter and stations

下载: 全尺寸图片幻灯片

1.2 数据处理

1）强震动观测记录的三分量方向（EW、NS、UD方向）与SH波、SV波及P波的实际振动方向并不完全一致. 以接收台站为参考点，采用反向方位角方法，将观测记录转换为与SH波和SV波方向对应的波形记录.

2）利用STA/LAT-AIC方法从每条强震动波形中识别P波到时，截取P波前与S波时间窗等长的部分确定噪声水平. 对于S波时间窗，使用Husid图识别S波到时，使用累积均方根函数S波结束的时间^[14]，保证S波的时间窗至少为4 s，以确保频谱分辨率. 图3（a）为该识别过程的示例.

图 3 强震动观测记录的处理示例

Figure 3. Processing of strong motion observation records

下载: 全尺寸图片幻灯片

3）截取S波时间窗后，首先，对于时域S波波形进行基线校正处理；接着，采用四阶Butterworth非因果滤波器进行双向滤波处理，设置滤波频率范围为0.4～20.0 Hz. 为了避免信号泄露，滤波完成后，在波形信号前后添加5%的汉宁窗函数^[15].

4）经上述处理后，对S波波形快速傅里叶变换，并使用带宽为20的Konno-Ohmachi滤波器进行平滑处理^[16]，以消除毛刺的影响. 最后，对0.4～20.0 Hz频域之间40个均匀对数间隔内的数据进行插值采样，处理效果如图3（b）所示.

2. 研究方法

2.1 线性反演

远场地震动S波的加速度记录O_ij在频域内可以表示为震源S_i、传播路径P_ij、场地G_i的乘积，如式（1）.

${O_{ij}}(f) = {S_i}(f) {P_{ij}}(f) {G_i}(f) \text{，}$

(1)

式中：i、j分别代表事件与台站，f为频率.

传播路径项P常仅考虑几何衰减 $\gamma ({R_{ij}})$ 和地震波经过地壳介质的非弹性衰减（通常表现为品质因子Q_s(f)的函数），如式（2）.

${P_{ij}} = \exp ( - {\text{π}} {R_{ij}}f/{Q_{\text{s}}}(f){V_{\text{s}}}) \gamma ({R_{ij}}) \text{，}$

(2)

式中： $R_{ij}$ 为震中距，V_s为3.6 km/s区域的地壳平均S波波速.

采用铰链三段式模型来表示几何衰减 $\gamma ({R_{ij}})$ ^[12]，该模型解释了地壳与地幔的边界对地震波传播的影响^[17]，如式（3）.

$\gamma ({R_{ij}})=\left\{ {\begin{array}{*{20}{l}} {R_{ij}^{ - {b_1}}},\quad {R \leqslant {R_{01}}}, \\ {R_{01}^{ - {b_1}} R_{01}^{ - {b_2}} R_{ij}^{ - b_2}},\quad{{R_{01}} < R \leqslant {R_{02}}}, \\ {R_{01}^{ - {b_1}} R_{01}^{{b_2}} R_{02}^{ - {b_2}} R_{02}^{{b_3}} R_{ij}^{ - b_3}},\quad{R > {R_{02}}} , \end{array}} \right.$

(3)

式中： R₀₁=1.5D，R₀₂=2.5D，b₁=1.0，b₂=0，b₃=0.5，D为地壳厚度，R为震中距离.

在确定传播路径P_ij之后，选取记录同一地震事件的露头岩石场地的观测数据作为约束. 通过土层台站与基岩台站的频谱比，直接消除同一地震事件的震源项S_i. 随后，对等式两侧同时取对数，最终表示为

${g_j} + {a_{ij}}Q_{\mathrm{s}}^{ - 1} = {d_{ij}} \text{，}$

(4)

式中： ${g_j} = \ln \dfrac{{{G_j}}}{{{G_{\mathrm{r}}}}}$ ； ${a_{ij}} = \dfrac{{{\text{π}} f({R_{i{\mathrm{r}}}} - {R_{ij}})}}{{{V_{\mathrm{s}}}}}$ ； ${d_{ij}} = \ln \dfrac{{{O_{ij}}\gamma ({R_{ij}})}}{{{O_{i{\mathrm{r}}}}\gamma ({R_{i{\mathrm{r}}}})}}$ ，下标r为参考台站.

将式（4）用矩阵的形式表示为

${\boldsymbol{Ax}} = {\boldsymbol{b}}\text{，}$

(5)

式中：x为模型参数矩阵，由式（4）中的场地效应g_j和Q_s组成.

线性反演的过程就是寻找式（5）中x的解，可通过最小二乘QR分解法一步求解^[18].

2.2 强震记录的HVSR法

HVSR求解场地效应可以直接表示^[6]为

${G_{{\mathrm{HVSR}}}}(f) = \frac{1}{N}\sum\limits_{{i} = 1}^N {\frac{{{{(|S_{\mathrm{H}}(f)||S_{\mathrm{V}}(f)|)}^{\frac{1}{2}}}}}{{U_{\mathrm{D}}(f)}}},$

(6)

式中： ${G_{{\mathrm{HVSR}}}}(f)$ 为强震记录HVSR的场地效应结果，N为记录地震事件的数量，S_H（f）、S_V（f）、U_D（f）分别为强震记录SH波、SV波以及UD分量的傅里叶幅值谱.

3. 几何扩散函数与Q_s值

线性反演法中，衰减模型的合理性对反演的可靠程度具有重要影响. 四川西部地区地壳结构复杂，区域平均地壳厚度在川西北地区的64 km逐渐减少到四川盆地的43 km，而在龙门山断裂带附近，地壳厚度从西北侧的52.5 km减少到东南侧的41.5 km^{[ 19]}. 为了构建区域几何衰减模型，该研究选取龙门山断裂两侧平均地壳厚度47 km，以此确定R₀₁为70 km、R₀₂为120 km（根据式（3）），进而得到式（7），该公式描述了考虑莫霍反射影响的龙门山区域三段式几何衰减模型.

$\gamma(R_{ij})=\left\{\begin{array}{*{20}{l}}R^{-1},\quad R\leqslant70\; \mathrm{km}， \\ 1/70,\quad70\; \mathrm{km} < R\leqslant120\; \mathrm{km}， \\ \left(1/70\right)\times(120/R)^{0.5},\quad R > 120\; \mathrm{km}.\end{array}\right.$

(7)

根据三段式几何衰减函数，反演了区域内S波的品质因子Q_s，随后使用正交最小二乘回归得到Q_s与频率的相关性，如式（8）.

${O_{\mathrm{s}}}(f) = 199.2{f^{0.8}}.$

(8)

常用Q_s(f)和Q_s(f)⁻¹来讨论地震波的衰减，此处采用与衰减成正比的Q_s（f）⁻¹进行讨论. 附近区域的其他研究结果与本研究的对比见图4^[20-22]，可以发现：不同区域内差异性显著. 表征品质因子的参数方程Q_s（f）=Q₀f ⁿ中，n值（频率相关程度）越大，区域内介质衰减速率越快. Q₀ （频率为1 Hz的Q_s值）大小可能与研究区域的震源距有关，震源距越大，地震波更容易在相对较深的介质中传播，通常地壳深部介质的Q_s值较大^[20].

图 4 品质因子Q_s值的对比

Figure 4. Comparison of quality factor Q_s values

下载: 全尺寸图片幻灯片

4. 场地效应分析

4.1 HVSR的理论基础

强震记录HVSR方法比地脉动的H/V谱比方法更为复杂，震源机制、震中距（以及场地各向异性、非线性等因素）都可能会对场地效应谱比结果造成影响^[23].

从机理的角度分析，对于HVSR结果近似为水平向S波场地效应的结论，Nakamura^[6]提出了2个重要的假设：1）首先，假设在坚硬的基底，地震波的传播方式具有各向同性的特点，地震基岩处地震波的水平和垂直分量大致相等. 该假设已通过Nakamura进行的井下微震测量实验得到验证；同时，对于强震动记录，Rong等^[24]也证实了地下基岩处水平分量与竖向分量的比值近似等于常数1. 2）假设垂直分量不受局部场地效应影响，即土层对地震波的垂直分量没有放大作用. 因此HVSR理论可以通过式（9）进行描述.

${G_{{\mathrm{HVSR}}}}(f) = \frac{{H(f)}}{{{H_{\mathrm{b}}}(f)}} \frac{{{H_{\mathrm{b}}}(f)}}{{{V_{\mathrm{b}}}(f)}} \frac{{{V_{\mathrm{b}}}(f)}}{{V(f)}}\text{，}$

(9)

式中： $H(f)$ 与 $V(f)$ 分别为地表强震记录的水平与垂直分量， ${H_{\mathrm{b}}}(f)$ 与 ${V_{\mathrm{b}}}(f)$ 分别为地下基岩强震记录的水平与垂直分量.

当 ${H_{\mathrm{b}}}(f) = {V_{{\mathrm{b}}} }(f)$ , ${V_{\mathrm{b}}}(f) = V(f)$ 时， ${G_{{\mathrm{HVSR}}}}$ 等于地表水平分量傅里叶谱与地下水平分量傅里叶谱的比值，即为场地效应.

4.2 方位角对HVSR的影响

HVSR理论简化了地震波的传播过程. 该方法基于基岩地震动的各向同性（根据Nakamura的假设1）消除了控制地震地面运动三要素的震源特性与传播路径因素. 然而，大小地震的破裂过程通常具有复杂的方向性效应. HVSR方法中，这种不确定的方向性作用均被归为场地效应. 一些研究发现，将强震记录旋转至地震与台站的大圆弧方向（SH波与SV波振动方向）进行H/V谱比似乎更有效^[23]，该研究数据处理阶段已经考虑到了这一点. 为了深入研究方向性效应对HVSR结果的影响，该研究展开了以下讨论：

1）对SH波、SV波以及水平分量的幅值差异进行了检验，以台站51HSD记录的地震事件为例，通过傅里叶谱比方法进行分析，如图5. 分析表明，水平分量的傅里叶幅值大于SH波和SV波. 此外，SH波、SV波的幅值同样具有差异，尤其是在低频部分，而中高频较为稳定. 低频长周期信号的异常通常由记录质量或背景噪声引起，导致SH、SV波在低频段展现出固有的差异. 该差异与方位角的影响较小，在后续的讨论中应进行区分.

图 5 SH、SV波、H分量的对比

Figure 5. Comparison of SH, SV wave, and H component

下载: 全尺寸图片幻灯片

2）该研究对来自不同方向的多次地震事件进行HVSR场地效应分析（多向HVSR），揭示了场地效应具有明显的优势方向. 以51HSD台站为例进行的多向HVSR显示，在震中相对于台站的SE方向（约135°）时，HVSR的场地放大系数最大，随方位角的增减而减小（图6（a））. 此外，水平分量的谱比结果最为显著，而SV波的谱比结果相对最小. 值得注意的是，不同方位角对场地卓越频率估计的影响较小（图6（b）），解释了使用HVSR方法估计场地卓越频率显示出稳定性的原因.

图 6 多向HVSR的结果

Figure 6. Results of multidirectional HVSR

下载: 全尺寸图片幻灯片

3）使用旋转HVSR的方法，进一步讨论方位角的影响. 以51HSD台站为例，以震中至台站的大圆弧方向作为0° 基准，对一次强震动波形数据在0～180° 的范围内，以5° 为间隔进行旋转，并计算了每个旋转间隔下的HVSR场地效应. 如图7所示，不同旋转方位角对HVSR结果有显著影响，并且在特定的旋转角度，HVSR的场地放大结果出现显著的“突变”，该现象验证了HVSR结果具有优势方向的观点. SH波与SV波的旋转HVSR结果存在差异，尽管“突变”位置有角度差异，场地放大值的差异却相对较小. 这种现象初步被认为是由于SH波和SV波固有的方向性以及其在低频的差异所导致的. 需要指出的是，采用单台单次记录进行旋转分析，结果可能具有一定的随机性，但多次强震记录求平均值的方式从机理上混淆了方位角，因此排除了使用该方案.

图 7 旋转HVSR的结果

Figure 7. Results of rotational HVSR

下载: 全尺寸图片幻灯片

4.3 线性反演与HVSR对比

为了确保安全性，主要基于幅值较大的水平分量记录进行场地效应评估. 比较基于S波水平分量的HVSR和线性反演法得出的场地效应值发现，这2种方法的结果存在显著差异（图8）.

图 8 HVSR与线性反演法场地效应的对比

Figure 8. Comparison of site effects between HVSR and linear inversion

下载: 全尺寸图片幻灯片

基于HVSR的2个关键理论基础，本文首先考察了地震波在基岩中传播的理论差异，如方位角的影响. 在场地效应评估中，通过采用多次记录的平均值，可以有效减轻方位角的影响. 其次，土层的竖向放大作用可能是导致HVSR与线性反演法在场地效应评估中存在差异的一个关键因素. 因此，部分研究重点强调竖向场地放大的作用^[25].

本文重新对汶川垂直分量波形记录进行了线性反演分析，以评估区域内土层的竖向放大效应. 如图9所示，大部分观测台站中，土层的竖向放大效应显著. 该现象将对HVSR的场地效应结果产生重要影响. 此时，式（9）中竖向放大被低估（高估），竖向放大与P波的传播相关，P波主要来自SV波的转换，而非P波的特定时间窗口部分. SV波将能量传递至垂直方向，进而导致了HVSR失效^[25]. 因此，为了补偿垂直放大对HVSR评估的影响，HVSR场地效应的修正方法可以表示为

图 9 线性反演计算的17个台站竖向场地效应的结果

Figure 9. Results of vertical site effects calculated by linear inversion for 17 stations

下载: 全尺寸图片幻灯片

${G_{{{\mathrm{RHVSR}}} }}(f) = {G_{{\mathrm{HVSR}}}}(f) {G_{\mathrm{V}}}(f),$

(10)

式中： ${G_{{{\mathrm{RHVSR}}} }}(f)$ 为竖向校正后的HVSR场地效应结果， ${G_{\mathrm{V}}}(f)$ 为线性反演法计算的土层竖向放大.

${G_{{{\mathrm{RHVSR}}} }}$ 与线性反演的水平向的对比结果如所示. 结果表明：大多数台站都存在着场地放大现象，且能清楚识别出卓越频率； ${G_{{{\mathrm{RHVSR}}} }}$ 的结果大于线性反演法，与HVSR结果偏小的普遍认识相悖. 该现象一方面表现出场地竖向校正产生了效果，使HVSR的结果接近于真实值；另一方面说明线性反演法的结果也存在着误差，具体原因可能是露头岩石台站51MXT存在场地效应. 参考台站场地效应的低估导致了线性反演结果偏小.

图 10 修正后的HVSR与线性反演法对比

Figure 10. Comparison of modified HVSR and linear inversion

下载: 全尺寸图片幻灯片

通过Delaunay三角测量方法将17个台站的线性反演结果以地理分布图的形式展示（图11，三角形标记是台站位置），场地效应与场地卓越频率呈现出一定的‘反比’关系，即场地效应随着卓越频率的增大而减小. 等值线地理图为研究区域内的场地效应、场地卓越频率提供参考价值.

图 11 区域场地放大系数及场地卓越频率的空间分布

Figure 11. Spatial distribution of regional site amplification factors and site predominant frequencies

下载: 全尺寸图片幻灯片

5. 结　论

本研究选取由17个强震动台站获取的21个汶川地震余震中233组强震动记录，通过线性反演法、强震记录HVSR法对比研究余震涉及区域内场地效应，主要工作及结论为

1）衰减因子Q_s值的估计. 结合龙门山地区实际情况，建立三段式衰减模型，本研究得到区域内地球介质衰减过程中S波的品质因子Q_s（f）=199.2f ^0.8.

2）方位角对HVSR的影响. 使用旋转HVSR、多向HVSR方案，研究了HVSR结果存在的方向性. 随着角度的变化，HVSR的放大值可能会出现突然放大的现象，但不同方向角对卓越频率的影响不大，该现象需要进一步的研究.

3） HVSR方法的经验修正. 竖向地震动的放大，即场地土层、工程基岩或中深部硬岩层耦合形成的竖向场地效应是HVSR与广义线性反演法产生差异的主要原因.

致谢：感谢中国地震局工程力学研究所强震动观测中心为本研究提供的数据支持；感谢研究所专项基金（DQJB24Z01、DQJB24R31、DQJB22Z03、DQJB21Z15）为本研究提供联合资助.

图 1 试验结果异常点检测

Figure 1. Outlier detection for test results

下载: 全尺寸图片幻灯片

图 2 相似材料配比流程

Figure 2. Flow chart of determining material proportion

下载: 全尺寸图片幻灯片

表 1 原型砂土与模型相似材料物理力学参数

Table 1. Mechanical parameters of prototype sand and model materials

项目	重度/（kN•m⁻³）	压缩模量/MPa	黏聚力/kPa	内摩擦角/（°）	渗透系数/（× 10⁻⁶cm•s⁻¹）
原型	17.8～20.1	35.30～66.80	13.00～25.10	25.3～32.1	100～600
模型	17.8～20.1	1.30～2.55	0.48～0.93	25.3～32.1	19～115

下载: 导出CSV

表 2 相似材料配比方案

Table 2. Test schemes of material proportions

试验工况	重晶石粉/ 标准砂	黏结剂浓度/%	凡硅比	石膏含量/%
1	0	2.0	1.0∶1.0	1.5
2	0	3.5	1.5∶1.0	3.0
3	0	5.0	2.0∶1.0	4.5
4	0	6.5	2.5∶1.0	6.0
5	0	8.0	3.0∶1.0	7.5
6	0.25	2.0	1.5∶1.0	4.5
7	0.25	3.5	2.0∶1.0	6.0
8	0.25	5.0	2.5∶1.0	7.5
9	0.25	6.5	3.0∶1.0	1.5
10	0.25	8.0	1.0∶1.0	3.0
11	0.50	2.0	2.0∶1.0	7.5
12	0.50	3.5	2.5∶1.0	1.5
13	0.50	5.0	3.0∶1.0	3.0
14	0.50	6.5	1.0∶1.0	4.5
15	0.50	8.0	1.5∶1.0	6.0
16	0.75	2.0	2.5∶1.0	3.0
17	0.75	3.5	3.0∶1.0	4.5
18	0.75	5.0	1.0∶1.0	6.0
19	0.75	6.5	1.5∶1.0	7.5
20	0.75	8.0	2.0∶1.0	1.5
21	1.00	2.0	3.0∶1.0	6.0
22	1.00	3.5	1.0∶1.0	7.5
23	1.00	5.0	1.5∶1.0	1.5
24	1.00	6.5	2.0∶1.0	3.0
25	1.00	8.0	2.5∶1.0	4.5

下载: 导出CSV

表 3 相似材料配比试验结果

Table 3. Experimental results of material proportions

试验工况	重度/ （kN•m⁻³）	压缩模量/MPa	内摩擦角/（°）	黏聚力/kPa	k₂₀/ （× 10⁻⁶ cm•s⁻¹）
1	17.5	2.22	28.3	0.26	701.00
2	17.8	2.93	26.7	14.66	24.80
3	17.7	2.52	26.0	14.14	15.10
4	16.4	1.91	20.7	25.14	525.00
5	16.6	6.25	25.6	0.52	439.00
6	19.0	1.82	29.6	0.74	2 330.00
7	18.8	8.28	34.9	0.83	123.00
8	19.1	4.08	35.4	0.01	69.20
9	19.5	2.28	35.8	2.30	16.20
10	19.0	1.82	31.4	4.73	1.87
11	20.2	5.15	27.2	18.25	45.40
12	20.5	2.53	36.1	0.21	1.18
13	19.4	2.29	34.3	0.87	24.90
14	19.1	5.44	33.7	0.41	4.50
15	18.7	3.34	25.7	1.01	96.40
16	21.0	1.85	37.2	0.63	3.90
17	20.3	5.11	37.3	1.13	109.00
18	20.5	1.90	33.2	0.98	6.20
19	19.3	3.24	31.3	0.61	0.85
20	20.2	1.03	30.2	0.67	0.17
21	21.3	14.64	27.2	3.50	140.00
22	21.0	7.43	23.1	34.30	0.57
23	20.8	1.00	32.3	1.57	0.47
24	20.6	2.16	34.4	1.92	0.46
25	19.9	1.74	30.4	2.73	0.42

下载: 导出CSV

表 4 重度极差分析

Table 4. Range analysis of weight density kN/m³

水平	重晶石粉/ 标准砂	黏结剂浓度	凡硅比	石膏含量
1	17.2	19.8	19.4	19.7
2	19.1	19.7	19.1	19.6
3	19.6	19.5	19.5	19.2
4	20.3	19.0	19.4	19.1
5	20.7	18.9	19.4	19.2
极差	3.5	0.9	0.4	0.5

下载: 导出CSV

表 5 压缩模量极差分析

Table 5. Range analysis of compression modulus MPa

水平	重晶石粉/ 标准砂	黏结剂浓度	凡硅比	石膏含量
1	3.17	5.14	3.76	1.81
2	3.66	5.25	2.46	2.21
3	3.75	2.36	3.83	3.32
4	2.62	3.00	2.42	6.01
5	5.39	2.84	6.11	5.23
极差	2.77	2.89	3.69	4.20

下载: 导出CSV

表 6 内摩擦角极差分析

Table 6. Range analysis of internal friction angle （°）

水平	重晶石粉/ 标准砂	黏结剂浓度	凡硅比	石膏含量
1	25.45	29.90	29.94	32.54
2	33.42	31.62	29.12	32.80
3	31.40	32.24	30.54	31.40
4	33.84	31.18	31.96	28.34
5	29.48	28.66	32.04	28.52
极差	7.97	3.58	2.92	4.46

下载: 导出CSV

表 7 黏聚力极差分析

Table 7. Range analysis of cohesion force kPa

水平	重晶石粉/ 标准砂	黏结剂浓度	凡硅比	石膏含量
1	10.34	4.68	8.14	1.00
2	1.72	10.23	3.72	4.56
3	4.15	3.51	7.16	3.83
4	0.81	6.07	5.74	6.29
5	8.80	1.93	1.67	10.74
极差	9.53	8.30	6.47	9.74

下载: 导出CSV

表 8 渗透系数极差分析

Table 8. Range analysis of permeability coefficient × 10⁻⁶cm/s

水平	重晶石粉/ 标准砂	黏结剂浓度	凡硅比	石膏含量
1	341	224	143	144
2	88	109	71	11
3	34	108	37	72
4	24	53	120	178
5	28	23	146	111
极差	317	201	109	167

下载: 导出CSV

表 9 多因素显著性分析

Table 9. Multivariate significance analysis

参数	重晶石粉/ 标准砂	黏结剂浓度	凡硅比	石膏含量
重度	0.875	−0.270	0.28	−0.143
压缩模量	0.165	−0.330	0.224	0.512
内摩擦角	0.266	−0.092	0.222	−0.393
黏聚力	−0.184	−0.156	−0.176	0.242
渗透系数	−0.332	−0.304	−0.109	0.030

下载: 导出CSV

表 10 多元线性回归分析系数表

Table 10. Coefficients of multivariate linear regression analysis

i	${\beta _{0,i}}$	${\beta _{1,i}}$	${\beta _{2,i}}$	${\beta _{3,i}}$	${\beta _{4,i}}$
1	18.87	3.29	−16.93	0.05	−8.93
2	0.26	1.37	−45.67	0.93	70.93
3	30.94	3.38	−19.47	1.41	−83.33
4	7.54	−2.08	−64.25	−2.18	141.35
5	246.91	−275.73	−1 168.62	11.12	675.61

下载: 导出CSV

表 11 多元线性回归分析系数（不含异常点）

Table 11. Multiple linear regression analysis coefficients （excluding abnormal points）

i	${\beta _{0,i}}$	${\beta _{1,i}}$	${\beta _{2,i}}$	${\beta _{3,i}}$	${\beta _{4,i}}$
1	19.66	3.05	−20.88	−0.07	−12.88
2	0.15	1.11	−43.49	1.06	66.57
3	30.94	1.92	−7.29	1.77	−71.15
4	7.29	−3.68	−50.94	−0.56	32.21
5	−133.31	−161.60	732.50	68.10	2 576.70

下载: 导出CSV

表 12 参考组配比

Table 12. Material proportion of reference group

参考组	重晶石粉/ 标准砂	黏结剂浓度/%	凡硅比	石膏含量/%
1	0.50	5.0	3∶1	3
2	0.45	5.5	3∶1	3
3	0.40	4.5	3∶1	3

下载: 导出CSV

表 13 参考组试验结果

Table 13. Test results of reference group

参考组	重度/ （kN•m⁻³）	压缩模量/MPa	内摩擦角/（°）	黏聚力/ kPa	k₂₀/ ( ×10⁻⁶ cm•s⁻¹）
1	19.4	2.29	34.3	0.87	24.9
2	19.2	1.89	32.1	0.81	36.8
3	18.9	2.01	31.7	0.63	52.1

下载: 导出CSV

参考文献(15)

左保成,陈从新,刘才华,等. 相似材料试验研究[J]. 岩土力学,2004,25(11): 1805-1808. doi: 10.3969/j.issn.1000-7598.2004.11.027

ZUO Baocheng, CHEN Congxin, LIU Caihua, et al. Research on similar material of slope simulation experiment[J]. Rock and Soil Mechanics, 2004, 25(11): 1805-1808. doi: 10.3969/j.issn.1000-7598.2004.11.027

郭信君. 超大直径越江隧道泥水盾构掘进模型试验研究[J]. 铁道标准设计,2008(10): 71-75. doi: 10.3969/j.issn.1004-2954.2008.10.022

GUO Xinjun. Model test study on slurry shield tunneling of super-large diameter river tunnel[J]. Railway Standard Design, 2008(10): 71-75. doi: 10.3969/j.issn.1004-2954.2008.10.022

何川,汪洋,方勇,等. 土压平衡式盾构掘进过程的相似模型试验[J]. 土木工程学报,2012,45(2): 162-169.

HE Chuan, WANG Yang, FANG Yong, et al. Similarity model test of earth-pressure-balanced shield tunneling process[J]. China Civil Engineering Journal, 2012, 45(2): 162-169.

徐前卫. 盾构施工参数的地层适应性模型试验及其理论研究[D]. 上海: 同济大学, 2006.

韩伯鲤,陈霞龄,宋一乐,等. 岩体相似材料的研究[J]. 武汉水利电力大学学报,1997,4(2): 6-9.

HAN Boli, CHEN Xialing, SONG Yile. Research on similar material of rockmass[J]. Journal of Wuhan University of Hydraulic and Electric Engineering, 1997, 4(2): 6-9.

张强勇,李术才,郭小红,等. 铁晶砂胶结新型岩土相似材料的研制及其应用[J]. 岩土力学,2008,29(8): 2126-2130. doi: 10.3969/j.issn.1000-7598.2008.08.021

ZHANG Qiangyong, LI Shucai, GUO Xiaohong, et al. Research and development of new typed cementatious geotechnical similar material for iron crystal sand and its application[J]. Rock and Soil Mechanics, 2008, 29(8): 2126-2130. doi: 10.3969/j.issn.1000-7598.2008.08.021

孙吉主,肖文辉. 盾构模型试验研究现状及新设想[J]. 土工基础,2011,25(6): 49-52. doi: 10.3969/j.issn.1004-3152.2011.06.015

SUN Jizhu, XIAO Wenhui. State of the art of shield tunnel model tests[J]. Soil Engineering and Foundation, 2011, 25(6): 49-52. doi: 10.3969/j.issn.1004-3152.2011.06.015

王凯,李术才,张庆松,等. 流-固耦合模型试验用的新型相似材料研制及应用[J]. 岩土力学,2016,37(9): 2521-2533.

WANG Kai, LI Shucai, ZHANG Qingsong, et al. Development and application of new similar materials of surrounding rock for a fluid-solid coupling model test[J]. Rock and Soil Mechanics, 2016, 37(9): 2521-2533.

李术才,周毅,李利平,等. 地下工程流-固耦合模型试验新型相似材料的研制及应用[J]. 岩石力学与工程学报,2012,31(6): 1128-1137. doi: 10.3969/j.issn.1000-6915.2012.06.006

LI Shucai, ZHOU Yi, LI Liping, et al. Development and application of a new similar material for underground engineering fluid-solid coupling model test[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(6): 1128-1137. doi: 10.3969/j.issn.1000-6915.2012.06.006

李树刚,赵鹏翔,林海飞,等. 煤岩瓦斯“固-气”耦合物理模拟相似材料特性实验研究[J]. 煤炭学报,2015,40(1): 80-86.

LI Shugang, ZHAO Pengxiang, LIN Haifeng, et al. Study on character of physical simulation similar material of coal-rock and gas solid-gas coupling[J]. Journal of the China Coal Society, 2015, 40(1): 80-86.

武伯弢, 朱合华, 徐前卫, 等. Ⅳ级软弱围岩相似材料的试验研究[J]. 岩土力学, 2013, 34（增刊1）: 109-116.

WU Botao, ZHU Hehua, XU Qianwei, et al. Experimental study of similar material for weak surrounding rock mass of class Ⅳ[J]. Rock and Soil Mechanics, 2013, 34（S1）: 109-116.

史小萌,刘保国,肖杰. 水泥和石膏胶结相似材料配比的确定方法[J]. 岩土力学,2015,36(5): 1357-1362.

SHI Xiaomeng, LIU Baoguo, XIAO Jie. A method for determining the ratio of similar materials with cement and plaster as bonding agents[J]. Rock and Soil Mechanics, 2015, 36(5): 1357-1362.

董金玉,杨继红,杨国香,等. 基于正交设计的模型试验相似材料的配比试验研究[J]. 煤炭学报,2012,37(1): 44-49.

DONG Jinyu, YANG Jihong, YANG Guoxiang, et al. Research on similar material proportioning test of model test based on orthogonal design[J]. Journal of the China Coal Society, 2012, 37(1): 44-49.

窦远明,王建宁,田贵州,等. 基于正交试验的软弱土质相似材料配比研究[J]. 铁道科学与工程学报,2017,14(3): 480-487. doi: 10.3969/j.issn.1672-7029.2017.03.007

DOU Yuanming, WANG Jianning, TIAN Guizhou, et al. Orthogonal test study on the mixture ratio of soft soil similar material[J]. Journal of Railway Science and Engineering, 2017, 14(3): 480-487. doi: 10.3969/j.issn.1672-7029.2017.03.007

曹玉茹. 虚拟变量回归在SPSS中的分析与实现[J]. 统计与决策,2018,34(10): 66-69.

CAO Yuru. Analysis and implementation of virtual variable regression in SPSS[J]. Statistics and Decision, 2018, 34(10): 66-69.

施引文献

附加材料(0)

访问统计

点击查看大图

图(2) / 表(13)

计量

文章访问数: 463
HTML全文浏览量: 338
PDF下载量: 15
被引次数: 0

1. 强震数据
1.1 数据选取
1.2 数据处理
2. 研究方法
2.1 线性反演
2.2 强震记录的HVSR法
3. 几何扩散函数与Qs值
4. 场地效应分析
4.1 HVSR的理论基础
4.2 方位角对HVSR的影响
4.3 线性反演与HVSR对比
5. 结　论

模拟深海环境砂土地层的材料配比试验研究

doi: 10.3969/j.issn.0258-2724.20180285

作者简介: 沈翔（1991—），男，博士研究生，研究方向为隧道与地下工程，E-mail：xshenbjtu@126.com

通讯作者: 袁大军（1961—），男，教授，博士，研究方向为隧道与地下工程，E-mail：djyuan603@163.com

计量

出版历程

Experiments on Material Proportions for Simulating Sandy Layer in Deep Sea

1. 强震数据

1.1 数据选取

1.2 数据处理

2. 研究方法

2.1 线性反演

2.2 强震记录的HVSR法

3. 几何扩散函数与Qs值

4. 场地效应分析

4.1 HVSR的理论基础

4.2 方位角对HVSR的影响

4.3 线性反演与HVSR对比

5. 结 论

计量

出版历程

目录

1. 强震数据

1.1 数据选取

1.2 数据处理

2. 研究方法

2.1 线性反演

2.2 强震记录的HVSR法

3. 几何扩散函数与Qs值

4. 场地效应分析

4.1 HVSR的理论基础

4.2 方位角对HVSR的影响

4.3 线性反演与HVSR对比

5. 结 论

作者简介:
沈翔（1991—），男，博士研究生，研究方向为隧道与地下工程，E-mail：xshenbjtu@126.com

通讯作者:
袁大军（1961—），男，教授，博士，研究方向为隧道与地下工程，E-mail：djyuan603@163.com

3. 几何扩散函数与Q_s值

5. 结　论

3. 几何扩散函数与Q_s值

5. 结　论