CN106094022B - 一种土层等效剪切波速和土层反射界面埋深的测量方法 - Google Patents
一种土层等效剪切波速和土层反射界面埋深的测量方法 Download PDFInfo
- Publication number
- CN106094022B CN106094022B CN201610395478.8A CN201610395478A CN106094022B CN 106094022 B CN106094022 B CN 106094022B CN 201610395478 A CN201610395478 A CN 201610395478A CN 106094022 B CN106094022 B CN 106094022B
- Authority
- CN
- China
- Prior art keywords
- soil layer
- formula
- shear wave
- buried depth
- wave velocity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000002689 soil Substances 0.000 title claims abstract description 295
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000005259 measurement Methods 0.000 claims abstract description 28
- 238000010008 shearing Methods 0.000 claims description 43
- 238000005553 drilling Methods 0.000 claims description 25
- 238000012360 testing method Methods 0.000 claims description 22
- 239000003673 groundwater Substances 0.000 claims description 17
- 238000012937 correction Methods 0.000 claims description 14
- 238000001028 reflection method Methods 0.000 claims description 11
- 238000012886 linear function Methods 0.000 claims description 10
- 238000012888 cubic function Methods 0.000 claims description 8
- 239000004927 clay Substances 0.000 claims description 7
- 238000011835 investigation Methods 0.000 claims description 7
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000013461 design Methods 0.000 claims description 5
- 239000011521 glass Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims 2
- 238000012876 topography Methods 0.000 claims 2
- 238000010304 firing Methods 0.000 claims 1
- 238000006467 substitution reaction Methods 0.000 claims 1
- 239000004576 sand Substances 0.000 description 4
- 238000009933 burial Methods 0.000 description 2
- 235000013312 flour Nutrition 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009189 diving Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明涉及土层等效剪切波速和土层反射界面埋深的测量方法，解决直接利用土层剪切波速和埋深关系的数学模型及其拟合参数测量土层等效剪切波速和土层反射界面埋深的方法，根据场地地震工程地质条件划分地震工程地质单元；在同一个地震工程地质单元内，对土层剪切波速数据统计，建立土层剪切波速和埋深关系的数学模型及其拟合参数；将拟合参数代入与土层剪切波速和埋深关系的数学模型相对应的土层等效剪切波速计算式，得该地震工程地质单元的土层等效剪切波速；将拟合参数代入与土层剪切波速和埋深关系的数学模型相对应的土层反射界面埋深计算式，得土层反射界面埋深，本发明提高了土层等效剪切波速和土层反射界面埋深测量的便捷性和测量精度。
Description
技术领域
本发明涉及地震工程领域，特别是一种土层等效剪切波速和土层反射界面埋深的测量方法。
背景技术
土层等效剪切波速是地震工程中的重要参数之一，主要用于建筑的场地类别划分，从而为建筑抗震设计提供依据。土层反射界面埋深是地震勘探的重要结果之一，主要为查明地下地质构造(含活动断层)和矿产等服务。
依据我国现行建筑抗震设计规范(GB 50011-2010)，2010)，建筑的场地类别，应根据土层等效剪切波速和场地覆盖层厚度划分。其中，土层等效剪切波速的计算式为：
式(1)中，vse为土层等效剪切波速(单位m·s-1)，d0为计算深度(单位m)，取场地覆盖层厚度和20m两者的小值，t为剪切波在地面至计算深度d0之间的传播时间(单位s)；即土层等效剪切波速是地面至计算深度d0处的土层剪切波速平均值。
钻孔剪切波速测试是取得钻孔内各个埋深处的土层剪切波速数据的重要手段。目前，在岩土工程勘察、地震安全性评价以及地震小区划等项目中，在场地钻孔中开展了大量钻孔剪切波速测试工作，积累了大批土层剪切波速数据。许多学者(陈国兴等，1998；刘红帅等，2010；邱志刚等，2011)在分析研究钻孔中土层剪切波速测试数据的基础上，对土层剪切波速和埋深关系进行了统计分析，常用的数学模型有：线性函数模型、一元二次函数模型和幂函数模型等，获得的模型参数即拟合参数带有明显的当地(统计区域)土层特征；一般认为这些土层剪切波速和埋深关系的数学模型及其拟合参数主要用于对当地(统计区域)土层剪切波速数值的初步估计，而难以有其它的用途，其主要原因在于：(1)在进行土层剪切波速和埋深关系的统计时，统计区域范围偏大，在统计区域范围内的场地的地震工程地质条件差异较大，因而获得的土层剪切波速和埋深关系的数学模型及其拟合参数的代表性较差；(2)通过大量统计工作后仅得到了土层剪切波速和埋深关系的数学模型及其拟合参数，缺乏能够利用这些数学模型及其拟合参数直接测量土层等效剪切波速的方法。
另一方面，在横波(剪切波)反射法地震勘探中，通过反射法地震勘探数据处理后首先得到地震反射时间剖面，在之后的地质解释中，需要把地震反射时间剖面转换为地质剖面，即所谓的时深转换，其计算式为：
式(2)中，H0为土层剪切波反射波组反射界面埋深，简称土层反射界面埋深(m),v0为剪切波自地面至土层剪切波反射波组反射界面处的土层剪切波速平均值(m·s-1)，目前常常由反射法地震勘探数据计算给出或者土层测井数据给出；t0为剪切波反射波组反射界面的双程垂直到时(s)，由地震反射时间剖面给出。
目前，式(2)中v0的取值直接来源于钻孔剪切波速测试数据的情况不多，其原因在于目前缺乏直接应用钻孔剪切波速测试数据测量土层剪切波反射波组反射界面埋深的具体方法。
发明内容
针对上述情况，为克服现有技术缺陷，本发明之目的就是提供一种土层等效剪切波速和土层反射界面埋深的测量方法，填补直接利用土层剪切波速和埋深关系的数学模型及其拟合参数测量土层等效剪切波速和土层反射界面埋深的方法空白，丰富钻孔剪切波速测试数据的处理方法。
本发明解决的技术方案是，1、根据场地地震工程地质条件划分地震工程地质单元；2、在同一个地震工程地质单元内，对土层剪切波速数据进行统计，建立土层剪切波速和埋深关系的数学模型及其拟合参数；3、将拟合参数代入与土层剪切波速和埋深关系的数学模型相对应的土层等效剪切波速计算式，即得到场地内该地震工程地质单元的土层等效剪切波速；4、将拟合参数代入与土层剪切波速和埋深关系的数学模型相对应的土层反射界面埋深计算式，即得到土层反射界面埋深。
本发明中土层等效剪切波速是地震工程中的重要参数之一，主要用于建筑的场地类别划分，从而为建筑抗震设计提供依据。土层反射界面埋深是地震勘探的重要结果之一，主要为查明地下地质构造(含活动断层)和矿产等服务。
本发明中土层等效剪切波速计算式为代数表达式，反射界面埋深计算式为代数表达式或代数方程，具有普遍意义；基于该方法，提高了土层等效剪切波速和土层反射界面埋深测量的便捷性和测量精度。
附图说明
图1为本发明测量方法的流程示意图。
图2为本发明的场地平面示意图。
图3为本发明的一个地震工程地质单元的土层剪切波速和埋深关系散点图。
其中，1—场地，2—地震工程地质单元，3—钻孔，4—剪切波反射法地震勘探测线。
具体实施方式
以下结合附图对本发明的具体实施方式作详细说明。
由图1-图3给出，本发明的测量方法，具体步骤如下：
(1)、根据场地地震工程地质条件划分场地的地震工程地质单元，方法是：
A、收集场地已有地震工程地质条件资料，或对场地进行地震工程地质条件调查、勘察或钻孔剪切波速测试，获得场地地震工程地质条件资料，所述的场地地震工程地质条件资料为场地地貌单元类型、场地地下水类型与单孔地下水水位埋深、单孔场地覆盖层厚度、不同深度的水平方向土层地质成因类型、水平方向土层整体结构(指土层岩性名称、土层分布、土层层面埋深)和土层剪切波速数据；所述的地震工程地质条件调查为携带场地地形图、地质罗盘、放大镜、GPS定位仪和测距仪，在场地及其邻区开展地形、地貌、地表地层(含土层)岩性、地质成因类型、地质构造、水井、泉的观察、测量和记录；所述的地震工程地质条件勘察为在场地上进行钻孔、采集钻孔内土层样品，对土层样品进行观察、记录和化验(土工试验)，并对该钻孔内单孔地下水水位埋深数值进行测量和记录；所述的钻孔剪切波速测试为在场地内进行钻孔后，将波速测试仪的检波器安置在钻孔内各个埋深处即各个测点处并紧贴孔壁进行剪切波速测试，钻孔内相邻测点的垂直间距为1-2m从而获得钻孔内各个埋深处的土层剪切波速数据；将地面至土层剪切波速大于500m·s-1且其下方土层剪切波速均不小于500m·s-1的土层顶面的距离确定为单孔场地覆盖层厚度，而场地覆盖层厚度由场地内多个单孔场地覆盖层厚度确定，在场地内同一个地震工程地质单元范围的场地覆盖层厚度取该地震工程地质单元内单孔场地覆盖层厚度的平均值；
B、将场地地貌单元类型相同、场地地下水类型相同、水平方向土层地质成因类型相同、不同深度的水平方向土层整体结构相近、单孔地下水水位埋深数值相差不大于3m、单孔场地覆盖层厚度数值相差不大于5m的区域划分为同一个地震工程地质单元；将一个场地按照其地震工程地质条件划分为一个地震工程地质单元或多个地震工程地质单元；所述水平方向土层整体结构相近指水平方向土层岩性名称相同、水平方向土层分布连续、同一土层层面埋深相差不大于5m；
(2)、在同一个地震工程地质单元内，对土层剪切波速数据进行统计，建立土层剪切波速和埋深关系的数学模型及其拟合参数，方法是：
A、当测量土层等效剪切波速时，埋深区间为地面至计算深度d0(m)，其中，计算深度d0取该地震工程地质单元内的场地覆盖层厚度和20m两者的小值(即场地覆盖层厚度数值小于20m时，取场地覆盖层厚度数值，当场地覆盖层厚度数值不小于20m时，取20m)；当测量土层反射界面埋深时,埋深区间为地面至钻孔内进行土层剪切波速测试的土层最大埋深，即埋深区间为地面至具有土层剪切波速数据的最大土层埋深处，是为了尽量扩大埋深区间的范围，也是为了充分利用土层剪切波速数据，埋深区间的范围越大，计算土层反射界面埋深的有效范围就越大；
B、采用下列5个土层剪切波速和埋深关系的数学模型，利用Origin软件对该地震工程地质单元内的土层剪切波速数据进行统计，分别得到下列5个数学模型的拟合参数及其校正决定系数(adjusted R2)、标准差(SD)：
土层剪切波速和埋深关系的数学模型为:
v＝ah+b h∈[h1,h2] 式(3)；
v＝ch2+eh+f h∈[h1,h2] 式(4)；
v＝jh3+kh2+lh+p h∈[h1,h2] 式(5)；
v＝qh4+rh3+uh2+wh+x h∈[h1,h2] 式(6)；
v＝ghi h∈[h1,h2] 式(7)；
式(3)为线性函数模型的表达式，式(4)为一元二次函数模型的表达式，式(5)为一元三次函数模型的表达式、式(6)为一元四次函数模型的表达式、式(7)为幂函数模型的表达式；
式(3)-式(7)中：h为土层埋深(单位m)，v为土层埋深h处的土层剪切波速(单位m·s-1)，a、b、c、e、f、j、k、l、p、q、r、u、w、x、g、i为拟合参数，由土层剪切波速数据统计后获得，h1、h2分别为埋深区间[h1,h2]的两个端点即统计的土层埋深的上界限值(单位m)和下限界值(单位m)；
上述数学模型是以地面为土层埋深h的零点，以垂直向下方向为土层埋深h的正方向(增加方向)；
上述校正决定系数和标准差是用来评价统计结果优劣的主要指标；
将该地震工程地质单元内的土层剪切波速数据及其对应的土层埋深数据作为输入数据，分别输入到Origin软件的上述5个数学模型的统计计算程序中，即可分别获得上述5个数学模型的拟合参数及其校正决定系数和标准差；
C、根据5个数学模型统计结果中的校正决定系数、标准差，找出5个数学模型中校正决定系数值最大和标准差最小的一个数学模型，建立该地震工程地质单元的土层剪切波速和埋深关系的数学模型，并同时获得该数学模型的拟合参数；即5个数学模型中的哪一个数学模型中的校正决定系数值最大和标准差最小，该数学模型及其拟合参数越好，即越与该地震工程地质单元内的土层剪切波速数据相吻合；
(3)、将上述(2)中的C步骤建立的该地震工程地质单元的土层剪切波速和埋深关系的数学模型的拟合参数代入与该建立的数学模型相对应的土层等效剪切波速计算式中(例如，当建立的数学模型为一元二次函数模型时，将建立的一元二次函数模型的拟合参数代入与一元二次函数模型相对应的土层等效剪切波速计算式(即式(15)中)，得该地震工程地质单元的土层等效剪切波速；
与土层剪切波速和埋深关系的数学模型相对应的土层等效剪切波速计算式推导：
根据高等数学中的积分中值定理，当函数v(h)在闭区间[h1,h2]上连续，则在积分区间[h1,h2]上至少存在一个点ζ，使下式成立：
式(8)中，h为土层埋深(单位m)，v为土层埋深h处的土层剪切波速(单位m·s-1)，v(h)为土层剪切波速和埋深关系的数学模型即函数表达式，为土层剪切波速在埋深区间[h1,h2]上的平均值，h1、h2分别为埋深区间的两个端点，ζ为埋深区间[h1,h2]上的一点即h1≤ζ≤h2；
将式(3)、式(4)、式(5)、式(6)、式(7)分别代入式(8)，化简后，得埋深区间[h1,h2]上的任意深度段即任意埋深区间[H1,H2]上土层剪切波速平均值(单位m·s-1)的计算式分别为:
式(9)-式(13)中，为在埋深区间[h1,h2]上的土层剪切波速平均值(单位m·s-1)，a、b、c、e、f、j、k、l、p、q、r、u、w、x、g、i分别为式(3)-式(7)的拟合参数，h1、h2分别为式(3)-式(7)的埋深区间[h1,h2]的两个端点(单位m)，H1、H2分别为埋深区间[h1,h2]上的任意深度段即任意埋深区间[H1,H2]的两个端点(单位m)，H1≥h1，H1＜H2≤h2，
式(9)、式(10)、式(11)、式(12)、式(13)分别为与土层剪切波速和埋深关系的线性函数模型、一元二次函数模型、一元三次函数模型、一元四次函数模型、幂函数模型相对应的埋深区间[h1,h2]上的任意埋深区间[H1,H2]上土层剪切波速平均值计算式；
依据我国现行建筑抗震设计规范(GB 50011-2010)，2010)，土层等效剪切波速是地面至计算深度d0处的土层剪切波速平均值，由于地面的土层埋深为0.0m，计算深度d0处的土层埋深为计算深度d0，因此，取H1为0.0m，H2为计算深度d0，分别代入式(9)、式(10)、式(11)、式(12)、式(13)，化简后，得与土层剪切波速和埋深关系的线性函数模型、一元二次函数模型、一元三次函数模型、一元四次函数模型、幂函数模型相对应的土层等效剪切波速vse(单位m·s-1)计算式，分别为:
式(14)-式(18)中，vse为土层等效剪切波速(m·s-1)，d0为计算深度(m)，取该地震工程地质单元内的场地覆盖层厚度值和20m两者的小值，a、b、c、e、f、j、k、l、p、q、r、u、w、x、g、i分别为式(3)-式(7)的拟合参数；
采用式(14)-式(18)计算得到的土层等效剪切波速为场地中该地震工程地质单元的土层等效剪切波速，代表性强，稳定性高，精度高，可以消除由传统方法通过单个钻孔剪切波速测试数据计算土层等效剪切波速带来的偶然误差；
一个场地按照地震工程地质条件划分为一个地震工程地质单元或多个地震工程地质单元；当一个场地划分为一个地震工程地质单元时，场地的区域范围与地震工程地质单元的区域范围相同，场地的土层等效剪切波速与地震工程地质单元的土层等效剪切波速相同；当一个场地划分为N个(N为正整数)地震工程地质单元时，场地的区域范围是N个地震工程地质单元的区域范围之和，这时，重复上述步骤(2)、步骤(3)，分别得N个地震工程地质单元的土层等效剪切波速；而场地的土层等效剪切波速采用分区域表示，场地在第n个(n＝1，2，3，…，N)地震工程地质单元区域范围的土层等效剪切波速为第n个地震工程地质单元的土层等效剪切波速；
(4)、将上述(2)中的C步骤建立的该地震工程地质单元的土层剪切波速和埋深关系的数学模型的拟合参数代入与该建立的数学模型相对应的土层反射界面埋深计算式中(例如，当建立的数学模型为一元二次函数模型时，将建立的一元二次函数模型的拟合参数代入与一元二次函数模型相对应的土层反射界面埋深计算式(即式(25))中)，得到土层反射界面埋深；
与土层剪切波速和埋深关系的数学模型相对应的土层反射界面埋深计算式推导：
在横波(剪切波)反射法地震勘探中，通过反射法地震勘探数据处理后，得到地震反射时间剖面，把地震反射时间剖面转换为地质剖面时，要进行时深转换，其计算式为：
式(2)中，H0为土层剪切波反射波组反射界面埋深,简称土层反射界面埋深(单位m),v0为剪切波自地面至土层剪切波反射波组反射界面处的土层剪切波速平均值(单位m·s-1)，t0为剪切波反射波组反射界面的双程垂直到时(单位s)；
在横波(剪切波)反射法地震勘探中，设土层反射界面埋深H0位于上述土层剪切波速和埋深关系的埋深区间[h1,h2]上的任意埋深区间[H1,H2]上，由于地面的土层埋深为0.0m，土层剪切波反射波组反射界面处的土层埋深为土层反射界面埋深H0，因此，取H1为0.0m，取H2为H0，分别代入式(9)、式(10)、式(11)、式(12)、式(13)，化简后，得与土层剪切波速和埋深关系的线性函数模型、一元二次函数模型、一元三次函数模型、一元四次函数模型、幂函数模型相对应的剪切波自地面至土层剪切波反射波组反射界面处的土层剪切波速平均值v0，计算式分别为：
式(19)-式(23)中，v0为剪切波自地面至土层剪切波反射波组反射界面埋深H0处的土层剪切波速平均值(单位m·s-1)，H0为土层反射界面埋深(单位m)，a、b、c、e、f、j、k、l、p、q、r、u、w、x、g、i分别为式(3)-式(7)的拟合参数；
将式(19)、式(20)、式(21)、式(22)、式(23)分别代入式(2)，化简后,得与土层剪切波速和埋深关系的线性函数模型、一元二次函数模型、一元三次函数模型、一元四次函数模型、幂函数模型相对应的土层反射界面埋深H0的计算式分别为：
式(24)、式(25)、式(28)、式(26)、式(27)中，H0为土层反射界面埋深(单位m)；t0为剪切波反射波组反射界面的双程垂直到时(单位s)，a、b、c、e、f、j、k、l、p、q、r、u、w、x、g、i分别为式(3)-式(7)的拟合参数；
式(26)、式(27)采用迭代法求解，式(25)的计算结果以及式(26)、式(27)的根的取舍结合土层剪切波速和埋深关系统计时所采用的埋深区间约束条件进行判别。
采用式(24)、式(25)、式(28)、式(26)、式(27)得到的土层剪切波反射波组反射界面埋深，由于是基于当地钻孔剪切波速测试数据，因此，计算结果更符合当地实际地质情况。
使用本发明的测量方法反复进行了测试，均取得了良好的效果，具体如下：
某县城场地1面积为52.0km2,求该场地的土层等效剪切波速和剪切波反射法地震勘探测线4在该场地内某一土层剪切波反射波组反射界面埋深。
首先，按照上述测量方法的(1)中的A步骤，在收集该县城场地地震工程地质条件资料的基础上，得到该县城场地地震工程地质条件为，场地地貌单元类型为黄河冲积平原，场地地下水类型为潜水，单孔地下水水位埋深为8.0—10.0m，水平方向土层地质成因类型为冲积层，土层整体结构为层状结构，在埋深0—27m范围内，土层岩性为黄褐-灰黄色稍密-中密粉土、黄褐色可塑粉质粘土，夹黄褐色稍密粉砂，在埋深27—75m范围内，土层岩性为浅灰黄色中密粉砂、黄褐色中密-密实粉土、浅灰黄色中密细砂、黄褐色硬塑粉质粘土，在埋深75—100m范围内，土层岩性为浅灰黄色密实细砂和黄褐色硬塑粉质粘土；该场地水平方向土层岩性名称相同、水平方向土层分布连续、同一土层层面埋深相差不大于5m；得到了该场地内20个100m深度的钻孔3的土层剪切波速数据(参见附图3)；该场地上述20个100m深度的钻孔的单孔场地覆盖层厚度最小值为89m，最大值为94m，平均值为91.6m，该场地的单孔场地覆盖层厚度数值相差不大于5m，按照上述测量方法的(1)中的B步骤，将整个场地划分为同一个地震工程地质单元2；由于该场地仅划分为同一个地震工程地质单元，该场地的区域范围与该地震工程地质单元的区域范围相同，该场地的土层等效剪切波速与该地震工程地质单元的土层等效剪切波速相同；
然后，按照上述测量方法的(2)中的A步骤，为了测量该场地的土层等效剪切波速，确定埋深区间的两个端点分别为0m、20m，即埋深区间为[0，20]，这是因为计算深度d0(m)应取场地覆盖层厚度91.6m与20m二者的小值即取20m；为了测量土层反射界面埋深,埋深区间的两个端点分别为0m、98m，即埋深区间为[0，98]，这是因为在100m深度的钻孔内，由于钻孔底部存在泥渣的原因，能够进行土层剪切波速测试的最大深度为98m；按照上述测量方法的(2)中的B步骤，采用上述5个土层剪切波速和埋深关系的数学模型，利用Origin软件对该地震工程地质单元内的土层剪切波速数据在埋深区间[0，20]、[0，98]上分别进行统计，分别得到上述5个数学模型在埋深区间[0，20]、[0，98]上的拟合参数及其校正决定系数、标准差，见表1；
表1某县城场地土层剪切波速和埋深关系统计结果
按照上述测量方法的(2)中的C步骤，由表1可知，在上述5个数学模型的统计结果中，在埋深区间[0，20]、[0，98]上一元四次函数模型的校正决定系数值最大和标准差最小，因此，建立的该地震工程地质单元的土层剪切波速和埋深关系的数学模型为5个数学模型中的一元四次函数模型，采用在埋深区间[0，20]上一元四次函数模型及其拟合参数用来测量土层等效剪切波速；采用在埋深区间[0，98]上一元四次函数模型及其拟合参数用来测量土层反射界面埋深；按照上述测量方法的步骤(3)，将表1中一元四次函数模型在埋深区间[0，20]上的拟合参数q＝-0.00326、r＝0.1416、u＝-1.85548、w＝11.71934、x＝139.59112代入与土层剪切波速和埋深关系的一元四次函数模型相对应的土层等效剪切波速vse(单位m·s-1)计算式(17)，即得该地震工程地质单元的土层等效剪切波速为188.1m·s-1；由于该场地的区域范围与该地震工程地质单元的区域范围相同，该场地的土层等效剪切波速也为188.1m·s-1；
在该地震工程地质单元内采用同样的土层等效剪切波速数据，采用传统的土层等效剪切波速测量方法的结果为该县城场地(该场地1的范围与该地震工程地质单元2的范围相同)内20个钻孔的单孔土层等效剪切波速的最小值为178m·s-1，最大值为192m·s-1，这20个钻孔的单孔土层等效剪切波速的算术平均值为184.2m·s-1；该场地的单孔土层等效剪切波速的算术平均值184.2m·s-1作为该县城场地的土层等效剪切波速；
采用本发明测量方法获得的场地的土层等效剪切波速188.1m·s-1与传统测量方法的184.2m·s-1之间的差值仅为3.9m·s-1，经反复测试验证，采用本发明测量方法获得的场地的土层等效剪切波速误差小，稳定性高，精度高。
之后，测量位于该地震工程地质单元内反射法地震勘探中剪切波反射波组反射界面的双程垂直到时t0＝0.475s的土层反射界面埋深H0，按照上述测量方法的步骤(4)，将表1中一元四次函数模型在埋深区间[0，98]上的拟合参数q＝0.0000327、r＝-0.00513、u＝0.19832、w＝3.69908、x＝137.60362代入与土层剪切波速和埋深关系的一元四次函数模型相对应的土层反射界面埋深H0(m)式(27)，解之，即得双程垂直到时t0为0.475s的土层反射界面埋深H0为74.5m，其值与该地震工程地质单元内由钻孔揭示的土层界面75m仅相差0.5m，即计算结果具有较高的精度。
采用式(27)得到的土层反射界面埋深，由于是基于当地钻孔剪切波速测试数据，其剪切波速数据精度较高，从而其测量结果具有较高的精度。
本发明与现有技术相比，具有如下的特点：
1)、本发明填补了直接利用土层剪切波速和埋深关系的数学模型及其拟合参数测量土层等效剪切波速和土层反射界面埋深的方法空白，丰富了钻孔土层剪切波速测试数据的处理方法，提高了钻孔剪切波速测试数据的直接应用价值；
2)、本发明测量得到的场地的土层等效剪切波速可以消除或减少由单个或少量钻孔土层剪切波速测试数据测量场地的土层等效剪切波速带来的偶然误差，测量结果稳定性高，精度高，满足工程需要；
3)、本发明测量得到的土层剪切波反射波组反射界面埋深由于是基于当地钻孔剪切波速测试数据，其测量结果具有较高的精度；
4)、本发明提供的土层等效剪切波速计算式为代数表达式，土层反射界面埋深计算式为代数表达式或代数方程，具有普遍意义，使用便捷。
Claims (2)
1.一种土层等效剪切波速的测量方法，其特征在于，具体步骤如下：
(1)、根据场地地震工程地质条件划分场地的地震工程地质单元，方法是：
A、收集场地已有地震工程地质条件资料，或对场地进行地震工程地质条件调查、勘察或钻孔剪切波速测试，获得场地地震工程地质条件资料，所述的场地地震工程地质条件资料为场地地貌单元类型、场地地下水类型与单孔地下水水位埋深、单孔场地覆盖层厚度、不同深度的水平方向土层地质成因类型、水平方向土层整体结构和土层剪切波速数据；所述的地震工程地质条件调查为携带场地地形图、地质罗盘、放大镜、GPS定位仪和测距仪，在场地及其邻区开展地形、地貌、地表地层岩性、地质成因类型、地质构造、水井、泉的观察、测量和记录；所述的地震工程地质条件勘察为在场地上进行钻孔、采集钻孔内土层样品，对土层样品进行观察、记录和化验，并对该钻孔内单孔地下水水位埋深数值进行测量和记录；所述的钻孔剪切波速测试为在场地内进行钻孔后，将波速测试仪的检波器安置在钻孔内各个埋深处即各个测点处并紧贴孔壁进行剪切波速测试，钻孔内相邻测点的垂直间距为1-2m从而获得钻孔内各个埋深处的土层剪切波速数据；将地面至土层剪切波速大于500m·s-1且其下方土层剪切波速均不小于500m·s-1的土层顶面的距离确定为单孔场地覆盖层厚度，而场地覆盖层厚度由场地内多个单孔场地覆盖层厚度确定，在场地内同一个地震工程地质单元范围的场地覆盖层厚度取该地震工程地质单元内单孔场地覆盖层厚度的平均值；
B、将场地地貌单元类型相同、场地地下水类型相同、水平方向土层地质成因类型相同、不同深度的水平方向土层整体结构相近、单孔地下水水位埋深数值相差不大于3m、单孔场地覆盖层厚度数值相差不大于5m的区域划分为同一个地震工程地质单元；将一个场地按照其地震工程地质条件划分为一个地震工程地质单元或多个地震工程地质单元；所述水平方向土层整体结构相近指水平方向土层岩性名称相同、水平方向土层分布连续、同一土层层面埋深相差不大于5m；
(2)、在同一个地震工程地质单元内，对土层剪切波速数据进行统计，建立土层剪切波速和埋深关系的数学模型及其拟合参数，方法是：
A、当测量土层等效剪切波速时，埋深区间为地面至计算深度d0，其中，计算深度d0取该地震工程地质单元内的场地覆盖层厚度和20m两者的小值；当测量土层反射界面埋深时,埋深区间为地面至钻孔内进行土层剪切波速测试的土层最大埋深；
B、采用下列5个土层剪切波速和埋深关系的数学模型，利用Origin软件对该地震工程地质单元内的土层剪切波速数据进行统计，分别得到下列5个数学模型的拟合参数及其校正决定系数、标准差：
土层剪切波速和埋深关系的数学模型为:
v＝ah+b h∈[h1,h2] 式(3)；
v＝ch2+eh+f h∈[h1,h2] 式(4)；
v＝jh3+kh2+lh+p h∈[h1,h2] 式(5)；
v＝qh4+rh3+uh2+wh+x h∈[h1,h2] 式(6)；
v＝ghi h∈[h1,h2] 式(7)；
式(3)为线性函数模型的表达式，式(4)为一元二次函数模型的表达式，式(5)为一元三次函数模型的表达式、式(6)为一元四次函数模型的表达式、式(7)为幂函数模型的表达式；
式(3)-式(7)中：h为土层埋深，单位m，v为土层埋深h处的土层剪切波速，单位m·s-1，a、b、c、e、f、j、k、l、p、q、r、u、w、x、g、i为拟合参数，由土层剪切波速数据统计后获得，h1、h2分别为埋深区间[h1,h2]的两个端点即统计的土层埋深的上界限值和下限界值，单位m；
上述数学模型是以地面为土层埋深h的零点，以垂直向下方向为土层埋深h的正方向；
上述校正决定系数和标准差是用来评价统计结果优劣的主要指标；
将该地震工程地质单元内的土层剪切波速数据及其对应的土层埋深数据作为输入数据，分别输入到Origin软件的上述5个数学模型的统计计算程序中，即可分别获得上述5个数学模型的拟合参数及其校正决定系数和标准差；
C、根据5个数学模型统计结果中的校正决定系数、标准差，找出5个数学模型中校正决定系数值最大和标准差最小的一个数学模型，建立该地震工程地质单元的土层剪切波速和埋深关系的数学模型，并同时获得该数学模型的拟合参数；即5个数学模型中的哪一个数学模型中的校正决定系数值最大和标准差最小，该数学模型及其拟合参数越好，即越与该地震工程地质单元内的土层剪切波速数据相吻合；
(3)、将上述(2)中的C步骤建立的该地震工程地质单元的土层剪切波速和埋深关系的数学模型的拟合参数代入与该建立的数学模型相对应的土层等效剪切波速计算式中，得该地震工程地质单元的土层等效剪切波速；
与土层剪切波速和埋深关系的数学模型相对应的土层等效剪切波速计算式推导：
根据高等数学中的积分中值定理，当函数v(h)在闭区间[h1,h2]上连续，则在积分区间[h1,h2]上至少存在一个点ζ，使下式成立：
式(8)中，h为土层埋深，v为土层埋深h处的土层剪切波速，v(h)为土层剪切波速和埋深关系的数学模型即函数表达式，为土层剪切波速在埋深区间[h1,h2]上的平均值，h1、h2分别为埋深区间的两个端点，ζ为埋深区间[h1,h2]上的一点即h1≤ζ≤h2；
将式(3)、式(4)、式(5)、式(6)、式(7)分别代入式(8)，化简后，得埋深区间[h1,h2]上的任意深度段即任意埋深区间[H1,H2]上土层剪切波速平均值的计算式分别为:
<mrow>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mi>a</mi>
<mn>2</mn>
</mfrac>
<mi>A</mi>
<mo>+</mo>
<mi>b</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mi>c</mi>
<mn>3</mn>
</mfrac>
<mi>B</mi>
<mo>+</mo>
<mfrac>
<mi>e</mi>
<mn>2</mn>
</mfrac>
<mi>A</mi>
<mo>+</mo>
<mi>f</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mi>j</mi>
<mn>4</mn>
</mfrac>
<mi>C</mi>
<mo>+</mo>
<mfrac>
<mi>k</mi>
<mn>3</mn>
</mfrac>
<mi>B</mi>
<mo>+</mo>
<mfrac>
<mi>l</mi>
<mn>2</mn>
</mfrac>
<mi>A</mi>
<mo>+</mo>
<mi>p</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mi>q</mi>
<mn>5</mn>
</mfrac>
<mi>D</mi>
<mo>+</mo>
<mfrac>
<mi>r</mi>
<mn>4</mn>
</mfrac>
<mi>C</mi>
<mo>+</mo>
<mfrac>
<mi>u</mi>
<mn>3</mn>
</mfrac>
<mi>B</mi>
<mo>+</mo>
<mfrac>
<mi>w</mi>
<mn>2</mn>
</mfrac>
<mi>A</mi>
<mo>+</mo>
<mi>x</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>H</mi>
<mn>2</mn>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>H</mi>
<mn>1</mn>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>H</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>H</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
式(9)-式(13)中，为在埋深区间[h1,h2]上的土层剪切波速平均值，a、b、c、e、f、j、k、l、p、q、r、u、w、x、g、i分别为式(3)-式(7)的拟合参数，h1、h2分别为式(3)-式(7)的埋深区间[h1,h2]的两个端点，单位m，H1、H2分别为埋深区间[h1,h2]上的任意深度段即任意埋深区间[H1,H2]的两个端点，单位m，H1≥h1，H1＜H2≤h2，A＝H1+H2，
式(9)、式(10)、式(11)、式(12)、式(13)分别为与土层剪切波速和埋深关系的线性函数模型、一元二次函数模型、一元三次函数模型、一元四次函数模型、幂函数模型相对应的埋深区间[h1,h2]上的任意埋深区间[H1,H2]上土层剪切波速平均值计算式；
依据我国现行建筑抗震设计规范GB 50011-2010，土层等效剪切波速是地面至计算深度d0处的土层剪切波速平均值，由于地面的土层埋深为0.0m，计算深度d0处的土层埋深为计算深度d0，因此，取H1为0.0m，H2为计算深度d0，分别代入式(9)、式(10)、式(11)、式(12)、式(13)，化简后，得与土层剪切波速和埋深关系的线性函数模型、一元二次函数模型、一元三次函数模型、一元四次函数模型、幂函数模型相对应的土层等效剪切波速vse计算式，分别为:
式(14)-式(18)中，vse为土层等效剪切波速，d0为计算深度，取该地震工程地质单元内的场地覆盖层厚度值和20m两者的小值，a、b、c、e、f、j、k、l、p、q、r、u、w、x、g、i分别为式(3)-式(7)的拟合参数。
2.一种土层反射界面埋深的测量方法，其特征在于，具体步骤如下：
(1)、根据场地地震工程地质条件划分场地的地震工程地质单元，方法是：
A、收集场地已有地震工程地质条件资料，或对场地进行地震工程地质条件调查、勘察或钻孔剪切波速测试，获得场地地震工程地质条件资料，所述的场地地震工程地质条件资料为场地地貌单元类型、场地地下水类型与单孔地下水水位埋深、单孔场地覆盖层厚度、不同深度的水平方向土层地质成因类型、水平方向土层整体结构和土层剪切波速数据；所述的地震工程地质条件调查为携带场地地形图、地质罗盘、放大镜、GPS定位仪和测距仪，在场地及其邻区开展地形、地貌、地表地层岩性、地质成因类型、地质构造、水井、泉的观察、测量和记录；所述的地震工程地质条件勘察为在场地上进行钻孔、采集钻孔内土层样品，对土层样品进行观察、记录和化验，并对该钻孔内单孔地下水水位埋深数值进行测量和记录；所述的钻孔剪切波速测试为在场地内进行钻孔后，将波速测试仪的检波器安置在钻孔内各个埋深处即各个测点处并紧贴孔壁进行剪切波速测试，钻孔内相邻测点的垂直间距为1-2m从而获得钻孔内各个埋深处的土层剪切波速数据；将地面至土层剪切波速大于500m·s-1且其下方土层剪切波速均不小于500m·s-1的土层顶面的距离确定为单孔场地覆盖层厚度，而场地覆盖层厚度由场地内多个单孔场地覆盖层厚度确定，在场地内同一个地震工程地质单元范围的场地覆盖层厚度取该地震工程地质单元内单孔场地覆盖层厚度的平均值；
B、将场地地貌单元类型相同、场地地下水类型相同、水平方向土层地质成因类型相同、不同深度的水平方向土层整体结构相近、单孔地下水水位埋深数值相差不大于3m、单孔场地覆盖层厚度数值相差不大于5m的区域划分为同一个地震工程地质单元；将一个场地按照其地震工程地质条件划分为一个地震工程地质单元或多个地震工程地质单元；所述水平方向土层整体结构相近指水平方向土层岩性名称相同、水平方向土层分布连续、同一土层层面埋深相差不大于5m；
(2)、在同一个地震工程地质单元内，对土层剪切波速数据进行统计，建立土层剪切波速和埋深关系的数学模型及其拟合参数，方法是：
A、当测量土层等效剪切波速时，埋深区间为地面至计算深度d0，其中，计算深度d0取该地震工程地质单元内的场地覆盖层厚度和20m两者的小值；当测量土层反射界面埋深时,埋深区间为地面至钻孔内进行土层剪切波速测试的土层最大埋深；
B、采用下列5个土层剪切波速和埋深关系的数学模型，利用Origin软件对该地震工程地质单元内的土层剪切波速数据进行统计，分别得到下列5个数学模型的拟合参数及其校正决定系数、标准差：
土层剪切波速和埋深关系的数学模型为:
v＝ah+b h∈[h1,h2] 式(3)；
v＝ch2+eh+f h∈[h1,h2] 式(4)；
v＝jh3+kh2+lh+p h∈[h1,h2] 式(5)；
v＝qh4+rh3+uh2+wh+x h∈[h1,h2] 式(6)；
v＝ghi h∈[h1,h2] 式(7)；
式(3)为线性函数模型的表达式，式(4)为一元二次函数模型的表达式，式(5)为一元三次函数模型的表达式、式(6)为一元四次函数模型的表达式、式(7)为幂函数模型的表达式；
式(3)-式(7)中：h为土层埋深，v为土层埋深h处的土层剪切波速，a、b、c、e、f、j、k、l、p、q、r、u、w、x、g、i为拟合参数，由土层剪切波速数据统计后获得，h1、h2分别为埋深区间[h1,h2]的两个端点即统计的土层埋深的上界限值和下限界值；
上述数学模型是以地面为土层埋深h的零点，以垂直向下方向为土层埋深h的正方向；
上述校正决定系数和标准差是用来评价统计结果优劣的主要指标；
将该地震工程地质单元内的土层剪切波速数据及其对应的土层埋深数据作为输入数据，分别输入到Origin软件的上述5个数学模型的统计计算程序中，即可分别获得上述5个数学模型的拟合参数及其校正决定系数和标准差；
C、根据5个数学模型统计结果中的校正决定系数、标准差，找出5个数学模型中校正决定系数值最大和标准差最小的一个数学模型，建立该地震工程地质单元的土层剪切波速和埋深关系的数学模型，并同时获得该数学模型的拟合参数；即5个数学模型中的哪一个数学模型中的校正决定系数值最大和标准差最小，该数学模型及其拟合参数越好，即越与该地震工程地质单元内的土层剪切波速数据相吻合；
(4)、将上述(2)中的C步骤建立的该地震工程地质单元的土层剪切波速和埋深关系的数学模型的拟合参数代入与该建立的数学模型相对应的土层反射界面埋深计算式中，得到土层反射界面埋深；
与土层剪切波速和埋深关系的数学模型相对应的土层反射界面埋深计算式推导：
在横波反射法地震勘探中，通过反射法地震勘探数据处理后，得到地震反射时间剖面，把地震反射时间剖面转换为地质剖面时，要进行时深转换，其计算式为：
式(2)中，H0为土层剪切波反射波组反射界面埋深,简称土层反射界面埋深,v0为剪切波自地面至土层剪切波反射波组反射界面处的土层剪切波速平均值，t0为剪切波反射波组反射界面的双程垂直到时；
在横波反射法地震勘探中，设土层反射界面埋深H0位于上述土层剪切波速和埋深关系的埋深区间[h1,h2]上的任意埋深区间[H1,H2]上，由于地面的土层埋深为0.0m，土层剪切波反射波组反射界面处的土层埋深为土层反射界面埋深H0，因此，取H1为0.0m，取H2为H0，分别代入式(9)、式(10)、式(11)、式(12)、式(13)，化简后，得与土层剪切波速和埋深关系的线性函数模型、一元二次函数模型、一元三次函数模型、一元四次函数模型、幂函数模型相对应的剪切波自地面至土层剪切波反射波组反射界面处的土层剪切波速平均值v0，计算式分别为：
式(19)-式(23)中，v0为剪切波自地面至土层剪切波反射波组反射界面埋深H0处的土层剪切波速平均值，H0为土层反射界面埋深，a、b、c、e、f、j、k、l、p、q、r、u、w、x、g、i分别为式(3)-式(7)的拟合参数；
将式(19)、式(20)、式(21)、式(22)、式(23)分别代入式(2)，化简后,得与土层剪切波速和埋深关系的线性函数模型、一元二次函数模型、一元三次函数模型、一元四次函数模型、幂函数模型相对应的土层反射界面埋深H0的计算式分别为：
式(24)、式(25)、式(28)、式(26)、式(27)中，H0为土层反射界面埋深；t0为剪切波反射波组反射界面的双程垂直到时，a、b、c、e、f、j、k、l、p、q、r、u、w、x、g、i分别为式(3)-式(7)的拟合参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610395478.8A CN106094022B (zh) | 2016-06-06 | 2016-06-06 | 一种土层等效剪切波速和土层反射界面埋深的测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610395478.8A CN106094022B (zh) | 2016-06-06 | 2016-06-06 | 一种土层等效剪切波速和土层反射界面埋深的测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106094022A CN106094022A (zh) | 2016-11-09 |
| CN106094022B true CN106094022B (zh) | 2018-03-27 |
Family
ID=57448339
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610395478.8A Active CN106094022B (zh) | 2016-06-06 | 2016-06-06 | 一种土层等效剪切波速和土层反射界面埋深的测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106094022B (zh) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107561590B (zh) * | 2017-08-11 | 2019-04-16 | 武汉中岩科技股份有限公司 | 一种土层剪切波测试装置及其方法 |
| CN107513991A (zh) * | 2017-09-01 | 2017-12-26 | 苏州市轨道交通集团有限公司 | 预测振点液化半径及三角形振点布置下振点间距的方法 |
| CN107957593B (zh) * | 2017-12-19 | 2019-07-02 | 中国民航大学 | 一种厚层地下冰退化监测系统及控制评价方法 |
| CN109763474B (zh) * | 2018-12-07 | 2021-06-15 | 东南大学 | 一种基于标准贯入试验的剪切波速估算方法 |
| CN110119587B (zh) * | 2019-05-21 | 2023-07-07 | 深圳市秉睦科技有限公司 | 一种地层分层方法及装置 |
| CN111622201B (zh) * | 2020-05-09 | 2021-07-02 | 大连理工大学 | 一种基于面波频散特征的土层等效剪切波速直接计算方法 |
| CN113267814B (zh) * | 2021-07-05 | 2022-07-01 | 中国地震局地球物理研究所 | 场地剪切波速测算方法及装置 |
| CN115186037B (zh) * | 2022-06-09 | 2023-06-13 | 中国地震局地球物理研究所 | 地图显示方法、装置、电子设备和计算机可读介质 |
| CN115147022B (zh) * | 2022-09-05 | 2022-12-02 | 中国科学院地质与地球物理研究所 | 一种复杂山区长大线状工程的工程地质区划方法及系统 |
| CN115437009B (zh) * | 2022-11-04 | 2023-03-24 | 中国地震局地球物理研究所 | 钻孔数据处理方法、装置、电子设备和计算机可读介质 |
| CN116972954B (zh) * | 2023-09-25 | 2023-12-19 | 山东省地震工程研究院 | 基于原位激振的岩土波速测量方法及装置 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2002491A1 (en) * | 1989-07-19 | 1991-01-19 | Kenneth E. Hanson | Method of shear wave velocity estimation |
| US6820010B1 (en) * | 2003-04-30 | 2004-11-16 | Conocophillips Company | Method for determining shear-wave velocity model for depth migration of mode-converted data |
| CN103424772B (zh) * | 2012-05-24 | 2016-09-14 | 中国石油化工股份有限公司 | 一种基于岩石物理的碳酸盐储层横波速度预测方法 |
-
2016
- 2016-06-06 CN CN201610395478.8A patent/CN106094022B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106094022A (zh) | 2016-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106094022B (zh) | 一种土层等效剪切波速和土层反射界面埋深的测量方法 | |
| KR101642951B1 (ko) | Ｇｉｓ 기반 실시간 지진피해 예측 방법 | |
| US10401489B2 (en) | Method for drilling and characterizing subsurface porosity | |
| Galloway et al. | Analysis and simulation of regional subsidence accompanying groundwater abstraction and compaction of susceptible aquifer systems in the USA | |
| Gomberg et al. | Lithology and shear-wave velocity in Memphis, Tennessee | |
| Sundell et al. | A probabilistic approach to soil layer and bedrock-level modeling for risk assessment of groundwater drawdown induced land subsidence | |
| CN103116186B (zh) | 一种小尺度非均质储集体的容积确定方法 | |
| Akingboye et al. | Characteristics and rippability conditions of near-surface lithologic units (Penang Island, Malaysia) derived from multimethod geotomographic models and geostatistics | |
| Tong et al. | Determination of shear wave velocity of Yangtze Delta sediments using seismic piezocone tests | |
| Howle et al. | Determination of specific yield and water-table changes using temporal microgravity surveys collected during the second injection, storage, and recovery test at Lancaster, Antelope Valley, California, November 1996 through April 1997 | |
| Holmøy | Significance of geologicalparameters for predicting waterleakage in hard rock tunnels | |
| Zainal et al. | Characterization of Landslide geometry using Seismic Refraction Tomography in the GayoLues, Indonesia | |
| CN114839678A (zh) | 一种基于地震波静力触探试验的复杂土层剪切波速计算方法 | |
| Barrows et al. | A high-precision gravity survey in the Delaware Basin of southeastern New Mexico | |
| Zhu et al. | How far will uncertainty of the subsurface limit the sustainability planning of the subsurface | |
| Andersson et al. | Testing the methodology for site descriptive modelling. Application for the Laxemar area | |
| Brown et al. | Contemporary uplift of the Diablo Plateau, west Texas, from leveling measurements | |
| Kargaranbafghi et al. | Application of seismic refraction and tomography methods to the assessment of the Sivaki dam site | |
| Singeisen | Initiation and mechanisms of rock slope failures during the 2016 Mw 7.8 Kaikōura earthquake, New Zealand | |
| Zhao et al. | Geological survey and drilling technology of karst land | |
| Carvalho et al. | A soil classification for seismic hazard assessment and mitigation of the Algarve | |
| Akingboye et al. | Supervised Statistical Optimization of Interfaced Geophysical–Geotechnical Datasets for Holistic Surface–Subsurface Modeling | |
| Shamet | Subsurface Characterization of Internally Eroded Soils in Karst Using the Cone Penetration Test | |
| Li et al. | The application of tunnel reflection tomography in tunnel geological advanced prediction based on regression analysis | |
| Mousa et al. | Microgravity Survey to Evaluate Earthquake Effects on the Derbendikhan Dam Site, Iraq |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |