A nonlinear multiparameter prestack seismic inversion method based on hybrid optimization approach

Multiparameter prestack seismic inversion is one of the most powerful techniques in quantitatively estimating subsurface petrophysical properties. However, it remains a challenging problem due to the nonlinearity and ill-posedness of the inversion process. Traditional regularization approach can stabilize the solution but at the cost of smoothing valuable geological boundaries. In addition, compared with linearized optimization methods, global optimization techniques can obtain better results regardless of initial models, especially for multiparameter prestack inversion. However, when solving multiparameter prestack inversion problems, the application of standard global optimization algorithms maybe limited due to the issue of high computational cost (e.g., simulating annealing) or premature convergence (e.g., particle swarm optimization). In this paper, we propose a hybrid optimization-based multiparameter prestack inversion method. In this method, we introduce a prior constraint term featured by multiple regularization functions, intended to preserve layered boundaries of geological formations; in particular, to address the problem of premature convergence existing in standard particle swarm optimization algorithm, we propose a hybrid optimization strategy by hybridizing particle swarm optimization and very fast simulating annealing to solve the nonlinear optimization problem. We demonstrate the effectiveness of the proposed inversion method by conducting synthetic test and field data application, both of which show encouraging results.     

巴音戈壁盆地南部塔木素铀矿床水岩作用特征及其与铀成矿关系研究

在较封闭的环境中,塔木素铀矿床高矿化度(平均矿化度为35.4 g/L)NaCl型地下水中的Na+替换了斜长石中的Ca2+,使斜长石由更长石(An≈13)转变为钠长石(An≈1),被替换出的Ca_2+(以CaO计,约占斜长石总质量的4.5%)与地下水中的HCO_3~-、CO_3~(2-)及Mg~(2+)等结合形成了白云石等碳酸盐胶结物。在此过程中发生了脱碳作用,促使地下水中以[UO_2(CO_3)_3]~(4-)、[UO_2(CO_3)_3]~(2-)等碳酸铀酰络合离子及Mg CO_3·Na_2UO_2(CO_3)_2复盐发生分离而形成了铀的沉淀,并受扩散作用影响,水中的铀趋向于向水岩作用相对强烈地段迁移,从而促使铀在特定的层位集中、富集。同时,斜长石因水岩作用(溶解、溶蚀等),在解理面及表面形成了次生的缝隙及孔洞等,为铀沉淀提供了空间。成岩后,含SO_4~(2-)等的酸性地表水沿层间下渗,溶解了砂岩中碳酸盐胶结物而形成了溶洞,为后期再次迁移的铀提供了沉淀空间,并形成了铀的进一步叠加、富集。     

车排子凸起与沙湾凹陷过渡带的构造模型

车排子凸起位于准噶尔盆地扎伊尔山前,是夹持在南部四棵树凹陷和东部沙湾凹陷之间的三角构造带。然而对这个三角构造带是原地掩伏构造,山前推覆构造,还是走滑构造控制形成的认识存在争议。目前对于车排子凸起与沙湾凹陷过渡带的构造变形特征以及变形机制的认识尚不清楚。而过渡带是研究该区构造变形与演化的良好区域,精细研究该区域构造几何学与运动学对于揭示其构造变形机制以及认识盆山之间的关系具有积极作用,并对于该区油气勘探具有重要意义。本文以研究区的连片三维地震资料为基础,运用不整合面识别、轴面分析及平衡剖面等方法,对车排子凸起与沙湾凹陷过渡带进行精细的构造解析、建模与构造复原。研究表明,过渡带发育下二叠统（P₁j-P₁f）,中二叠统（P₂x-P₂w）,上二叠统（P₃w）,三叠系（T）,侏罗系（J）,白垩系（K）与古近系-第四系（E-Q）7个构造变形层序。发育佳木河组底部及石炭系内部滑脱层。深部（前二叠系）构造楔发育于整个过渡带,其个数和几何学形态的变化直接影响了上覆背斜形态。构造楔在早二叠世开始楔入,挤压作用持续至三叠纪末期。晚三叠世发育生长地层,显示了P-T构造层褶皱变形的时间,随后三叠纪末期,过渡带受到来自西北缘造山带方向的逆冲推覆构造的叠加作用,表现在后缘发育一系列叠瓦构造,推覆距离约为8 km,南部逆冲断层上盘二叠系-三叠系被剥蚀完全,北部残余少量下二叠统佳木河组。     

