Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.     

川陕丁家林金矿区探矿权的评估

为准确、客观估价丁家林金矿区探矿权价值，根据矿区地质报告中提交的储量级别、矿体控制程度和工程质量，确定了基础和可采储量、生产规模、矿山服务年限、贴现率、有关实物工作量等6项参数。综合运用贴现现金流量和重置成本法对探矿权进行评估，对达到详查程度的基础储量采用贴现现金流量法，控制程度低的资源量采用重置成本法，二者之和综合反映探矿权的价值。此方法适用于地质工作程度不均衡矿区的探矿权评估，贴现现金流量法要充分考虑储量的可信度和获利临界品位．重置成本法关键是准确确定有关的实物工作量。     

Depositional processes on oceanic island shelves – Evidence from storm‐generated Neogene deposits from the mid‐North Atlantic

Oceanic islands – such as the Azores in the mid‐North Atlantic – are periodically exposed to large storms that often remobilize and transport marine sediments along coastlines, and into deeper environments. Such disruptive events create deposits – denominated tempestites – whose characteristics reflect the highly dynamic environment in which they were formed. Tempestites from oceanic islands, however, are seldom described in the literature and little is known about storm‐related sediment dynamics affecting oceanic island shelves. Therefore, the geological record of tempestite deposits at oceanic islands can provide invaluable information on the processes of sediment remobilization, transport and deposition taking place on insular shelves during and after major storms. In Santa Maria Island (Azores), a sequence of Neogene tempestite deposits was incorporated in the island edifice by the ongoing volcanic activity (thus preserved) and later exposed through uplift and erosion. Because it was overlain by a contemporary coastal lava delta, the water depth at the time of deposition could be inferred, constituting an excellent case‐study to gain insight on the still enigmatic processes of insular shelf deposition. Sedimentological, palaeontological, petrographic and palaeo‐water depth information allowed the reconstruction of the depositional environment of these sediments. The sequence typifies the characteristics of a tempestite (or successive tempestites) formed at ca 50 m depth, in a steep, energetic open insular shelf, and with evidence for massive sediment remobilization from the nearshore to the middle or outer shelf. The authors claim that cross‐shelf transport induced by storm events is the main process of sediment deposition acting on steep and narrow shelves subjected to high‐energetic environments, such as the insular shelves of open‐sea volcanic islands.     

地下资源开采与诱发地震：现状与进展

页岩气、地热等清洁能源正成为全球能源开发日益重要的组成部分,也是我国未来能源产量增长的主体,对保障能源资源安全、支撑双碳目标意义重大。然而随着页岩气、地热等地下工业开采活动规模的不断扩大,由此引发的地震灾害问题正日益受到人们的关注。强化工业活动与诱发地震关系研究、管控相关诱发地震灾害风险成为全球面临的重要课题和挑战。基于此,笔者等调研了工业活动诱发地震的全球典型案例,系统介绍了不同国家(美国、加拿大和韩国)工业活动诱发地震的现状及其差异,并对流体注入诱发地震的主要途径、物理机制、相关作业参数和诱发地震关系等做了全面的综述。发现诱发地震不仅与流体注入参数(体积、速率等)相关,还受控于特定的区域地质背景和断裂条件。合理管控地震灾害风险、促进相关能源的开采,应该统筹开采区地质地球物理的综合调查和开采过程的动态监测。     

煤层采空区“三相”的激发极化特征

激发极化法对煤层采空区、采空积水区、裂隙富水区等情况的探测,可以取得不同的响应特征,具有良好的区分效果。将该方法的物理-电学特性与土力学理论相结合,从土力学中的\     

圆回线和矩形回线源瞬变电磁响应差异性计算

在回线源瞬变电磁发展的早期,往往使用圆回线产生的场来近似模拟多边形回线的场。为了对这种近似的影响进行分析,本文以矩形回线为例,对圆回线近似模拟多边形回线的误差进行计算。首先,根据给出的圆回线和矩形回线内外任一点的感生电动势新的计算方法,对两种回线装置内的感生电动势分布进行对比。研究结果表明:在扩散场早期,圆回线和矩形回线场的等值线分布形态存在较大的差别,而在扩散场中晚期,两种回线感生电动势相对误差较小。两种回线感生电动势相对误差的变化特征可以用感应电流烟圈扩散的理论进行解释。通过理论模型的分析可知,圆回线近似矩形回线会对浅层电性结构分析带来较大影响。     

