一种基于纵波波速预估高温后花岗岩围岩孔隙度的方法技术

本发明专利技术公开了一种基于纵波波速预估高温后花岗岩围岩孔隙度的方法，首先基于Boltzmann函数来建立孔隙度和弹性波纵波波速之间的关系，然后选取有代表性的测试点制作花岗岩岩样经行室内试验，测试弹性波纵波波速和孔隙度，最后通过现场测试花岗岩围岩的弹性波纵波波速，代入花岗岩围岩孔隙度和纵波波速的关系式中，即可预估围岩的孔隙度。本发明专利技术基于Boltzmann函数建立了花岗岩围岩孔隙度和弹性波纵波波速的关系公式，创造了通过测试花岗岩围岩弹性波纵波波速来预估围岩孔隙度的方法；在实际工程中，缩短了进行围岩孔隙度室内试验的过程，极大地节省了围岩孔隙度测试的时间。

【技术实现步骤摘要】
一种基于纵波波速预估高温后花岗岩围岩孔隙度的方法
本专利技术属于岩石物理研究
，涉及一种花岗岩围岩孔隙度预估方法，尤其是涉及一种基于纵波波速预估高温后花岗岩围岩孔隙度的方法。
技术介绍
当岩体工程，例如隧道工程，遭受火灾后，围岩的强度及稳定性大大下降，其原因是火灾引发的高温会诱导围岩开裂，造成岩石性能劣化。了解高温后围岩的开裂程度对于灾后岩体工程的安全性评估，加固和修复至观重要。孔隙度是描述围岩开裂程度的一个重要指标。目前，测量围岩孔隙度的方法主要有毛管压力曲线法(包括压汞法、半渗透隔板法、离心机法、气体吸附法等)、图像分析法(包括扫描电镜法、CT扫描法、铸体薄片法等)、三维孔隙结构模拟法(包括切片组合法、X射线衍射成像、薄片图像重建法等)、核磁共振法和测井法等。这些方法均是提取围岩岩心直接进行室内试验。在实际火灾后岩体工程安全性评估工作中，测试试样数量多，价格昂贵，测量周期长，不方便用于围岩整体孔隙度状态评估当中。围岩的弹性波纵波波速可以通过现场测试技术方便地测出。超声波波速测试已经是一种相当成熟的技术了，通常采用一发双收单孔测试法进行现场测试，原理及步骤为：首先在围岩断面上确定测试点，在测试点处钻孔；然后将圆管状声波探头置入钻孔，孔内注水以使探头与孔壁有良好的声耦合；采用一发双收的装置，包括一个超声波发射器，两个超声波接收器；发射器在钻孔中发射声波，在孔壁周围产生滑行波沿着钻孔壁传播。当首波传播到第一个接收器时，开始计时；当滑行波继续传播到第二个接收器时计时停止。声波在两个接收器之间的时间为Δt，两个接收器之间的距离为ΔL，根据公式：计算声波波速，然后逐点进行测试，直至各点测试完毕。通过建立一种高温后围岩的孔隙度与弹性波纵波波速的关系，用简单易测的纵波波速去预估现场不方便大量测量的孔隙度，对于提高灾后工程安全性评估工作效率，节约成本具有重要意义。
技术实现思路
本专利技术提出一种基于纵波波速预估高温后花岗岩围岩孔隙度的方法，采用数学模型来建立花岗岩围岩孔隙度和纵波波速的关系，通过测量火灾后花岗岩围岩的弹性波纵波波速，来预估其孔隙度。本方法主要用于火灾后岩体工程的安全性和稳定性评估工作中，有利于提高工作效率，节约成本。为了达到上述目的，本专利技术提供了一种基于纵波波速预估高温后花岗岩围岩孔隙度的方法，首先基于Boltzmann函数来建立孔隙度和弹性波纵波波速之间的关系，然后选取有代表性的测试点制作花岗岩岩样经行室内试验，测试弹性波纵波波速和孔隙度，最后通过现场测试花岗岩围岩的弹性波纵波波速，代入花岗岩围岩孔隙度和纵波波速的关系式中，即可预估围岩的孔隙度。该方法的具体实现过程如下：S1、选用Boltzmann函数建立高温后花岗岩围岩孔隙度与弹性波纵波波速的关系式；采用Boltzmann函数表达高温后花岗岩围岩孔隙度和弹性波纵波波速的关系，公式为其中，P为孔隙度，V为纵波波速；A1、A2、A3、A4分别为拟合系数；A1为拟合曲线下渐进线的纵坐标值；A2为拟合曲线上渐进线的纵坐标值；A3为拟合曲线拐点的横坐标值；A4为拟合曲线中间部分的斜率参数。S2试验测量点实施；S2.1.选取有代表性的测试点测试其弹性波纵波波速和孔隙度；花岗岩围岩测试点选取数量要大于五个，测试点之间纵波波速的差异性要较大；弹性波纵波波速通过超声波试验测得；孔隙度通过毛管压力曲线(包括压汞法、半渗透隔板法、离心机法、气体吸附法等)、图像分析(包括扫描电镜法、CT扫描法、铸体薄片法等)、三维孔隙结构模拟(包括切片组合法、X射线衍射成像、薄片图像重建法等)或核磁共振测得；S2.2.对室内试验测试得到的试验数据采用Boltzmann函数进行拟合，求得函数系数；将室内试验得到的测试点的数据进行拟合，得到拟合曲线；由拟合曲线的特征得出公式中的参数A1、A2、A3、A4；A1为拟合曲线下渐进线的纵坐标值；A2为拟合曲线上渐进线的纵坐标值；A3为拟合曲线拐点的横坐标值；A4为拟合曲线中间部分的斜率参数。确定函数系数后得到完整的花岗岩围岩孔隙度和弹性波纵波波速的关系公式，用于高温后的花岗岩围岩孔隙度的预估。S3.测试其他位置测试点的弹性波纵波波速，采用花岗岩围岩孔隙度和弹性波纵波波速关系公式对高温后的花岗岩围岩孔隙度预估；S3.1.测量火灾后花岗岩围岩的弹性波纵波波速；在测试点处钻孔，然后将超声波一发双收装置放入钻孔，钻孔内注水以使探头与孔壁进行声耦合，采用的一发双收装置包括一个超声波发射器和两个超声波接收器；超声波发射器在钻孔中发射声波，在围岩的孔壁周围产生滑行波沿着钻孔壁传播；当首波传播到第一个超声波接收器时，开始计时；当滑行波继续传播到第二个超声波接收器时，计时停止；两个超声波接收器接收超声波的时间差为Δt，两个超声波接收器之间的距离为ΔL；计算岩样弹性波纵波波速：然后逐个对各个测试点进行测试，直至各测试点测试完毕。S3.2.预估高温后的花岗岩围岩测试点的孔隙度；将测试点测得花岗岩围岩弹性波纵波波速V代入S2得出的花岗岩围岩孔隙度和弹性波纵波波速的关系式中，计算得出花岗岩围岩孔隙度P；由此得到的孔隙度数据为预估的高温后的花岗岩围岩孔隙度。与现有技术相比较，本专利技术基于Boltzmann函数建立了花岗岩围岩孔隙度和弹性波纵波波速的关系公式，创造了通过测试花岗岩围岩弹性波纵波波速来预估围岩孔隙度的方法。在实际工程中，通过测量围岩的弹性波纵波波速来预估围岩孔隙度，改变了围岩孔隙度不方便进行现场大量测量的局限性，缩短了进行围岩孔隙度室内试验的过程，极大地节省了围岩孔隙度测试的时间。通过进行现场超声波测试，得到花岗岩围岩的弹性波纵波波速，代入公式得到花岗岩围岩的孔隙度，利用花岗岩围岩的孔隙特征分析围岩的稳定性和安全性，进而为灾后工程的加固和修复提供参考依据。附图说明图1方法流程图图2现有的实验数据拟合后的结果。图3案例一结果对比图。图4案例二结果对比图。具体实施方式以下结合附图和实施例对本专利技术方法进行详细说明。一种基于纵波波速预估高温后花岗岩围岩孔隙度的方法，该方法的实施步骤如下：S1、选用Boltzmann函数采用Boltzmann函数可以很好地表达高温后花岗岩弹性波纵波波速和孔隙度的关系，附图1为对现有的实验数据采用Boltzmann函数拟合后的结果，表1为相关拟合系数，相关系数R2反应拟合的准确程度，由R2>0.9可以看出拟合程度高，采用Boltzmann函数是可行的。表1相关拟合系数基于Boltzmann函数的高温后花岗岩孔隙度与岩石弹性波纵波波速关系公式为其中P为孔隙度，V为纵波波速，A1、A2、A3、A4分别为拟合系数；A1为拟合曲线下渐进线的纵坐标值；A2为拟合曲线上渐进线的纵坐本文档来自技高网...

