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.     

Accurate local upscaling with variable compact multipoint transmissibility calculations

We propose a new single-phase local upscaling method that uses spatially varying multipoint transmissibility calculations. The method is demonstrated on two-dimensional Cartesian and adaptive Cartesian grids. For each cell face in the coarse upscaled grid, we create a local fine grid region surrounding the face on which we solve two generic local flow problems. The multipoint stencils used to calculate the fluxes across coarse grid cell faces involve the six neighboring pressure values. They are required to honor the two generic flow problems. The remaining degrees of freedom are used to maximize compactness and to ensure that the flux approximation is as close as possible to being two-point. The resulting multipoint flux approximations are spatially varying (a subset of the six neighbors is adaptively chosen) and reduce to two-point expressions in cases without full-tensor anisotropy. Numerical tests show that the method significantly improves upscaling accuracy as compared to commonly used local methods and also compares favorably with a local–global upscaling method.     

葛店矿双庙扩大区三维地震勘探

葛店煤矿双庙扩大区位于豫、皖两省交界处。勘探区为一轴向北北东向的不对称向斜构造,地层倾角26°～30°,测区深部有岩浆岩侵入体,严重影响煤层赋存形态;地表条件复杂,村庄密集、河流较多,测区东侧塌陷区水塘大面积分布．三维地震勘探野外施工难度大。本区地震勘探观测系统设计为12线18炮制,为避免偏移归位及叠加次数降低等问题,在向斜轴部两侧、村子周边及河岸处适当加密炮点,并在塌陷区利用钢钎插置检波器以保证CDP点完整性。数据处理采取了叠前反褶积、叠前时间偏移等技术,数据处理效果显著。遵循多波组对比、多个数据体综合解释的原则,以纵向、横向和任意时间剖面相结合,充分利用三维可视化技术进行地震资料解释。本次勘探成果查明了主采煤层三2^2、二2煤层的构造形态及埋藏深度;解释落差5m以上断层34条,其中修改原有断层5条;另外新发现小于5m的异常断点22个。经巷道验证煤层底板标高误差0．43％～1．96％,巷道已揭露断层7条,与解释成果基本一致。     

第二次土地调查若干问题探讨

第二次土地调查的目的是全面查清全国土地利用状况,掌握真实的土地基础数据,如何保证调查的质量?文章以江阴市为例进行探讨。     

基坑潜水疏干的简化计算在昊华水泥厂中的应用

以流线、流面、汇点的概念为基础,对稳定流双井干扰和直线隔水边界附近涌水量理论公式进行对比分析,提出了二个虚拟界面,其中虚拟界面Ⅰ,运用流线、流面的性质,流线方程等给出证明;虚拟界面Ⅱ则通过半无限条形降落漏斗的分析,应用元流和总流的能量方程得到流量为零,流线为零的平面。在同样条件下,条形无限涌水量是半无限潜含水层涌水量的二倍。应用总流能量方程对三种情况水头损失的分析,解释了这种关系存在的合理性,得出虚拟界面Ⅱ,并以此得出该界面内的最大残余水头计算公式。将基坑降水运用虚拟界面简化为扇形,条形半无限含水层,从而实现单井预测,该方法应用到昊华水泥厂基坑降水中,预测效果理想。     

Group Time Delay in the Reflection of P-Polarized Electromagnetic Plane Waves at a Stratigraphic Interface

The group time delay (GTD) formulae of quasi-total reflection of inhomogeneous P-polarized electromagnetic waves (P-PEW, whose electric field parallels to the plane of incidence), caused by the Goos–Hänchen effect, is derived by the phase shift of the wave. A numerical example where the frequency equals to 1 GHz is given. The curves of calculating results show that there are three discontinuous points at the critical angle of phase shift vector, the critical angle of attenuation vector and 90^○. When the angle of incidence is equal to one of these three angles, the GTD will become infinite, and the electromagnetic wave will propagate along the interface. The GTD will be very large, when the angle of incidence is close to one of these three angles. The results indicate that there is an important relationship between the two critical angles and the conductivity and permittivity of the two strata, and that the GTD has an important relationship with the angle of incidence. These results can be used to determine the lithology of the strata and to delimit the interfaces of strata more effectively. It is suggested that this research may prove useful in electromagnetic logging analysis and, perhaps, in the design of logging instruments.     

