CN112363209B - 一种断层三维数字化检测方法 - Google Patents
一种断层三维数字化检测方法 Download PDFInfo
- Publication number
- CN112363209B CN112363209B CN202011241490.6A CN202011241490A CN112363209B CN 112363209 B CN112363209 B CN 112363209B CN 202011241490 A CN202011241490 A CN 202011241490A CN 112363209 B CN112363209 B CN 112363209B
- Authority
- CN
- China
- Prior art keywords
- sphere
- center position
- fault
- standard
- correction
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 54
- 238000002474 experimental method Methods 0.000 claims abstract description 6
- 238000012937 correction Methods 0.000 claims description 51
- 238000000034 method Methods 0.000 claims description 37
- 238000005259 measurement Methods 0.000 claims description 30
- 238000012545 processing Methods 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 18
- 239000000523 sample Substances 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 7
- 239000010979 ruby Substances 0.000 claims description 6
- 229910001750 ruby Inorganic materials 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 3
- 238000004904 shortening Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000004441 surface measurement Methods 0.000 description 1
- 230000000007 visual effect Effects 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
-
- 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
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
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)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种断层三维数字化检测方法,利用DMIS技术建立地震超前观测模型,并采集和处理地震数据,开始利用探测仪器进行断层探测实验,通过PET彩色图像技术向断层表面发射彩色光源,将差值最小值所对应的起始轮廓区作为PET发射图像上构成断面的轮廓区并以此锁定断面的方位,与传统的断面定位方法相比,此方法能够利用PET技术对不同结构断层进行色度编码,以此对断层进行标记和与预设阀值比对,更加精确地知晓不同断层所在方位以及其结构轮廓,观察断层时也更加直观和清晰并判断断层位置,利用键入四点边界方法对断层进行自动测量,并加以记录,对测量结果进行数据处理并对结果进行自行比对。
Description
技术领域
本发明涉及地震勘探中地震数据解释技术领域,特别涉及一种断层三维数字化检测方法。
背景技术
断层表面由自由曲面构成,它的表面形态直接关系到面波勘探调查和近断层地震作用下的土质边坡动力响应等问题,所以精度要求高且数据信息量大,传统的三坐标测量方法精度低、耗时长,不能满足要求,因此断层表面形态的检测更趋向于逆向工程。对于断层表面形态的自由曲面测量是将三坐标测量技术移植到地震物理模拟中的关键。地壳岩层因受力达到一定强度而发生破裂,并沿破裂面有明显相对移动的构造称断层。
地震往往是由断层活动引起的,是断层活动的一种变现,而隧道穿越或接近断层场地会加剧地震动力反应,导致围岩稳定性降低,引起衬砌的整体或局部的破环,断层与隧道的相对位置关系不同,引起的地震应力也不同,因此有必要对断层进行有效的检测,减小地震造成的危害性。
现有的断层检测时无法精确地对不同结构断层的方位和轮廓做进一步地标记和比对,同时无法对微小断层或其他不同范围的断层进行有效检测。
针对上述问题,本发明提出了一种断层三维数字化检测方法,它具有能精确地对不同结构断层的方位和轮廓做进一步地标记和比对,同时能够对微小断层或其他不同范围的断层进行有效检测等优点。
发明内容
本发明的目的在于提供一种断层三维数字化检测方法,利用DMIS技术建立地震超前观测模型,并采集和处理地震数据,开始利用探测仪器进行断层探测实验,通过PET彩色图像技术向断层表面发射彩色光源,并判断断层位置,利用键入四点边界方法对断层进行自动测量,并加以记录,对测量结果进行数据处理并对结果进行自行比对,它具有能精确地对不同结构断层的方位和轮廓做进一步地标记和比对,同时能够对微小断层或其他不同范围的断层进行有效检测等优点,可以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:一种断层三维数字化检测方法,包括以下步骤:
S101:利用DMIS技术建立地震超前观测模型,并采集和处理地震数据;
S102:开始利用探测仪器进行断层探测实验;
S103:通过PET彩色图像技术向断层表面发射彩色光源,并判断断层位置;
S104:利用键入四点边界方法对断层进行自动测量,并加以记录;
S105:对测量结果进行数据处理并对结果进行自行比对。
进一步地,针对步骤S102,还包括以下步骤:
S1021:探测仪器选用PH10M测头,且测头包括三个电子接触器,前端红宝石探针以及自动旋转测座;
S1022:利用探针接触被测物体并使探针偏转一定角度,并至少断开一个接触器;
S1023:利用机器读出此时的X、Y、Z光栅,并记录此时探针所检测到的标准球的球心位置;
S1024:从标准球上重取五点进行重新校正,以此消除误差,同时更换探测角度再重复步骤S1022和S1023。
进一步地,针对步骤S1021,PH10M测头采用偏心转动远离,且测量软件显示的坐标数据不是测量灯光栅数据,而是红宝石球中心数据。
进一步地,针对步骤S1024,既可利用重取五点进行校正也可以在测量前将全部使用的角度预先校准完全后再进行测量,无需再重新进行校正。
进一步地,针对步骤S1024,还包括以下步骤:
S10241:在断层处更换测头角度,并将其初步确定为40度;
S10242:记录更换角度时的两点坐标,并将断层部分作为一个独立五边形进行测量;
S10243:待角度更换完成后键入所记录的坐标,并再次记录其中两点坐标;
S10244:利用步骤S10243所记录的两点坐标作为下一个被测五边形的一条边,以延伸出下一个被测多边形;
S10245:更换测头角度,并将其确定为50度,再重复步骤S10241-S10244。
进一步地,针对步骤S103,还包括以下步骤:
S1031:利用图像采集仪器获取PET发射图像上像素点的像素值,且此像素值包括最大像素值;
S1032:将大于或等于预设的起始阈值的像素值所对应的PET发射图像上的像素点作为起始轮廓区,且起始阈值为小于或等于最大像素值的像素值;
S1033:调整起始阀值,并统计和获取调整后的起始阀值所对应的起始轮廓区上和其内部所包含的像素点个数;
S1034:获取均匀断面上对应的PET发射图像对应位置的像素点个数;
S1035:计算调整后的起始轮廓区内的像素点数与均匀断面上对应位置像素点之间的差值;
S1036:将差值最小值所对应的起始轮廓区作为PET发射图像上构成断面的轮廓区并以此锁定断面的方位。
进一步地,针对步骤S104,还包括以下步骤:
S1041:在测量时,先确定四点边界的一条边并将其设定为起始边长,同时利用超声波探测技术测定断面到测点的距离;
S1042:若需衔接的四边形个数较多,则加长起始边长,并将重新加长后的边长重新作为待测四点边界的边长,再对断层进行测量;
S1043:若需衔接的四边形个数较少,则适当缩短起始边长,并再次对断层进行测量;
S1044:在测量的过程中更换任意角度,且角度的最大值为180度。
进一步地,针对步骤S1043,分割后的四边形边界可以完全衔接,且相邻的四边形衔接处无间隙,缩短时对四边形的多条边长同步进行缩短。
进一步地,针对步骤S1041,在利用超声波探测技术测定断面到测点的距离,具体工作步骤还包括:
S10411:在测距开始时,通过控制单元向超声测距模块中的发射端以及接收端发送控制启动指令;
S10412:所述发射端与所述接收端在接收到所述控制启动指令时,控制计时器开始进行计时;
S10413:所述接收端在接收到所述发射端传输的回波信号后,通过所述控制单元控制所述计时器停止计时,同时,记录所述回波信号到达所述接收端的时间点;
S10414:根据所述回波信号到达所述接收端的时间点与所述断面到测点可探测的距离范围,将所述回波信号划分为N个子信号;
S10415:获取N个子信号中每个子信号的能量值,同时,根据每个子信号的能量值从N个子信号中选取M个候选子信号;
S10416:设定权值函数,并根据所述权值函数获取M个候选子信号中每个候选子信号的选波概率,并将M个候选子信号的选波概率和能量值生成选波参数;
S10417:根据M个候选子信号的选波参数确定超声测距模块的测距目标回波信号;
S10418:对接收到所述测距目标回波信号进行时间信息数据化处理,并将处理结果发送至控制单元;
S10419:所述控制单元接收到所述处理结果后,再次交替向不同超声波测距模块中的发射端和接收端发送启动指令,其中,所述处理结果中包括:所述断面到测点的距离权值;
S10420:重复步骤S10412—S10419,交替启动不同超声波测距模块中的一一对应的发射端和接收端进行测距,并获取所述断面到测点的多组距离权值;
S10421:通过控制单元将获取的所有所述断面到测点的距离权值进行综合分析计算,获得所述断面到测点的最终距离。
进一步地,针对步骤S1023,通过所述标准球的球心位置,同时,获取所述标准球的球心位置的误差系数,并根据误差系数,对所述标准球位置进行重新校正,具体工作过程包括:
通过所述X、Y、Z光栅,获取所述标准球的球心位置,同时,将所述球心位置与所述控制模块中预先设定好的球心位置进行预估;
若所述球心位置与所述控制模块中预先设定好的球心位置发生偏差,控制数据计算模块计算所述球心位置与预先设定好的球心位置之间的偏差系数;
其中,η表示所述球心位置与预先设定好的球心位置之间的偏差系数,A{x,y,z}表示所述标准球球心的实际坐标点,A{X,Y,Z}表示预先设定好的所述标准球球心的目标坐标点,V表示所述标准球的体积,r表示所述标准球的半径,α表示所述标准球的球心位置与预先设定好的所述标准球的球心位置连线与水平线之间的夹角,k表示所述标准球的张力参数;
基于所述球心位置与预先设定好的球心位置之间的偏差系数,通过数据计算模块,计算所述标准球位置需要校正的校正值;
其中,Q表示所述标准球位置需要校正的校正值,η表示所述球心位置与预先设定好的球心位置之间的偏差系数,m表示所述标准球的质量,r表示所述标准球的半径,ρ表示所述标准球的密度,ζ表示校正系数,w表示位置校正的校正频率,W表示所述数据计算模块的工作频率,t表示所述数据计算模块的工作时间,v表示所述数据计算模块的工作速度,l表示校正幅值;
通过控制模块并基于所述校正值控制校正模块对所述标准球的球心位置进行校正;
根据X、Y、Z光栅,重新获取校正后的所述标准球的球心位置,并继续与预先设定好的球心位置进行预估;
若所述球心位置与所述控制模块中预先设定好的球心位置的偏差系数为零,控制校正模块停止工作,完成对所述标准球的球心位置的校正。
与现有技术相比,本发明的有益效果是:
1.本发明提出的一种断层三维数字化检测方法,利用图像采集仪器获取PET发射图像上像素点的像素值,且此像素值包括最大像素值,将大于或等于预设的起始阈值的像素值所对应的PET发射图像上的像素点作为起始轮廓区,且起始阈值为小于或等于最大像素值的像素值,调整起始阀值,并统计和获取调整后的起始阀值所对应的起始轮廓区上和其内部所包含的像素点个数,再获取均匀断面上对应的PET发射图像对应位置的像素点个数,将计算调整后的起始轮廓区内的像素点数与均匀断面上对应位置像素点之间的差值,最后将差值最小值所对应的起始轮廓区作为PET发射图像上构成断面的轮廓区并以此锁定断面的方位,与传统的断面定位方法相比,此方法能够利用PET技术对不同结构断层进行色度编码,以此对断层进行标记和与预设阀值比对,更加精确地知晓不同断层所在方位以及其结构轮廓,观察断层时也更加直观和清晰。
2.本发明提出的一种断层三维数字化检测方法,在测量时,先确定四点边界的一条边并将其设定为起始边长,同时利用超声波探测技术测定断面到测点的距离,若需衔接的四边形个数较多,则加长起始边长,并将重新加长后的边长重新作为待测四点边界的边长,再对断层进行测量,若需衔接的四边形个数较少,则适当缩短起始边长,并再次对断层进行测量,且在测量的过程中可更换任意角度,与一般的键入四点边界方法不同的是其能够先利用超声波探测技术确定所测断面的大致方位,在测量时能够根据所测断面的范围加长或缩短四点边界的边长,并根据不同范围的断面选用不同边长的四点边界将其进行分隔和衔接,以此达到拼接更完全和优化四点检测的目的,也能够对一些微小断层进行检测,解决了传统断层识别技术无法识别微小断层特性的技术问题。
3.本发明提出的一种断层三维数字化检测方法,在断层处更换测头角度,并将其初步确定为40度,记录更换角度时的两点坐标,并将断层部分作为一个独立五边形进行测量,待角度更换完成后键入所记录的坐标,并再次记录其中两点坐标,利用所记录的两点坐标作为下一个被测五边形的一条边,以延伸出下一个被测多边形,更换测头角度,并将其确定为50度,再重复多次,依此类推,测量任务可一次完成,既满足了多边形完全拼接的目的,且测头角度可多次选择,避免断层处的数据不完整或边界部分缺失甚至造成数据出现重叠,在一定程度上也提高了测量精度。
4.在测距开始时,通过控制单元向超声测距模块中的发射端以及接收端发送控制启动指令,从而有利于准确控制计时器进行计时,当接收端在接收到发射端传输的回波信号后,准确控制计时器停止计时,从而有效记录回波信号到达接收端的时间点,根据时间点与断面到测点可探测的距离范围,可以将回波信号划分为N个子信号,通过获取N个子信号中每个子信号的能量值,从而准确快速从N个子信号中选取M个候选子信号,通过设定权值函数,精确获取M个候选子信号中每个候选子信号的选波概率,并将M个候选子信号的选波概率和能量值生成选波参数,从而有利于确定超声测距模块的测距目标回波信号,并对目标回波信号进行时间信息数据化处理,同时将处理结果发送至控制单元,之后再次交替向不同超声波测距模块中的发射端和接收端发送启动指令,并重复上述步骤,获取断面到测点的多组距离权值,同时进行综合分析计算,使得获得的断面到测点的最终距离更加准确。
5.通过X、Y、Z光栅准确获取标准球的球心位置,并将球心位置与控制模块中预先设定好的球心位置进行预估,从而有利于快速判定球心位置是否发生偏差,如果发生偏差,通过球心的实际坐标与目标坐标之间的差异计算该球心位置与预先设定好的球心位置之间的偏差系数,并基于球心位置与预先设定好的球心位置之间的偏差系数,通过数据计算模块,精准计算标准球位置需要校正的校正值,通过控制模块并基于校正值控制校正模块对标准球的球心位置进行校正,同时,在校正后重新获取球心位置,从而继续计算偏差系数,当偏差系数为0时,控制校正模块停止工作,从而完成对标准球的球心位置的校正,该过程大大提升了该方法的智能性与高效性。
附图说明
图1为本发明断层三维检测方法的总系统流程图;
图2为本发明断层三维检测方法的步骤S103的系统流程图;
图3为本发明断层三维检测方法的步骤S104的系统流程图;
图4为本发明断层三维检测方法的步骤S102的系统流程图;
图5为本发明断层三维检测方法的步骤S1024的系统流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参阅图1和图4,一种断层三维数字化检测方法,包括以下步骤:
S101:利用DMIS技术建立地震超前观测模型,并采集和处理地震数据;
S102:开始利用探测仪器进行断层探测实验;
S103:通过PET彩色图像技术向断层表面发射彩色光源,并判断断层位置;
S104:利用键入四点边界方法对断层进行自动测量,并加以记录;
S105:对测量结果进行数据处理并对结果进行自行比对。
针对步骤S102,还包括以下步骤:
S1021:探测仪器选用PH10M测头,且测头包括三个电子接触器,前端红宝石探针以及自动旋转测座;
S1022:利用探针接触被测物体并使探针偏转一定角度,并至少断开一个接触器;
S1023:利用机器读出此时的X、Y、Z光栅,并记录此时探针所检测到的标准球的球心位置;
S1024:从标准球上重取五点进行重新校正,以此消除误差,同时更换探测角度再重复步骤S1022和S1023。
针对步骤S1021,PH10M测头采用偏心转动远离,且测量软件显示的坐标数据不是测量灯光栅数据,而是红宝石球中心数据。
针对步骤S1024,既可利用重取五点进行校正也可以在测量前将全部使用的角度预先校准完全后再进行测量,无需再重新进行校正。
参阅图5,一种断层三维数字化检测方法,针对步骤S1024,还包括以下步骤:
S10241:在断层处更换测头角度,并将其初步确定为40度;
S10242:记录更换角度时的两点坐标,并将断层部分作为一个独立五边形进行测量;
S10243:待角度更换完成后键入所记录的坐标,并再次记录其中两点坐标;
S10244:利用步骤S10243所记录的两点坐标作为下一个被测五边形的一条边,以延伸出下一个被测多边形;
S10245:更换测头角度,并将其确定为50度,再重复步骤S10241-S10244。
依此类推,测量任务可一次完成,既满足了多边形完全拼接的目的,且测头角度可多次选择,避免断层处的数据不完整或边界部分缺失甚至造成数据出现重叠,在一定程度上也提高了测量精度。
参阅图2-3,一种断层三维数字化检测方法,针对步骤S103,还包括以
下步骤:
S1031:利用图像采集仪器获取PET发射图像上像素点的像素值,且此像素值包括最大像素值;
S1032:将大于或等于预设的起始阈值的像素值所对应的PET发射图像上的像素点作为起始轮廓区,且起始阈值为小于或等于最大像素值的像素值;
S1033:调整起始阀值,并统计和获取调整后的起始阀值所对应的起始轮廓区上和其内部所包含的像素点个数;
S1034:获取均匀断面上对应的PET发射图像对应位置的像素点个数;
S1035:计算调整后的起始轮廓区内的像素点数与均匀断面上对应位置像素点之间的差值;
S1036:将差值最小值所对应的起始轮廓区作为PET发射图像上构成断面的轮廓区并以此锁定断面的方位。
与传统的断面定位方法相比,此方法能够利用PET技术对不同结构断层进行色度编码,以此对断层进行标记和与预设阀值比对,更加精确地知晓不同断层所在方位以及其结构轮廓,观察断层时也更加直观和清晰。
针对步骤S104,还包括以下步骤:
S1041:在测量时,先确定四点边界的一条边并将其设定为起始边长,同时利用超声波探测技术测定断面到测点的距离;
S1042:若需衔接的四边形个数较多,则加长起始边长,并将重新加长后的边长重新作为待测四点边界的边长,再对断层进行测量;
S1043:若需衔接的四边形个数较少,则适当缩短起始边长,并再次对断层进行测量;
S1044:在测量的过程中更换任意角度,且角度的最大值为180度。
与一般的键入四点边界方法不同的是其能够先利用超声波探测技术确定所测断面的大致方位,在测量时能够根据所测断面的范围加长或缩短四点边界的边长,并根据不同范围的断面选用不同边长的四点边界将其进行分隔和衔接,以此达到拼接更完全和优化四点检测的目的,也能够对一些微小断层进行检测,解决了传统断层识别技术无法识别微小断层特性的技术问题。
进一步地,针对步骤S1043,分割后的四边形边界可以完全衔接,且相邻的四边形衔接处无间隙,缩短时对四边形的多条边长同步进行缩短。
综上所述:本发明提出的一种断层三维数字化检测方法,利用DMIS技术建立地震超前观测模型,并采集和处理地震数据,开始利用探测仪器进行断层探测实验,在断层处更换测头角度,并将其初步确定为40度,记录更换角度时的两点坐标,并将断层部分作为一个独立五边形进行测量,待角度更换完成后键入所记录的坐标,并再次记录其中两点坐标,利用所记录的两点坐标作为下一个被测五边形的一条边,以延伸出下一个被测多边形,更换测头角度,并将其确定为50度,再重复多次,依此类推,测量任务可一次完成,既满足了多边形完全拼接的目的,且测头角度可多次选择,避免断层处的数据不完整或边界部分缺失甚至造成数据出现重叠,在一定程度上也提高了测量精度,通过PET彩色图像技术向断层表面发射彩色光源,并判断断层位置,利用图像采集仪器获取PET发射图像上像素点的像素值,且此像素值包括最大像素值,将大于或等于预设的起始阈值的像素值所对应的PET发射图像上的像素点作为起始轮廓区,且起始阈值为小于或等于最大像素值的像素值,调整起始阀值,并统计和获取调整后的起始阀值所对应的起始轮廓区上和其内部所包含的像素点个数,再获取均匀断面上对应的PET发射图像对应位置的像素点个数,将计算调整后的起始轮廓区内的像素点数与均匀断面上对应位置像素点之间的差值,最后将差值最小值所对应的起始轮廓区作为PET发射图像上构成断面的轮廓区并以此锁定断面的方位,与传统的断面定位方法相比,此方法能够利用PET技术对不同结构断层进行色度编码,以此对断层进行标记和与预设阀值比对,更加精确地知晓不同断层所在方位以及其结构轮廓,观察断层时也更加直观和清晰,利用键入四点边界方法对断层进行自动测量,并加以记录,在测量时,先确定四点边界的一条边并将其设定为起始边长,同时利用超声波探测技术测定断面到测点的距离,若需衔接的四边形个数较多,则加长起始边长,并将重新加长后的边长重新作为待测四点边界的边长,再对断层进行测量,若需衔接的四边形个数较少,则适当缩短起始边长,并再次对断层进行测量,且在测量的过程中可更换任意角度,与一般的键入四点边界方法不同的是其能够先利用超声波探测技术确定所测断面的大致方位,在测量时能够根据所测断面的范围加长或缩短四点边界的边长,并根据不同范围的断面选用不同边长的四点边界将其进行分隔和衔接,以此达到拼接更完全和优化四点检测的目的,也能够对一些微小断层进行检测,解决了传统断层识别技术无法识别微小断层特性的技术问题,对测量结果进行数据处理并对结果进行自行比对。
本发明提供了一种断层三维数字化检测方法,针对步骤S1041,在利用超声波探测技术测定断面到测点的距离,具体工作步骤还包括:
S10411:在测距开始时,通过控制单元向超声测距模块中的发射端以及接收端发送控制启动指令;
S10412:所述发射端与所述接收端在接收到所述控制启动指令时,控制计时器开始进行计时;
S10413:所述接收端在接收到所述发射端传输的回波信号后,通过所述控制单元控制所述计时器停止计时,同时,记录所述回波信号到达所述接收端的时间点;
S10414:根据所述回波信号到达所述接收端的时间点与所述断面到测点可探测的距离范围,将所述回波信号划分为N个子信号;
S10415:获取N个子信号中每个子信号的能量值,同时,根据每个子信号的能量值从N个子信号中选取M个候选子信号;
S10416:设定权值函数,并根据所述权值函数获取M个候选子信号中每个候选子信号的选波概率,并将M个候选子信号的选波概率和能量值生成选波参数;
S10417:根据M个候选子信号的选波参数确定超声测距模块的测距目标回波信号;
S10418:对接收到所述测距目标回波信号进行时间信息数据化处理,并将处理结果发送至控制单元;
S10419:所述控制单元接收到所述处理结果后,再次交替向不同超声波测距模块中的发射端和接收端发送启动指令,其中,所述处理结果中包括:所述断面到测点的距离权值;
S10420:重复步骤S10412—S10419,交替启动不同超声波测距模块中的一一对应的发射端和接收端进行测距,并获取所述断面到测点的多组距离权值;
S10421:通过控制单元将获取的所有所述断面到测点的距离权值进行综合分析计算,获得所述断面到测点的最终距离。
该实施例中,发送控制启动指令可以是人工操作发送,也可以是系统自动发送控制启动指令。
该实施例中,回波信号可以是当接收端不能完全吸收波的能量时,波的一部分能量会被反射回来,就形成了回波信号。
该实施例中,权值函数可以是根据候选子信号的信号能量均值所构成的函数。
该实施例中,选波概率可以是与回波信号相匹配的候选子信号的概率。
该实施例中,选波参数可以是衡量超声波测量距离的参数。
该实施例中,测距目标回波信号,可以根据当超声波测距反射回的能量符合选波参数的标准时的信号。
该实施例中,距离权值可以是通过不同超声波模块获取的断面到测点的距离均值的频数。
上述技术方案的工作原理及有益效果是:在测距开始时,通过控制单元向超声测距模块中的发射端以及接收端发送控制启动指令,从而有利于准确控制计时器进行计时,当接收端在接收到发射端传输的回波信号后,准确控制计时器停止计时,从而有效记录回波信号到达接收端的时间点,根据时间点与断面到测点可探测的距离范围,可以将回波信号划分为N个子信号,通过获取N个子信号中每个子信号的能量值,从而准确快速从N个子信号中选取M个候选子信号,通过设定权值函数,精确获取M个候选子信号中每个候选子信号的选波概率,并将M个候选子信号的选波概率和能量值生成选波参数,从而有利于确定超声测距模块的测距目标回波信号,并对目标回波信号进行时间信息数据化处理,同时将处理结果发送至控制单元,之后再次交替向不同超声波测距模块中的发射端和接收端发送启动指令,并重复上述步骤,获取断面到测点的多组距离权值,同时进行综合分析计算,使得获得的断面到测点的最终距离更加准确。
本发明提供了一种断层三维数字化检测方法,针对步骤S1023,通过所述标准球的球心位置,同时,获取所述标准球的球心位置的误差系数,并根据误差系数,对所述标准球位置进行重新校正,具体工作过程包括:
通过所述X、Y、Z光栅,获取所述标准球的球心位置,同时,将所述球心位置与所述控制模块中预先设定好的球心位置进行预估;
若所述球心位置与所述控制模块中预先设定好的球心位置发生偏差,控制数据计算模块计算所述球心位置与预先设定好的球心位置之间的偏差系数;
其中,η表示所述球心位置与预先设定好的球心位置之间的偏差系数,A{x,y,z}表示所述标准球球心的实际坐标点,A{X,Y,Z}表示预先设定好的所述标准球球心的目标坐标点,V表示所述标准球的体积,r表示所述标准球的半径,α表示所述标准球的球心位置与预先设定好的所述标准球的球心位置连线与水平线之间的夹角,k表示所述标准球的张力参数;
基于所述球心位置与预先设定好的球心位置之间的偏差系数,通过数据计算模块,计算所述标准球位置需要校正的校正值;
其中,Q表示所述标准球位置需要校正的校正值,η表示所述球心位置与预先设定好的球心位置之间的偏差系数,m表示所述标准球的质量,r表示所述标准球的半径,ρ表示所述标准球的密度,ζ表示校正系数,w表示位置校正的校正频率,W表示所述数据计算模块的工作频率,t表示所述数据计算模块的工作时间,v表示所述数据计算模块的工作速度,l表示校正幅值;
通过控制模块并基于所述校正值控制校正模块对所述标准球的球心位置进行校正;
根据X、Y、Z光栅,重新获取校正后的所述标准球的球心位置,并继续与预先设定好的球心位置进行预估;
若所述球心位置与所述控制模块中预先设定好的球心位置的偏差系数为零,控制校正模块停止工作,完成对所述标准球的球心位置的校正。
上述技术方案的工作原理及有益效果是:通过X、Y、Z光栅准确获取标准球的球心位置,并将球心位置与控制模块中预先设定好的球心位置进行预估,从而有利于快速判定球心位置是否发生偏差,如果发生偏差,通过球心的实际坐标与目标坐标之间的差异计算该球心位置与预先设定好的球心位置之间的偏差系数,并基于球心位置与预先设定好的球心位置之间的偏差系数,通过数据计算模块,精准计算标准球位置需要校正的校正值,通过控制模块并基于校正值控制校正模块对标准球的球心位置进行校正,同时,在校正后重新获取球心位置,从而继续计算偏差系数,当偏差系数为0时,控制校正模块停止工作,从而完成对标准球的球心位置的校正,该过程大大提升了该方法的智能性与高效性。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (7)
1.一种断层三维数字化检测方法,其特征在于,包括以下步骤:
S101:利用DMIS技术建立地震超前观测模型,并采集和处理地震数据;
S102:开始利用探测仪器进行断层探测实验;
S103:通过PET彩色图像技术向断层表面发射彩色光源,并判断断层位置;
S104:利用键入四点边界方法对断层进行自动测量,并加以记录;
S105:对测量结果进行数据处理并对结果进行自行比对;
针对步骤S102,还包括以下步骤:
S1021:探测仪器选用PH10M测头,且测头包括三个电子接触器,前端红宝石探针以及自动旋转测座;
S1022:利用探针接触被测物体并使探针偏转一定角度,并至少断开一个接触器;
S1023:利用机器读出此时的X、Y、Z光栅,并记录此时探针所检测到的标准球的球心位置;
S1024:从标准球上重取五点进行重新校正,以此消除误差,同时更换探测角度再重复步骤S1022和S1023;
针对步骤S103,还包括以下步骤:
S1031:利用图像采集仪器获取PET发射图像上像素点的像素值,且此像素值包括最大像素值;
S1032:将大于或等于预设的起始阈值的像素值所对应的PET发射图像上的像素点作为起始轮廓区,且起始阈值为小于或等于最大像素值的像素值;
S1033:调整起始阈值,并统计和获取调整后的起始阈值所对应的起始轮廓区上和其内部所包含的像素点个数;
S1034:获取均匀断面上对应的PET发射图像对应位置的像素点个数;
S1035:计算调整后的起始轮廓区内的像素点数与均匀断面上对应位置像素点数之间的差值;
S1036:将差值最小值所对应的起始轮廓区作为PET发射图像上构成断面的轮廓区并以此锁定断面的方位;
针对步骤S104,还包括以下步骤:
S1041:在测量时,先确定四点边界的一条边并将其设定为起始边长,同时利用超声波探测技术测定断面到测点的距离;
S1042:若需衔接的四边形个数较多,则加长起始边长,并将重新加长后的边长重新作为待测四点边界的边长,再对断层进行测量;
S1043:若需衔接的四边形个数较少,则适当缩短起始边长,并再次对断层进行测量;
S1044:在测量的过程中更换任意角度,且角度的最大值为180度。
2.如权利要求1所述的一种断层三维数字化检测方法,其特征在于,针对步骤S1021,PH10M测头采用偏心转动远离,且测量软件显示的坐标数据不是测量灯光栅数据,而是红宝石球中心数据。
3.如权利要求1所述的一种断层三维数字化检测方法,其特征在于,针对步骤S1024,既可利用重取五点进行校正也可以在测量前将全部使用的角度预先校准完全后再进行测量,无需再重新进行校正。
4.如权利要求1所述的一种断层三维数字化检测方法,其特征在于,针对步骤S1024,还包括以下步骤:
S10241:在断层处更换测头角度,并将其初步确定为40度;
S10242:记录更换角度时的两点坐标,并将断层部分作为一个独立五边形进行测量;
S10243:待角度更换完成后键入所记录的坐标,并再次记录其中两点坐标;
S10244:利用步骤S10243所记录的两点坐标作为下一个被测五边形的一条边,以延伸出下一个被测多边形;
S10245:更换测头角度,并将其确定为50度,再重复步骤S10241-S10244。
5.如权利要求1所述的一种断层三维数字化检测方法,其特征在于,针对步骤S1043,分割后的四边形边界可以完全衔接,且相邻的四边形衔接处无间隙,缩短时对四边形的多条边长同步进行缩短。
6.如权利要求1所述的一种断层三维数字化检测方法,其特征在于,针对步骤S1041,在利用超声波探测技术测定断面到测点的距离,具体工作步骤还包括:
S10411:在测距开始时,通过控制单元向超声测距模块中的发射端以及接收端发送控制启动指令;
S10412:所述发射端与所述接收端在接收到所述控制启动指令时,控制计时器开始进行计时;
S10413:所述接收端在接收到所述发射端传输的回波信号后,通过所述控制单元控制所述计时器停止计时,同时,记录所述回波信号到达所述接收端的时间点;
S10414:根据所述回波信号到达所述接收端的时间点与所述断面到测点可探测的距离范围,将所述回波信号划分为N个子信号;
S10415:获取N个子信号中每个子信号的能量值,同时,根据每个子信号的能量值从N个子信号中选取M个候选子信号;
S10416:设定权值函数,并根据所述权值函数获取M个候选子信号中每个候选子信号的选波概率,并将M个候选子信号的选波概率和能量值生成选波参数;
S10417:根据M个候选子信号的选波参数确定超声测距模块的测距目标回波信号;
S10418:对接收到所述测距目标回波信号进行时间信息数据化处理,并将处理结果发送至控制单元;
S10419:所述控制单元接收到所述处理结果后,再次交替向不同超声波测距模块中的发射端和接收端发送启动指令,其中,所述处理结果中包括:所述断面到测点的距离权值;
S10420:重复步骤S10412—S10419,交替启动不同超声波测距模块中的一一对应的发射端和接收端进行测距,并获取所述断面到测点的多组距离权值;
S10421:通过控制单元将获取的所有所述断面到测点的距离权值进行综合分析计算,获得所述断面到测点的最终距离。
7.如权利要求1所述的一种断层三维数字化检测方法,其特征在于,针对步骤S1023,通过所述标准球的球心位置,同时,获取所述标准球的球心位置的误差系数,并根据误差系数,对所述标准球位置进行重新校正,具体工作过程包括:
通过所述X、Y、Z光栅,获取所述标准球的球心位置,同时,将所述球心位置与控制单元中预先设定好的球心位置进行预估;
若所述球心位置与所述控制单元中预先设定好的球心位置发生偏差,控制数据计算模块计算所述球心位置与预先设定好的球心位置之间的偏差系数;
其中,η表示所述球心位置与预先设定好的球心位置之间的偏差系数,A{x,y,z}表示所述标准球球心的实际坐标点,A{X,Y,Z}表示预先设定好的所述标准球球心的目标坐标点,V表示所述标准球的体积,r表示所述标准球的半径,α表示所述标准球的球心位置与预先设定好的所述标准球的球心位置连线与水平线之间的夹角,k表示所述标准球的张力参数;
基于所述球心位置与预先设定好的球心位置之间的偏差系数,通过数据计算模块,计算所述标准球位置需要校正的校正值;
其中,Q表示所述标准球位置需要校正的校正值,η表示所述球心位置与预先设定好的球心位置之间的偏差系数,m表示所述标准球的质量,r表示所述标准球的半径,ρ表示所述标准球的密度,ζ表示校正系数,w表示位置校正的校正频率,W表示所述数据计算模块的工作频率,t表示所述数据计算模块的工作时间,v表示所述数据计算模块的工作速度,l表示校正幅值;
通过控制单元并基于所述校正值控制校正模块对所述标准球的球心位置进行校正;
根据X、Y、Z光栅,重新获取校正后的所述标准球的球心位置,并继续与预先设定好的球心位置进行预估;
若所述球心位置与所述控制单元中预先设定好的球心位置的偏差系数为零,控制校正模块停止工作,完成对所述标准球的球心位置的校正。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011241490.6A CN112363209B (zh) | 2020-11-09 | 2020-11-09 | 一种断层三维数字化检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011241490.6A CN112363209B (zh) | 2020-11-09 | 2020-11-09 | 一种断层三维数字化检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112363209A CN112363209A (zh) | 2021-02-12 |
CN112363209B true CN112363209B (zh) | 2023-07-25 |
Family
ID=74509154
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011241490.6A Active CN112363209B (zh) | 2020-11-09 | 2020-11-09 | 一种断层三维数字化检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112363209B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101506686A (zh) * | 2006-06-21 | 2009-08-12 | 特拉斯帕克地球科学公司 | 地质沉积体系的解释 |
AU2013200609A1 (en) * | 2008-04-11 | 2013-02-21 | Terraspark Geosciences, Llc | Visulation of geologic features using data representations thereof |
CN104459768A (zh) * | 2013-09-22 | 2015-03-25 | 中国石油化工股份有限公司 | 一种基于可视化的三维空间目标地质体追踪方法 |
CN106443667A (zh) * | 2010-02-14 | 2017-02-22 | 弗米尔制造公司 | 用于地下对象检测的导数成像 |
CN107589449A (zh) * | 2017-08-29 | 2018-01-16 | 电子科技大学 | 基于曲线Gabor滤波的三维数据断层增强方法 |
CN109709300A (zh) * | 2018-11-16 | 2019-05-03 | 高军 | 一种隧道断层泥注浆加固正交实验分析方法及实验系统 |
CN210198334U (zh) * | 2019-08-12 | 2020-03-27 | 苏州市计量测试院 | 一种几何量测量校准的标准装置 |
WO2020127302A1 (en) * | 2018-12-21 | 2020-06-25 | Joswig & Walter Patent Gbr | Arrangement and process for locating and/or characterising fracture events in the earth crust, particularly suitable to monitor fracking |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8081796B2 (en) * | 2008-11-24 | 2011-12-20 | Ingrain, Inc. | Method for determining properties of fractured rock formations using computer tomograpic images thereof |
US9618639B2 (en) * | 2012-03-01 | 2017-04-11 | Drilling Info, Inc. | Method and system for image-guided fault extraction from a fault-enhanced seismic image |
-
2020
- 2020-11-09 CN CN202011241490.6A patent/CN112363209B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101506686A (zh) * | 2006-06-21 | 2009-08-12 | 特拉斯帕克地球科学公司 | 地质沉积体系的解释 |
AU2013200609A1 (en) * | 2008-04-11 | 2013-02-21 | Terraspark Geosciences, Llc | Visulation of geologic features using data representations thereof |
CN106443667A (zh) * | 2010-02-14 | 2017-02-22 | 弗米尔制造公司 | 用于地下对象检测的导数成像 |
CN104459768A (zh) * | 2013-09-22 | 2015-03-25 | 中国石油化工股份有限公司 | 一种基于可视化的三维空间目标地质体追踪方法 |
CN107589449A (zh) * | 2017-08-29 | 2018-01-16 | 电子科技大学 | 基于曲线Gabor滤波的三维数据断层增强方法 |
CN109709300A (zh) * | 2018-11-16 | 2019-05-03 | 高军 | 一种隧道断层泥注浆加固正交实验分析方法及实验系统 |
WO2020127302A1 (en) * | 2018-12-21 | 2020-06-25 | Joswig & Walter Patent Gbr | Arrangement and process for locating and/or characterising fracture events in the earth crust, particularly suitable to monitor fracking |
CN210198334U (zh) * | 2019-08-12 | 2020-03-27 | 苏州市计量测试院 | 一种几何量测量校准的标准装置 |
Non-Patent Citations (2)
Title |
---|
"基于DMIS的断层三维表面形态自动检测方法研究";席嘉文;《地震工程学报》;全文 * |
"谱分解技术在断层描述及储层预测中的应用";黄龙泽;《海洋石油》;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112363209A (zh) | 2021-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105203023B (zh) | 一种车载三维激光扫描系统安置参数的一站式标定方法 | |
EP4335739A2 (en) | Method and system for determining the draft of a vessel | |
CN102636785B (zh) | 一种水下目标三维定位方法 | |
CN110133639B (zh) | 一种传力杆施工质量检测方法 | |
CN108311545B (zh) | 一种y型轧机连轧对中及孔型检测系统及方法 | |
CN113553763A (zh) | 一种隧道围岩快速分级方法和系统 | |
JP2017207457A (ja) | 領域変位算出システム、領域変位算出方法、及び領域変位算出プログラム | |
CN114994061A (zh) | 一种基于机器视觉的钢轨智能化检测方法及系统 | |
CN106526598A (zh) | 侧扫声纳水下高精度定位系统及方法 | |
CN109001694B (zh) | 一种动态自适应天线扫描特性模拟方法及系统 | |
CN113487722A (zh) | 一种基于三维激光扫描方法的混凝土构件自动化检测方法 | |
CN114037706A (zh) | 一种基于三维点云模型的预制梁段表面平整度检测方法 | |
CN108827200A (zh) | 一种船体分段智能检测系统及方法 | |
CN206321792U (zh) | 侧扫声纳水下高精度定位系统 | |
CN102224523A (zh) | 立体匹配处理系统、立体匹配处理方法和记录媒介 | |
CN112363209B (zh) | 一种断层三维数字化检测方法 | |
CN116433885A (zh) | 一种基于亚像素边缘的多开口销定位方法 | |
NO332515B1 (no) | Seismisk undersokelse og analyse | |
CN111854705A (zh) | 一种水下测绘方法 | |
NO300190B1 (no) | Fremgangsmåte ved innsamling av seismiske data | |
JP2020038218A5 (ja) | 処理装置、検査システム、処理方法、プログラム、及び記憶媒体 | |
CN113075732B (zh) | 一种消除高分辨率小多道地震地层异常起伏的方法 | |
CN113568041B (zh) | 时移地震三维拖缆采集数据的可重复性分析方法及系统 | |
KR102606726B1 (ko) | 선박의 블록 분석시스템 | |
CN116576806B (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 |