【技术保护点】
1.一种基于纵波波速预估高温后花岗岩围岩孔隙度的方法，其特征在于：首先基于Boltzmann函数来建立孔隙度和弹性波纵波波速之间的关系，然后选取有代表性的测试点制作花岗岩岩样经行室内试验，测试弹性波纵波波速和孔隙度，最后通过现场测试花岗岩围岩的弹性波纵波波速，代入花岗岩围岩孔隙度和纵波波速的关系式中，即可预估围岩的孔隙度；/n该方法的具体实现过程如下：/nS1、选用Boltzmann函数建立高温后花岗岩围岩孔隙度与弹性波纵波波速的关系式；/n采用Boltzmann函数表达高温后花岗岩围岩孔隙度和弹性波纵波波速的关系，公式为/n

【技术特征摘要】
1.一种基于纵波波速预估高温后花岗岩围岩孔隙度的方法，其特征在于：首先基于Boltzmann函数来建立孔隙度和弹性波纵波波速之间的关系，然后选取有代表性的测试点制作花岗岩岩样经行室内试验，测试弹性波纵波波速和孔隙度，最后通过现场测试花岗岩围岩的弹性波纵波波速，代入花岗岩围岩孔隙度和纵波波速的关系式中，即可预估围岩的孔隙度；
该方法的具体实现过程如下：
S1、选用Boltzmann函数建立高温后花岗岩围岩孔隙度与弹性波纵波波速的关系式；
采用Boltzmann函数表达高温后花岗岩围岩孔隙度和弹性波纵波波速的关系，公式为



