四川盆地五峰-龙马溪组页岩储层非均质性特征

通过对比焦石坝、长宁和威远气田五峰-龙马溪组页岩产层段的岩矿、物性、含气性以及可压裂性等地质参数及其纵向变化规律,综合分析了页岩储层的平面和纵向非均质性及其对“甜点”优选和储层改造的影响。研究表明,纵向非均质特征在不同地区具有总体一致性：石英含量、有机碳含量和含气量均呈现下高上低的基本趋势;底部深水陆棚沉积层段普遍具有高有机碳含量、高脆性、高孔隙度和高含气量的特征,一系列属性参数的相互耦合和匹配决定了五峰组-鲁丹阶成为“甜点”的主要发育层段。平面非均质特征主要表现为储层品质以及“甜点”发育层段和厚度的区域差异：长宁和威远地区的钙质含量明显高于焦石坝地区,发育钙质硅质混合页岩相;焦石坝地区的水平应力差明显低于长宁和威远地区,且天然裂缝更为发育;长宁和焦石坝地区的“甜点层”发育于五峰-鲁丹阶下部,而威远地区的“甜点层”则发育于鲁丹阶-埃隆阶中部。页岩储层的非均质性决定了“甜点”的分布范围、发育层段和厚度;矿物组成、天然裂缝、岩石力学参数和地应力场的空间变化是进行压裂设计的重要依据。     

天津滨海新区地热资源循环利用研究——馆陶组热储回灌技术研究与示范

可持续回灌是地热资源循环采用的基本。本文通过热储敏感性、堵塞原因、成井结构以及地面装置分析了回灌衰减的原因。认为物理堵塞是回灌堵塞的主要原因,成井结构是维持可持续回灌的关键,标准的地面设备可避免由于物理和化学原因引起的堵塞。本文采用成功的成井结构和标准的地面装置开展馆陶组热储回灌地热资源循环采用的示范工程建设,回灌率达到100~123 m3/h,实现了天津滨海新区馆陶组热储回灌突破,保证了滨海新区馆陶组地热资源循环利用。     

Operations optimization of multireservoir systems using storage moments equations

Fletcher–Ponnambalam presented a new model for considering the balance equation of the storage volume of the reservoir using indicator functions. For stochastic inflows, the two storage moments of this balance equation, namely, the mean and variance, calculated using a random release policy were found to be quite accurate unlike any known models. Significantly, this model required no discretization of storage volumes or releases. In this paper, this work has been extended to two new cases: for multireservoir systems which require further consideration of the stochastic releases and for arbitrary distribution of inflows using the Beta-equivalent Kumaraswamy distribution which has a simpler form than Beta. The randomized release policies are easy to use even in a multireservoir problem. The Parambikulm-Aliyar Project from India is used as a case study. Results show accurate predictions of mean of storages in the multireservoir case but show the need for further improvement in the standard deviations of storages. The optimal benefits and the policies obtained are shown to be at least as good as obtained using Monte Carlo-based methods.     

基于双相介质理论的储层参数反演方法

传统基于单相介质理论的储层参数反演方法将孔隙流体与固体骨架等效为单一固体,弱化了孔隙流体的影响,反演结果精度不高.本文提出根据双相介质理论反演储层参数的方法.首先,在前人研究的基础上,利用岩石物理模型建立弹性参数与孔隙度、饱和度、泥质含量等储层参数间的关系,进而将双相介质反射系数推导为储层参数的函数;其次,根据贝叶斯反演理论,在高斯噪声假设的基础上,采用更加符合实际情况的修正柯西分布函数描述反射系数的稀疏性,推导出储层物性参数目标反演函数;最后,应用差分进化非线性全局寻优算法来求解目标反演函数,使得反演结果与实际资料间误差最小.新方法旨在突出流体对介质反射系数的影响,以期得到较高的储层参数反演精度.模型与实际资料测试均表明该方法可行、有效且反演精度较高.     

