CN114166399B - 一种用于反演坝面水下爆炸压力的方法 - Google Patents
一种用于反演坝面水下爆炸压力的方法 Download PDFInfo
- Publication number
- CN114166399B CN114166399B CN202111307888.XA CN202111307888A CN114166399B CN 114166399 B CN114166399 B CN 114166399B CN 202111307888 A CN202111307888 A CN 202111307888A CN 114166399 B CN114166399 B CN 114166399B
- Authority
- CN
- China
- Prior art keywords
- explosion
- dam
- pressure
- underwater
- inverting
- 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
- 238000004880 explosion Methods 0.000 title claims abstract description 94
- 238000000034 method Methods 0.000 title claims abstract description 19
- 230000035939 shock Effects 0.000 claims description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 7
- 238000005259 measurement Methods 0.000 abstract description 9
- 238000002474 experimental method Methods 0.000 description 4
- 238000009434 installation Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
- G01L5/14—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for measuring the force of explosions; for measuring the energy of projectiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L25/00—Testing or calibrating of apparatus for measuring force, torque, work, mechanical power, or mechanical efficiency
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
本发明提供一种用于反演坝面水下爆炸压力的方法，包括：将水下爆炸压力测量装置布置于水下测点及坝面测点处，以获取现场爆炸不同测点的爆炸压力；对测得的数据进行处理以获得水下测点爆炸压力时程曲线和爆炸压力峰值，并进行拟合获得衰减公式；根据衰减公式以及坝面测点的实测爆炸压力相关数据推算获得爆源与传感器的理论距离L1和传感器与坝面的理论距离L2；将L1及L2进行叠加作为爆源与坝面的理论距离，并根据衰减公式反演坝面入射压力理论值；假定反射波与入射波的衰减规律相同，根据实测反射超压并结合衰减公式反演坝面反射超压理论值。本发明能根据现场水下爆炸实测数据确定传感器布置的真实位置，校正测量误差，实现坝面水下爆炸压力的准确合理反演。
Description
技术领域
本发明属于水下爆炸的技术领域，具体涉及一种用于反演坝面水下爆炸压力的方法。
背景技术
坝面爆炸压力是估算坝面混凝土模态参数的重要参考依据，获取可靠的坝面爆炸压力的是准确估算混凝土模态参数的前提。现有坝面水下爆炸压力的确定方法主要是利用传感器等水下爆炸压力测量装置直接获取实测数据，然而现有水下爆炸冲击波传感器工作状态始终与坝面有一定距离，所获取的数据并非坝面的真实压力。同时由于现场实验条件限制和各种不确定因素，坝面传感器的布置难免产生安装误差。例如大坝迎水壁面与竖直方向存在一定小角度时，水下传感器产生的测量误差可能进一步扩大，这对准确确定坝面水下爆炸压力和坝面混凝土模态参数提出了挑战。因此，通过一定方法利用压力传感器实测数据反演坝面水下爆炸压力对准确评估坝面水下爆炸压力及估算混凝土模态参数具有重要意义。
发明内容
本发明的目的在于针对现有技术的不足之处，提供一种用于反演坝面水下爆炸压力的方法，该方法能根据现场水下爆炸实验实测数据，确定传感器布置的真实位置，校正其测量误差，实现坝面水下爆炸压力的准确合理反演。
为解决上述技术问题，本发明采用如下技术方案：
一种用于反演坝面水下爆炸压力的方法，包括如下步骤：
步骤1，将水下爆炸压力测量装置布置于水下测点及坝面测点处，以获取现场爆炸不同测点的爆炸压力；
步骤2，对步骤1测得的水下测点数据进行处理以获得水下测点爆炸压力时程曲线和爆炸压力峰值，并根据不同测点的爆炸压力峰值、装药当量和爆心距拟合获得水下爆炸冲击波的衰减公式；
步骤3，根据步骤2拟合的衰减公式以及坝面测点的实测爆炸压力相关数据推算获得爆源与传感器的理论距离L1和传感器与坝面的理论距离L2；
步骤4，将步骤3中爆源与传感器的距离L1及传感器与坝面的距离L2叠加作为爆源与坝面的理论距离，并根据步骤2拟合的衰减公式反演坝面入射压力理论值；
步骤5，假定反射波与入射波的衰减规律相同，根据实测反射超压并结合步骤2拟合的衰减公式反演坝面反射超压理论值。
进一步地，步骤2中，根据Cole公式：P＝K(Q1/3/R)α拟合水下爆炸冲击波的衰减规律，其中P为爆炸压力峰值，Q为装药当量，R为爆心距，K、α为相关系数，在拟合时，根据不同测点的爆炸压力峰值、装药当量和爆心距确定K、α值从而获得衰减公式。
进一步地，爆源与传感器的理论距离其中P1为实测入射压力峰值；
传感器与坝面的理论距离其中c为冲击波在水中的传播速度，t为实测入、反射压力峰值的时间差。
进一步地，根据步骤2拟合的衰减公式反演坝面入射压力理论值
进一步地，根据步骤2拟合的衰减公式反演坝面反射超压理论值 其中，P2为实测反射超压。
与现有技术相比，本发明的有益效果为：本发明能够在水下爆炸压力测量装置实测数据的基础上，确定传感器的真实安装位置，对实测数据进行校正，可反演出更合理可靠的坝面水下爆炸压力；此外，本发明实施简便，不用对测量装置的安装精度提出较高要求，为受现场实验条件的限制时反演合理可靠的坝面水下爆炸压力提供保障。
附图说明
图1为本发明实施例一种用于反演坝面水下爆炸压力的方法流程图；
图2为本发明实施例中某水下爆炸实验现场布置的示意图；
图3为本发明实施例中坝面水下爆炸压力时程曲线示意图。
具体实施方式
下面将结合本发明实施例对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
下面结合具体实施例对本发明作进一步说明，但不作为本发明的限定。
如图1所示，本发明提供一种用于反演坝面水下爆炸压力的方法，包括以下步骤：
步骤1，将水下爆炸压力测量装置布置于水下测点及坝面测点以获取现场爆炸实验不同测点的爆炸压力；以某现场实验为例，其中6个水下测点垂直于坝面布置在水深10m处，4个水下测点平行于坝面布置在不同水深；7个坝面测点沿水平方向布置于水深10m处，4个坝面测点沿垂直方向布置于不同水深；其中部分布置见图2，图2示意了沿垂直方向布置的3个坝面测点；
步骤2，对步骤1测得的水下测点数据进行处理以获得水下测点爆炸压力时程曲线和爆炸压力峰值，并根据不同测点的爆炸压力峰值、装药当量和爆心距拟合获得水下爆炸冲击波的衰减公式；
在该步骤中，筛选并整理步骤1得到的水下测点数据以获得水下测点爆炸压力时程曲线，在本实施例中得到如图3所示的曲线，并从爆炸压力时程曲线上获取爆炸压力峰值，且根据Cole公式：P＝K(Q1/3/R)α拟合水下爆炸冲击波的衰减规律，其中，P为爆炸压力峰值，Q为装药当量，R为爆心距，K、α为相关系数，在确定爆炸压力峰值、装药当量和爆心距后即可确定K、α值，从而获得水下爆炸冲击波的衰减公式；在本实施例中，拟合得到的水下爆炸冲击波衰减公式为P＝53.41(Q1/3/R)1.14；
步骤3，根据步骤2拟合的衰减公式以及坝面测点的实测爆炸压力相关数据推算获得爆源与传感器的理论距离L1和传感器与坝面的理论距离L2；
在本实施例中，首先整理分析坝面测点爆炸压力时程曲线的特点，爆炸压力时程曲线呈现如图3所示双峰值的特点，结合现场实验条件，坝面与竖直方向存在9°夹角，推测出现有传感器无法准确贴合坝面以及其安装误差是产生该现象的原因；根据误差产生的原因利用传感器实测入射压力峰值并结合步骤2拟合的衰减公式推算爆源与传感器的理论距离L1。
其中，Q为装药当量，P1为实测入射压力峰值。
再利用实测入、反射压力峰值的时间差t推算传感器与坝面的理论距离L2为：
其中c为冲击波在水中的传播速度，取1480m/s。
在本实施例中，如图2所示，以BP13测点为例，装药当量Q＝20kg，入射压力实测值P1为7.30MPa，实测入、反射压力峰值的时间差t为0.61ms，可求得爆源与传感器的理论距离L1＝15.6m，传感器与坝面的理论距离L2＝0.5m，爆源与坝面的实际L＝18.3m。
步骤4，将步骤3中爆源与传感器的距离L1及传感器与坝面的距离L2叠加作为爆源与坝面的理论距离，并根据步骤2拟合的衰减公式反演坝面入射压力理论值；
在该步骤中，设定爆源与坝面的理论距离为步骤3中爆源与传感器的距离L1及传感器与坝面的距离L2的叠加，根据拟合的衰减公式反演坝面入射压力理论值为：
在本实施例中，根据上述数据可求得测点BP13的坝面入射压力理论值P1’＝7.07MPa。
步骤5，假定反射波与入射波的衰减规律相同，根据实测反射超压并结合步骤2拟合的衰减公式反演坝面反射超压理论值。
在该步骤中，假定反射波与入射波衰减规律相同，即其K、α值相等，利用实测反射超压，根据衰减公式反演坝面反射超压理论值为
测点BP13的反射超压实测值P2＝3.60MPa，可求得坝面反射超压理论值P2’＝3.66MPa。
由于传感器的真实位置与坝面存在较大距离，现根据实验实测数据求得坝面反射系数为49.3％，而根据本发明的方法反演的水下爆炸压力数据求得坝面的反射系数为51.8％，与数值模拟结果更为接近，这表明本发明提出的水下爆炸压力反演方法可以有效反演合理可靠的坝面水下爆炸压力数据。
以上仅为本发明较佳的实施例，并非因此限制本发明的实施方式及保护范围，对于本领域技术人员而言，应当能够意识到凡运用本发明说明书内容所作出的等同替换和显而易见的变化所得到的方案，均应当包含在本发明的保护范围内。
Claims (5)
1.一种用于反演坝面水下爆炸压力的方法，其特征在于，包括如下步骤：
步骤1，将水下爆炸压力测量装置布置于水下测点及坝面测点处，以获取现场爆炸不同测点的爆炸压力；
步骤2，对步骤1测得的水下测点数据进行处理以获得水下测点爆炸压力时程曲线和爆炸压力峰值，并根据不同测点的爆炸压力峰值、装药当量和爆心距拟合获得水下爆炸冲击波的衰减公式；
步骤3，根据步骤2拟合的衰减公式以及坝面测点的实测爆炸压力相关数据推算获得爆源与传感器的理论距离L1和传感器与坝面的理论距离L2；
步骤4，将步骤3中爆源与传感器的距离L1及传感器与坝面的距离L2叠加作为爆源与坝面的理论距离，并根据步骤2拟合的衰减公式反演坝面入射压力理论值；
步骤5，假定反射波与入射波的衰减规律相同，根据实测反射超压并结合步骤2拟合的衰减公式反演坝面反射超压理论值。
2.根据权利要求1所述的一种用于反演坝面水下爆炸压力的方法，其特征在于，步骤2中，根据Cole公式：P＝K(Q1/3/R)α拟合水下爆炸冲击波的衰减规律，其中P为爆炸压力峰值，Q为装药当量，R为爆心距，K、α为相关系数，在拟合时，根据不同测点的爆炸压力峰值、装药当量和爆心距确定K、α值从而获得衰减公式。
3.根据权利要求2所述的一种用于反演坝面水下爆炸压力的方法，其特征在于，步骤3中，爆源与传感器的理论距离其中P1为实测入射压力峰值；
传感器与坝面的理论距离其中c为冲击波在水中的传播速度，t为实测入、反射压力峰值的时间差。
4.根据权利要求2所述的一种用于反演坝面水下爆炸压力的方法，其特征在于，步骤4中，根据步骤2拟合的衰减公式反演坝面入射压力理论值
5.根据权利要求2所述的一种用于反演坝面水下爆炸压力的方法，其特征在于，步骤5中，根据步骤2拟合的衰减公式反演坝面反射超压理论值其中，P2为实测反射超压。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111307888.XA CN114166399B (zh) | 2021-11-05 | 2021-11-05 | 一种用于反演坝面水下爆炸压力的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111307888.XA CN114166399B (zh) | 2021-11-05 | 2021-11-05 | 一种用于反演坝面水下爆炸压力的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114166399A CN114166399A (zh) | 2022-03-11 |
| CN114166399B true CN114166399B (zh) | 2023-11-10 |
Family
ID=80478106
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111307888.XA Active CN114166399B (zh) | 2021-11-05 | 2021-11-05 | 一种用于反演坝面水下爆炸压力的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114166399B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114741863B (zh) * | 2022-03-31 | 2024-08-16 | 西北核技术研究所 | 一种城市街道中偶然爆炸冲击波超压峰值的快速计算方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2009143947A (ru) * | 2009-11-27 | 2011-06-10 | Российская Федерация, от имени которой выступает Министерство промышленности и торговли Российской Федерации (Минпромторг России) (RU) | Способ измерения давления |
| CN102967189A (zh) * | 2012-11-22 | 2013-03-13 | 中北大学 | 爆炸冲击波超压时空场重建方法 |
| CN108828008A (zh) * | 2018-04-19 | 2018-11-16 | 西安近代化学研究所 | 一种基于冲击波超压和准静压的炸药内爆炸威力评价方法 |
| CN109767471A (zh) * | 2019-01-15 | 2019-05-17 | 北京理工大学 | 一种动态爆心定位方法及系统 |
| WO2020206865A1 (zh) * | 2019-04-11 | 2020-10-15 | 中国矿业大学(北京) | 一种水中爆炸压力测试实验系统及方法 |
| CN112270140A (zh) * | 2020-09-29 | 2021-01-26 | 北京遥测技术研究所 | 一种用于爆炸冲击波超压预测的远场反演近场分布的方法 |
| CN112989726A (zh) * | 2021-04-26 | 2021-06-18 | 中国船舶科学研究中心 | 一种不同能量结构装药的水下爆炸载荷模型构造方法 |
-
2021
- 2021-11-05 CN CN202111307888.XA patent/CN114166399B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2009143947A (ru) * | 2009-11-27 | 2011-06-10 | Российская Федерация, от имени которой выступает Министерство промышленности и торговли Российской Федерации (Минпромторг России) (RU) | Способ измерения давления |
| CN102967189A (zh) * | 2012-11-22 | 2013-03-13 | 中北大学 | 爆炸冲击波超压时空场重建方法 |
| CN108828008A (zh) * | 2018-04-19 | 2018-11-16 | 西安近代化学研究所 | 一种基于冲击波超压和准静压的炸药内爆炸威力评价方法 |
| CN109767471A (zh) * | 2019-01-15 | 2019-05-17 | 北京理工大学 | 一种动态爆心定位方法及系统 |
| WO2020206865A1 (zh) * | 2019-04-11 | 2020-10-15 | 中国矿业大学(北京) | 一种水中爆炸压力测试实验系统及方法 |
| CN112270140A (zh) * | 2020-09-29 | 2021-01-26 | 北京遥测技术研究所 | 一种用于爆炸冲击波超压预测的远场反演近场分布的方法 |
| CN112989726A (zh) * | 2021-04-26 | 2021-06-18 | 中国船舶科学研究中心 | 一种不同能量结构装药的水下爆炸载荷模型构造方法 |
Non-Patent Citations (1)
| Title |
|---|
| 爆轰产物状态方程的水下爆炸反演理论研究;杨晨琛;李晓杰;闫鸿浩;王小红;王宇新;爆炸与冲击(009);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114166399A (zh) | 2022-03-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107132520B (zh) | 一种基于水声超短基线定位系统的声线修正方法及系统 | |
| Rennie et al. | Measurement of bed load velocity using an acoustic Doppler current profiler | |
| CN114397643B (zh) | 一种基于超短基线水声定位系统的声线修正方法 | |
| CN114166399B (zh) | 一种用于反演坝面水下爆炸压力的方法 | |
| CN109696234B (zh) | 一种发射点与接收点之间水平距离的确定方法及系统 | |
| CN111142071B (zh) | 一种结合半经验公式的单阵元匹配场爆炸声源定位方法 | |
| CN105444939A (zh) | 一种波浪砰击载荷测量装置及方法 | |
| CN105004413A (zh) | 用于水下目标定位的声传播路径综合速度测定方法与装置 | |
| CN110780263A (zh) | 一种基于卡西尼卵形线的多基地声纳系统定位精度分析方法 | |
| CN110132281A (zh) | 一种基于询问应答模式的水下高速目标高精度自主声学导航方法 | |
| Zheng et al. | Study on the calibration method of USBL system based on ray tracing | |
| RU2550576C1 (ru) | Способ измерения дистанции до шумящего объекта | |
| CN105718666B (zh) | 一种浪-流耦合海面的雷达后向散射截面数值模拟方法 | |
| CN114563778A (zh) | 基于时延冗余测量的阵形校准方法 | |
| CN109767471B (zh) | 一种动态爆心定位方法及系统 | |
| CN112763582B (zh) | 一种软岩崩解实时监测装置及监测方法 | |
| CN110412312B (zh) | 涌潮流速垂向分布的一种计算方法 | |
| CN114608567B (zh) | 一种小俯仰角条件下的usbl定位方法 | |
| Fayek et al. | A least square optimization approach for determining the soil boundary and absolute volume of unsaturated soils | |
| CN210268565U (zh) | 一种管道蠕变测量系统 | |
| CN112945034A (zh) | 基于水下冲击波实测波形的电子雷管延时精度测试方法 | |
| CN114578332B (zh) | 一种主动声呐目标径向尺度宽容性估计方法 | |
| CN112098953A (zh) | 一种大气折射修正量计算快速迭代方法及装置 | |
| CN110346783A (zh) | 基于幂级数等效测距模型的声线修正方法 | |
| CN117092588B (zh) | 一种估算水声定位系统时钟偏差方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |