CN113640869A - 基于三维真地表的观测系统自动排障量化方法 - Google Patents
基于三维真地表的观测系统自动排障量化方法 Download PDFInfo
- Publication number
- CN113640869A CN113640869A CN202010390756.7A CN202010390756A CN113640869A CN 113640869 A CN113640869 A CN 113640869A CN 202010390756 A CN202010390756 A CN 202010390756A CN 113640869 A CN113640869 A CN 113640869A
- Authority
- CN
- China
- Prior art keywords
- observation system
- shot
- point
- points
- earth surface
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000011002 quantification Methods 0.000 title claims abstract description 19
- 238000010276 construction Methods 0.000 claims abstract description 111
- 238000001514 detection method Methods 0.000 claims abstract description 63
- 230000005284 excitation Effects 0.000 claims abstract description 41
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 15
- 230000004888 barrier function Effects 0.000 claims description 35
- 238000003384 imaging method Methods 0.000 claims description 22
- 238000011158 quantitative evaluation Methods 0.000 claims description 19
- 238000007689 inspection Methods 0.000 claims description 12
- 238000009827 uniform distribution Methods 0.000 claims description 12
- 238000010521 absorption reaction Methods 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 7
- 238000005516 engineering process Methods 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 6
- 239000000523 sample Substances 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000035772 mutation Effects 0.000 claims description 4
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 230000008030 elimination Effects 0.000 claims 2
- 238000011156 evaluation Methods 0.000 claims 2
- 238000013139 quantization Methods 0.000 claims 2
- 238000004445 quantitative analysis Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 17
- 238000013461 design Methods 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000003079 shale oil Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
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. analysis, for interpretation, for correction
- G01V1/282—Application of seismic models, synthetic seismograms
Abstract
本发明提供一种基于三维真地表的观测系统自动排障量化方法,包括:步骤1,设计规则的观测系统,建立三维真地表构造模型;步骤2,进行野外施工障碍物的自动识别和圈定,将圈定的障碍物进行分类,设定激发点施工可行距离;步骤3,根据圈定的地表障碍物位置,自动选取需要避障的炮点和检波点;步骤4,利用自动避障算法,将选取的炮点和检波点进行移动;步骤5,输出自动避障后的炮点和检波点位置信息,量化评价自动避障后的观测系统。该基于三维真地表的观测系统自动排障量化方法针对三维真地表构造模型进行观测系统自动避障,在大幅提高施工效率的同时,确保了观测系统符合施工要求,从而为采集高品质数据奠定了基础。
Description
技术领域
本发明涉及油气田勘探开发技术领域,特别是涉及到一种基于三维真地表的观测系统自动排障量化方法。
背景技术
随着复杂构造油气藏、薄层岩性油气藏、页岩油气资源等勘探开发程度的加深,在地震数据野外采集中,需要设计全方位、高密度的观测系统,观测系统具有宽方位、宽频带、高密度、高覆盖次数、单点接收、小炮间距、小道间距的特点。设计“两宽一高”或高密度全方位的野外观测系统,对于改善地震资料成像效果、提高储层描述精度具有重要作用,但是由于复杂地表障碍物的不利影响,使得实际施工观测系统存在严重的变观问题,从而影响了后续的地震资料处理和综合解释。国内外专家学者对采集观测系统的设计准则、施工参数选取以及野外施工步骤进行了详细的研究,基于波动方程的地下储层地震照明分析,可以优化复杂地表区观测系统的野外施工方案,并可以针对地下目的层能量屏蔽区与弱能量反射区,通过炮点变观、加密炮点、增加排列等手段,对观测系统进行变观,从而保证施工质量。
常规方法通常采用二维卫星图片平面图,进行手工的拾取和标定障碍物位置,然后将障碍物内的炮点和检波点移动出去,存在施工效率低、变观没有依据、缺乏监控手段的问题。
在申请号:201010165002.8的中国专利申请中,涉及到基于地理信息的三维地震勘探采集观测系统自动避障优化设计技术,采用以下步骤实现:(1)首先建立适合地震勘探采集设计分析应用的地理信息集合数据结构,完成地理信息勘探属性的分类和定义。(2)其次根据地理信息的勘探属性和避障属性值自动扩展生成点、线、面地物禁炮区域,以及无避障属性地物的禁接收区域、可激发区域、避高区域,建立准确的地表勘探地理信息矢量化数字模型。(3)然后进行拟定的理论观测系统采集覆盖周期内各CMP固有采集属性统计计算,包括:覆盖次数、炮检距分布、炮检方位分布、纵向炮检距分布和非纵距分布。(4)完成排列随调整激发点同向移动、不移动和反向移动采集观测方式任意激发点与排列对应关系自动定义。(5)以拟定的理论观测系统固有采集属性为基准;以建立的地表勘探地理信息矢量化数字模型为基础;以勘探目标为主体完成所有可布设勘探物理点的采集属性分析、与CMP固有采集属性吻合度分析,进行勘探采集观测系统自动避障优化设计,实现:①基于地理信息的采集观测系统自动避障优化设计;②基于已知可用激发点空间坐标的采集观测系统自动优化设计;③基于地理信息的跨束非纵采集观测系统自动避障优化设计;④基于地理信息或基于已知可用激发点空间坐标的本束横向对称偏移激发点采集观测系统自动优化设计;⑤基于地理信息或基于已知可用激发点空间坐标的挂炮线采集观测系统自动优化设计。(6)局部大障碍区目标采集观测系统定义及目标采集观测系统与束状采集观测系统衔接采集观测模拟分析计算。地表高陡地形、地下疏松条带、地下空腔和地下巨厚低速区严重影响采集质量,该专利建立的地表勘探地理信息矢量化数字模型只包含地表信息,未考虑地表以下低降速带速度与厚度、Q体的影响,地表以下100米范围内的近地表对激发和接收影响剧烈,炮点在不同地层激发产生的能量、频率差异明显,这是影响采集质量的关键因素,所以在自动避障中必须考虑。另外,该专利无法自动圈定障碍物,采用数学集合进行避障,对避障后的观测系统进行定性分析,缺乏定量化的评价准则。为此我们发明了一种新的基于三维真地表的观测系统自动排障量化方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种可以快速高效地得到最优的野外采集观测系统的基于三维真地表的观测系统自动排障量化方法。
本发明的目的可通过如下技术措施来实现:基于三维真地表的观测系统自动排障量化方法,该基于三维真地表的观测系统自动排障量化方法包括:步骤1,设计理论的规则观测系统,建立三维真地表构造模型;步骤2,进行野外施工障碍物的自动识别和圈定,将圈定的障碍物进行分类,设定激发点施工可行距离;步骤3,根据圈定的地表障碍物位置,自动选取需要避障的炮点和检波点;步骤4,利用自动避障算法,将选取的炮点和检波点进行移动;步骤5,输出自动避障后的炮点和检波点位置信息,量化评价自动避障后的观测系统,定量评价准则指导炮点和检波点位置的再优化。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,根据施工工区的地表条件、构造特征和地质任务,设计规则的观测系统,并形成规则观测系统文本资料,导入施工工区的卫星图片、近地表速度、Q体和地表高程,结合野外踏勘情况,建立三维真地表构造模型,真地表包含了地表、低降速带、高速层顶界面等地层信息。
在步骤1中,所述三维真地表构造模型是由卫星图片、近地表速度、Q体、地表高程和野外踏勘照片建立的,通过经纬度坐标、低降速带厚度和地表高程真实的反映了地表障碍物的变化情况,同时包含不同位置的吸收衰减特性,实现了由二维平面标量图向三位立体向量图的转变。
在步骤2中,在三维构造模型中进行野外施工障碍物的自动识别和圈定,结合障碍物特点及安全施工要求,将圈定的障碍物进行分类,并根据分类情况设定激发点施工可行距离。
在步骤2中,根据四色图像数字化处理技术,在三维真地表构造模型中自动圈定地表障碍物,地表障碍物信息包括:经纬度坐标、低降速带厚度、吸收衰减特性、地表高程;结合障碍物特点及安全施工要求,将圈定的障碍物进行分类,障碍物包括:村镇、工厂、河流、公路、养殖场、疏松沉积带、巨厚低速区、井场和高陡山丘,其中:养殖场、井场、高陡山丘属于I类障碍物;村镇、工厂属于II类障碍物;公路、河流、疏松沉积带、巨厚低速区属于III类障碍物;I类障碍物激发点施工可行距离20m,II类障碍物激发点施工可行距离10m,III类障碍物激发点施工可行距离5m。
在步骤3中,将规则的观测系统投影到真地表三维构造模型中,根据圈定的地表障碍物位置,结合障碍物分类情况,自动选取需要避障的炮点和检波点。
在步骤3中,将理论观测系统资料导入,炮点和检波点投影到三维真地表构造模型中,根据圈定的地表障碍物位置,结合障碍物分类情况,自动选取需要避障的炮点和检波点;所述自动选取需要避障的炮点和检波点方法为,在三维向量图中,假定炮点或检波点为P,障碍物范围为一多边形区域Area,同时考虑地表高程的影响,从P点向左侧或右侧水平方向试射射线,射线与多边形区域Area会产生交点数,根据交点个数判断P点是否在多边形区域Area内;如果交点数为奇数或高程突变超出预设阈值,则P点在多边形区域Area内;若交点数为偶数,则P点不在多边形区域Area内;遍历搜索所有炮点和检波点,从而得到落入障碍物中的炮检点,并形成无需避障的观测系统资料、需要避障的观测系统资料。
在步骤4中,利用自动避障算法,将选取的炮点和检波点进行移动,检波点无需考虑施工可行距离,直接完成自动避障,在进行炮点移动时,根据激发点施工可行距离,将炮点移动到施工可行位置,从而完成炮点的自动避障。
在步骤4中,利用自动避障算法,导入需要避障的观测系统资料,将选取的炮点和检波点进行移动;检波点无需考虑施工可行距离,直接完成自动避障;在进行炮点移动时,结合激发点施工可行距离,将炮点移动到施工可行位置,从而完成炮点的自动避障;所述自动避障算法为,读取某一障碍物范围内的炮检点,根据炮检点落入障碍物的位置,按照从中间到两边、从下到上的原则,确定障碍物内每一个炮检点的移动顺序,然后按照规划的移动顺序,上下交替搜索障碍物之外的可放置空点,将炮检点移动到障碍物外侧;多个障碍物的情况下,障碍物计算顺序从上到下依次进行;在进行炮点移动时,结合三类障碍物的激发点施工可行距离,进行炮点的自动避障。
在步骤5中,输出自动避障后的炮点和检波点位置信息,计算观测系统覆盖次数、炮道密度、方位角、炮检距均匀分布和成像算子属性图,量化评价自动避障后的观测系统。
在步骤5中,输出自动避障后的炮点和检波点位置信息,计算观测系统覆盖次数、炮道密度、方位角、炮检距均匀分布和成像算子属性图;提出的量化评价公式为:θ=α×Fold+ρ×ST+β×Zip+φ×RS+λ×IM,其中:α为覆盖次数权重系数,Fold为覆盖次数,ρ为炮道密度权重系数,ST为炮道密度,β为方位角权重系数,Zip为方位角,φ为炮检距均匀分布权重系数,RS为炮检距均匀分布,λ为成像算子权重系数,IM为成像算子;综合考虑观测系统覆盖次数、炮道密度、方位角、炮检距均匀分布和成像算子属性,评价自动避障后的观测系统,根据量化评价结果和定量评价准则决定是否需要优化加密炮检点或输出最终的观测系统资料。
该基于三维真地表的观测系统自动排障量化方法还包括,在步骤5之后,根据量化评价结果和定量评价准则,若不满足观测系统施工要求,则进行炮点和检波点的加密;若满足观测系统施工要求,则输出最终的观测系统资料,完成观测系统的自动避障。
本发明中的基于三维真地表的观测系统自动排障量化方法,从建立三维真地表构造模型入手,通过障碍物的自动识别、圈定和分类,利用自动避障算法,实现炮检点的自动避障,并通过覆盖次数、炮道密度、方位角、炮检距均匀分布和成像算子属性图,量化评价自动避障后的观测系统。在大幅提高施工效率的同时,确保了观测系统符合施工要求,从而为采集高品质数据奠定了基础。该方法根据卫星图片、近地表速度与厚度、Q体、地表高程和野外踏勘照片建立三维真地表构造模型,综合评估经纬度坐标、低降速带厚度、吸收衰减特性、地表高程等障碍物信息,通过障碍物的自动识别、圈定和分类,同时充分考虑了激发点施工可行距离,利用自动避障算法,实现了炮检点的自动避障,并通过覆盖次数、炮道密度、方位角、炮检距均匀分布和成像算子属性,量化评价自动避障后的观测系统,定量评价准则能够确保观测系统避障的科学性。针对三维真地表进行观测系统自动排障量化,在大幅提高施工效率的同时,确保了观测系统符合施工要求,从而为采集高品质数据奠定了基础,在单点高密度地震勘探中具有重要作用。
附图说明
图1为本发明的一具体实施例中设计规则观测系统示意图;
图2为本发明的一具体实施例中设计观测系统在障碍物投影的分布图;
图3为本发明的一具体实施例中观测系统自动避障后的炮检点分布图;
图4为本发明的一具体实施例中自动避障前后观测系统属性对比图;
图5为本发明的基于三维真地表的观测系统自动排障量化方法的一具体实施例的流程图。
具体实施方式
应该指出,以下详细说明都是示例性的,旨在对本发明提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本发明的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作和/或它们的组合。
如图5所示,图5为本发明的基于三维真地表的观测系统自动排障优化方法的流程图。
步骤101、根据施工工区的地表条件、构造特征和地质任务,设计规则的观测系统,并形成规则观测系统SPS(Shell Procession Support Format)文件,导入施工工区的卫星图片、近地表速度、Q体和地表高程,结合野外踏勘情况,建立三维真地表构造模型;
所述三维真地表构造模型是由卫星图片、近地表速度与厚度、Q体、地表高程和野外踏勘照片建立的,通过经纬度坐标、低降速带厚度和地表高程真实的反映了地表障碍物的变化情况,同时包含不同位置的吸收衰减特性,实现了由二维平面标量图向三位立体向量图的转变。
步骤102、在三维构造模型中进行野外施工障碍物的自动识别和圈定,结合障碍物特点及安全施工要求,将圈定的障碍物进行分类,并根据分类情况设定激发点施工可行距离;
根据四色图像数字化处理技术,在三维真地表构造模型中自动圈定地表障碍物,地表障碍物信息包括:经纬度坐标、低降速带厚度、吸收衰减特性、地表高程等。结合障碍物特点及安全施工要求,将圈定的障碍物进行分类,障碍物包括:村镇、工厂、河流、公路、养殖场、疏松沉积带、巨厚低速区、井场和高陡山丘等,其中:养殖场、井场、高陡山丘属于I类障碍物;村镇、工厂属于II类障碍物;公路、河流、疏松沉积带、巨厚低速区属于III类障碍物。I类障碍物激发点施工可行距离20m,II类障碍物激发点施工可行距离10m,III类障碍物激发点施工可行距离5m。
步骤103、将理论观测系统投影到真地表三维构造模型中,根据圈定的地表障碍物位置,结合障碍物分类情况,自动选取需要避障的炮点和检波点;
将理论观测系统SPS文件导入,炮点和检波点投影到三维真地表构造模型中,根据圈定的地表障碍物位置,结合障碍物分类情况,自动选取需要避障的炮点和检波点。所述自动选取需要避障的炮点和检波点方法为,在三维向量图中,假定炮点或检波点为P,障碍物范围为一多边形区域Area,同时考虑地表高程的影响,从P点向左侧或右侧水平方向试射射线,射线与多边形区域Area会产生交点数,根据交点个数判断P点是否在多边形区域Area内,大幅提高了准确度和生产效率。若果交点数为奇数或高程突变超出预设阈值,则P点在多边形区域Area内;若交点数为偶数,则P点不在多边形区域Area内。遍历搜索所有炮点和检波点,从而得到落入障碍物中的炮检点,并形成无需避障的观测系统SPS文件、需要避障的观测系统SPS文件。
步骤104、利用自动避障算法,将选取的炮点和检波点进行移动,检波点无需考虑施工可行距离,直接完成自动避障,在进行炮点移动时,根据激发点施工可行距离,将炮点移动到施工可行位置,从而完成炮点的自动避障;
利用自动避障算法,导入需要避障的观测系统SPS文件,将选取的炮点和检波点进行移动。检波点无需考虑施工可行距离,直接完成自动避障。在进行炮点移动时,结合激发点施工可行距离,将炮点移动到施工可行位置,从而完成炮点的自动避障。所述自动避障算法为,读取某一障碍物范围内的炮检点,根据炮检点落入障碍物的位置,按照从中间到两边、从下到上的原则,确定障碍物内每一个炮检点的移动顺序,然后按照规划的移动顺序,上下交替搜索障碍物之外的可放置空点,将炮检点移动到障碍物外侧。多个障碍物的情况下,障碍物计算顺序从上到下依次进行。在进行炮点移动时,结合三类障碍物的激发点施工可行距离,进行炮点的自动避障。
步骤105、输出自动避障后的炮点和检波点位置信息,计算观测系统覆盖次数、炮道密度、方位角、炮检距均匀分布和成像算子属性图,量化评价自动避障后的观测系统;
输出自动避障后的炮点和检波点位置信息,计算观测系统覆盖次数、炮道密度、方位角、炮检距均匀分布和成像算子属性图。提出的量化评价公式为:θ=α×Fold+ρ×ST+β×Zip+φ×RS+λ×IM,其中:α为覆盖次数权重系数,Fold为覆盖次数,ρ为炮道密度权重系数,ST为炮道密度,β为方位角权重系数,Zip为方位角,φ为炮检距均匀分布权重系数,RS为炮检距均匀分布,λ为成像算子权重系数,IM为成像算子。综合考虑观测系统覆盖次数、炮道密度、方位角、炮检距均匀分布和成像算子属性,评价自动避障后的观测系统,根据量化评价结果决定是否需要优化加密炮检点或输出最终SPS文件。
步骤106、根据量化评价结果和定量评价准则,若不满足观测系统施工要求,则进行炮点和检波点的位置优化或加密;若满足观测系统施工要求,则输出最终的观测系统SPS文件,完成了观测系统的自动避障。
在应用本发明的一具体实施例中,包括了以下步骤:
(1)根据施工工区的地表条件、构造特征和地质任务,设计理论的规则观测系统,并形成规则观测系统SPS文件,设计规则观测系统示意图如图1所示。
(2)导入施工工区的卫星图片、近地表速度与厚度、Q体和地表高程,结合野外踏勘情况,建立三维真地表构造模型,该模型是由卫星图片、近地表速度、Q体、地表高程和野外踏勘照片建立的,通过经纬度坐标、低降速带厚度和地表高程真实的反映了地表障碍物的变化情况,同时包含不同位置的吸收衰减特性,实现了由二维平面标量图向三位立体向量图的转变。在三维立体向量图中不仅可以标定地表分布的障碍物,而且可以标定地表高陡地形、地下疏松条带、地下空腔和地下巨厚低速区等特殊障碍物。
(3)在三维真地表构造模型中进行野外施工障碍物的自动识别和圈定,根据四色图像数字化处理技术,在三维真地表构造模型中自动圈定地表障碍物。所述四色图像数字处理技术为,以红色、绿色、黄色为三元基色,对三维真地表构造模型进行图像变化,并以蓝色表示地表高程、低降速带厚度和速度变化情况,从而实现了三维模型的图像数字化处理。在考虑了近地表变化的图像中进行障碍物自动识别和圈定,圈定出来的地表障碍物信息包括:经纬度坐标、低降速带厚度、吸收衰减特性、地表高程等。结合障碍物特点及安全施工要求,将圈定的障碍物进行分类,障碍物包括:村镇、工厂、河流、公路、养殖场、疏松沉积带、巨厚低速区、井场和高陡山丘等。
(4)结合障碍物特点及安全施工要求,将圈定的障碍物进行分类,其中:养殖场、井场、高陡山丘属于I类障碍物;村镇、工厂属于II类障碍物;公路、河流、疏松沉积带、巨厚低速区属于III类障碍物。根据分类情况设定激发点施工可行距离,I类障碍物激发点施工可行距离20m,II类障碍物激发点施工可行距离10m,III类障碍物激发点施工可行距离5m。激发点施工可行距离是通过理论公式模拟和野外试验数据分析得到的,激发点施工可行距离充分考虑了障碍物的特点和施工安全要求,在不同地区的激发点施工可行距离是不同的。
(5)将理论观测系统SPS文件导入,将理论观测系统的炮点和检波点投影到真地表三维构造模型中,根据圈定的地表障碍物位置,结合障碍物分类情况,自动选取需要避障的炮点和检波点,设计观测系统在障碍物投影的分布图如图2所示。
(6)读取SPS文件中的炮点和检波点坐标,在三维向量图中,假定炮点或检波点为P,障碍物范围为一多边形区域Area,同时考虑地表高程的影响,从P点向左侧或右侧水平方向试射射线,射线与多边形区域Area会产生交点数,根据交点个数判断P点是否在多边形区域Area内,大幅提高了准确度和生产效率。若果交点数为奇数或高程突变超出预设阈值,则P点在多边形区域Area内;若交点数为偶数,则P点不在多边形区域Area内。遍历搜索所有炮点和检波点,从而得到落入障碍物中的炮检点,并形成无需避障的观测系统SPS文件、需要避障的观测系统SPS文件。
(7)利用自动避障算法,导入需要避障的观测系统SPS文件,将选取的炮点和检波点进行移动,检波点无需考虑施工可行距离,直接完成自动避障。读取某一障碍物范围内的检波点,根据检波点落入障碍物的位置,按照从中间到两边、从下到上的原则,确定障碍物内每一个检波点的移动顺序,然后按照规划的移动顺序,上下交替搜索障碍物之外的可放置空点,将检波点移动到障碍物外侧。多个障碍物的情况下,障碍物计算顺序从上到下依次进行。
(8)在进行炮点移动时,根据激发点施工可行距离,将炮点移动到施工可行位置,从而完成炮点的自动避障。读取某一障碍物范围内的炮点,根据炮点落入障碍物的位置,按照从中间到两边、从下到上的原则,确定障碍物内每一个炮点的移动顺序,然后按照规划的移动顺序,上下交替搜索障碍物之外的可放置空点,将炮点移动到障碍物外侧。多个障碍物的情况下,障碍物计算顺序从上到下依次进行。在进行炮点移动时,结合三类障碍物的激发点施工可行距离,进行炮点的自动避障,观测系统自动避障后的炮检点分布图如图3所示。
(9)输出自动避障后的炮点和检波点位置信息,计算观测系统覆盖次数、炮道密度、方位角、炮检距均匀分布和成像算子属性图,量化评价自动避障后的观测系统,自动避障前后观测系统属性图对比如图4所示。
(10)观测系统自动避障后的量化评价数学公式为:θ=α×Fold+ρ×ST+β×Zip+φ×RS+λ×IM,其中:α为覆盖次数权重系数,Fold为覆盖次数,ρ为炮道密度权重系数,ST为炮道密度,β为方位角权重系数,Zip为方位角,φ为炮检距均匀分布权重系数,RS为炮检距均匀分布,λ为成像算子权重系数,IM为成像算子。综合考虑观测系统覆盖次数、炮道密度、方位角、炮检距均匀分布和成像算子属性,权重系数根据地下构造和三维速度模型进行权值大小的匹配,利用综合属性评价自动避障后的观测系统,根据量化评价结果决定是否重新避障或输出最终SPS文件。
(11)根据量化评价结果,若不满足观测系统施工和地质任务要求,则则进行炮点和检波点的加密;若满足观测系统施工和地质任务要求,则输出最终的观测系统SPS文件,完成了基于三维真地表的观测系统的自动避障。
本发明中的基于三维真地表的观测系统自动排障量化方法,针对三维真地表进行观测系统自动避障,创新提出了考虑地表高陡地形、地下疏松条带、地下空腔和地下巨厚低速区影响的观测优化方法,在大幅提高施工效率的同时,确保了观测系统符合施工要求,从而为采集高品质数据奠定了基础,在单点高密度地震勘探中具有重要作用。
最后说明的是,以上优选实施例仅用以说明本发明的技术方案而非限制,尽管通过上述优选实施例已经对本发明进行了详细的描述,但本领域技术人员应当理解,可以在形式上和细节上对其作出各种各样的改变,而不偏离本发明权利要求书所限定的范围。
Claims (12)
1.基于三维真地表的观测系统自动排障量化方法,其特征在于,该基于三维真地表的观测系统自动排障量化方法包括:
步骤1,设计规则的观测系统,建立三维真地表构造模型;
步骤2,进行野外施工障碍物的自动识别和圈定,将圈定的障碍物进行分类,设定激发点施工可行距离;
步骤3,根据圈定的地表障碍物位置,自动选取需要避障的炮点和检波点;
步骤4,利用自动避障算法,将选取的炮点和检波点进行移动;
步骤5,输出自动避障后的炮点和检波点位置信息,量化评价自动避障后的观测系统,定量评价准则指导炮点和检波点位置的再优化。
2.根据权利要求1所述的基于三维真地表的观测系统自动排障量化方法,其特征在于,在步骤1中,根据施工工区的地表条件、构造特征和地质任务,设计规则的观测系统,并形成规则观测系统文本资料,导入施工工区的卫星图片、近地表速度与厚度、Q体和地表高程,结合野外踏勘情况,建立三维真地表构造模型,真地表包含了地表、低降速带、高速层顶界面这些地层信息。
3.根据权利要求2所述的基于三维真地表的观测系统自动排障量化方法,其特征在于,在步骤1中,所述三维真地表构造模型是由卫星图片、近地表速度与厚度、Q体、地表高程和野外踏勘照片建立的,通过经纬度坐标、低降速带厚度和地表高程真实的反映了地表障碍物的变化情况,同时包含不同位置近地表地层的吸收衰减特性,实现了由二维平面标量图向三维立体向量图的转变。
4.根据权利要求1所述的基于三维真地表的观测系统自动排障量化方法,其特征在于,在步骤2中,在三维构造模型中进行野外施工障碍物的自动识别和圈定,结合障碍物特点及安全施工要求,将圈定的障碍物进行分类,并根据分类情况设定激发点施工可行距离。
5.根据权利要求4所述的基于三维真地表的观测系统自动排障量化方法,其特征在于,在步骤2中,根据四色图像数字化处理技术,在三维真地表构造模型中自动圈定地表障碍物,地表障碍物信息包括:经纬度坐标、低降速带厚度、吸收衰减特性、地表高程;结合障碍物特点及安全施工要求,将圈定的障碍物进行分类,障碍物包括:村镇、工厂、河流、公路、养殖场、疏松沉积带、巨厚低速区、井场和高陡山丘,其中:养殖场、井场、高陡山丘属于I类障碍物;村镇、工厂属于II类障碍物;公路、河流、疏松沉积带、巨厚低速区属于III类障碍物;I类障碍物激发点施工可行距离20m,II类障碍物激发点施工可行距离10m,III类障碍物激发点施工可行距离5m。
6.根据权利要求1所述的基于三维真地表的观测系统自动排障量化方法,其特征在于,在步骤3中,将规则观测系统投影到真地表三维构造模型中,根据圈定的地表障碍物位置,结合障碍物分类情况,自动选取需要避障的炮点和检波点。
7.根据权利要求6所述的基于三维真地表的观测系统自动排障量化方法,其特征在于,在步骤3中,将规则观测系统文本资料导入,炮点和检波点投影到三维真地表构造模型中,根据圈定的地表障碍物位置,结合障碍物分类情况,自动选取需要避障的炮点和检波点;所述自动选取需要避障的炮点和检波点方法为,在三维立体向量图中,假定炮点或检波点为P,障碍物范围为一多边形区域Area,同时考虑地表高程的影响,从P点向左侧或右侧水平方向试射射线,射线与多边形区域Area会产生交点数,根据交点个数判断P点是否在多边形区域Area内;如果交点数为奇数或高程突变超出预设阈值,则P点在多边形区域Area内;若交点数为偶数,则P点不在多边形区域Area内;遍历搜索所有炮点和检波点,从而得到落入障碍物中的炮检点,并形成无需避障的观测系统资料、需要避障的观测系统资料。
8.根据权利要求1所述的基于三维真地表的观测系统自动排障量化方法,其特征在于,在步骤4中,利用自动避障算法,将选取的炮点和检波点进行移动,检波点无需考虑施工可行距离,直接完成自动避障,在进行炮点移动时,根据激发点施工可行距离,将炮点移动到施工可行位置,从而完成炮点的自动避障。
9.根据权利要求8所述的基于三维真地表的观测系统自动排障量化方法,其特征在于,在步骤4中,利用自动避障算法,导入需要避障的观测系统资料,将选取的炮点和检波点进行移动;检波点无需考虑施工可行距离,直接完成自动避障;在进行炮点移动时,结合激发点施工可行距离,将炮点移动到施工可行位置,从而完成炮点的自动避障;所述自动避障算法为,读取某一障碍物范围内的炮检点,根据炮检点落入障碍物的位置,按照从中间到两边、从下到上的原则,确定障碍物内每一个炮检点的移动顺序,然后按照规划的移动顺序,上下交替搜索障碍物之外的可放置空点,将炮检点移动到障碍物外侧;多个障碍物的情况下,障碍物计算顺序从上到下依次进行;在进行炮点移动时,结合三类障碍物的激发点施工可行距离,进行炮点的自动避障。
10.根据权利要求1所述的基于三维真地表的观测系统自动排障量化方法,其特征在于,在步骤5中,输出自动避障后的炮点和检波点位置信息,计算观测系统覆盖次数、炮道密度、方位角、炮检距均匀分布和成像算子属性图,量化评价自动避障后的观测系统。
11.根据权利要求10所述的基于三维真地表的观测系统自动排障量化方法,其特征在于,在步骤5中,输出自动避障后的炮点和检波点位置信息,计算观测系统覆盖次数、炮道密度、方位角、炮检距均匀分布和成像算子属性图;提出的量化评价公式为:θ=α×Fold+ρ×ST+β×Zip+φ×RS+λ×IM,其中:α为覆盖次数权重系数,Fold为覆盖次数,ρ为炮道密度权重系数,ST为炮道密度,β为方位角权重系数,Zip为方位角,φ为炮检距均匀分布权重系数,RS为炮检距均匀分布,λ为成像算子权重系数,IM为成像算子;综合考虑观测系统覆盖次数、炮道密度、方位角、炮检距均匀分布和成像算子属性,评价自动避障后的观测系统,实现了变观后观测系统的定量分析,根据量化评价结果和量化评价准则决定是否需要优化加密炮检点或输出最终的观测系统资料。
12.根据权利要求1所述的基于三维真地表的观测系统自动排障量化方法,其特征在于,该基于三维真地表的观测系统自动排障量化方法还包括,在步骤5之后,根据量化评价结果和定量评价准则,若不满足观测系统施工要求,则进行炮点和检波点的位置优化或加密;若满足观测系统施工要求,则输出最终的观测系统文本资料,完成观测系统的自动避障。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010390756.7A CN113640869B (zh) | 2020-05-11 | 2020-05-11 | 基于三维真地表的观测系统自动排障量化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010390756.7A CN113640869B (zh) | 2020-05-11 | 2020-05-11 | 基于三维真地表的观测系统自动排障量化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113640869A true CN113640869A (zh) | 2021-11-12 |
CN113640869B CN113640869B (zh) | 2023-12-01 |
Family
ID=78415265
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010390756.7A Active CN113640869B (zh) | 2020-05-11 | 2020-05-11 | 基于三维真地表的观测系统自动排障量化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113640869B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102023310A (zh) * | 2009-09-19 | 2011-04-20 | 中国石油化工集团公司 | 三维观测系统均匀性分析方法 |
CN102236103A (zh) * | 2010-05-07 | 2011-11-09 | 杨世奇 | 基于地理信息的三维地震勘探采集观测系统自动避障优化设计技术 |
CN103605152A (zh) * | 2013-10-30 | 2014-02-26 | 中国石油天然气集团公司 | 一种基于障碍物安全区的自动炮点偏移方法 |
CN104142514A (zh) * | 2013-10-29 | 2014-11-12 | 中国石油化工股份有限公司 | 一种三维地震观测系统定量设计方法 |
WO2015005889A1 (ru) * | 2013-07-08 | 2015-01-15 | Mostoviy Vasyl | Активный мониторинг, способ диагностики, использование |
CN104597493A (zh) * | 2013-10-30 | 2015-05-06 | 中国石油天然气集团公司 | 基于实际地震数据的观测系统变观设计评价方法及装置 |
CN106501841A (zh) * | 2016-09-19 | 2017-03-15 | 中国石油天然气集团公司 | 一种二维宽线地震观测系统的优化设计方法及装置 |
CN107490808A (zh) * | 2017-06-14 | 2017-12-19 | 中国石油化工股份有限公司 | 一种高可靠性地震勘探观测系统的建立方法 |
CN108169818A (zh) * | 2017-12-15 | 2018-06-15 | 西南石油大学 | 一种新型观测系统质量评价方法 |
US20180356544A1 (en) * | 2017-06-08 | 2018-12-13 | Total Sa | Method for evaluating a geophysical survey acquisition geometry over a region of interest, related process, system and computer program product |
CN110927777A (zh) * | 2018-09-19 | 2020-03-27 | 中国石油化工股份有限公司 | 一种三维地震采集观测系统障碍区炮点的移动方法及装置 |
-
2020
- 2020-05-11 CN CN202010390756.7A patent/CN113640869B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102023310A (zh) * | 2009-09-19 | 2011-04-20 | 中国石油化工集团公司 | 三维观测系统均匀性分析方法 |
CN102236103A (zh) * | 2010-05-07 | 2011-11-09 | 杨世奇 | 基于地理信息的三维地震勘探采集观测系统自动避障优化设计技术 |
WO2015005889A1 (ru) * | 2013-07-08 | 2015-01-15 | Mostoviy Vasyl | Активный мониторинг, способ диагностики, использование |
CN104142514A (zh) * | 2013-10-29 | 2014-11-12 | 中国石油化工股份有限公司 | 一种三维地震观测系统定量设计方法 |
CN103605152A (zh) * | 2013-10-30 | 2014-02-26 | 中国石油天然气集团公司 | 一种基于障碍物安全区的自动炮点偏移方法 |
CN104597493A (zh) * | 2013-10-30 | 2015-05-06 | 中国石油天然气集团公司 | 基于实际地震数据的观测系统变观设计评价方法及装置 |
CN106501841A (zh) * | 2016-09-19 | 2017-03-15 | 中国石油天然气集团公司 | 一种二维宽线地震观测系统的优化设计方法及装置 |
US20180356544A1 (en) * | 2017-06-08 | 2018-12-13 | Total Sa | Method for evaluating a geophysical survey acquisition geometry over a region of interest, related process, system and computer program product |
CN107490808A (zh) * | 2017-06-14 | 2017-12-19 | 中国石油化工股份有限公司 | 一种高可靠性地震勘探观测系统的建立方法 |
CN108169818A (zh) * | 2017-12-15 | 2018-06-15 | 西南石油大学 | 一种新型观测系统质量评价方法 |
CN110927777A (zh) * | 2018-09-19 | 2020-03-27 | 中国石油化工股份有限公司 | 一种三维地震采集观测系统障碍区炮点的移动方法及装置 |
Non-Patent Citations (3)
Title |
---|
晁如佑;付英露;石一青;魏燕;成云;李金莲;: "复杂障碍区三维地震观测系统变观设计方法及应用", 复杂油气藏, no. 04 * |
李成博等: "CSI:基于压缩感知的高精度高效率地震资料采集技术", 石油物探 * |
白雪等: "基于单点高密度的宽频宽方位技术在野外采集中的应用", 中国石油学会2019年物探技术研讨会 * |
Also Published As
Publication number | Publication date |
---|---|
CN113640869B (zh) | 2023-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111927552B (zh) | 基于天空地勘察技术的复杂山区长大深埋隧道勘察方法 | |
CN102419453A (zh) | 长导线源瞬变电磁地空探测方法 | |
CN109188506A (zh) | 一种适用于高铁隧底地震ct的纯地表三维观测系统 | |
CN102213769A (zh) | 一种利用三维垂直地震剖面资料确定各向异性参数的方法 | |
CN109031462A (zh) | 用于评估感兴趣的区域上的地球物理勘测采集几何结构的方法 | |
CN106443765A (zh) | 基于多分量观测系统的城市工程地震探测综合成像方法 | |
CN102236103A (zh) | 基于地理信息的三维地震勘探采集观测系统自动避障优化设计技术 | |
CN105093299A (zh) | 一种基于炮检距向量片技术优化观测系统的方法及装置 | |
CN102841379A (zh) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 | |
CN111415413B (zh) | 基于野外露头三维模型的训练图像建立方法 | |
Ismail et al. | Application of combined terrestrial laser scanning and unmanned aerial vehicle digital photogrammetry method in high rock slope stability analysis: A case study | |
US20130083967A1 (en) | System and Method for Extracting Features in a Medium from Data Having Spatial Coordinates | |
CN105866833A (zh) | Vsp-cdp叠加方法及三维vsp覆盖次数计算方法 | |
Hansman et al. | Workflow: From photo-based 3-D reconstruction of remotely piloted aircraft images to a 3-D geological model | |
CN106154315A (zh) | 一种对称均匀的高密度地震采集方法 | |
CN103116183B (zh) | 一种石油地震采集面元覆盖次数属性体切片成图方法 | |
CN113640869B (zh) | 基于三维真地表的观测系统自动排障量化方法 | |
MX2014005103A (es) | Extraccion de datos de cizallamiento sv a partir de datos sismicos de onda p. | |
CN113643430B (zh) | 一种结合地表和地下双因素的观测系统智能变观方法 | |
Gong et al. | Active suspension control based on estimated road class for off-road vehicle | |
Jia et al. | Artificial seismic source field research on the impact of the number and layout of stations on the microseismic location error of mines | |
RU2603856C1 (ru) | Способ проведения прогнозно-поисковых работ месторождений полезных ископаемых на исследуемой площади | |
CN104155690A (zh) | 基于椭球展开的三维地震数据叠加速度求取方法 | |
Aabø et al. | Establishing a high resolution 3D fracture dataset in chalk: possibilities and obstacles working with outcrop data | |
Nurpeisova et al. | Innovative methods of the rock massif fractures survay and treatment of its results |
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 |