CN113408038A - 一种基于数值模拟的地形插值方法及系统 - Google Patents

一种基于数值模拟的地形插值方法及系统 Download PDF

Info

Publication number
CN113408038A
CN113408038A CN202110809141.8A CN202110809141A CN113408038A CN 113408038 A CN113408038 A CN 113408038A CN 202110809141 A CN202110809141 A CN 202110809141A CN 113408038 A CN113408038 A CN 113408038A
Authority
CN
China
Prior art keywords
terrain
numerical simulation
sample point
finite element
interpolation
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
Application number
CN202110809141.8A
Other languages
English (en)
Other versions
CN113408038B (zh
Inventor
郑凯凯
史俊
袁健
于素慧
孙国瑞
王炜
张毅
高成强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Central South University
Rocket Force University of Engineering of PLA
Original Assignee
Harbin Institute of Technology
Central South University
Rocket Force University of Engineering of PLA
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology, Central South University, Rocket Force University of Engineering of PLA filed Critical Harbin Institute of Technology
Priority to CN202110809141.8A priority Critical patent/CN113408038B/zh
Publication of CN113408038A publication Critical patent/CN113408038A/zh
Application granted granted Critical
Publication of CN113408038B publication Critical patent/CN113408038B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Complex Calculations (AREA)

Abstract

一种基于数值模拟的地形插值方法及系统,涉及地球物理模型建立技术领域,所述方法从待表征的地形中随机选取若干样本点,利用有限元方法建立各样本点的空间相关性,进而表征实际地形全貌。本发明提供的方法没有迭代计算过程和解方程问题,实施过程简单清晰,有效提高空间插值的计算效率,并且应用实际采集的样本点三维地理数据构造权重函数能显著提高插值精度。

Description

一种基于数值模拟的地形插值方法及系统
技术领域
本发明涉及地球物理模型建立技术领域,尤其涉及一种基于数值模拟的地形插值方法及系统。
背景技术
空间插值方法广泛应用于从分散采样点获取空间连续数据,并在计算机图形学、环境科学和数值分析等多个学科中得到了应用。实际上,准确估计未测位置的实验数据,对于由工程师做出有效决策并由科学家做出合理的解释具有实际意义。到目前为止,已有30多个空间插值方法得到了发展,包括从全局到局部、从确定性到统计、从单一方法到多种方法组合。但在实际应用中,很难从现有的空间插值方法中选择合适的方法。此外,它们估计的准确性较差,难以满足实际需要。
空间插值方法的估计准确性主要取决于是否充分且合理地利用了所有样本信息。然而,通常的空间插值方法仅关注样本的空间分布和数值大小而忽略了其背后的隐藏信息,导致插值精度不足。例如,常规方法在处理变形和温度的样本数据之间没有区别。此外,常见的空间插值方法通常假定一些数学模型来补充空间相关信息,但他们的确定是经验的和模糊的,插值精度难以保证。因此,很有必要研究一种具有高精度和明确物理意义的新插值方法。实际上,样本背后的特定物理模型包含许多隐藏信息,并且这些模型信息可以通过数值模拟来部分反映。因此,数值模拟和插值的结合可成为上述问题的有效解决方案。
发明内容
有鉴于此,本发明提出了一种基于数值模拟的地形插值方法及系统,没有迭代计算过程和解方程问题,实施过程简单清晰,有效提高空间插值的计算效率,并且应用实际采集的样本点三维地理数据构造权重函数能显著提高插值精度。
为了实现上述目的,本发明采用的技术方案是:
根据本发明的第一方面,提供了一种基于数值模拟的地形插值方法,所述方法从待表征的地形中随机选取若干样本点,利用有限元方法建立各样本点的空间相关性,进而表征实际地形全貌。
进一步的,所述方法包括以下步骤:
S1:在待表征的地形中随机选取若干样本点,并采集所述若干样本点的三维地理数据;
S2:依据所述三维地理数据建立有限元模型,并划分网格;
S3:分别构造各样本点在有限元模型中的离散权重函数,表征各样本点之间的空间相关性;
S4:结合各样本点的离散权重函数与三维地理数据,即可得到待表征地形的差值场,从而表征实际地形。
进一步的,若干样本点均位于所述有限元模型的网格节点上。
进一步的,所述S3具体包括为:
S31:在所述有限元模型中的样本点i处施加单位一的垂直于地面方向的位移,并限制其他样本点的刚体位移;i为正整数,用于表示样本点序号;
S32:通过数值模拟方式得到有限元模型中全部网格节点的垂直于地面方向位移的模拟值,即为样本点i的离散权重函数;
S33:重复上述步骤S31~S32,直至获得全部样本点的离散权重函数。
进一步的,所述离散权重函数具体为:
Ni={u1,u2,...,un};
其中,Ni为样本点i的离散权重函数;u1,u2,...,un为有限元模型中第1个节点、第2个节点至第n个节点处的垂直于地面方向的位移模拟值;n为正整数,取有限元模型中的网格节点总数。
进一步的,样本点的离散权重函数随空间位置的变化而改变。
进一步的,所述S4具体包括:
将各样本点的三维地理数据Ui代入公式
Figure BDA0003167515920000021
即可计算得到待表征地形的插值场,从而表征实际地形;
其中,m为样本点总数;i为正整数,用于表示样本点序号;Ni为样本点i的离散权重函数;Ui为样本点i的应变;U为待表征地形的插值场。
进一步的,所述样本点处的有限元网格比非样本点处的有限元网格密集。
进一步的,所述三维地理数据包括样本点在地心坐标系或参心坐标系下的坐标数据。
根据本发明的第二方面,提供了一种基于数值模拟的地形插值系统,包括:
处理器和用于存储可执行指令的存储器;
其中,所述处理器被配置为执行所述可执行指令,以执行上述的基于数值模拟的地形插值方法。
根据本发明的第三方面,提供了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现上述的基于数值模拟的地形插值方法。
相对于现有技术,本发明所述的一种基于数值模拟的地形插值方法及系统,具有如下优势:
(1)本发明提出的基于数值模拟方法构建的权重函数,构造过程简单,对地形地貌上的样本点数量和分布没有限制,避免传统方法中繁琐的地形数据收集过程;
(2)本发明实施过程中没有解方程和优化问题,有效减少了插值的繁杂程度,才做简单,方便推广应用;
(3)本发明提出的地形插值方法具有精确性、全局性、非凸性、可外插性等优异特性;
(4)基于数值模拟的插值方法赋予了插值明确的物理意义,即满足地形物理模型上的最小势能原理,能显著减低地形插值畸变,并提高插值的精度。
附图说明
说明书附图用来提供对本发明的进一步理解,构成本发明的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。
图1是本发明实施例所述方法的流程图;
图2是本发明实施例地形插值中的有限元模型和离散权重函数示意图;
图3是本发明实施例地形插值的最终效果图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本发明的一些方面相一致的装置和方法的例子。
本发明的说明书和权利要求书中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例,例如能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
多个,包括两个或者两个以上。
和/或,应当理解,对于本发明中使用的术语“和/或”,其仅仅是一种描述关联对象的关联关系,表示可以存在三种关系。例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。
一种基于数值模拟的地形插值方法,包括:
a、建立插值空间的有限元模型,并将其划分为适当多的单元;
b、施加特定的边界条件,用数值模拟方法来获得各样本点的离散权重函数Ni={u1,u2,...,uk,...,un}来表征测点之间的空间相关性,其中uk为第k个节点处的模拟值,n为节点总数;
c、将t时刻各样本点的数值Ui代入公式
Figure BDA0003167515920000041
即可计算得到插值场U={U1,U2,...,Uk,...,Un},其中m为样本点总数。
优选地,步骤a中所有样本点应正好处于单元网格节点上,测点处的网格应划分较密。
优选地,步骤b中所述的特定边界条件是指在计算第i个样本点处的权重函数Ni时,在该点处加单位一的位移,而在其余样本点处设为零。
优选地,步骤b所述的数值模拟是用有限元软件对模型进行分析,不考虑材料的非线性,模拟所得各节点的物理量数值即组成样本点i的离散权重函数Ni
优选地,步骤c中,U为插值场的离散表达,Ni是关于空间坐标的离散函数,与时间无关,即不同时刻同一个样本点的权重函数相同。
优选地,用ANSYS-APDL语言编制自动循环加约束的程序,并批量输出离散权重函数向量。
优选地,本发明特别适用于具有传递特性的物理量,如变形、速度、温度、浓度等,其当物系偏离平衡状态时,就会发生转移现象使物系趋向平衡状态。
实施例:
由于地形地势是地壳受构造应力的变形结果,因此可用板壳的位移来表征地形的起伏,物理意义比较接近。在某地形地貌中随机选取十个位置作为样本点,形成插值空间,且高程样本数据也是一列0到1的随机浮点数。如图1所示,按照以下步骤构建这十个测点的权重函数并计算地形插值场:
a、首先用大型通用软件ANSYS建立插值空间大小的薄板有限元模型,尺寸为2×1km2,如图1所示。将其划分为40×20个Shell 181单元,任意设定材料的弹性模量为1000MPa,不考虑大变形和弹塑性。薄板的厚度设为10m,其厚度越薄则插值效果越平滑。
b、先以构造样本点3处的权重函数为例。根据权重函数的构造性质,在样本点3处加z向1m的位移,在其他样本点处加z向的固定约束。另外,为了限制板其它方向的刚体位移,在样本点1处加x向和y向的约束,在样本点2处加y向的约束。然后进行静力分析,得到的z向位移场即为离散权重函数N1={u1,u2,...,uk,...,un},如图2中的a所示。
c、类似地,可以得到其它的离散权重函数Ni(i=1,2,…,10)。其中N8如图2中的b所示。
d、在得到各样本点离散权重函数后,然后将各样本点的数值Ui代入公式
Figure BDA0003167515920000051
可得到地形的插值场,从而展示完整地貌,如图3所示。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护之内。

Claims (10)

1.一种基于数值模拟的地形插值方法,其特征在于,所述方法从待表征的地形中随机选取若干样本点,利用有限元方法建立各样本点的空间相关性,进而表征实际地形全貌。
2.根据权利要求1所述的一种基于数值模拟的地形插值方法,其特征在于,所述方法具体包括:
S1:在待表征的地形中随机选取若干样本点,并采集所述若干样本点的三维地理数据;
S2:依据所述三维地理数据建立有限元模型,并划分网格;
S3:分别构造各样本点在有限元模型中的离散权重函数,表征各样本点之间的空间相关性;
S4:结合各样本点的离散权重函数与三维地理数据,即可得到待表征地形的差值场,从而表征实际地形全貌。
3.根据权利要求2所述的一种基于数值模拟的地形插值方法,其特征在于,若干样本点均位于所述有限元模型的网格节点上。
4.根据权利要求2所述的一种基于数值模拟的地形插值方法,其特征在于,所述S3具体包括为:
S31:在所述有限元模型中的样本点i处施加单位一的垂直于地面方向的位移,并限制其他样本点的刚体位移;i为正整数,用于表示样本点序号;
S32:通过数值模拟方式得到有限元模型中全部网格节点的垂直于地面方向位移的模拟值,即为样本点i的离散权重函数;
S33:重复上述步骤S31~S32,直至获得全部样本点的离散权重函数。
5.根据权利要求4所述的一种基于数值模拟的地形插值方法,其特征在于,所述离散权重函数具体为:
Ni={u1,u2,...,un};
其中,Ni为样本点i的离散权重函数;u1,u2,...,un为有限元模型中第1个节点、第2个节点至第n个节点处的垂直于地面方向的位移模拟值;n为正整数,取有限元模型中的网格节点总数。
6.根据权利要求2所述的一种基于数值模拟的地形插值方法,其特征在于,所述S4具体包括:
将各样本点的三维地理数据Ui代入公式
Figure FDA0003167515910000021
即可计算得到待表征地形的插值场,从而表征实际地形;
其中,m为样本点总数;i为正整数,用于表示样本点序号;Ni为样本点i的离散权重函数;Ui为样本点i的应变;U为待表征地形的插值场。
7.根据权利要求2所述的一种基于数值模拟的地形插值方法,其特征在于,所述样本点处的有限元网格比非样本点处的有限元网格密集。
8.根据权利要求2所述的一种基于数值模拟的地形插值方法,其特征在于,所述三维地理数据包括样本点在地心坐标系或参心坐标系下的坐标数据。
9.一种基于数值模拟的地形插值系统,其特征在于,包括:
处理器和用于存储可执行指令的存储器;
其中,所述处理器被配置为执行所述可执行指令,以执行如权利要求1至8中任一项所述的基于数值模拟的地形插值方法。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至8中任一项所述的基于数值模拟的地形插值方法。
CN202110809141.8A 2021-07-16 2021-07-16 一种基于数值模拟的地形插值方法及系统 Active CN113408038B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110809141.8A CN113408038B (zh) 2021-07-16 2021-07-16 一种基于数值模拟的地形插值方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110809141.8A CN113408038B (zh) 2021-07-16 2021-07-16 一种基于数值模拟的地形插值方法及系统

Publications (2)

Publication Number Publication Date
CN113408038A true CN113408038A (zh) 2021-09-17
CN113408038B CN113408038B (zh) 2022-07-05

Family

ID=77686735

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110809141.8A Active CN113408038B (zh) 2021-07-16 2021-07-16 一种基于数值模拟的地形插值方法及系统

Country Status (1)

Country Link
CN (1) CN113408038B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110106507A1 (en) * 2009-10-30 2011-05-05 Igeoss Method for building a depositional space corresponding to a geological domain
US20130096890A1 (en) * 2011-10-13 2013-04-18 William Brian Vanderheyden Material point method modeling in oil and gas reservoirs
US20150066460A1 (en) * 2013-08-30 2015-03-05 Jimmy Klinger Stratigraphic function
CN112307640A (zh) * 2020-11-11 2021-02-02 中南大学 基于自然单元-无限元的三维多源自然电位数值模拟方法
CN112800657A (zh) * 2021-04-15 2021-05-14 中南大学 基于复杂地形的重力场数值模拟方法、装置和计算机设备

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110106507A1 (en) * 2009-10-30 2011-05-05 Igeoss Method for building a depositional space corresponding to a geological domain
US20130096890A1 (en) * 2011-10-13 2013-04-18 William Brian Vanderheyden Material point method modeling in oil and gas reservoirs
US20150066460A1 (en) * 2013-08-30 2015-03-05 Jimmy Klinger Stratigraphic function
CN112307640A (zh) * 2020-11-11 2021-02-02 中南大学 基于自然单元-无限元的三维多源自然电位数值模拟方法
CN112800657A (zh) * 2021-04-15 2021-05-14 中南大学 基于复杂地形的重力场数值模拟方法、装置和计算机设备

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
于国兴 等: "基于有限元的混凝土面板堆石坝三维渗流场研究", 《决策探索(中)》, 30 November 2018 (2018-11-30), pages 40 - 42 *
张毅 等: "西昌市耕地土壤有效钼空间变异特征及影响因素分析", 《湖北农业科学》, vol. 51, no. 13, 31 July 2012 (2012-07-31), pages 2701 - 2705 *

Also Published As

Publication number Publication date
CN113408038B (zh) 2022-07-05

Similar Documents

Publication Publication Date Title
Hötzer et al. The parallel multi-physics phase-field framework Pace3D
US6208982B1 (en) Method and apparatus for solving complex and computationally intensive inverse problems in real-time
Nightingale et al. Adaptive semi-linear inversion of strong gravitational lens imaging
JP6784780B2 (ja) 大規模再生可能エネルギーのデータについて確率モデルを構築する方法
CN108563837A (zh) 一种冲积河流水沙模型的模型参数实时校正方法和系统
CN106683185A (zh) 一种基于大数据的高精度曲面建模方法
CN112257021A (zh) 一种可选的克里金空间插值降雨量估算方法
CN109033181B (zh) 一种复杂地形地区风场地理数值模拟方法
CN110110406B (zh) 一种基于Excel计算平台实现LS-SVM模型的边坡稳定性预测方法
CN114692523A (zh) 基于图卷积的自适应高维流体动力方程的流速预测方法
CN113408038B (zh) 一种基于数值模拟的地形插值方法及系统
Lingsch et al. Vandermonde neural operators
CN116861822A (zh) 基于笛卡尔网格的物面边界的处理方法及装置
CN115795402B (zh) 一种基于变分法的多源降水数据融合方法和系统
CN104376120A (zh) 一种信息检索方法及系统
CN107247828A (zh) 一种基于逆kriging函数的结构有限元模型修正方法
Ding et al. Accelerating multi‐dimensional interpolation using moving least‐squares on the GPU
Lingsch et al. A structured matrix method for nonequispaced neural operators
CN111210877A (zh) 一种推断物性参数的方法及装置
CN113408040B (zh) 一种土木工程中的模拟数据修正方法及系统
CN115470675A (zh) 一种重载压力容器裂纹扩展路径的预测方法及系统
CN105426626B (zh) 基于相似数据样式聚类的多点地质统计学建模方法
DiPrete et al. Reinforcement Learning for Block Decomposition of CAD Models
CN109960776B (zh) 一种用于水力走时和水力信号衰减反演计算的改进方法
Cobos et al. Rangemap: an R package to explore species' geographic ranges

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 710025 No. 2 Tongxin Road, Baqiao District, Xi'an City, Shaanxi Province

Applicant after: Rocket Force University of Engineering of PLA

Applicant after: HARBIN INSTITUTE OF TECHNOLOGY

Applicant after: CENTRAL SOUTH University

Address before: 410083 Hunan province Changsha Lushan Road No. 932

Applicant before: CENTRAL SOUTH University

Applicant before: HARBIN INSTITUTE OF TECHNOLOGY

Applicant before: Rocket Force University of Engineering of PLA

CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Gao Chengqiang

Inventor after: Yuan Jian

Inventor after: Wang Wei

Inventor after: Yu Suhui

Inventor after: Zhang Yi

Inventor after: Shi Jun

Inventor after: Zheng Kaikai

Inventor before: Zheng Kaikai

Inventor before: Shi Jun

Inventor before: Yuan Jian

Inventor before: Yu Suhui

Inventor before: Sun Guorui

Inventor before: Wang Wei

Inventor before: Zhang Yi

Inventor before: Gao Chengqiang

GR01 Patent grant
GR01 Patent grant