基于垂向一阶导数与解析信号比值的欧拉反演方法

常规欧拉反褶积法中构造指数的选取以及分散解存在较多的问题,提出了基于联立垂向一阶导数与解析信号的欧拉齐次方程的RDAS-Euler反演方法。该方法可以更为精确的估计场源的范围及埋深,且不需考虑构造指数N的影响,避免了因构造指数不当而引起的反演误差。通过对单一地质体及组合地质体模型的实验证明本文方法能有效地完成目标体的反演工作,反演结果与理论值之间的误差小于10%,且相对于常规欧拉反褶积法更加稳定准确,能够更好的得到地质体边界及深度信息。将RDAS-Euler法应用于黑龙江省虎林盆地实测布格重力异常数据,获得了丰富的断裂信息,说明RDASEuler法增强了对断裂平面位置的识别能力。     

风沙运动的欧拉双流体模型参数研究

利用数值方法研究了双流体模型中动力黏度μ_s,kin、动量交换系数β和碰撞恢复系数e_ss对风沙运动的影响.结果显示:动力黏度μ_s,kin影响沙床面起沙率,动量交换系数β影响输沙率沿流程的变化,碰撞恢复系数e_ss影响输沙量的垂向分布.利用研究得到的参数,模拟风沙流运动,并与已有的诸多实验和数值模拟结果对比,结果吻合较好,说明该参数的选取可以有效地模拟风沙流运动.     

Multi‐transient electromagnetic repeatability experiment over the North Sea Harding field‡

We present results of synthetic time‐lapse and real repeatability multi‐transient electromagnetic surveys over the North Sea Harding field. Using Archie's law to convert porosity and fluid saturation to resistivity we created 3D isotropic models of the reservoir resistivity at different stages of production from the initial state in 1996 through to complete hydrocarbon production by 2016 and, for each stage, we simulated an east‐west transient electromagnetic survey line across Harding. Unconstrained 1D full‐waveform Occam inversions of these synthetic data show that Harding should be detectable and its lateral extent reasonably well‐defined. Resistivity changes caused by hydrocarbon production from initial pre‐production state to production of the oil rim in 2011 are discernible as are significant changes from 2011–2016 during the modelled gas blowdown phase. The 2D repeatability surveys of 2007 and 2008 tied two wells: one on and the other off the structure. Between the two surveys the segment of the field under investigation produced 3.9 million barrels of oil – not enough to generate an observable time‐lapse electromagnetic anomaly with a signal‐to‐noise ratio of 40 dB. Processing of the 2007 and 2008 data included deconvolution for the measured source current and removal of spatially‐correlated noise, which increased the signal‐to‐noise ratio of the recovered impulse responses by about 20 dB and resulted in a normalized root‐mean‐square difference of 3.9% between the data sets. 1D full‐waveform Occam inversions of the real data showed that Harding was detectable and its lateral extent was also reasonably well‐defined. The results indicate that the multi‐transient electromagnetic method is suitable for exploration, appraisal and monitoring hydrocarbon production.     

由精密控制人工震源观测到的绵竹5.6级地震前后波速变化

地震危险区断层的应力状态变化及其精确测量一直是地球物理学关注的热点问题之一,主动震源探测技术的发展为解决上述问题带来了可能.汶川大地震后,结合汶川地震断裂带科学钻探,在四川省绵竹市九龙镇--跨越龙门山断裂带前山断裂,利用10 t精密控制人工震源进行了连续监测试验,在试验观测期间的2009年6月30日凌晨2时3分,当地发生了M_s5.6级地震,地震前后观测数据的对比分析表明,该地震发生前后,穿过断裂带的直达S波走时发生了时延为5~9 ms的微弱变化,射线路径上的平均波速相对下降了约0.3%,该时延主要是由地震同震效应引起的断裂带内应力调整和物性变化造成的.所使用的精密控制人工震源性能稳定,适于开展断层监测.这一结果对于发展断裂带主动震源动态监测技术具有重要意义.     

中国东北地区重力场研究——利用改进的三方向小子域滤波划分主构造线及大地构造单元

