CN103207410A - 一种针对崎岖海底的混合网格模型建立方法 - Google Patents
一种针对崎岖海底的混合网格模型建立方法 Download PDFInfo
- Publication number
- CN103207410A CN103207410A CN2013101340398A CN201310134039A CN103207410A CN 103207410 A CN103207410 A CN 103207410A CN 2013101340398 A CN2013101340398 A CN 2013101340398A CN 201310134039 A CN201310134039 A CN 201310134039A CN 103207410 A CN103207410 A CN 103207410A
- Authority
- CN
- China
- Prior art keywords
- grid
- rugged
- seabed
- rugged seabed
- physical property
- 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
Images
Abstract
本发明涉及一种针对崎岖海底的混合网格模型建立方法,包括以下步骤:1)读入崎岖海底的水深、计算区域的尺度、介质物性参数分布信息;2)根据崎岖海底的水深、计算区域的尺度、介质物性参数分布信息构建出整个模型框架;3)采用正方形网格剖分步骤2)构建出的整个模型框架;4)采用三角网格剖分崎岖海底附近区域;5)根据网格节点邻近网格单元介质的物性参数分布、几何面积信息,采用面积加权平均法计算网格节点上的物性参数;6)将网格节点上的物性参数、计算参数信息存储并输出。本发明能够更准确地刻画崎岖海底的形态,提高计算效率,保证算法的稳定性,同时可以避免边界绕射和角点散射等问题。本发明可以广泛用于崎岖海底条件下的地震正演模拟、层析成像、偏移成像、速度建模等领域中。
Description
技术领域
本发明涉及一种地震勘探模型建立方法,特别是关于一种针对崎岖海底的混合网格模型建立方法。
背景技术
崎岖海底是海洋地震勘探普遍会遇到的问题。针对崎岖海底的准确精细的模型建立方法,是研究崎岖海底条件下的地震正演模拟、数据处理等技术的基础性工具。因此很有必要专门研究针对崎岖海底的模型建立方法。而模型建立方法需要解决的核心问题是模型的网格剖分问题。在数值计算问题上,对崎岖海底模型进行网格剖分涉及的关键技术问题是对不规则计算边界的描述。在处理类似于崎岖海底引起的不规则计算边界的网格剖分问题时,现有技术主要有正方形阶梯网格和曲线网格两类。但是:
如图5所示,采用正方形阶梯网格剖分不规则计算边界。虽然能让整个计算区域均为正方形网格剖分,但是从图中可以看出,阶梯状的正方形网格对模型的刻画非常粗糙,对崎岖海底形状的刻画也非常不准确。同时从数值计算的角度看,因为正方形阶梯状网格存在很多界面处的角点,并且整个界面不连续且跳变很严重,所以易引起数值计算时的边界绕射和角点散射等问题;采用正方形阶梯状网格剖分崎岖海底模型,将对不规则计算边界的形态刻画过于粗糙,并且阶梯状的网格还会引起数值计算时的边界绕射和角点散射等问题。
如图6所示,采用曲网格中的正交贴体网格剖分不规则计算边界。虽然不对不规则界面作任何近似,但是生成这样的网格剖分需要求解线性方程组,并且在求解三维问题时这样的线性方程组非常庞大,由此产生的网格生成消耗的计算量就更大,即需要付出更多的额外的网格生成成本。同时从数值计算的角度看,曲网格的非均匀性很强,直接在此网格中进行数值计算必然会遇到一些额外的困难。当然曲网格技术的数值求解基本均在曲坐标系下进行,不过曲坐标系下的多数物性方程均比笛卡尔坐标系下的方程复杂很多,更大的问题是这些方程(比如说Eikonal方程,此方程为非线性方程)在数值实现时,由于受网格生成质量的影响而很难保证数值的稳定性。
发明内容
针对上述问题,本发明的目的是提供一种既能准确精细刻画崎岖海底形状,又能保证数值计算的绝大多数计算是在规则正方形网格中进行的针对崎岖海底的混合网格模型建立方法。
为实现上述目的,本发明采取以下技术方案:一种针对崎岖海底的混合网格模型建立方法,包括以下步骤:1)读入崎岖海底的水深、计算区域的尺度、介质物性参数分布信息;2)根据崎岖海底的水深、计算区域的尺度、介质物性参数分布信息构建出整个模型框架;3)采用正方形网格剖分步骤2)构建出的整个模型框架;4)采用三角网格剖分崎岖海底附近区域;5)根据网格节点邻近网格单元介质的物性参数分布、几何面积信息,采用面积加权平均法计算网格节点上的物性参数;6)将网格节点上的物性参数、计算参数信息存储并输出。
本发明由于采取以上技术方案,其具有以下优点:1、本发明由于采用混合网格方式对崎岖海底模型进行剖分,因此绝大部分区域是采用正方形网格剖分,仅在崎岖海底附近采用精细三角网格剖分,因此能够更准确地刻画崎岖海底的形态。2、本发明由于采用混合网格对崎岖海底模型进行剖分,因此既能够使数值计算的绝大部分工作在规则的正方形网格中进行,提高计算效率,保证算法的稳定性。3、本发明由于在崎岖海底附近采用了三角网络,因此能够建立起三角网络的界面点与正方形网格节点的位置关系,进而可以有效地避免边界绕射和角点散射等问题。本发明可以被广泛地应用于崎岖海底条件下的地震正演模拟、层析成像、偏移成像、速度建模等领域中。
附图说明
图1是本发明基于混合网格的崎岖海底模型建立流程示意图
图2是本发明混合网格法剖分琼东南盆地北部宝岛凹陷坡折水深剧变模型
图3是本发明混合网格法剖分白云凹陷垂直于滑坡走向的滑坡根部断块凹谷模型
图4是采用本发明混合网格对不规则计算边界的剖分示意图
图5是已有技术正方形阶梯网格剖分不规则计算边界示意图
图6是已有技术曲网格中的正交贴体网格剖分不规则计算边界示意图
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明针对崎岖海底的混合网格模型建立方法,主要包括以下步骤:
1)读入崎岖海底的水深、计算区域的尺度、介质的物性参数分布等信息;
2)根据步骤1)读入的信息构建出整个模型框架;
3)采用正方形网格剖分步骤1)构建出的整个模型框架;
4)采用三角网格剖分崎岖海底附近区域;
5)根据网格节点邻近网格单元介质的物性参数分布、几何面积等信息,采用面积加权平均法计算网格节点上的物性参数;
6)根据各数值模拟、数据处理等技术的需要,将网格节点上的物性参数、计算参数等信息存入对应数组中,并输出数据以备数值模拟、数据处理等技术使用。
上述步骤5)中的物性参数是指介质的密度、地震波速度等。
为了更好的说明本发明的特点,列举以下具体实施例:
实施例1
如图2所示,根据已发表资料《琼东南盆地深水区地震资料品质影响因素探析》(邓勇,李列,柴继堂等,中国海上油气,2010,22(6):381-386),本发明建立了一个琼东南盆地北部宝岛凹陷坡折水深剧变模型。图中显示了模型介质的速度参数分布,基于本发明的模型建立方法,采用22×10(实际数值计算时远远比这个网格数多,此处为了显示方便采用了较少的网格,以下相同,不再赘述)的网格对模型进行了剖分。
实施例2:
如图3所示,根据已发表资料《南海北部白云大型海底滑坡的几何形态与变形特征》(孙运宝,吴时国,王志君等,海洋地质与第四纪地质,2008,28(6):69-77)本发明建立了一个白云凹陷垂直于滑坡走向的滑坡根部断块凹谷模型。图中的显示了模型介质的密度参数分布,基于本发明的模型建立方法,采用22×12的网格对模型进行了剖分。
如图4所示,本发明采用混合网格方法对不规则计算边界的剖分,首先通过与正方形阶梯网格剖分(如图5所示)对比,可以发现:正方形阶梯网格对崎岖海底的形状刻画得非常粗糙,并且崎岖海底的形状也由光滑的曲线变为了具有很多角点的阶梯状折线;其次通过与曲网格剖分(如图6所示)对比,曲网格在计算区域的很大一部分存在明显的不均匀性(即网格尺度大小不一),并且网格的很多地方存在明显的不完全正交性(即网格线在相交处不垂直)。与上述两类网格不同,如图4所示,本发明采用混合网格的绝大部分区域为正方形网格剖分,并且崎岖海底附近的精细三角网格能够更准确地刻画崎岖海底的形态。
此外,在已有技术的两类网格和本发明混合网格中,从数值计算的角度上可以看出:首先,因为正方形阶梯状网格存在很多界面处的角点,并且整个界面不连续且跳变很严重,所以易引起数值计算时的边界绕射和角点散射等问题;其次,曲网格存在不完全正交性和非均匀性等问题,直接在此网格中进行数值计算必然会遇到一些额外的困难,当然可以在曲坐标系下进行相应曲网格的数值计算,但是曲坐标系下的基本物性方程均比笛卡尔坐标系下的方程复杂很多,更大的问题是这些方程(比如说Eikonal方程,此方程为非线性方程)在数值实现时,由于受网格生成质量的影响而很难保证数值稳定性。最后,与以上两种方法不同,本发明的混合网格首先能保证数值计算的绝大部分工作在规则的正方形网格中进行,这样能很好地保证算法的稳定性,并且三角网络还建立起三角网络的界面点和正方形网格节点的位置关系,这样不易产生边界绕射和角点散射等问题。
上述各实施例仅用于说明本发明,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (2)
1.一种针对崎岖海底的混合网格模型建立方法,包括以下步骤:
1)读入崎岖海底的水深、计算区域的尺度、介质物性参数分布信息;
2)根据崎岖海底的水深、计算区域的尺度、介质物性参数分布信息构建出整个模型框架;
3)采用正方形网格剖分步骤2)构建出的整个模型框架;
4)采用三角网格剖分崎岖海底附近区域;
5)根据网格节点邻近网格单元介质的物性参数分布、几何面积信息,采用面积加权平均法计算网格节点上的物性参数;
6)将网格节点上的物性参数、计算参数信息存储并输出。
2.如权利要求1所述的一种针对崎岖海底的混合网格模型建立方法,其特征在于:所述步骤5)中的所述物性参数分布包括密度和地震波速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310134039.8A CN103207410B (zh) | 2013-04-17 | 2013-04-17 | 一种针对崎岖海底的混合网格模型建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310134039.8A CN103207410B (zh) | 2013-04-17 | 2013-04-17 | 一种针对崎岖海底的混合网格模型建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103207410A true CN103207410A (zh) | 2013-07-17 |
CN103207410B CN103207410B (zh) | 2015-12-09 |
Family
ID=48754685
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310134039.8A Active CN103207410B (zh) | 2013-04-17 | 2013-04-17 | 一种针对崎岖海底的混合网格模型建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103207410B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104535391A (zh) * | 2014-12-30 | 2015-04-22 | 中国科学院地质与地球物理研究所 | 一种基于层状地质模型的地球物理数据处理方法 |
CN105513131A (zh) * | 2015-12-16 | 2016-04-20 | 中国地质大学(武汉) | 一种顾及褶皱的三维地质结构模型角点网格剖分方法 |
CN111538083A (zh) * | 2020-06-08 | 2020-08-14 | 吉林大学 | 基于速度梯度的崎岖海底界面的光滑处理方法 |
CN112649859A (zh) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | 一种地震波速度自适应无网格场节点建立方法及系统 |
CN112684504A (zh) * | 2020-12-03 | 2021-04-20 | 中国矿业大学 | 一种基于全散射模型的城市地下空洞快速探测方法 |
CN112684499A (zh) * | 2020-12-03 | 2021-04-20 | 南方科技大学 | 一种地震波仿真分析方法、存储介质及设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6954698B2 (en) * | 2002-06-28 | 2005-10-11 | Gedex Inc. | System and method for surveying underground density distributions |
CN102183790A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 基于时空双变网格的弹性波正演模拟技术 |
-
2013
- 2013-04-17 CN CN201310134039.8A patent/CN103207410B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6954698B2 (en) * | 2002-06-28 | 2005-10-11 | Gedex Inc. | System and method for surveying underground density distributions |
CN102183790A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 基于时空双变网格的弹性波正演模拟技术 |
Non-Patent Citations (4)
Title |
---|
孙卫涛 等: "地形构造中地震波传播的非对称交错网格模拟", 《应用数学和力学》 * |
孙建国: "复杂地表条件下地球物理场数值模拟方法评述", 《世界地质》 * |
张永刚: "地震波场数值模拟方法", 《石油物探》 * |
罗延钟 等: "二维地形不平条件下均匀外电场的有限差分模拟", 《物化探计算技术》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104535391A (zh) * | 2014-12-30 | 2015-04-22 | 中国科学院地质与地球物理研究所 | 一种基于层状地质模型的地球物理数据处理方法 |
CN104535391B (zh) * | 2014-12-30 | 2017-01-04 | 中国科学院地质与地球物理研究所 | 一种基于层状地质模型的地球物理数据处理方法 |
CN105513131A (zh) * | 2015-12-16 | 2016-04-20 | 中国地质大学(武汉) | 一种顾及褶皱的三维地质结构模型角点网格剖分方法 |
CN105513131B (zh) * | 2015-12-16 | 2016-12-14 | 中国地质大学(武汉) | 一种顾及褶皱的三维地质结构模型角点网格剖分方法 |
CN112649859A (zh) * | 2019-10-12 | 2021-04-13 | 中国石油化工股份有限公司 | 一种地震波速度自适应无网格场节点建立方法及系统 |
CN112649859B (zh) * | 2019-10-12 | 2024-03-22 | 中国石油化工股份有限公司 | 一种地震波速度自适应无网格场节点建立方法及系统 |
CN111538083A (zh) * | 2020-06-08 | 2020-08-14 | 吉林大学 | 基于速度梯度的崎岖海底界面的光滑处理方法 |
CN112684504A (zh) * | 2020-12-03 | 2021-04-20 | 中国矿业大学 | 一种基于全散射模型的城市地下空洞快速探测方法 |
CN112684499A (zh) * | 2020-12-03 | 2021-04-20 | 南方科技大学 | 一种地震波仿真分析方法、存储介质及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN103207410B (zh) | 2015-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103207410A (zh) | 一种针对崎岖海底的混合网格模型建立方法 | |
Mandli et al. | Adaptive mesh refinement for storm surge | |
CN102147479B (zh) | 一种储层空间物性参数的建模方法 | |
CN107180450A (zh) | 一种基于dem的河谷横断面形态的算法 | |
CN103886135B (zh) | 基于二维非恒定流数值模型的电力工程选址方法 | |
EP3189357A1 (en) | Sediment transport simulation with parameterized templates for depth profiling | |
CN114218840A (zh) | 河口航道水沙运动及其地形演变整体建模及可视化系统 | |
Deng et al. | Morphogenetic modelling of coastal and estuarine evolution | |
CN104091065A (zh) | 一种求解浅水问题模拟间断水流数值的方法 | |
CN110400371A (zh) | 一种水平构造地貌实体的三维模型构建方法 | |
CN110990926B (zh) | 一种基于面积修正率的城市地表建筑水动力学仿真方法 | |
Raabe et al. | MEPBAY and SMC: software tools to support different operational levels of headland-bay beach in coastal engineering projects | |
CN116774292A (zh) | 一种地震波走时确定方法、系统、电子设备及存储介质 | |
Yang et al. | The impact of coastal reclamation on tidal and storm surge level in Sanmen Bay, China | |
CN114692471A (zh) | 一种岩溶地下水系统流网模拟方法 | |
CN103376463A (zh) | 一种基于断层控制的反演建模方法 | |
CN102567594A (zh) | 一种近岸岛礁型人工鱼礁群流场仿真建模的方法 | |
CN103217715B (zh) | 多尺度规则网格层析反演静校正方法 | |
Pender et al. | 12 Flood Inundation Modelling to Support Flood Risk Management | |
Ma et al. | A high-precision hydrodynamic model coupled with the hydrological habitat suitability model to reveal estuarine vegetation distribution | |
CN105204090A (zh) | 复杂形态砂体的布尔模拟方法 | |
CN113327323B (zh) | 基于散点数据的水体环境地形构建方法 | |
CN114925624A (zh) | 一种天然河道三维水流数值模拟方法 | |
Nicholas | Roughness parameterization in CFD modelling of gravel‐bed rivers | |
CN114544286A (zh) | 一种基于粒径分布的颗粒随机生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder | ||
CP01 | Change in the name or title of a patent holder |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee after: CNOOC research institute limited liability company Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee before: CNOOC Research Institute Patentee before: China National Offshore Oil Corporation |