Geostatistical inversion of coupled problems: dealing with computational burden and different types of data

Practical application of geostatistical inversion to coupled problems is hampered by a number of difficulties. In this paper, we address two of them: first, the computational cost of sensitivity (Jacobian) matrices and, second, the evaluation of the relative weights of different types of data. Regarding the first, we revise the adjoint state equations to propose a form whose cost is independent of the number of unknown parameters and only grows with the number of observation wells. Regarding the second, we derive expressions for the relative weights of different types of data. These expressions are based on minimizing the expected likelihood, rather than the likelihood itself. The efficiency of both improvements is tested on a synthetic example. The example analyzes a wide range of groundwater flow and solute transport conditions. Yet, the expected likelihood consistently yields the optimal weights. The proposed form of the adjoint state equations leads to one order of magnitude reduction in CPU time with respect to the conventional sensitivity equations.     

吐哈盆地中央构造带正反转演化特征

吐哈盆地中央构造带由火焰山构造和七克台构造组成。中央构造带形成于三叠纪晚期至侏罗纪早期,表现为伸展构造特征,生长断层上盘地层厚度明显大于下盘,并于断层上盘所在的台北凹陷形成沉降中心。晚侏罗世,由于拉萨陆块与欧亚大陆的碰撞作用导致吐哈盆地由伸展盆地转变为挤压盆地,中央构造带也于此时发生构造反转,由早期的伸展正断层转变为挤压逆断层。发生于55Ma的喜山构造事件对天山地区产生了深刻的影响,但影响时间略有滞后,大致发生在晚渐新世至早中新世,中央构造带即在此次构造事件中强烈变形,逆冲出露于地表。     

POSITIVE INVERSION STRUCTURE OF THE CENTRAL STRUCTURE BELT IN TURPAN-HAMI BASIN

The central structure belt in Turpan-Hami basin is composed of the Huoyanshan structure and Qiketai structure formed in late Triassic-early Jurassic, and is characterized by extensional tectonics. The thickness of strata in the hanging wall of the growth fault is obviously larger than that in the footwall, and a deposition center was evolved in the Taibei sag where the hanging wall of the fault is located. In late Jurassic the collision between Lhasa block and Eurasia continent resulted in the transformation of the Turpan-Hami basin from an extensional structure into a compressional structure, and consequently in the tectonic inversion of the central structure belt of the Turpan-Hami basin from the extensional normal fault in the earlier stage to the compressive thrust fault in the later stage. The Tertiary collision between the Indian plate and the Eurasian plate occurred around 55Ma, and this Himalayan orogenic event has played a profound role in shaping the Tianshan area, only the effect of the collision to this area was delayed since it culminated here approximately in late Oligocene-early Miocene. The central structure belt was strongly deformed and thrusted above the ground as a result of this tectonic event.     

Geomechanical model for the post-Variscan evolution of the Permocarboniferous–Mesozoic basins in Northeast Germany

The Permocarboniferous basins in Northeast Germany formed on the heterogeneous and eroded parts of the Variscan orogene and its deformed northern foreland. Transtensional tectonic movements and thermal re-equilibration lead to medium-scale crustal fragmentation, fast subsidence rates and regional emplacement of large amounts of mostly acidic volcanics. The later basin formation and differentiation was triggered by reversals of the large-scale stress field and reactivation of prominent zones of weakness like the Elbe Fault System and the Rhenohercynian/Saxothuringian boundary that separate different Variscan basement domains in the area. The geomechanical behaviour of the latter plays an important role for the geodynamic evolution of the medium to large-scale structural units, which we can observe today in three dimensions on structural maps, geophysical recordings and digital models. This study concentrates on an area that comprises the southern Northeast German Basin, the Saale Basin, the Flechtingen High, the Harz Mountains High and the Subhercynian Basin. The presented data include re-evaluations of special geological and structural maps, the most recent interpretation of the DEKORP BASIN 9601 seismic profile and observations of exposed rock sections in Northeast Germany. On the basis of different structural inventories and different basement properties, we distinguish two structural units to the south and one structural unit to the north of the Elbe Fault System. For each unit, we propose a geomechanical model of basin formation and basin inversion, and show that the Rhenohercynian Fold and Thrust Belt domain is deformed in a thin-skinned manner, while the Mid-German Crystalline Rise Domain, which is the western part of the Saxothuringian Zone, rather shows a thick-skinned deformation pattern. The geomechanical model for the unit north to the Elbe Fault System takes account to the fact that the base of the Zechstein beneath the present Northeast German basin shows hardly any evidence for brittle deformation, which indicates a relative stable basement. Our geomechanical model suggests that the Permocarboniferous deposits may have contributed to the structural stiffness by covering small to medium scale structures of the upper parts of the brittle basement. It is further suggested that the pre-Zechstein successions underneath the present Northeast German basin were possibly strengthening during the Cretaceous basin inversion, which resulted in stress transfer to the long-lived master faults, as indicated for example by the shape of the salt domes in the vicinity of the latter faults. Contrary to this, post-Zechstein successions deformed in a different and rather complex way that was strongly biased by intensive salt tectonic movements.     

