收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

三维数值流形方法研究及其在地学中的初步应用

武艳强  
【摘要】:伴随地震孕育、发生与震后调整过程的应力、应变问题在地震研究中占有重要地位。以GPS为代表的空间观测技术为大尺度地壳运动监测提供了高效、稳定、精确的测量结果,并为研究上述应力、应变问题提供了有效数据约束。GPS速度场和应变率场结果既可用于变形分析又可作为数值模拟的有效约束和检验,为了获取研究区域可靠的应变率场结果,有必要开展GPS应变率场计算方法的对比研究。另一方面,由于在整个地震过程中应力、应变存在动态演化现象,因此需要研究适合于动态资料约束、大变形模拟、连续与非连续变形耦合的三维数值流形方法。基于上述分析,本文包括如下三个方面研究内容。 第一、通过模拟与实际数据分析,讨论了四种连续应变率场计算方法的解算精度与抗差性等问题,发展了能够较客观反映实际变形分布的最小二乘配置球面应变率解算方法,并与其它几种常用方法进行了对比分析。 首先,对大空间尺度(东经75°~135°,北纬20°~50°)模拟数据计算结果的分析表明,采用1°×1°采样数据及其50%限定(对1°×1°数据进行50%数据量的离散化采样,并剔除2个5°×10°区域)采样数据作为输入的情况下,Delaunay三角形方法因噪声对解算结果影响太大不可取,其它三种整体应变率场解算结果具有一致性,但抗差性有所差别。通过分析理论结果与附加了不同误差的计算结果的相关系数表明,抗差性由好到坏排列如下:最小二乘配置方法、球谐函数方法、多面函数方法和Delaunay三角形方法。从输入数据的稀疏程度对不同应变率计算方法影响程度看,在数据采样率介于2°~1°网格之间时最小二乘配置方法受数据稀疏的影响最小,球谐函数和多面函数对输入数据的密集性要求较高。 其次,中等空间尺度(东经90°~120°,北纬25°~40°)模拟数据(1°~0.5°网格)结果表明,三种整体方法在测点分布足够密的情况下均能满足实际计算需求。此时,三种方法对附加的误差敏感程度有一定差别但不显著,对比而言最小二乘配置稍强于其他两种方法。通过对不同空间采样数据的应变率计算结果的分析表明,随着输入数据越来越稀疏,多面函数和球谐函数方法计算结果与理论值的相关性减弱的幅度要快于最小二乘配置方法,表明此两种方法对数据的分布密度要求较高。 最后,通过对中国大陆1999~2004期应变率计算结果的分析,发现最小二乘配置方法的稳定性最高,即使仅用50%的数据作为输入得到的应变率结果与全部数据输入基本一致,误差水平也没有明显的变化。球谐函数方法受点位分布稀疏程度的影响较大,表现为边缘效应的量值有所增大,影响范围也有所扩大。多面函数方法受测点稀疏程度的影响较大,表现为计算的不稳定性。另外,多面函数方法和球谐函数方法受点位分布的几何形状影响较大,比如中国东北地区应变率计算结果的失真。 因此,总体而言,从抗差性、边缘效应、误差分布、稳定性角度来看最小二乘配置方法最佳。究其原因在于最小二乘配置的协方差函数是对实际观测数据统计计算得到的,能够反映数据的真实分布特征,而球谐函数的展开阶数,多面函数的光滑因子、核函数、平差结点等的选择都是通过反复试算得到的,很难做到对输入数据的最优描述。在输入数据分布密度满足的情况下,最小二乘配置方法的计算过程无需人工进行参数选择,即使不同人员进行解算也能得到一致的结果。 第二、通过单纯形积分特性研究、具体的矩阵元素表达式推导、模拟结果的理论测试等过程,从多角度对三维数值流形方法开展了研究。 首先,通过对单纯形积分公式的分析,给出了二维和三维单纯形积分的C++算法实现,并结合具体实例描述了单纯形积分的求解过程。在对中空规则形状积分区域的面积(体积)和重心的计算值与理论值对比分析的基础上,对单纯形积分的精度进行了讨论。通过对比不同边长比(10~(-6)~10~6)情况下计算值与理论值的差异,讨论了积分区域图形条件对积分结果的影响;通过对中空非规则形状积分区域计算结果的分析,讨论了单纯形积分的普适性特征。最后,分析了高阶多项式在中空非规则形状积分区域的积分精度。结果表明,积分项为低价或高价多项式情况下计算值与理论值的相对误差均约为10~(-15)~10~(-14),图形条件对积分结果的影响极小且不存在系统性。综合分析认为,计算结果与理论结果的差异是由计算机的计算误差造成的,单纯形积分结果具有高度的稳定性和精确性。该积分算法为本文三维数值流形方法的研究奠定了基础。 其次,在真实坐标下推导了三维数值流形方法的所有矩阵的元素表达式,具体包括单元刚度矩阵、初应力矩阵、点荷载矩阵、体荷载矩阵、惯性矩阵、速度矩阵、点位移矩阵、接触矩阵、摩擦矩阵。将所有矩阵装配起来即可实现连续、非连续三维数值流形方法模拟。由于数值流形方法不做等参变换,因此该方法的点荷载、点位移等约束可以施加到单元内任意一点。同时,在三维数值流形方法实现过程中,添加了多点位移组合约束算法,可实现相对运动约束,其求解策略包括最小二乘方法和强制约束等。 最后,通过连续变形和非连续变形具体计算实例分析,验证了固定点矩阵、点荷载矩阵、体荷载矩阵、惯性矩阵、速度矩阵、接触矩阵、摩擦矩阵等的正确性。16单元悬臂梁模拟结果与理论结果的对比分析表明程序具有高效性,同时还验证了应力逐步累积模块的有效性;多接触面剪切破坏实例模拟了弹性回跳现象,验证了开闭迭代算法的有效性,证明现有程序可以在极限平衡条件下开展静-动单次转化模拟。 通过对单纯形积分算法、矩阵元素表达式推导、关键算法描述、实例测试等的研究,对三维数值流形方法的算法实现及模拟效果进行了讨论,验证了程序的有效性和模拟结果的高精度特性。 第三、汶川地震的发生与青藏高原的运动与变形密切相关,由于巴颜喀拉地块东向运动受阻,致使龙门山断裂带积累了大量的应变能。因此,开展震前区域变形特征研究有利于更加清楚的认识汶川地震孕育过程。该部分从主应力与主应变率方向差异、不同空间尺度的GPS应变率场分布特征和多步三维数值流形方法模拟等三个方面对汶川地震前地壳变形特征开展了针对性研究。 首先,利用最小二乘配置方法对中国大陆及周边World Stress Map(WSM)计划应力方向数据和GPS速度场数据进行了处理,经过数据预处理、经验协方差参数估计和拟合推估等过程得到了中国大陆范围内应力方向和应变率场结果。数据处理过程中,对应力方向数据和GPS速度场拟合残差的无偏性和正态性进行了检验,并对应力方向和主压应变率网格结果的误差分布特征及其成因进行了分析。在充分考虑数据结果精度的基础上,讨论了中国大陆主压应力和主压应变率方向特征的异同,结果表明二者分布总体一致。在青藏地块东部存在显著差异,应力方向表现为近NE向而GPS主压应变率表现为近EW向,该特征表明青藏地块东部地区近年来东西向挤压应变积累显著;另一显著差异区域位于西域地块西部,应力方向表现为自西向东由NW向NE的转变过程,GPS主压应变率以近NS向为主。 其次,利用前文GPS应变率计算方法对比分析择优推荐的最小二乘配置球面应变率解算方法,基于1999~2007年中国大陆西部约700个GPS测点的速度场结果,分析了汶川8.0级地震前的区域应变率场分布特征。GPS主应变率场结果显示青藏高原主要表现为南北—北北东向挤压、东西—北西西向拉伸的特性。青藏中东部地区呈现由南到北、自西向东主压应变率方向逐步向东偏转的有序分布特征。青藏块体东西向GPS应变率场结果表明,块体西部(92°E以西)以东西向拉张变形为主;青藏块体东部(92°E以东、100°E以西)以东西向挤压为主,该结果反映了青藏高原物质东向流动受到了华北-华南地块的阻挡的动力学特征。因此,应变能在块体边界处积累,有利于汶川地震的发生。虽然汶川8.0级地震震源区不处于剪应变率的高值区,但是处于面应变率挤压区和东西向挤压应变率高值区;发震断裂在震前存在挤压兼右旋剪切变形背景,且震源区西侧的挤压变形幅度明显大于东侧;在断裂带尺度上,应变积累较为缓慢、存在变形趋于极限现象。 最后,通过多步三维数值流形模拟方法研究了汶川地震前挤压应力特征和变形特征。结果表明,最小二乘配置GPS速度场、应变率场结果和三维数值流形方法模拟结果一致性较高,验证了三维数值流形方法模拟结果的可靠性。GPS主应变率场和三维NMM模拟结果均显示在1999~2007年间龙门山断裂带南段的应变积累速率大于中北段,中北段的挤压应变积累已趋于极限,南段依然可以继续进行应变积累。 汶川地震前地壳变形特征表明孕震后期龙门山断裂带及其附近区域的弹性变形已接近极限,与邻近区域相比处于显著弱变形状态。龙门山断裂带南段的应变积累速率大于中北段的结果表明,龙门山断裂带中北段的闭锁程度要大于南段,更有利于地震破裂发生。因此,在研究GPS应变率场分布特征时应该结合构造变形的背景从不同空间尺度进行分析。 总体而言,本文侧重于方法研究,其中GPS应变率场既可用于地壳变形分析又可用于检验模拟结果的正确性;三维数值流形方法的研究侧重公式推导、结果测试等,以保证算法的正确性为首要任务。文章最后给出的三维数值流形模拟实例较为简单,在地学领域进行更广泛的应用尚需开展更加深入的研究工作。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 周硕愚,吴云,秦小军,帅平,施顺英,邓干金,周昕,刘序俨,林继华;基于多种GPS数据研究福建及其邻近海域1994-1997年地壳水平运动[J];地球物理学报;2000年04期
2 许才军,李志才,王华;华北地区活动地块运动时空变化特征[J];大地测量与地球动力学;2002年02期
3 沈正康,王敏,甘卫军,张祖胜;中国大陆现今构造应变率场及其动力学成因研究[J];地学前缘;2003年S1期
4 张效亮;谢富仁;史保平;;利用GPS数据估算川滇南部地震(M_S≥6.5)平均复发间隔[J];地震学报;2010年01期
5 段维波;吴云;陈慧杰;;利用GPS资料求解中国大陆应变场的新方法[J];大地测量与地球动力学;2011年S1期
6 武艳强;江在森;杨国华;刘晓霞;赵静;;汶川地震前GPS资料反映的应变率场演化特征[J];大地测量与地球动力学;2011年05期
7 A.Copley;吕春来;;青藏高原东南边缘的运动学和动力学[J];世界地震译丛;2008年05期
8 朱守彪,石耀霖;广义应变花方法及其在川滇地区地应变率场计算中的应用[J];中国科学院研究生院学报;2005年03期
9 朱守彪;赵晓燕;刘杨;缪淼;;1999年集集地震前后台湾地区地应变率场的分布及其动力学成因[J];大地测量与地球动力学;2011年02期
10 吕国皓;官忠信;;理想可压弹塑性介质中Ⅰ型准静态扩展裂纹端的局部场渐近分析[J];沈阳理工大学学报;1990年03期
11 陈小斌;;中国陆地现今水平形变状况及其驱动机制[J];中国科学(D辑:地球科学);2007年08期
12 江在森;武艳强;方颖;李鹏;王武星;;汶川8.0级地震前区域地壳运动与应变场动态特征[J];地震;2009年01期
13 蒋锋云;张希;张晓亮;王双绪;耿晓民;;用二维有限元方法研究青藏地块地壳水平运动[J];地震研究;2009年03期
14 P.S.Symonds,华云龙 ,余同希;非弹性动力学20年的进展[J];力学进展;1986年03期
15 毛亚明;叶银灿;;横向荷载作用下桩侧极限土抗力的上、下限分析[J];海洋学研究;2010年03期
16 李金平;;基于GPS位移资料的地壳现今形变特征研究[J];测绘通报;2008年01期
17 周硕愚,吴云,施顺英,杨福平;中国大陆东南边缘海现时地壳运动与地震动力学综合研究[J];地壳形变与地震;2001年01期
18 任金卫;中国大陆壳幔运动耦合作用分析[J];地震地质;2004年04期
19 荆燕,任金卫,王辉;用活断层资料分析中国大陆及其邻区晚更新世以来的构造应变特征[J];地震研究;2005年04期
20 江在森;牛安福;王敏;黎凯武;方颖;张希;张晓亮;;活动断裂带构造变形定量分析[J];地震学报;2005年06期
中国重要会议论文全文数据库 前9条
1 朱守彪;石耀霖;蔡永恩;;中国大陆及邻区现今地应变率场的分布及其误差分析[A];中国地球物理学会第二十三届年会论文集[C];2007年
2 朱守彪;石耀霖;蔡永恩;;中国大陆及邻区现今地应变率场的分布及其误差分析[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集(下)[C];2007年
3 孟国杰;申旭辉;Vladimir Smirnov;Eugene A Rogozhin;伍吉仓;;勘察加地区现今地壳运动与变形特征研究[A];中国地震学会空间对地观测专业委员会2009年学术研讨会论文摘要集[C];2009年
4 武艳强;江在森;杨国华;魏文薪;刘晓霞;;GPS应变率场解算方法的对比分析及其可靠性研究[A];中国地震学会第14次学术大会专题[C];2012年
5 江在森;武艳强;;汶川地震前后区域地壳形变场及孕震特征研究[A];中国地震学会第14次学术大会专题[C];2012年
6 杨国华;杨博;韩月萍;张风霜;;华北地区最近几年形变场的动态变化[A];中国地震学会成立三十年学术研讨会论文摘要集[C];2009年
7 于雪梅;程伟;唐立强;吴爽;;压力敏感性材料宏观屈服准则的构建[A];塑性力学新进展——2011年全国塑性力学会议论文集[C];2011年
8 朱守彪;石耀霖;;利用GPS资料计算地应变率场不同方法的比较研究[A];中国地球物理.2003——中国地球物理学会第十九届年会论文集[C];2003年
9 周硕愚;;中国大陆东南边缘海现代地壳运动与区域地震动力学——多种GPS测量、多种地面形变测量和震源机制结果的综合研究[A];2000年中国地球物理学会年刊——中国地球物理学会第十六届年会论文集[C];2000年
中国博士学位论文全文数据库 前3条
1 武艳强;三维数值流形方法研究及其在地学中的初步应用[D];中国地震局地质研究所;2012年
2 秦训鹏;烧结材料三维变形广义上限元理论与应用研究[D];武汉理工大学;2008年
3 侯波;面向制造的铝合金材料失效准则建模与应用研究[D];上海交通大学;2012年
中国硕士学位论文全文数据库 前10条
1 荆燕;利用第四纪活断裂资料模拟中国大陆地壳变形速度场及应变率场[D];中国地震局地质研究所;2002年
2 党学会;用中国大陆构造环境监测网络GPS数据研究中国大陆现今应变场特征[D];中国地震局地震研究所;2012年
3 吕健;基于大地测量数据对青藏高原地区现今地壳形变的分析[D];中国地震局地震研究所;2012年
4 曹鑫;基于非线性破坏准则下的岩石地基承载力上限分析[D];中南大学;2007年
5 杜伟吉;基于GPS和跨断层水准的汾渭断陷带地壳形变特征研究[D];长安大学;2013年
6 张效亮;利用GPS观测数据评估川滇南部地区活动断裂地震危险性[D];中国地震局地壳应力研究所;2009年
7 赵静;块体模型与负位错模型分析川滇地区块体应变与断裂带变形特征[D];中国地震局地震预测研究所;2012年
8 潘博;河套断陷带北缘主要活动断裂带地震危险性研究[D];中国地质大学(北京);2012年
9 李强;汶川地震后青藏高原东北缘GPS地壳形变场特征研究[D];中国地震局地震预测研究所;2013年
10 刘同文;青藏东北缘三维地壳形变特征研究[D];长安大学;2012年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978