Generalized Born inversion of seismic reflection data

Summary . Born inverse methods give accurate and stable results when the source wavelet is impulsive. However, in many practical applications (reflection seismology) an impulsive source cannot be realized and the inversion needs to be generalized to include an arbitrary source function. In this paper, we present a Born solution to the seismic inverse problem which can accommodate an arbitrary source function and give accurate and stable results. It is shown that the form of the generalized inversion algorithm reduces to a Wiener shaping ***filter, which is solved efficiently using a Levinson recursion algorithm. Numerical examples of synthetic and real field data illustrate the validity of our method.     

On the spatial distribution of seismicity and the 3D tectonic stress field in western Greece

We analyzed a large number of focal mechanisms and relocated earthquake hypocenters to investigate the geodynamics of western Greece, the most seismically active part of the Aegean plate-boundary zone. This region was seismically activated multiple times during the last decade, providing a large amount of enhanced quality new information that was obtained by the Hellenic Unified Seismological Network (HUSN). Relocated seismicity using a double-difference method appears to be concentrated above ∼35 km depth, exhibiting spatial continuity along the convergence boundary and being clustered elsewhere. Earthquakes are confined within the accreted sediments escarpment of the down-going African plate against the un-deformed Eurasian hinterland. The data arrangement shows that Pindos constitutes a seismic boundary along which large stress heterogeneities occur. In Cephalonia no seismicity is found to be related with the offshore Cephalonia Transform Fault (CTF). Onshore, NS crustal extension dominates, while in central and south Peloponnesus the stress field appears rotated by 90°. Shearing-stress obliquity by 30° is indicated along the major strike-slip faults, consistent with clockwise crustal rotation. Within the lower crust, the stress field appears affected by plate kinematics and distributed deformation of the lower crust and upper mantle, which guide the regional geodynamics.     

Rockburst research at Falconbridge's Strathcona Mine,Sudbury, Canada

Microseismic systems at five mines in the Sudbury Basin provide the basic data for Falconbridge Limited's rockburst research. Daily and long-term analysis of this data as well as underground observations have confirmed the fault-slip mechanism at three mines. A detailed analysis of the complete history of Falconbridge Mine is being conducted and Distinct Element numerical models are being used to simulate both the stick-slip behaviour of faults and the dynamic effects of the induced vibrations on rock and backfill. University research includes acoustic tomographic imaging of the rock mass based on seismic wave propagation and collection of full microseismic waveforms to allow application of advanced seismic and statistical analysis techniques.Formerly with Mines Technical Services, Falconbridge Limited, Sudbury Operations, Falconbridge , Ont., P0M 1S0.Presented at the Fred Leighton Memorial Workshop on Mining Induced Seismicity, Montreal, Canada, August 30, 1987.     

青藏高原及邻近地区地壳与上地幔剪切波三维速度结构

本文利用30个基准台所记录的238条长周期面波资料,经过适配滤波和分格频散反演,得到中国大陆及邻区147个分格10—105s的纯路径频散,进而反演出青藏高原及邻近地区深至170km的剪切波三维速度结构.研究表明,青藏高原中西部地区和东部地区的地壳平均厚度分别为70±7km和65±7km,地壳平均剪切波速度分别为3.55和3.62km/s,上地幔顶盖平均速度分别为4.63和4.61km/s; 岩石层厚度均为120±10km;东部地区下地壳内30—40km深度处普遍存在低速层;青藏高原及其东侧的上地幔低速层内有横贯东西且明显向上隆起的低速腔.滇西缅北地区的地壳厚45±5km,上地壳及下地壳内都有低速层;上地幔顶盖的速度为4.42km/s,比青藏高原本体及恒河平原都低.恒河平原地壳厚34±2km,速度平均为3.45km/s;上地幔顶盖厚86±10km,速度平均为4.63km/s,顶盖内55—83km深处有一个低速夹层.     

三河-平谷8.0级地震区浅层三维密度结构反演研究