其中，P为孔隙度，V为纵波波速；A1、A2、A3、A4分别为拟合系数；A1为拟合曲线下渐进线的纵坐标值；A2为拟合曲线上渐进线的纵坐标值；A3为拟合曲线拐点的横坐标值；A4为拟合曲线中间部分的斜率参数；
S2试验测量点实施；
S2.1.选取有代表性的测试点测试其弹性波纵波波速和孔隙度；
花岗岩围岩测试点选取数量要大于五个，测试点之间纵波波速的差异性要较大；
弹性波纵波波速通过超声波试验测得；
孔隙度通过毛管压力曲线、图像分析、三维孔隙结构模拟或核磁共振测得；
S2.2.对室内试验测试得到的试验数据采用Boltzmann函数进行拟合，求得函数系数；
将室内试验得到的测试点的数据进行拟合，得到拟合曲线；由拟合曲线的特征得出公式中的参数A1、A2、A3、A4；A1为拟合曲线下渐进线的纵坐标值；A2为拟合曲线上渐进线的纵坐标值；A3为拟合曲线拐点的横坐标值；A4为拟合曲线中间部分的斜率参数；确定函数系数后得到完整的花岗岩围岩孔隙度和弹性波纵波波速的关系公式，用于高温后的花岗岩围岩孔隙度的预估；
S3.测试其他...

【专利技术属性】
技术研发人员：范立峰，樊一丁，高经纬，杜修力，席岩，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：北京;11