青藏高原北缘前寒武纪地质演化：进展与讨论

青藏高原北缘前寒武纪块体主要出露于塔里木盆地南缘的铁克里克和阿尔金、祁连山及柴达木盆地周缘等地区。这些古老块体分布于相对年轻的古生代秦祁昆造山带之中,大多经历了古生代造山事件的改造叠加。近10年来,对该地区前寒武纪地层时代格架、物质组成和构造属性研究取得了许多进展：1）揭示铁克里克、北阿尔金、全吉和阿拉善地块具有相似的2.5～2.2 Ga、2.1～1.7 Ga的基底演化和1.7 Ga之后的盖层沉积记录,反映塔里木克拉通、全吉地块和华北克拉通在Columbia 超大陆中具有相关和相似性。2）一系列地层时代被重新厘定,为建立新的前寒武纪地层划分方案提供了依据;3）厘定出铁克里克新元古代裂谷火山—沉积盆地,其可以与华北克拉通南缘1 000～830 Ma的裂解事件对比,它们在新元古代早期可能构成同一个克拉通裂解边缘;4）在南阿尔金与祁连、柴达木等地块识别出性质相似的新元古代岩浆作用和变质作用,它们可能共同代表一个新元古代早期活动大陆边缘的构造杂岩,为古生代汇聚、增生于塔里木克拉通或华北克拉通南缘的外来地块,存在于秦祁昆造山带之中。目前关于青藏高原北缘前寒武纪块体的物质组成、时代格架、构造属性和古环境变化研究还存在许多问题,仍需要开展进一步的大量工作。     

Rainfall threshold determination for flash flood warning in mountainous catchments with consideration of antecedent soil moisture and rainfall pattern

Flash flood disaster is a prominent issue threatening public safety and social development throughout the world, especially in mountainous regions. Rainfall threshold is a widely accepted alternative to hydrological forecasting for flash flood warning due to the short response time and limited observations of flash flood events. However, determination of rainfall threshold is still very complicated due to multiple impact factors, particular for antecedent soil moisture and rainfall patterns. In this study, hydrological simulation approach (i.e., China Flash Flood-Hydrological Modeling System: CNFF-HMS) was adopted to capture the flash flood processes. Multiple scenarios were further designed with consideration of antecedent soil moisture and rainfall temporal patterns to determine the possible assemble of rainfall thresholds by driving the CNFF-HMS. Moreover, their effects on rainfall thresholds were investigated. Three mountainous catchments (Zhong, Balisi and Yu villages) in southern China were selected for case study. Results showed that the model performance of CNFF-HMS was very satisfactory for flash flood simulations in all these catchments, especially for multimodal flood events. Specifically, the relative errors of runoff and peak flow were within?±?20%, the error of time to peak flow was within?±?2 h and the Nash–Sutcliffe efficiency was greater than 0.90 for over 90% of the flash flood events. The rainfall thresholds varied between 93 and 334 mm at Zhong village, between 77 and 246 mm at Balisi village and between 111 and 420 mm at Yu village. Both antecedent soil moistures and rainfall temporal pattern significantly affected the variations of rainfall threshold. Rainfall threshold decreased by 8–38 and 0–42% as soil saturation increased from 0.20 to 0.50 and from 0.20 to 0.80, respectively. The effect of rainfall threshold was the minimum for the decreasing hyetograph (advanced pattern) and the maximum for the increasing hyetograph (delayed pattern), while it was similar for the design hyetograph and triangular hyetograph (intermediate patterns). Moreover, rainfall thresholds with short time spans were more suitable for early flood warning, especially in small rural catchments with humid climatic characteristics. This study was expected to provide insights into flash flood disaster forecasting and early warning in mountainous regions, and scientific references for the implementation of flash flood disaster prevention in China.     

