用有限元方法计算通海地震断层参数

在用地形变资料研究地震的震源过程以及地壳形变与地震的关系上,不少学者进行了研究。一般利用实测的地面位移场与根据弹性位错理论推导出理论位移场,采用试算法反复进行对比得到。此外,由于在三角平差计算中会得到不同的位移场,文献[2]应     

基于不等式约束的三维电阻率探测混合反演方法

三维电阻率探测的线性反演和非线性反演中均存在着多解性的固有难题.电阻率线性反演方法的效率较高,但反演结果对初始模型的依赖性较强,易陷入局部极小;而非线性反演方法不依赖初始模型,但搜索效率极低,尚未见到关于三维电阻率非线性反演的文献.针对上述问题,融合线性与非线性反演方法的互补优势,提出了最小二乘法（线性方法）与改进遗传算法（非线性方法）相结合的混合反演方法的概念和思想.首先,提出了将介质电阻率变化范围作为不等式约束引入反演方程的思路,以实现压制多解性、提高可靠性的目标.提出了宽松不等式约束和基于钻孔推断的局部严格不等式约束的获取及定义方法.在此基础上,分别提出了基于不等式约束的最小二乘线性反演方法和遗传算法非线性反演方法.其次,对于遗传算法在变异搜索方向控制、初始群体产生等方面进行了改进,优化了其搜索方向和初始群体多样性.然后,提出了混合反演方法及其实现方案,利用改进遗传算法进行第一阶段反演,发挥其对初始模型的依赖程度低的优势,搜索到最优解附近的空间,输出当前最优个体;利用最小二乘法进行第二阶段反演,将遗传算法得到的当前最优个体作为初始模型,在最优解附近空间执行高效率的局部线性搜索,最终实现地电结构的三维成像.最后,开展了合成数据与实际工程算例验证,与传统最小二乘方法进行了对比,发现混合反演方法在压制多解性、摆脱初始模型依赖和提高反演效果方面有较好效果.

     

跨孔电阻率法有限元数值模拟及反演

油气资源和矿产资源勘探中,常需要较大的勘探深度和较高的深部分辨能力.充分利用现有钻孔的井地或跨孔电阻率成像技术适应这种要求,但引起的视电阻率响应复杂,从而给资料的正确处理与解释带来一定的困难,为此有必要通过正演数值模拟和反演相结合的手段,对理论模型井地或跨孔装置的响应异常规律进行研究.文章通过设计一个典型地质模型,利用Ansys有限元软件对模型进行2D跨孔电阻率法有限元正演计算,并将模拟的响应视电阻率作为反演的输入,结合平滑约束最小二乘反演,分析了异常特征和分布规律.     

遗传算法中的光滑约束反演及其在青藏高原面波研究中的应用

光滑约束技术在线性反演中具有重要的作用,但在遗传算法的反演中则很难直接施加于模型参数,其原因是采用光滑处理后的模型参与迭代后,模型的多样性受到很强的压制,并在少量的迭代过程中使种群的各个模型趋向一致,从而得不到满足条件的最优解．本文给出了一种可用于遗传算法反演的间接光滑约束方法．该方法将遗传算法迭代过程中产生的模型经处理后得到的光滑模型,作为误差函数计算的输入模型．迭代过程仍采用原模型,避免了模型的多样性损失,在面波反演和接收函数反演的试验中取得了良好的效果．我们利用该方法对青藏高原地区的瑞利波相速度资料进行了反演,揭示了青藏高原中部地区S波速度结构的横向变化特征．结果表明,青藏高原北部地区地壳Ｓ波速度较南部地区低;大多数路径在15～40ｋｍ 深度范围内,存在12～25ｋｍ 厚的地壳低速层;上地幔低速层位于100ｋｍ 深度以下,厚度主要在40～80ｋｍ 范围内变化,个别路径可达100ｋｍ 以上．安多台以北、玛沁和玉树以西之间,在上地幔90～230ｋｍ 深度范围存在明显的低速层,最低速度约4.2～ 4.3km／s．根据不同路径的S波速度结构和前人的资料,我们认为印度板块的俯冲可能以雅鲁藏布缝合带附近为界．     

地震波频散效应与反Q滤波相位补偿

