珠江口盆地西江主洼油气差异分布机制

西江主洼是珠江口盆地一个低勘探程度洼陷，油气差异聚集特征明显，但其机理不清。本文在断陷盆地油气成藏理论指导下，利用研究区地质、地球物理和地球化学等资料，开展了烃源岩、断盖组合、储集体系和运聚模式等方面的研究。结果表明，裂陷期沉降、沉积中心有序迁移，造成主力烃源岩自东向西由文四段迁移至文三段、文一+二段，控制东、西部油气差异分布。恩平组区域性泥岩发育且晚期断裂缺乏，导致油气纵向上更易聚集于下构造层；区域泥岩减薄尖灭或晚期断裂切开盖层的区域，上构造层有一定油气分布。"源-汇"类型从宏观上控制储层优劣，影响下构造层油气富集程度；"仓储"运移是控制上构造层油气规模聚集的主要模式。下构造层古近系应作为重点部署方向，东部围绕文四段、西部围绕文三段、文一+二段烃源岩构成的含油气系统展开；上构造层勘探需关注珠海组，在隆起周边寻找具备"仓储"运移模式的有利区带。     

地铁施工对济南白泉泉群流量的影响

为研究地铁建设对济南白泉泉群的影响，在综合分析白泉泉域地质、水文地质条件的基础上，假定研究区岩溶强径流带位置及水力性质，利用FEFLOW软件建立地下水流数值模型。以规划地铁M1号线为研究对象，分析了济南东站、梁王站、梁王东站分别施工及3个站同时施工4种情景下，采用施工降水或施工降水+人工回灌两种施工方式对白泉泉群流量的影响。结果表明：单独采用施工降水的施工方式使得白泉泉群流量衰减，其中3个站同时施工对泉流量的影响最大，泉流量最大衰减达5.48%；各站分别施工时，济南东站对泉流量影响最大，泉流量较未施工时减少了0.043×10⁴ m³/d。采用施工降水+人工回灌的施工方式，能够有效缓解泉流量的衰减，各车站施工时的泉流量衰减由仅施工降水时的2.26%~5.48%降低至0.08%~1.21%。岩溶强径流带有利于地下水形成优势径流，促进白泉泉群补给，一定程度上缓解因地铁施工引起的泉流量衰减。     

雄安新区地热地质模型探究：来自地球物理的证据

雄安新区内地热资源丰富，区内有牛驼镇地热田、容城地热田和高阳地热田，地热资源开发利用较早，但是对其深部热源机制仍未形成统一观点。为了研究雄安新区内地热田深部热源机制，在新区及外围进行了深反射地震和长周期大地电磁探测，对取得的同剖面的深反射地震和大地电磁数据进行处理和综合解释，探明了研究区从地表至莫霍面范围内地质构造和电性结构。下地壳结构在深反射地震剖面与大地电磁剖面上有很好的对应关系。电阻率低值区对应着在深反射地震剖面上存在一系列反射同相轴，且同相轴可以延续到莫霍面，电阻率高值区对应着在深反射地震剖面上无明显连续反射同相轴，尤其是在莫霍面之上呈现地震反射近似"空白区"。结合区域地热资料构建了研究区深部地热地质模型，对新区内深部地热机制进行了解释。该模型为"二元"生热模型，其热源包含两个部分，深部地幔热源和地壳放射性元素衰变生热。放射性元素衰变生热占地表热流的接近30%，而幔源热流在地表热流中的占比可达约70%。在牛驼镇下方，莫霍面以上，由于地幔热物质上涌造成下地壳上隆，幔源岩浆底侵作用于下地壳形成了局部热异常，该热异常具有低速高导的地球物理特征，认为是牛驼镇地热田和容城地热田的深部热源；以区域断裂为热通道，大地热流由深部向上传导、扩散到牛驼镇凸起和容城凸起顶部，对碳酸盐岩储水层进行加热，形成地热储层；上覆新近系沉积地层是良好的热盖层。     

船式支护在基坑中的应用

在深厚软土地区的基坑,对基坑底软土进行搅拌加固形成强度较高的底板而保证坑底稳定,侧壁采用支护桩支护保证侧壁的稳定使其不倾倒,支护桩与加固底板两者共同作用形成的船状结构支护体系类似轮船的船侧壁及船底板,以保证基坑的稳定。结合具体的工程实例,介绍了船式支护的侧壁及底板的设计方法。工程案例成功实施,监测结果表明,船式支护既可避免支护桩嵌固深度过长,又利于软土的土方施工及工程桩的保护,是安全可靠的,对类似工程有一定的参考意义。     

A Block-Based Theoretical Model Suited to Gravitational Sliding

An original theoretical model has been devised to simulate mass flow over hill slopes due to gravitational sliding. The sliding mass is discretized into a sequence of contiguous blocks which are subjected to gravitational forces, to bottom friction and to surface resistance stresses that are generally negligible for subaerial flows, but are relevant for submarine slides. The blocks interact with each other while sliding down the hill flanks because of internal forces that dissipate mechanical energy and produce a momentum exchange between the individual blocks, yet conserving the total momentum of the mass. Internal forces are expressed in terms of interaction coefficients depending on the instantaneous distance between the block centers of mass, which is a measure of the deformation experienced by the blocks: the functional dependence includes three parameters, namely the interaction intensity ¯, the deformability parameter and the shape parameter , by means of which a wide range of interaction types can be fully accounted for. The time integration is performed numerically by solving the equations for the block velocities and positions at any time ti by means of the block accelerations at the previous time t_i-1, and by subsequently updating the block accelerations, which allows to proceed iteratively to the following times. The model has been tested against laboratory results available from literature and by means of several numerical experiments involving a simplified geometry both for the sliding body and the basal surface, with the purpose of clarifying the influence of the model parameters on the slide dynamics. The model improves the performance of the existing kinematic models for slides, moreover preserving an equivalent numerical simplicity. Future applications and possible improvements of this model are suggested.     