利用重力异常进行大地构造单元的划分是区域重力资料处理与解释的重要任务.为了精细地划分中国东北地区大地构造单元及其分界线位置,对三方向小子域滤波进行了改进,探讨了改进的三方向小子域滤波的合理性.利用小子域及改进的三方向小子域滤波对研究区布格重力异常进行处理,并与欧拉反褶积反演结果进行对比分析.结果表明,改进的三方向小子域滤波在识别平面断裂的应用中具有更好的效果.综合分析东北地区布格重力异常特征、小子域和改进的三方向小子域滤波结果,划分出延续性较好、规模较大的断裂18条.对研究区重力异常进行了分区,描述了一、二级异常的分区特点,给出了东北地区大地构造单元的基本格架.黑河-齐齐哈尔-白城等重要构造线平面位置的确定突破了传统的重力异常分区模式,指出东北地区NE向构造具有走滑特征,且走滑距离具有由西向东逐渐增大的特点.     

实测板块运动模型修正导航系统监测站坐标精度分析

为了解决卫星导航系统各类地面站的坐标修正问题,必需精确给出各个地面站受地壳板块运动引起的位置变化速度。采用中国地壳运动观测网GPS观测数据计算的实测速度场,利用滑动欧拉矢量建模法建立局部虚拟的板块运动模型,并对虚拟板块运动模型计算的测站速度场的精度进行了分析,结果表明,在东西方向的精度为2.092 mm/a,在南北方向的精度为2.072 mm/a,在点位分布较少的地区和我国边界地区,其精度并没有明显降低,据此指出,该模型可用于我国卫星导航系统各类地面站的坐标修正。     

移动荷载作用下双层Euler梁模型土动力响应分析

本文以弹性波动理论为基础,利用分层法计算了双层连续支撑Euler梁下地基土的位移,同时考虑了不同的车辆荷载模型,并将计算的结果和单层Euler梁下地基土的位移进行了比较,比较表明：轨道结构本身和车辆悬挂体系具有一定的减振功能。因此,该研究对由交通荷载引起的环境振动分析和评价具有一定的应用价值。     

Earthquake damage detection in the Imperial County Services Building III: Analysis of wave travel times via impulse response functions

The majority of structural health monitoring methods are based on detecting changes in the modal properties, which are global characteristics of the structure, and are not sensitive to local damage. Wave travel times between selected sections of a structure, on the other hand, are local characteristics, and are potentially more sensitive to local damage. In this paper, a structural health monitoring method based on changes in wave travel times is explored using strong motion data from the Imperial Valley Earthquake of 1979 recorded in the former Imperial County Services (ICS) Building, severely damaged by this earthquake. Wave travel times are measured from impulse response functions computed from the recorded horizontal seismic response in three time windows—before, during, and after the largest amplitude response, as determined from previous studies of this building, based on analysis of novelties in the recorded response. The results suggest initial spatial distribution of stiffness consistent with the design characteristics, and reduction of stiffness following the major damage consistent with the spatial distribution of the observed damage. The travel times were also used to estimate the fundamental fixed-base frequency of the structure f₁ (assuming the building deformed as a shear beam), and its changes during this earthquake. These estimates are consistent with previous estimates of the soil–structure system frequency, f_sys, during the earthquakes (f₁<f_sys as expected from soil–structure interaction studies), and with other estimates of frequency (f₁ from ETABS models, and f_sys from ambient vibration tests, and “instantaneous” f₁ from high-frequency pulse propagation).     

系泊船非线性波浪力时域计算:二维模型

为找到具有工程实用价值的港口系泊船波浪力的时域计算方法,建立了在港口中存在系泊船时非线性波浪力时域计算的垂直二维耦合模型:用Boussinesq方程计算船的两侧的外域,用欧拉方程计算船底面下的内域,两域在交界面处的连接条件是流量连续和压力相等.将复平面内的边界元方法应用于所研究问题,对耦合模型进行了验证.进行了相关模型实验,实验结果与数值计算结果比较表明这两种数值计算模型都具有满意的精度,但耦合模型的计算效率要远远高于边界元方法的计算效率.本耦合模型的数学处理简单,可适用于工程计算.