随钻测录井地质导向二维分解实时绘图方法

根据随钻测录井实时地质导向和大斜度、水平井评价成图技术需求，针对传统绘图方法存在的弊端，提出了将测录井信息、井眼轨迹和地质模型进行二维分解的实时绘图方法。针对二维分解绘图方法绘图时空复杂度较高的问题，给出了不同事件驱动下的局部实时计算和拷屏重绘算法，控制了对CPU和内存的消耗、提高了绘图效率，消除了实时绘图的闪烁和卡顿现象。应用实例表明，二维分解实时绘图方法能够实现大尺度随钻测录井地质导向图形的流畅、无卡顿实时绘图，可提高大斜度、水平井储层模型评价的刻画精度和时效。     

中国大陆第四纪地壳运动时程

本文从中国大陆有代表性的青藏高原、黄土高原、三峡地区着手,进行了新构造运动剖析。并结合全国火山活动期次,总结出中国大陆第四纪地壳运动的时程变化。从而为进一步探讨新构造运动特点和规律,为国民经济建设及人类生存环境提供了一定的地质背景资料。     

Treatment of high fluoride concentration wastewater by layered double hydroxides： Mechanism studies

An appropriate concentration of fluoride in drinking water is required to prevent dental cavities, but long-term ingestion of water that contains more than a suitable level of fluoride can cause bone disease and mottling of the teeth. Fluoride ions can be found in wastewater from the fluoride chemical industry, as well as the semiconductor, metal processing, fertilizer, and glass-manufacturing industries. The discharge standard for fluoride in industrial wastewater in China is 10 mg/L. Efficient treatment of fluoride-containing wastewater is therefore of major concern in China, following the rapid development of the fluoride chemical industry. Several methods have been used to remove fluoride from water, such as adsorption, chemical precipitation, electrodialysis, ion exchange and electrochemical processes. Layered double hydroxides （LDHs） are anionic clays with high anion exchange capacities which are effective adsorbents for removal of a variety of anionic pollutants. LDHs have been studied as potential adsorbents for removing toxic anionic species such as CrO4^2-, TcO4^-, SeO3^2-, pesticides, and anionic surfactants from aqueous systems. One of the main attractions of using LDHs for fluoride removal, is that unlike other chemical treatment methods, no chemical sludge should be produced. In the present study, an attempt was made to investigate the mechanism of fluoride removal by LDHs under different conditions using batch methods. In addition, the release of fluoride adsorbed on LDHs by treatment with an aqueous solution of Na2CO3 was studied. The residual fluoride was found to be 10 mg/L in a solution with an initial concentration of 1000 mg/L, which meets the discharge standard for fluoride in industrial wastewater in China.     

Structural Evolvement of Heating Treatment of Mg/AI-LDH and Preparation of Mineral Mesoporous Materials

Although hydrotalcite, or layered double hydroxides （LDHs）, is not a common mineral, it is an important material that can be easily synthesized in laboratory. In this study, structural evolvement and BET surface area changes of heat treated Mg/AI-LDH is evaluated by XRD, TEM and N2-BET analyses. The results indicate that the magnesium-aluminum LDH with carbonate as interlayer anion, periclase-like oxides was formed at temperatures of 400-800℃. Meanwhile, 2-3 nanometer mesoporous were formed during decomposition of LDH. However, the heat treated samples still preserve the morphology of the original LDH plates. Periclase-like formed from LDH heat treatment may re-hydrolyze and recover the structure of LDH. However, crystallinity of the recovered LDH is lower than that of the original LDH. This heat treatment will result in formation of （Mg, Al）-oxide nano-crystals and nanopores among the nano-crystals. When heating temperature exceeds 1000, the periclase-like （Mg, Al）-oxide is transformed into a composite with periclase （MgO） and spinel phases. The periclase can be re-hydrolyzed and dissolved in HCl solution. After acid treatment, the sample with a high surface area is composed of spinel nano-crystals and nanopores among them. Our results will provide a new and economic way to synthesize mesoporous materials through pathways of phase transformation of precursor materials with different composition.