用传输线模型法计算非均匀介质中的电场响应

介绍了一种非均介质中电场响应的数值模拟方法.传输线模型法是以惠更斯原理为基础,通过建立麦克斯韦方程中电场和磁场与传输线方程中电压和电流的等效关系,按时间步进的办法来模拟电磁场响应的一种电磁场数值计算方法.从一维介质中电场传播的物理过程出发,导出了点电流源情形下地质介质中电场响应迭代方程;编程计算了非均匀介质中的电场响应,其结果表明,用这种方法对非均匀介质中的电场响应进行数值模拟是行之有效的.     

1997年玛尼地震对青藏川滇地区构造块体系统稳定性影响的三维DDA+FEM方法数值模拟

本文用三维非连续变形与有限元相结合（DDA+FEM）的方法,在青藏川滇地区三维构造块体相互制约的大背景中,通过用GPS资料做位移速率边界约束和震源机制约束,计算得到研究区的初始位移场和应力场与该地区GPS测量结果和震源机制分布结果基本一致.在此基础上进一步数值模拟1997年玛尼7.9级大震的发生过程,研究大震引起研究区各块体边界断层应力状态变化的特征.（1）发震断层两侧发生左旋走滑错动,最大水平位错大约7 m;（2）深部位错面上位错分布与用地震波资料震源反演的结果类似;（3）最大差应力变化等值线图与由星载D\|INSAR技术获取的地表形变场图像相似;（4）地表垂直位移表明地震断层面略向北逆冲.计算模拟得到了玛尼地震发生引起青藏川滇地区构造块体系统各边界断层上库仑破裂应力变化的分布,表明玛尼大震的发生除了使其发震断层的两端库仑破裂应力增大,应力进一步集中外,位于上地壳层上东昆仑断裂中段的2001年昆仑山8.1大震（H=11 km）发震断层段的库仑破裂应力增加约2 MPa,位于中地壳层上喀拉昆仑断裂带中的2008年改则6.9级地震（H=30 km）发震断层段的库仑破裂应力也增加约0.7 MPa,可见这两个已接近破裂强度地段的失稳对发生大震起了一定促进作用.研究结果也表明：作者发展的三维DDA+FEM方法能有效地用于大震活动与各构造块体相互作用关系的研究.     

地震波数值模拟的非规则网格PML吸收边界

以格子法为基础,以声波方程为例研究非规则网格PML(Perfectly Matched Layer)方法.本方法的核心是建立局部坐标系下的分裂方程和基于积分近似的微分方程弱形式.该非规则网格模拟方法允许在计算域内设置任意形状的人工边界.对于二维半空间问题,与采用矩形人工边界相比,采用半圆形人工边界可减少计算量20%以上.采用光滑的曲边界,不仅可减少计算区域,还可避免常规的PML吸收边界在吸收带角点区域的特殊处理.本方法事先计算和存储边界单元的局部几何参数,在计算的每一时间步查表调用这些参数,与常规的直边界PML方法相比,不增加任何计算量.     

高层钢结构抗震pushover分析的侧向力分布模式及其影响

本文针对高层钢结构抗震Pushover分析的适用性和侧向力分布模式的选取等问题,首先,介绍和分析了Pushover分析常用的五种侧向力分布模式及其本质特征;然后,给出了判断侧向力分布模式合理性的准则;最后,通过对比一座15层的平面钢框架结构在Ⅳ类场地的典型地震动下的非线性时程分析和五种侧向力分布下的Pushover分析中的破坏形式和弹塑性反应,研究了高层钢结构抗震Pushover分析的适用性和各侧向力分布模式的合理性;在此基础上,本文给出了高层钢结构Pushover分析时侧向力分布模式的选取建议。     

汶川M_S 8.0地震InSAR同震形变场特征分析

采用7条地震前后日本ALOS/PALSAR整轨数据,利用D-InSAR技术提取了2008年5月12日四川汶川地震500km×450km的地表连续同震形变场。形变场覆盖了金川—石棉、黑水—乐山、松盘—彭山、南坪—简阳、康县—重庆所有区域,包括了理县、汶川、茂县、北川、青川等重灾区域。结果显示,整个形变场影响范围较大,四川盆地出现了不同程度的地表形变。发震断层附近的非相干区域显示此次地震地表主破裂带在北川-映秀断裂带上,可以追踪出地表破裂带从汶川县映秀镇西南震中附近一直到青川县苏河北侧,全长约为230km。在发震断层西南段汶川至茂县一带,非相干条带的宽度明显大于其它段落,这与彭县-灌县断裂(前山断裂)的都江堰—安县段地表破裂段密切相关,地表破裂带长约70km。在远离发震断层的区域,西北盘总体表现为抬升,东南盘表现为沉降。但在发震断层附近,断层两侧均表现为局部抬升,且沿断层形变量分布很不均匀,表现出较强的分段性,显示出发震断层以逆冲为主的断层性质。从美国哈佛大学(Harvard)、美国地质调查局(USGS)与美国国家地震信息中心(NEIC)、中国地震台网中心(CENC)所给出的震中位置和发震时刻的差异来看,也反映出     