龙滩水库蓄水后库区中小地震震源机制

本文在对龙滩库区2006年9月30日至2013年5月26日发生的3682次地震进行精定位的基础上,利用FOCMEC方法和改进的格点尝试法反演了ML2.0以上地震震源机制及区域构造应力场,并综合龙滩库区地层岩性、断裂构造和渗透条件等资料,探讨了地震活动类型与库区蓄水过程的关系,获得以下认识:(1)龙滩水库蓄水后的地震活动主要丛集在罗妥、八茂、拉浪、坝首和布柳河5个深水区,地震类型以逆断为主,正断和走滑也占有一定比例.蓄水初期,库区地震类型呈现多样性,蓄水约4年3个月后,地震主要发生在浅部地层中,并大多为逆断型地震.(2)5个地震丛构造应力场最大主应力方向以NW—SE为主,倾角均较小,中等和最小主应力分布较凌乱,在此构造应力环境中,龙滩库区主要断裂整体上呈现走滑运动的性质,局部伴生不同程度的逆倾滑或正倾滑运动分量.(3)库区深、浅部地震活动水平和地震性质之所以会随蓄水过程发生变化,可能与深、浅部构造应力环境、岩体力学性质和渗透性能的差异有关.     

Direct forecasting of subsurface flow response from non-linear dynamic data by linear least-squares in canonical functional principal component space

Inverse modeling is widely used to assist with forecasting problems in the subsurface. However, full inverse modeling can be time-consuming requiring iteration over a high dimensional parameter space with computationally expensive forward models and complex spatial priors. In this paper, we investigate a prediction-focused approach (PFA) that aims at building a statistical relationship between data variables and forecast variables, avoiding the inversion of model parameters altogether. The statistical relationship is built by first applying the forward model related to the data variables and the forward model related to the prediction variables on a limited set of spatial prior models realizations, typically generated through geostatistical methods. The relationship observed between data and prediction is highly non-linear for many forecasting problems in the subsurface. In this paper we propose a Canonical Functional Component Analysis (CFCA) to map the data and forecast variables into a low-dimensional space where, if successful, the relationship is linear. CFCA consists of (1) functional principal component analysis (FPCA) for dimension reduction of time-series data and (2) canonical correlation analysis (CCA); the latter aiming to establish a linear relationship between data and forecast components. If such mapping is successful, then we illustrate with several cases that (1) simple regression techniques with a multi-Gaussian framework can be used to directly quantify uncertainty on the forecast without any model inversion and that (2) such uncertainty is a good approximation of uncertainty obtained from full posterior sampling with rejection sampling.     

新丰江水库上游氮磷污染的时空变化

新丰江水库流域是广东省最大的饮用水源保护区,对东江下游和香港地区的用水安全有十分重要的作用,水质长期保持在国家地表水I类标准。但随着经济的发展,水库上游地区人类活动加剧,对水质和水量造成潜在威胁。通过2007、2008年对新丰江水库上游流域采样分析氮磷质量分数,发现总氮和总磷超过地表水Ⅱ类标准的采样点比例分别为48%和12%。与总磷相比,新丰江水库上游总氮的污染要更加严重,丰水期超标点明显少于枯水期,而枯水期氮-磷正相关关系较好（P＜0.000 1,R²＝0.868 3）,污染物可能为同一污染源或同期释放。空间上超标点多位于新丰、连平县城和大型城镇居民点如马头、忠信、灯塔镇附近,河流污染主要来自生活污水和养殖废水等人为点源污染。尽管目前新丰江水库库区的水质情况总体较好,但上游水质也存在恶化风险,可能导致水库自净能力降低,水体富营养化加剧,因此需要对新丰江水库上游的水质进行定期监测,通过控制内、外源性氮磷措施,以减少饮用水源安全风险和保证新丰江水库水质长期稳定。     

大容量气枪震源子波激发特性分析