Research on temporal and spatial distribution, evolutionary character and mechanism of crustal deformation field before and after the Tangshan earthquake

Ｒｅｓｅａｒｃｈｏｎｔｅｍｐｏｒａｌａｎｄｓｐａｔｉａｌｄｉｓｔｒｉｂｕ┐ｔｉｏｎ,ｅｖｏｌｕｔｉｏｎａｒｙｃｈａｒａｃｔｅｒａｎｄｍｅｃｈａｎｉｓｍｏｆｃｒｕｓｔａｌｄｅｆｏｒｍａｔｉｏｎｆｉｅｌｄｂｅｆｏｒｅａｎｄａｆｔｅｒｔｈｅＴａｎｇｓｈａｎ...     

The Ms = 6.2, June 15, 1995 Aigion earthquake (Greece): evidence for low angle normal faulting in the Corinth rift

We present the results of a multidisciplinary study of the M_s = 6.2, 1995, June 15, Aigion earthquake (Gulf of Corinth, Greece). In order to constrain the rupture geometry, we used all available data from seismology (local, regional and teleseismic records of the mainshock and of aftershocks), geodesy (GPS and SAR interferometry), and tectonics. Part of these data were obtained during a postseismic field study consisting of the surveying of 24 GPS points, the temporary installation of 20 digital seismometers, and a detailed field investigation for surface fault break. The Aigion fault was the only fault onland which showed detectable breaks (< 4 cm). We relocated the mainshock hypocenter at 10 km in depth, 38 ° 21.7 N, 22 ° 12.0 E, about 15 km NNE to the damaged city of Aigion. The modeling of teleseismic P and SH waves provides a seismic moment M_o = 3.4 10¹⁸ N.m, a well constrained focal mechanism (strike 277 °, dip 33 °, rake – 77°), at a centroidal depth of 7.2 km, consistent with the NEIC and the revised Harvard determinations. It thus involved almost pure normal faulting in agreement with the tectonics of the Gulf. The horizontal GPS displacements corrected for the opening of the gulf (1.5 cm/year) show a well-resolved 7 cm northward motion above the hypocenter, which eliminates the possibility of a steep, south-dipping fault plane. Fitting the S-wave polarization at SERG, 10 km from the epicenter, with a 33° northward dipping plane implies a hypocentral depth greater than 10 km. The north dipping fault plane provides a poor fit to the GPS data at the southern points when a homogeneous elastic half-space is considered: the best fit geodetic model is obtained for a fault shallower by 2 km, assuming the same dip. We show with a two-dimensional model that this depth difference is probably due to the distorting effect of the shallow, low-rigidity sediments of the gulf and of its edges. The best-fit fault model, with dimensions 9 km E–W and 15 km along dip, and a 0.87 m uniform slip, fits InSAR data covering the time of the earthquake. The fault is located about 10 km east-northeast to the Aigion fault, whose surface breaks thus appears as secondary features. The rupture lasted 4 to 5 s, propagating southward and upward on a fault probably outcropping offshore, near the southern edge of the gulf. In the shallowest 4 km, the slip – if any – has not exceeded about 30 cm. This geometry implies a large directivity effect in Aigion, in agreement with the accelerogram aig which shows a short duration (2 s) and a large amplitude (0.5 g) of the direct S acceleration. This unusual low-angle normal faulting may have been favoured by a low-friction, high pore pressure fault zone, or by a rotation of the stress directions due to the possible dip towards the south of the brittle-ductile transition zone. This fault cannot be responsible for the long term topography of the rift, which is controlled by larger normal faults with larger dip angles, implying either a seldom, or a more recently started activity of such low angle faults in the central part of the rift.     

参数的期望估计及其在形变分析中的应用

利用“参数的期望估计”能准确定位、定量粗差及参数的期望估计不受粗差影响折特殊性质,寻找地壳变形区域和不变形区域,以确定拟稳点,进而进行拟稳变换,分析地菜变,将是一种有效方法。     

四川木里耳泽岩溶型金矿床形成条件和成矿机制

四川本里耳泽金-菱铁矿建造矿床,赋存于上二叠统海相碳酸盐岩层中,矿化受古岩溶的严格控制,组成矿石的基本矿物为菱铁矿,金不均匀地分布于菱铁矿矿体内。成矿溶液来自加热的循环地下水。成矿过程分为早阶段和中-晚阶段。成矿温度为156-210℃,成矿深度小于1km.成矿物质主要来自地层。成矿热液呈碱性和弱还原性,贫硫,含中等盐度。矿床形成时的地质构造环境较稳定。这是我国仅见的一种特殊金矿类型。     

Estimation for partially observed Markov processes

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.