地下开采引发地面沉陷的未确知聚类预测方法

对未确知聚类预测法进行优化,并将其应用于开采地面沉陷的预测研究。采用开采地面沉陷的实测数据按最大沉陷量进行分类,利用各分类影响因素的均值表示各分类中心,并确定各影响因素的未确知测度函数。由待测对象指标的综合未确知测度与各分类指标的未确知测度间的距离来确定待预测对象所属等级,给出了预测值的计算公式。经计算验证,该方法的正确率为75％。但在实际应用中,为了保证地表建筑设施等更加安全,允许预测级高判,则正确率可达100％。针对某铁矿一观测点进行预测,并与实测数据比较,结果表明,未确知聚类预测的结果是令人满意的,为开采地面沉陷的预测提供了一种新思路。     

Characteristics of rapid giant landslides in China

Factual data for 70 rapid, giant landslides since 1900 show that the occurrence of these landslides was largely predisposed by tectonics, geological structures, lithology and topography, and often triggered by rainfall and earthquakes. In terms of mobile behavior, the giant landslides can be classified into three types: slides, slide-flows and flows. It is found that each type of landslide was constrained to certain geologic and topographic regimes. There are good correlations between kinematic parameters of landslides and slope geometries, which confirm the important role played by topographical condition in the mobile behavior of landslides. Moreover, it is also found that each type of landslide presents distinct geotechnical characteristics in terms of nature of the slip zone and properties of sliding mass. Brief analyses of five typical cases illustrate that landslide mechanisms can be conceptually depicted by failure mechanisms of their slip zones prior to onset of movement and following energy conversion during movement. Problems and questions related to experience in China suggest that comprehensive and systematic investigation and study on rapid giant landslides are urgently needed.     

Relative molar mass distributions of chromophoric colloidal organic matter in coastal seawater determined by Flow Field-Flow Fractionation with UV absorbance and fluorescence detection

A new method for the characterization of chromophoric colloidal organic matter in seawater has been applied to samples from the Baltic Sea, Kattegatt and Skagerrak seas. Size fractionation of the sample by Flow Field-Flow Fractionation and measurement of the fluorescent and UV absorbing properties of the individual size fractions result in a relative molar mass distribution (RMM) of the optical properties. The RMM distributions have been used to estimate number and weight average relative molar masses, and polydispersity indices. At least two sources of coloured organic matter were identified from the ratio of fluorescence to UV: the Baltic surface water and the Skagerrak deep water. The dominating processes were mixing and dilution, but processes such as photo bleaching of fluorescence are also believed to be important. The RMM distribution derived from UV detection (1150–1300 Dalton) increased with increasing salinity while that derived for fluorescence (1500–1250 Dalton) decreased with increasing salinity. The specific UV absorbance taken as a proxy of the aromaticity of the chromophoric organic material showed decreasing trend with both increasing salinity and increasing UV derived weight average relative molar mass. Increasing polydispersity of the colloidal material was also observed as a function of salinity.     

Mitigation of Ice-Induced Vibrations for Offshore Platforms Using Tuned Mass Damper

The problems of ice-induced vibration have been noticed and concerned since the 1960s, but it has not been well resolved. One reason is that the dynamic interaction between ice and structure is so complicated that practical ice force model has not been developed. The recent full-scale tests conducted on jacket platforms in the Bohai Sea presented that ice could cause intense vibrations which endanger the facilities on the deck and make discomfort for the crew. In this paper, the strategy of mitigation of ice-induced offshore structure vibration is discussed. Based on field observations and understanding of the interaction between ice and structure, the absorption mitigation method to suppress ice-induced vibration is presented. The numerical simulations were conducted for a simplified model of platform attached with a Tuned Mass Danlper （TMD） under ice force function and ice force time history. The simulation results show that TMD can fa- vorably reduce ice-induced vibrations, therefore, it can be considered to be an alternative approach to utilize. Finally, the application possibilities of utilizing TMDs on other miniature offshore structures in ice-covered areas of marginal oil fields are discussed.