CN105676272B - 炸药震源安全激发参数确定方法 - Google Patents
炸药震源安全激发参数确定方法 Download PDFInfo
- Publication number
- CN105676272B CN105676272B CN201610034486.XA CN201610034486A CN105676272B CN 105676272 B CN105676272 B CN 105676272B CN 201610034486 A CN201610034486 A CN 201610034486A CN 105676272 B CN105676272 B CN 105676272B
- Authority
- CN
- China
- Prior art keywords
- shock attenuation
- big gun
- safe
- shooting parameter
- parameter
- 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
- 238000000034 method Methods 0.000 title claims abstract description 44
- 239000010410 layer Substances 0.000 claims abstract description 41
- 238000012544 monitoring process Methods 0.000 claims abstract description 24
- 230000035939 shock Effects 0.000 claims abstract description 22
- 239000002689 soil Substances 0.000 claims abstract description 18
- 239000002344 surface layer Substances 0.000 claims abstract description 15
- 239000011159 matrix material Substances 0.000 claims abstract description 11
- 238000011835 investigation Methods 0.000 claims abstract description 8
- 239000000284 extract Substances 0.000 claims abstract description 5
- 239000002245 particle Substances 0.000 claims description 15
- 230000005284 excitation Effects 0.000 claims description 11
- 239000002360 explosive Substances 0.000 claims description 7
- 230000003068 static effect Effects 0.000 claims description 7
- 238000005422 blasting Methods 0.000 claims description 4
- 239000011435 rock Substances 0.000 claims description 4
- 239000003814 drug Substances 0.000 claims description 3
- 229940079593 drug Drugs 0.000 claims description 3
- 238000009835 boiling Methods 0.000 claims description 2
- 230000009514 concussion Effects 0.000 claims description 2
- 238000007689 inspection Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 238000000605 extraction Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 4
- 230000003993 interaction Effects 0.000 description 4
- 239000004927 clay Substances 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000005428 wave function Effects 0.000 description 2
- 239000000571 coke Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000013312 flour Nutrition 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000012360 testing method 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/02—Generating seismic energy
- G01V1/104—Generating seismic energy using explosive charges
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/02—Generating seismic energy
- G01V1/04—Details
-
- 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/30—Analysis
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)
Abstract
本发明提供了一种确定多层介质中,炸药震源安全激发参数的方法。该方法包括:(a)采用近地表调查方法初步划定地表层岩土分层层位;(b)针对保护的主要目标选取安全小药量试炮,试炮时采用多通道强震监测仪采集主要保护目标周围多点地表震动数据并提取最大震动速度值;(c)根据初步划定的地表岩土层位、试炮点与震动监测点的位置关系和监测提取的最大震动速度,建立地层震动衰减参数反演矩阵;(d)求解(c)建立的震动衰减参数矩阵,得到地表层各层震动衰减关系;(e)根据需保护的主要目标的震动安全阈值,基于地表层震动衰减关系,对不同药量在不同层位激发震动预测分析,确定较优的安全激发参数。本发明提供的炸药震源安全激发参数确定方法速度快、适用范围广,而且成本低、简便易实施,是一种有效的炸药震源安全激发参数确定方法。
Description
技术领域
本发明涉及地球物理勘探领域,更特别地是涉及一种炸药震源安全激发参数确定方法。
背景技术
地球物理勘探的方法很多,如地震勘探、电磁勘探、重力勘探等。现阶段最主要的勘探方式是地震勘探,与其他勘探方法相比,地震勘探具有高精度、高分辨率、大探测深度等优点。
地震勘探首先需要人工方法激发地震波,激发地震波的效果至关重要,直接影响地震勘探的成像深度和精度。目前地震波的激发源可分为炸药震源和非炸药震源两种,原理大致相同,非炸药震源主要有重锤震源、气动震源、空气枪震源等。从19世纪20年代至今,炸药震源一直是地震勘探的主要常规震源,炸药产生的能量大,很容易激发出地震波,实践证明,在地震勘查、石油勘探中使用炸药震源具有良好的效果。
炸药震源的主要缺点是不安全,而且现在的勘探范围已由原来的人口稀少区域深入到人口密集的区域,地形条件复杂和施工难度也进一步加大,除了要达到较好的勘探效果,还要重点考虑安全因素,因此炸药震源的地震安全性激发,已成为各方面密切关注的焦点。但是目前安全激发参数的获取方法精度较低,无法满足石油勘探过程中,在人口密集地区进行激发的精度要求。对此本发明充分考虑了近地表多层条件对激发参数的影响,并提出了多层介质中高精度炸药震源安全激发参数确定方法。
发明内容
目的:本发明提供的炸药震源安全激发参数确定方法速度快、适用范围广,而且成本低、简便易实施,是一种有效的炸药震源安全激发参数确定方法。
本发明提供了一种炸药震源安全激发参数确定方法,该方法包括:
(a)采用近地表调查方法初步划定地表层岩土层位;
(b)针对保护的主要目标选取安全小药量试炮,试炮时采用多通道强震监测仪采集主要保护目标周围多点地表震动数据并提取最大震动速度值;
(c)根据初步划定的地表岩土层位、试炮点与震动监测点的位置关系和监测提取的最大震动速度,建立地层震动衰减参数反演矩阵;
(d)求解(c)建立的震动衰减参数矩阵,得到地表层各层震动衰减关系;
(e)根据需保护的主要目标的震动安全阈值,基于地表层震动衰减关系,对不同药量在不同层位激发震动预测分析,确定较优的安全激发参数。
步骤(a)包括以下步骤:
(a1)采用近地表调查方法可采用触探方法,包括静力初探和动力触探,亦可采用小折射方或微测井等方法。
(a2)通过近地表调查方法,将距离地面20~30米范围分成3~5个层段,确定每个层的波速
步骤(b)包括以下步骤:
(b1)根据需要保护的目标,根据经验公式,取确定地层岩土层段中经验衰减最小的参数估算试炮估计药量,在实际试炮时取0.6~0.8倍估计药量为实际试炮药量;
(b2)试炮时,在距离炮点0~50米范围内布置强震速度传感器,记录这些位置点处地表震动历史,并提取各个监测点的峰值质点速度(PPV);
(b3)试炮时,激发炮数和监测点数量由(a2)确定的地表层段数量相关,(激发炮数*监测点数)>(2*地表层段数);
(b4)试炮激发深度一般取在阻抗较大的岩土层内,取1~2个深度;
(b5)试炮时,(b2)监测点位置在0~50米范围内均匀布置,一般取在1~2条测线方向上,每个测点距离炮点距离不同,监测点一般取5~10个点;
步骤(c)包括以下步骤:
(c1)根据峰值质点速度(PPV)在不同层段衰减关系,依据炮点和监测点位置关系,建立经多层段衰减后监测点峰值指点速度关系,建立地层震动衰减参数反演矩阵;
步骤(d)包括以下步骤:
(d1)根据最小二乘法,求解地层震动衰减参数,确定各层段衰减参数;
步骤(e)包括以下步骤:
(e1)根据所要保护的目标,由爆破震动安全标准确定安全阈值;
(e2)根据岩层段参数,确定炮孔位置和震源炸药阻抗匹配的激发岩土层深度;
(e3)改变激发药量,由地表各层震动衰减关系,计算地表震动峰值质点速度,计算激发药量取为0.25~10kg,间隔0.25kg;
(e4)根据计算得到保护目标位置处峰值质点速度,与安全阈值相对比,小于安全阈值判定为安全,实际安全激发药量取为能保证在安全阈值内最大激发药量。
附图说明
图1是根据本发明的炸药震源安全激发参数确定方法的流程图;
图2静力触探得到的图;
图3试炮布置图;
图4是炮点和监测点地震波传播图
图5地震震动衰减模型进行分析得到的衰减参数关系示意图;
图6采用本发明的应用的炸药震源安全激发参数确定方法应用效果对比
具体实施方式
实例1:三层土介质结构安全激发参数确定方法:
步骤1:对目标地区单孔静力触探以及岩土识别图的分析可知(如图2),在该点处0~30m范围内主要的土的类型有两种:粉土(及粉粘土)和粉砂。对于粉土可细分为粉土、粉质粘土以及粉土与粉质粘土层三种类型。
根据静力触探结果,将距离地面30m厚的岩土分为3个层段。
针对3个层段的土介质,利用环刀法检测其密度;
利用公式确定每层介质的波速,其中弹性模量E可以利用材料库进行查询。
步骤2:根据《爆破安全规程》(6722-2014)查得土层的衰减系数为:k=250~350,α=1.8~2.0,选择其中衰减最小的系数,即k=350,α=2.0。
根据《爆破安全规程》(6722-2014)查得露天浅孔爆破产生地震波的主频在10-100Hz之间,其中一般民用建筑在主频10~50Hz范围内地震波作用下允许的最大振动速度为2~2.5cm/s
根据以上查得的内容,可以计算出为了保护50m距离外的目标建筑不受损坏,试炮估计药量为54.1kg,因此实际估计药量为32kg。
将炮点确定后,沿测线方向0~50m距离内均匀布设若干强震速度传感器。
激发深度选择在第三层岩土层。
在激发一发的条件下,计算三层介质的衰减系数,需要沿测线方向布设6个以上检波器,具体布设如图3。
步骤3:根据峰值质点速度(PPV)在不同层段衰减关系,依据炮点和监测点位置关系,建立经多层段衰减后监测点峰值指点速度关系,建立地层震动衰减参数反演矩阵如下;
y=a·b·x,其中y=lnv,a=lnk,b=α,
步骤4:利用最小二乘法进行线性回归后,可以确定三个不同岩土层的衰减系数
步骤5:根据《爆破安全规程》(6722-2014)查得工业和商业建筑物允许的质点最大振动速度为2.5-5cm/s
根据岩层段参数,确定炮孔位置和震源炸药阻抗匹配的激发岩土层深度
将初始激发药量设定为0.25kg,通过反演出的衰减参数计算目标建筑处地表质点振动峰值,如果没有超过安全阈值则继续增大药量,直到振动速度超过安全阈值为止。
选择恰好使目标建筑附近地表振动在安全阈值以内的最大药量为激发药量。
实例2:五层土介质结构安全激发参数确定方法:
步骤1:对目标地区单孔静力触探以及岩土识别,对于粉土可细分为粉土、粉质粘土以及粉土与粉质粘土层三种类型。
根据静力触探结果,将距离地面30m厚的岩土分为5个层段。
针对3个层段的土介质,利用环刀法检测其密度;
利用公式确定每层介质的波速,其中弹性模量E可以利用材料库进行查询。
步骤2:根据《爆破安全规程》(6722-2014)查得土层的衰减系数为:k=250~350,α=1.8~2.0,选择其中衰减最小的系数,即k=350,α=2.0。
根据《爆破安全规程》(6722-2014)查得露天浅孔爆破产生地震波的主频在10-100Hz之间,其中一般民用建筑在主频10~50Hz范围内地震波作用下允许的最大振动速度为2~2.5cm/s
根据以上查得的内容,可以计算出为了保护50m距离外的目标建筑不受损坏,试炮估计药量为54.1kg,因此实际估计药量为32kg。
将炮点确定后,沿测线方向0~50m距离内均匀布设若干强震速度传感器。
激发深度选择在第五层岩土层。
在激发一发的条件下,计算五层介质的衰减系数,需要沿测线方向布设10个以上检波器。
步骤3:根据峰值质点速度(PPV)在不同层段衰减关系,依据炮点和监测点位置关系,建立经多层段衰减后监测点峰值指点速度关系,建立地层震动衰减参数反演矩阵如下;
y=a·b·x,其中y=lnv,a=lnk,b=α,
步骤4:利用最小二乘法进行线性回归后,可以确定三个不同岩土层的衰减系数
步骤5:根据《爆破安全规程》(6722-2014)查得工业和商业建筑物允许的质点最大振动速度为2.5-5cm/s
根据岩层段参数,确定炮孔位置和震源炸药阻抗匹配的激发岩土层深度
将初始激发药量设定为0.25kg,通过反演出的衰减参数计算目标建筑处地表质点振动峰值,如果没有超过安全阈值则继续增大药量,直到振动速度超过安全阈值为止。
选择恰好使目标建筑附近地表振动在安全阈值以内的最大药量为激发药量。
Claims (6)
1.一种多层介质中炸药震源安全激发参数确定方法,其特征是:该方法包括,
(a)采用近地表调查方法初步划定地表层岩土层位;
(b)针对保护的主要目标选取安全小药量试炮,试炮时采用多通道强震监测仪采集主要保护目标周围多点地表震动数据并提取最大震动速度值;
(c)根据初步划定的地表层岩土层位、试炮点与震动监测点的位置关系和监测提取的最大震动速度值,建立地层震动衰减参数反演矩阵;
(d)求解(c)建立的地层震动衰减参数反演矩阵,得到地表层不同层位的震动衰减关系;
(e)根据需保护的主要目标的震动安全阈值,基于地表层不同层位震动衰减关系,对不同药量在不同层位激发震动预测分析,确定安全激发参数。
2.根据权利要求1所述炸药震源安全激发参数确定方法,其特征是:步骤(a)包括以下步骤:
(a1)采用近地表调查方法中的触探方法,包括静力触探和动力触探,或采用小折射或微测井方法;
(a2)通过近地表调查方法,将距离地面20~30米范围分成3~5个层段,确定每个层段的波速。
3.根据权利要求2所述炸药震源安全激发参数确定方法,其特征是:步骤(b)包括以下步骤:
(b1)根据需要保护的主要目标,根据经验公式,取确定地层岩土层段中经验衰减最小的参数估算试炮估计药量,在实际试炮时取0.6~0.8倍试炮估计药量为实际试炮药量;
(b2)试炮时,在距离炮点0~50米范围内布置强震速度传感器,记录这些位置点处地表震动历史,并提取各个监测点的峰值质点速度;
(b3)试炮时,激发炮数和监测点数量与(a2)确定的地表层段数量相关,(激发炮数*监测点数)>(2*地表层段数);
(b4)试炮激发深度取在阻抗较大的岩土层内,取1~2个深度;
(b5)试炮时,(b2)监测点位置在0~50米范围内均匀布置,取在1~2条测线方向上,每个测点距离炮点距离不同,监测点取5~10个点。
4.根据权利要求1所述炸药震源安全激发参数确定方法,其特征是:步骤(c)包括以下步骤:
(c1)根据峰值质点速度在不同层段衰减关系,依据炮点和监测点位置关系,建立经多层段衰减后监测点峰值质点速度关系,建立地层震动衰减参数反演矩阵。
5.根据权利要求1所述炸药震源安全激发参数确定方法,其特征是:步骤(d)包括以下步骤:
(d1)根据最小二乘法,求解地层震动衰减参数反演矩阵,确定各层段震动衰减参数。
6.根据权利要求1所述炸药震源安全激发参数确定方法,其特征是:步骤(e)包括以下步骤:
(e1)根据所要保护的主要目标,由爆破震动安全标准确定震动安全阈值;
(e2)根据岩层段参数,确定炮孔位置和震源炸药阻抗匹配的激发岩土层深度;
(e3)改变激发药量,由地表层不同层位的震动衰减关系,计算地表震动峰值质点速度,计算激发药量取为0.25~10kg,间隔0.25kg;
(e4)根据计算得到保护的主要目标位置处峰值质点速度,与震动安全阈值相对比,小于震动安全阈值判定为安全,实际震动安全激发药量取为能保证在震动安全阈值内最大激发药量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610034486.XA CN105676272B (zh) | 2016-01-19 | 2016-01-19 | 炸药震源安全激发参数确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610034486.XA CN105676272B (zh) | 2016-01-19 | 2016-01-19 | 炸药震源安全激发参数确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105676272A CN105676272A (zh) | 2016-06-15 |
CN105676272B true CN105676272B (zh) | 2018-11-06 |
Family
ID=56301505
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610034486.XA Active CN105676272B (zh) | 2016-01-19 | 2016-01-19 | 炸药震源安全激发参数确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105676272B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108732611A (zh) * | 2017-04-20 | 2018-11-02 | 中国石油化工股份有限公司 | 地震采集激发药量计算方法及计算机可读存储介质 |
CN107831531B (zh) * | 2017-11-15 | 2020-02-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 地震勘探炸药震源激发点的安全布设方法及判断方法 |
CN110082820B (zh) * | 2018-01-26 | 2020-11-10 | 中石化石油工程技术服务有限公司 | 炸药震源混合分布式宽频激发的方法 |
CN110749376B (zh) * | 2019-11-14 | 2022-08-05 | 北京龙德时代技术服务有限公司 | 基于连续数字测振的爆破安全分析方法及分析系统 |
CN112505756B (zh) * | 2020-10-30 | 2024-03-26 | 中国石油天然气集团有限公司 | 地震勘探野外炮点安全距离确定方法及装置 |
CN117970488B (zh) * | 2024-04-02 | 2024-05-31 | 山东省地震局 | 一种三维复杂地学模型中地震精确定位方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2012114039A (ru) * | 2012-04-10 | 2013-10-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Кузбасский государственный технический университет имени Т.Ф. Горбачева" (КузГТУ) | Способ определения радиуса сейсмоопасной зоны |
CN102819041B (zh) * | 2012-08-28 | 2014-12-10 | 中国石油化工股份有限公司 | 一种利用土层阻力信息确定地震炸药激发深度的方法 |
-
2016
- 2016-01-19 CN CN201610034486.XA patent/CN105676272B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105676272A (zh) | 2016-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105676272B (zh) | 炸药震源安全激发参数确定方法 | |
US20120243368A1 (en) | Fracture Imaging Methods Employing Skeletonization of Seismic Emission Tomography Data | |
CN105425281B (zh) | 分布式炸药震源激发参数确定方法 | |
CN104114811A (zh) | 定位诱导的地下岩层裂缝中的示踪支撑剂的与岩性和井眼状况无关的方法 | |
Wang et al. | Experimental investigation on the fracture behaviour of black shale by acoustic emission monitoring and CT image analysis during uniaxial compression | |
CN110018528B (zh) | 幔源油气综合地质分析方法 | |
CN106767371A (zh) | 一种探地雷达快速探测复垦土壤覆盖表土厚度的方法 | |
CN103678778A (zh) | 一种放射性物化探信息综合方法 | |
Zhang et al. | Sensitivity and regression analysis of acoustic parameters for determining physical properties of frozen fine sand with ultrasonic test | |
Colombero et al. | Multiscale geophysical characterization of an unstable rock mass | |
US9695684B2 (en) | System and method for predicting the front arrival time in reservoir seismic monitoring | |
Leucci et al. | Ground Penetrating Radar (GPR): an application for evaluating the state of maintenance of the building coating | |
RU2390805C1 (ru) | Способ контроля геометрических и гидродинамических параметров гидроразрыва пласта | |
Movassagh et al. | A comparison of fractal methods for evaluation of hydraulic fracturing surface roughness | |
Łój et al. | Geophysical surveys and modelling for recognizing of gypsum karst | |
Orlando et al. | Characterization of a pre‐Trajan wall by integrated geophysical methods | |
Binder et al. | Joint microseismic event location with surface geophones and downhole DAS at the FORGE geothermal site | |
EP2577359B1 (en) | Method and system for determining the rigidity of a geological layer | |
Riahi et al. | Seismic refraction and downhole survey for characterization of shallow depth materials of Bam city, southeast of Iran | |
CN117631024A (zh) | 一种灰岩地区炸药震源定量化精确激发设计方法 | |
Maksimov et al. | Prediction and Early Detection of Karsts—An Overview of Methods and Technologies for Safer Drilling in Carbonates | |
Cumerlato et al. | Application of refraction tomography to map the extent of blast-induced fracturing | |
Liu et al. | [Retracted] Application of Internet of Things Technology in the Construction of Bored Pile Foundation in a High Altitude Environment | |
Larson et al. | Surface-wave-based inversions of shallow seismic structure | |
EP2751601A1 (en) | Measurement of permeability in borehole in the presence of mudcake |
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 |