本文以Futterman提出的地震波振幅衰减和相速度频散表达式为基础,从井震匹配的角度出发,推导了不进行反Q滤波、仅反Q滤波振幅补偿、仅反Q滤波相位补偿,以及反Q滤波相位与振幅同时补偿四种情况下,地震波速度频散与相速度及地震记录振幅谱间的表达式,并从理论上说明了反Q滤波相位补偿的必要性.通过相同观测系统,炸药震源和可控震源分别激发采集的零偏移距VSP资料实例,验证了本文所推导的速度频散表达式的合理性;通过井震标定实例,进一步说明了反Q滤波相位补偿,可有效消除地震子波速度频散,提升地面地震资料与零偏移距VSP走廊叠加剖面的匹配度,最终提高地震资料成果的可信度.     

半全局约束下的多基线立体影像MVLL匹配方法

多视铅垂线轨迹法（multi-view vertical line locus,MVLL）能够以物方地面铅垂线为几何约束,匹配获取物方空间点的最佳高程,是一种实用的多基线立体影像匹配方法。针对MVLL匹配方法中的物方空间点独立匹配、缺少整体性约束问题,将物方空间局部光滑特性应用至匹配过程,提出了一种半全局约束下的多基线立体影像MVLL匹配方法。首先,将物方空间点沿地面铅垂线上的准确高程搜索,等效至像方空间沿核线方向上的准确视差搜索;其次,采用MVLL匹配方法计算物方空间点在多张影像上的等效像方匹配测度;然后,采用半全局匹配法（semi-global matching, SGM）对匹配测度进行多路径聚合分析,得到物方空间局部光滑约束下的等效视差图;最后,将等效视差图转化为物方空间点匹配高程,并对多基线立体影像进行多分辨率匹配处理,实现整体最优条件下的MVLL匹配。采用多类型地形特征与局部区域影像进行匹配实验与对比分析,实验结果表明,所提方法能够对多类型地形特征的物方空间匹配测度进行优化,获取更加可靠的匹配结果,具有更高的影像匹配性能。     

基于CORS数据反演江苏上空电离层电子密度

本文利用IRI2007模型和江苏CORS网数据,结合附加约束的同时迭代重构算法,反演了地磁平静日江苏上空电子密度的结构分布。结果表明,整个研究区域电子密度在不同纬度和不同高度上存在明显的日变特征,电子密度最大值出现在北京时13～15时,随着时间推移,峰值幅度衰减;白天出现明显的E层,夜晚消失;反演结果说明了电离层三维层析技术为监测电离层时空结构提供了一种强有力的实验支持。     

CT技术用于激光等离子体成象研究

本文研究适用于ＩＣＦ实验的ＣＴ技术，并已用于“星光Ⅱ”装置的激光等离子体的实验Ｘ光测量，重建出Ｘ光的三维分布，获得了有意义的结果，表明该技术在ＩＣＦ实验中的应用的可能性。     

高度约束MHE算法在伪距单点定位中的应用

针对目前全球卫星导航系统(GNSS)中伪距单点定位(SPP)技术的定位精度已不能满足现代大多数应用场景的需求,提出了一种附加高度约束的滚动时域估计(MHE)算法,以此来改善SPP的定位性能. 附加高度约束的MHE算法是将接收机的位置高度作为非线性约束加入到SPP的估计参数中,并采用近似MHE算法来进一步提高定位精度的优化算法. 结果表明:高度约束的MHE滤波比传统最小二乘(LS)的滤波具有更好的平滑特性,同时随其视窗大小的增加,其定位精度得到了进一步改善. 验证了附加高度约束MHE方案的有效性、可行性,所得结果对SPP的实际应用具有重要的参考意义.      

室内楼层平面导航路网模型的自动生成

室内路网模型作为室内导航研究的基础,如何自动生成室内路网模型成为近年的研究热点.对于目前自动路网模型研究中出现的复杂环境适应度不够的问题,本文提出了一种室内楼层平面路网模型的自动提取方法.它将室内空间分为公共空间和专属空间两类,从而形成公共空间路径和专属空间路径,并用转换点将两类路径进行衔接;最终利用提取的公共空间路径...