Water system characteristics of Karst river basins in South China and their driving mechanisms of hydrological drought

Natural Hazards - It is common knowledge that drought is considered one of the most damaging natural disasters in terms of economic costs, societal problems, and ecological impacts. In this study,...     

应用陆坡形态-迁移轨迹组合法寻找大型富砂深水扇

在海外选区研究阶段,盆地勘探程度极低,资料贫乏,导致常规深水沉积分析工作很难开展。而陆坡形态与陆坡迁移轨迹富含丰富的地质信息,可借助二者的组合对深水沉积进行宏观分析。以二维地震资料为基础,应用陆坡形态-迁移轨迹组合法寻找大型富砂深水扇的有利发育区。A盆地发育上升型、平直型、下降型三种陆坡迁移轨迹类型,其中平直-下降型迁移轨迹对应的陆坡进积距离更大,更有利于发育大型深水扇沉积;同时A盆地发育上凸型、下凹型、直线型和坡坪型四种陆坡类型,上凸型与下凹型陆坡的坡度更缓,宽度更大,更有利于重力流的砂泥分异,从而形成富砂深水扇沉积。将同一时期的大型、富砂深水扇的有利分布区相叠合,重叠区域即为该时期大型富砂深水扇的有利发育区,主要位于盆地东部拐角带与盆地西部。同时应用组合法还可建立盆地漂移期的沉积演化模式,分析漂移期沉积演化历史。     

江南造山带黄金洞金矿床成矿机制:矿物形成环境与金成矿物理化学条件制约

黄金洞超大型金矿床位于江南造山带中段,赋存于新元古界浅变质岩系中,受控于NNE-NE向长平断裂带,金资源储量达100吨。该矿床可由南至北分为金福、金塘、杨山庄和曲溪矿段,主要矿化类型有石英-硫化物脉型、构造角砾岩型和黄铁毒砂绢英岩型。金属矿物主要发育有自然金、毒砂、黄铁矿、黄铜矿、闪锌矿、方铅矿、磁黄铁矿等;非金属矿物有石英、绢云母、方解石和菱铁矿等,其中金主要以自然金与不可见金形式存在。根据野外与镜下观察,金成矿作用分为Ⅰ石英-金-毒砂-黄铁矿、Ⅱ石英-金-多金属硫化物-白钨矿和Ⅲ石英-辉锑矿-绿泥石3个阶段,前二者为主要成矿阶段。曲溪矿段Ⅱ阶段毒砂相对不发育、而磁黄铁矿和自然金显著发育,绿泥石主要发育于Ⅲ阶段中,与辉锑矿及闪锌矿共生。根据不同矿段各阶段毒砂与Ⅲ阶段绿泥石成分,计算其温度、lgf(S_2)与lgf(O_2),可见Ⅰ阶段成矿温度与硫逸度高于Ⅱ阶段:杨山庄矿段两阶段成矿温度分别为300～378℃、260～300℃,lgf(S_2)分别为-11～-7.2、-11.9～-10.1;金塘两阶段成矿温度为240～311℃、245～298℃;金福Ⅱ阶段成矿温度上限为297℃;曲溪矿段成矿温度为268～368℃,Ⅱ阶段毒砂lgf(S_2)与Ⅲ阶段绿泥石lgf(O_2)分别为-13.2～-8.7、-50.9～-40.1。根据不同阶段矿物之间的相互关系及成矿温度与硫逸度演化特征,推断Ⅰ、Ⅱ成矿阶段伴随强烈的硫化作用,金以类质同象方式进入毒砂和黄铁矿中,形成不可见金;其中Ⅱ阶段由于成矿流体压力骤降,含金流体发生相分离作用,H2S等气体大量逃逸,导致成矿流体中硫含量骤降,加以硫化作用持续消耗流体中的硫,促进了含金络合物分解与自然金的沉淀。     