中越红河断裂带考察与研究（一）

基于对越南境内红河断裂带的野外考察,结合中国境内红河断裂带已有考察研究成果,从断裂几何结构、构造演化与地层发育、新生代地质地貌特征等多个方面,初步探讨了红河断裂带现今活动的基本特征.认为:(1)红河断裂带在断裂几何结构、断裂长期的发育演化过程、断裂活动的地质地貌表现等诸多方面都表现出异常突出的时空不均匀性和分段性;(2)红河断裂新构造运动以来的夷平面发育与解体、新生代盆地发育与堆积、水系及河流阶地发育与变形等方面的特征表明,中国境内红河断裂带第四纪以来总体活动水平由北向南逐渐减弱的趋势一直延续到了越南境内;(3)青藏高原侧向挤出作用的影响范围足有限的,这一认识对于探索川滇地区地壳运动与强震活动的动力学机制具有重要意义.     

复杂介质地震定位中震源轨迹的计算

交切法是最基本的地震定位方法之一.该方法利用震源轨迹进行定位,不需要求解方程,即使仅有少量的地震记录,利用交切法也能获得有价值的震中信息.但另一方面,交切法的定位精度较低,特别是震源深度的误差常常较大.     

强震前地震非均匀度异常与前兆群体异常关系的研究

通过对比研究华北地区部分中强以上地震震前地震非均匀度异常与区域地震前兆观测群体异常的变化过程, 发现两者在地震前的异常时空演化过程中, 具有一定的关联性。 在时间上, 部分地震震前的地震非均匀度异常出现后, 半年左右前兆观测异常在数量上将出现一次较大幅度的增加; 在空间上, 前兆异常展布范围与地震非均匀度异常范围大体一致。 研究结果表明, 地震非均匀度(GL值)在强地震前出现的异常变化从震例的角度, 表明可能与局部应力场变化有关。 地震非均匀度和前兆群体异常的这种准同步变化, 可能对判定地震孕育进入中短期阶段具有一定的指示意义。     

含流体致密砂岩的纵波频散及衰减:基于双重双重孔隙结构模型描述的特征分析

致密砂岩普遍具有低孔、低渗及微裂缝发育的地质特征,并且呈现出很强的非均匀性.致密砂岩储层与常规砂岩储层比较,具有明显的岩石物理性质、渗流力学性质方面的差异.致密砂岩内部的非均匀性对弹性波频散、衰减有显著影响,其中包括孔隙结构的非均匀性,即岩石内部孔隙参数的不均一性,以及孔隙内部不相混溶流体的非均匀分布;此外,非均匀性的尺度也决定了波出现显著频散与衰减的频段.综合考虑致密砂岩孔隙结构非均匀性及流体斑块状饱和的非均匀性,本文采用双双重孔隙介质结构模拟了致密砂岩的弹性波响应,分析了同时具备两类非均质性岩石中的波传播特征.调查分析了两组分别来自中国鄂尔多斯盆地苏里格气田及四川盆地广安气田的不同类型致密砂岩储层的岩芯超声波实验数据,给出了岩石样本的弹性波速度频散与衰减曲线.结果显示理论模型预测结果与完全饱和、部分饱和岩石的实验数据吻合良好.对两个地区致密砂岩岩芯数据进行对比分析,苏里格致密砂岩样本总体上比广安致密砂岩渗透率高,在各孔隙度范围内,特征模拟显示苏里格样本的裂隙尺寸明显大于广安样本.广安致密砂岩在低孔隙度范围内发育了更多、更小的颗粒裂隙/接触.致密砂岩的速度频散与衰减结果受流体黏度、晶体破裂及流体斑块状饱和的共同影响.此外,孔隙度越大,部分饱和岩石中斑块状饱和机制对总衰减的贡献越低,与之相对,结构非均质性所占的比重则有所增强.     

面向目标自适应有限元法的带地形三维大地电磁各向异性正演模拟

传统三维大地电磁各向异性模拟均是基于规则六面体网格,计算精度有限且较难拟合复杂地质条件.本文采用面向目标自适应非结构矢量有限元法,对三维大地电磁各向异性介质进行模拟.首先从电场双旋度方程出发,利用伽辽金方法建立变分方程;然后利用电流密度连续性条件构建适合大地电磁各向异性问题的加权后验误差估计方法,实现面向目标的网格自适应正演;最后通过典型算例分析各向异性对网格自适应和大地电磁响应的影响特征以及各向异性的识别方法.本文算法能够高精度地拟合起伏地表和任意各向异性介质,适用于分析复杂地电条件大地电磁响应特征,为提高大地电磁资料解释水平提供了理论基础.