Numerical models of tectonic deformation at the Baltica–Avalonia transition zone during the Paleocene phase of inversion

Predictions from dynamic modelling of the lithospheric deformation are presented for Northern Europe, where several basins underwent inversion during the Late Cretaceous and Early Cenozoic and contemporary uplift and erosion of sediments occurred. In order to analyse the evolution of the continental lithosphere, the equations for the deformation of a continuum are solved numerically under thin sheet assumption for the lithosphere. The most important stress sources are assumed to be the Late Cretaceous Alpine tectonics; localized rheological heterogeneities can also affect local deformation and stress patterns. Present-day observations available in the studied region and coming from seismic structural interpretations and stress measurements have been used to constrain the model. Our modelling results show that lateral variation in lithospheric strength below the basin systems in Central Europe strongly controls the regional deformation and the stress regime. Furthermore, we have demonstrated that the geometry of the boundary between Baltica and Avalonia, together with different rheological characteristics of the two plates, had a crucial role on local crustal deformation and faulting regime resulting in the Baltica–Avalonia transition zone from the S–N Alpine convergence.     

小湾电站高边坡系统锚固与排水的优化设计方法

据小湾电站地下厂房进水口开挖高边坡的岩体结构与工程地质特征,采用基于结构面网络模拟的随机楔体稳定分析方法对边坡的稳定性进行了三维分析。进一步运用风险分析理论对边坡系统锚杆与预应力锚索的布置进行了优化设计,最后运用结构面网络模拟理论对排水方案进行了优化设计。     

高耸塔架结构节点损伤基于神经网络的两步诊断法

本文针对高耸钢塔架结构的损伤特点，建立了具有节点损伤的有限元模型，提出了一种分层神经网络两步诊断法：第一步，由基于区域残余力理论的第一层神经网络进行结构损伤区域的初诊；第二步，由基于应变模态理论的第二层神经网络进行损伤区域内的具体损伤节点位置和程度的诊断。对一平面塔架结构的数值仿真分析表明：本文提出的损伤诊断方法的结果是令人满意的。     

比值参数在高密度电法中的应用

探讨了比值参数在高密度电阻率法中的应用；根据工程实例说明比值参数在压制局部干扰和邻近异常分离中有良好效果；将比值参数断面与电阻率反演断面比较研究，能取得较好的异常解释效果。     

重磁遗传算法三维反演中高速计算及有效存储方法技术

将地下场源区域规则划分成很多小长方体单元，并且通过反演确定这些单元的物性变 化，勾画出场源的分布图像，这种方式逐步成为重磁反演，特别是三维反演的重要方向；遗 传算法等非线性技术进行该类反演将逐步成为发展趋势. 本文指出，在应用遗传算法进行该 类反演过程中，隐含着数据量较大时超常规的计算量，它已成为制约该类反演充分发挥作用 的瓶颈问题；同时，本文提出了针对性的分离并存储几何格架的计算策略、以及独特的几何 格架等效压缩存储技术，可以从根本上提高非线性反演计算速度，为该类反演的有效应用奠 定了坚实的基础.     

二维波动方程约束反演的大范围收敛广义脉冲谱方法

针对二维波动方程，将大范围收敛的同伦 方法引入算子参数识别反演的过程，并结合求解不适定问题的Tikhonov正则化方法，设计了二维波动方程反演的大范围收敛广义脉冲谱方法，并利用测井资料进行了约束反演. 大量 的数值模拟结果以及抗噪实验均表明了方法的有效性.