南阿尔金泥盆纪中基性杂岩体的成因及其地质意义研究

出露于南阿尔金西段的塔什萨依辉长岩和闪长岩及其内花岗质细脉是古生代多期岩浆活动的产物,它们构成了南阿尔金俯冲碰撞杂岩带的一部分,是深入探讨该地区构造演化过程的良好材料。本文通过对这些地质体野外地质特征、岩石学、地球化学特征、锆石U-Pb年代学和锆石Hf同位素特征的系统研究,以期确定其形成时代、源区性质和岩石成因,为深入探讨南阿尔金俯冲碰撞杂岩带的演化过程提供进一步约束。本文研究表明:辉长岩和闪长岩具相对较高的SiO_2含量(50.4%~53.5%)和Mg#(57~80),低的K_2O/Na_2O比值(0.17~0.67),为钙碱性系列岩石;辉长岩的∑REE偏低(36.9×10~(-6)~83.1×10~(-6)),Eu正异常(δEu=1.17~1.85);闪长岩∑REE稍高(182×10~(-6)~190×10~(-6)),Eu负异常(δEu=0.85~0.86);二者相对富集LREE和LILE,亏损HFSE。其主量元素特征、微量元素(高Nb、Zr含量,低Zr/Y比值)、锆石Hf同位素和Sr-Nd同位素组成表明二者应为大陆板内伸展背景下软流圈交代的岩石圈地幔部分熔融的产物,并在其岩浆演化过程中经历了地壳混染和分离结晶作用。花岗质细脉具高的Si O2含量(72.4%～75.5%)和K_2O/Na_2O比值(2.33～2.56),属于钾玄岩系列岩石;与闪长岩相比,两者微量元素蛛网图相似,但花岗质细脉轻重稀土分馏程度较大,Eu负异常(δEu=0.55～0.61)和P、Ti的亏损更加显著。锆石U-Pb定年结果显示,辉长岩形成时代约为400～420Ma;闪长岩和其内花岗质细脉的形成时代基本一致,约为400Ma,并具有相似的锆石Hf同位素组成和Sr-Nd同位素组成,结合地球化学特征及野外产状,推测花岗质细脉为闪长岩分异的产物。结合南阿尔金区域地质背景综合分析,塔什萨依辉长岩和闪长岩起源于造山后伸展背景下软流圈交代的岩石圈地幔的部分熔融,并经历了地壳混染和分离结晶作用;而闪长岩进一步分异形成花岗质细脉,并在其形成过程中可能伴随少量壳源熔融物质的加入。     

Dynamic Response to Seismic Waves at the Shallow Slope Supported by Bolts

In the soil slope supported by bolts, longitude waves instead of transverse waves, generated by earthquakes, first reach the slope surface. With the dynamic response of the P (pressure) wave along the anchorage structure, first, a theoretical study was conducted to investigate the propagation characteristics of the interference superposition, generated by the SV (shear-vertical) and the P waves. The SV wave was formed by the wave, originating from the bottom and reflected from the free surface of the slope, whereas the latter was the incident P wave, propagating in the slope. In addition, the structural measures, restraining the seismic wave, and the characteristics of the restraint effect at the free segment of the bolt were investigated. According to the wave-way difference between the incident P wave and the reflected SV wave, the minimum critical slope angle, influenced by the interference at the shallow slope, and the maximum influencing depth of the dynamic response, acting vertically to the slope surface, were obtained. The results indicate that the maximum influencing depth linearly correlated with the slope angle. Furthermore, based on the propagation characteristics of the P wave along the bolt, and the coupled relation between the wave length and the anchorage design parameters, the axial acceleration of the wave propagating along the bolt axis was obtained. Then, the theoretical length of the anti-seismic bolt, subjected to seismic waves, and the compensation force of the anchorage structure were obtained. Finally, a numerical study, based on FlAC3D, properly verified the theoretical conclusions.