为研究三河-平谷地震区浅层构造背景、地震孕育机理及地震与构造的关系,利用高精度重力异常数据,采用基于块体生长模式的重力三维反演算法对地震区浅层三维密度结构进行反演,并通过模拟试算验证基于块体生长模式反演方法的有效性和稳定性。高精度布格重力异常显示,三河-平谷8.0级地震位于大兴重力局部高、三河-马坊重力局部高与大厂重力低之间的交汇过渡低值区域。研究区浅层三维密度结构反演结果表明,1679年三河-平谷M8.0地震明显受NE向夏垫断裂控制,断裂两侧密度差异明显且向下延伸约10 km,推测发震部位深约10 km。     

青藏高原西部及邻区构造应力场反演和断层滑动趋势分析

基于Global CMT提供的1976-01-01~2015-12-31期间1 255个M_W≥3.5地震的震源机制解,利用FMSI方法对青藏高原西部及邻区（新疆西部地区）应力场进行反演,得到区域应力场特征,并分析西昆仑断裂带及阿尔金断裂带的断层滑动趋势。结果表明,青藏高原西部及邻区的最大主压应力轴σ1整体呈NNE-SSW向,自西向东顺时针旋转;最小主压应力轴σ3整体呈NWW-SEE向。最大主压应力轴σ1倾角大小约为5°~9°,说明该地区的构造应力场以水平运动为主。西昆仑断裂带及阿尔金断裂带具有较强的滑动趋势（T_s≥0.5）,说明西昆仑断裂带及阿尔金断裂带未来失稳滑动的可能性较大。     

岩石脆性预测在吉木萨尔凹陷致密油勘探中的应用

岩石的脆性是非常规油气藏勘探开发中需考虑的关键岩石力学参数之一。岩石脆,压裂过程中易产生缝网,从而有效改造储层,达到提高单井产量的目的。基于岩石弹性参数的归一化杨氏模量和泊松比方法,是当前广泛采用的可用于地震预测的岩石脆性评价技术,但其结果受归一化参数取值的影响较大,适用性受到影响。利用杨氏模量和泊松比的商作为脆性评价指标,并用大量岩石力学实验刻度,大大提高了岩石脆性判别的可信度。利用地震叠前弹性参数反演,开展吉木萨尔凹陷二叠系芦草沟组致密油岩石脆性预测,指导水平井部署。为保证更好的储层改造效果,在水平井轨迹设计过程中不但考虑了储层物性,还充分考虑了岩石脆性。在优选的有利区钻探水平井,并实施了体积压裂,微地震监测结果显示压裂产生了缝状网,储层改造效果理想。该水平井试油获得了稳定高产工业油流。     

计算机互联网上时间传递的一种设计方案

目前,国内外时间传递方式主要有长波授时、短波地、电话授时、双星比对、ＧＰＳ时间信息等。信息产业的不断发展,使得互联网络已成为一种新的信息传递手段和资源共享方式 。通过计算机网络传输标准时间、频率将是授时领域的新课题,也是时间同步的新手段。讨论了计算机互联网上时间传递的一种设计方案,给出利用计算机网络系统实现标准时间的传递的具体方法,并介绍一种计算机互联网时间传递软件的使用方法等。     

基于非结构化有限元的三维井地电阻率法约束反演

电磁探测反演是典型的不适定问题,易造成反演结果的多解性,不适定性是反演自身固有的特征,没有求解的附加信息这一本质困难是很难克服的,解决该问题的有效方法是研究约束反演。本文采用目前较为主流的高斯牛顿—共轭梯度法(GN-CG),在反演目标函数中直接施加约束条件,将介质电阻率的取值范围作为先验信息和约束条件以外点罚函数法的方式引入到反演目标函数中,与常规三维电阻率反演目标函数相比,增加了不等式约束项的目标函数,理论上可以压制反演的多解性。通过多种理论模型的测试结果表明,本文基于不等式约束的三维井地电阻率反演算法有效地改善了反演结果的精度,以惩罚函数法施加不等式约束条件的方式是现实可行及有效的。