大容量气枪水库激发作为陆地震源的可行性与有效性已经得到成功验证.为进一步提高气枪震源激发效果,本文通过水库气枪激发试验对单枪容量为2000 in³的气枪震源激发子波特征及规律进行了研究.依据近场水听器和远场短周期地震仪记录数据,分析气枪震源沉放深度、工作压力等不同激发条件对压力脉冲和气泡脉冲的影响.有助于人们根据不同尺度地下结构探测对震源激发信号的要求,调整气枪激发参数和激发环境,获得最佳激发效果.试验结果表明：（1）沉放深度对压力脉冲波形影响较小,其优势频率不随沉放深度而改变;（2）随着沉放深度从5 m增加到11 m,气泡脉冲的优势频率由5 Hz增加至7 Hz,其最大振幅亦近线性递增;（3）工作压力越大,激发压力脉冲能量越强,而对气泡脉冲的影响主要体现在主频降低.适合远距离深穿透地下结构探测的低频信号主要来自大容量气枪所激发气泡的反复振荡,由于气枪振荡过程非常复杂,本文通过较为简洁的数学和物理模型进行了解释.     

Application of Multiple Point Geostatistics to Non-stationary Images

Simulation of flow and solute transport through aquifers or oil reservoirs requires a precise representation of subsurface heterogeneity that can be achieved by stochastic simulation approaches. Traditional geostatistical methods based on variograms, such as truncated Gaussian simulation or sequential indicator simulation, may fail to generate the complex, curvilinear, continuous and interconnected facies distributions that are often encountered in real geological media, due to their reliance on two-point statistics. Multiple Point Geostatistics (MPG) overcomes this constraint by using more complex point configurations whose statistics are retrieved from training images. Obtaining representative statistics requires stationary training images, but geological understanding often suggests a priori facies variability patterns. This research aims at extending MPG to non-stationary facies distributions. The proposed method subdivides the training images into different areas. The statistics for each area are stored in separate frequency search trees. Several training images are used to ensure that the obtained statistics are representative. The facies probability distribution for each cell during simulation is calculated by weighting the probabilities from the frequency trees. The method is tested on two different object-based training image sets. Results show that non-stationary training images can be used to generate suitable non-stationary facies distributions.     

The genetic relationship between mafic dike swarms and plutonic reservoirs in the mesozoic of central chile (30°–33°45′S): insights from AMS and geochemistry

Five mafic dike swarms between 30° and 33°45′S were studied for their geochemical signature and kinematics of magma flow directions by means of AMS data. In the Coastal Range of central Chile (33°−33°45′S), Middle Jurassic dike swarms (Concón and Cartagena dike swarms, CMDS and CrMDS, respectively) and an Early Cretaceous dike swarm (El Tabo Dike Swarm, ETDS) display the presence of dikes of geochemically enriched (high-Ti) and depleted (low-Ti) basaltic composition. These dikes show geochemical patterns that are different from the composition of mafic enclaves of the Middle Jurassic Papudo-Quintero Complex, and this suggests that the dikes were injected from reservoirs not related to the plutonic complex. The mantle source appears to be a depleted mantle for Jurassic dikes and a heterogeneous-enriched lithospheric mantle for Cretaceous dikes. In the ETDS, vertical and gently plunging magma flow vectors were estimated for enriched and depleted dikes, respectively, which suggest, together with variations in dike thickness, that reservoirs were located at different depths for each dike family. In the Elqui Dike Swarm (EDS) and Limarí Mafic Dike Swarm (LMDS), geochemical patterns are similar to those of the mafic enclaves of the Middle Jurassic plutons. In the LMDS, east to west magma flow vectors are coherent with injection from neighbouring pluton located to the east. In the EDS, some dikes show geochemical and magma flow patterns supporting the same hypothesis. Accordingly, dikes do not necessarily come from deep reservoir; they may propagate in the upper crust from coeval shallow pluton chamber. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.