首都圈地区现今断层活动的数学模拟及地震危险性讨论

以首都圈地区现今活动断层上近２０年的位移测量资料为依据，用二维线弹性有限元对该区断层的活动特征进行了拟合，并结合有关资料讨论了该区的地震危险性。结果表明：１９７７－１９８６年间该区主压应力优势方位与华北较一致，约为Ｎ４５°－８０°Ｅ；１９８６－１９９０年其主压应力方向向北偏转，约为Ｎ５°－６０°Ｅ；１９９０年至现在其主压应力又向东偏转，角度大于第一阶段，约为Ｎ８０°－９５°Ｅ。该区的张家口－延庆一带近期有发生中强震的可能；丰镇－阳高－大同地区和凉城－古营盘地区应力较高，也应引起注意。     

Directions and possibilities of mathematical geology

This paper considers the present state of mathematical geology. Three directions are recognized: applied, theoretical, and mathematical. Applied mathematical geology includes formal use of mathematics to solve problems and computer processing of data. Success is achieved by a correspondence of mathematical methods used to the nature of geological data. This correspondence can be demonstrated by purely mathematical means. Theoretical mathematical geology uses mathematics as a language of geology; however, a number of methodological problems must be solved: formalization of initial geological concepts and creation of a strict conceptual basis, substantiation of initial principles of mathematical simulation, creation of theoretical geological models, problems of elementary and coincidence in geology, and methodological substantiations of possibilities of any mathematical model to approximate geological models. The essense and significance of these problems are considered. The main task of mathematical geology is to prove its correspondence to the nature of the geological objects studied, geological data obtained, and geological problems solvable. Finally, the main problems of mathematical geology are not so much mathematical as geological and methodological.     

岩体三维节点偶对分析及其应用

本文讨论了节点偶对分析的三维扩展。节点偶对分析是针对摩擦滑动节理单元的约束边界进行的序列矢量判定。在三维节理单元的分析中,几何约束和相应的力学判据极为复杂,采用节点偶对等效结合力方法才能迎刃而解。文中对三维转向节理模型的分析得到合理的结果,它表明本文所提出的方法可成功地应用于节理岩体的有限元分析。     

地壳演化与成矿作用——以川滇地洼系“四层楼”铜矿床序列为例

提出川滇地洼系“四层楼”铜矿床序列的形成与陆壳演化的成生联系,是与本区陆壳由前地槽—地槽—地台—地洼演化各阶段与之相匹配的成矿作用的产物.与此同时,并总结了本区“四层楼”铜矿床序列的成矿作用具有明显的继承性、新生性、旋回性及层控性四大特点和多因复成矿床的成矿模式.     

Alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands)

Ancient fluvial successions often act as hydrocarbon reservoirs. Sub‐surface data on the alluvial architecture of fluvial successions are often incomplete and modelling is performed to reconstruct the stratigraphy. However, all alluvial architecture models suffer from the scarcity of field data to test and calibrate them. The purposes of this study were to quantify the alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands) and to determine spatio‐temporal trends in the architecture. Five north–south orientated cross‐sections, perpendicular to the general flow direction, were compiled for the fluvial‐dominated part of the delta. These sections were used to calculate the width/thickness ratios of fluvial sandbodies (SBW/SBT) and the proportions of channel‐belt deposits (CDP), clastic overbank deposits (ODP) and organic material (OP) in the succession. Furthermore, the connectedness ratio (CR) between channel belts was calculated for each cross‐section. Distinct spatial and temporal trends in the alluvial architecture were found. SBW/SBT ratios decrease by a factor of ca 4 in a downstream direction. CDP decreases from ca 0·7 (upstream) to ca 0·3 (downstream). OP increases from less than 0·05 in the upstream part of the delta to more than 0·25 in the downstream delta. ODP is approximately constant (0·4). CR is ca 0·25 upstream, which is approximately two times larger than in the downstream part of the delta. Furthermore, CDP in the downstream Rhine–Meuse delta increases after 3000 cal yr BP. These trends are attributed to variations in available accommodation space, floodplain geometry and channel‐belt size. For instance, channel belts tend to narrow in a downstream direction, which reduces SBW/SBT, CDP and CR. Tectonics cause local deviations in the general architectural trends. In addition, the positive correlation between avulsion frequency and the ratio of local to regional aggradation rate probably influenced alluvial architecture in the Rhine–Meuse delta. The Rhine–Meuse data set can be a great resource when developing more sophisticated models for alluvial architecture simulation, which eventually could lead to better characterizations of hydrocarbon reservoirs. To aid such usage of the Rhine–Meuse data set, constraints for relevant parameters are provided at the end of the paper.     

The Lower Permian Wasp Head Formation,Sydney Basin: high‐latitude,shallow marine sedimentation following the late Asselian to early Sakmarian glacial event in eastern Australia

The Lower Permian Wasp Head Formation (early to middle Sakmarian) is a ～95 m thick unit that was deposited during the transition to a non‐glacial period following the late Asselian to early Sakmarian glacial event in eastern Australia. This shallow marine, sandstone‐dominated unit can be subdivided into six facies associations. (i) The marine sediment gravity flow facies association consists of breccias and conglomerates deposited in upper shoreface water depths. (ii) Upper shoreface deposits consist of cross‐stratified, conglomeratic sandstones with an impoverished expression of the Skolithos Ichnofacies. (iii) Middle shoreface deposits consist of hummocky cross‐stratified sandstones with a trace fossil assemblage that represents the Skolithos Ichnofacies. (iv) Lower shoreface deposits are similar to middle shoreface deposits, but contain more pervasive bioturbation and a distal expression of the Skolithos Ichnofacies to a proximal expression of the Cruziana Ichnofacies. (v) Delta‐influenced, lower shoreface‐offshore transition deposits are distinguished by sparsely bioturbated carbonaceous mudstone drapes within a variety of shoreface and offshore deposits. Trace fossil assemblages represent distal expressions of the Skolithos Ichnofacies to stressed, proximal expressions of the Cruziana Ichnofacies. Impoverished trace fossil assemblages record variable and episodic environmental stresses possibly caused by fluctuations in sedimentation rates, substrate consistencies, salinity, oxygen levels, turbidity and other physio‐chemical stresses characteristic of deltaic conditions. (vi) The offshore transition‐offshore facies association consists of mudstone and admixed sandstone and mudstone with pervasive bioturbation and an archetypal to distal expression of the Cruziana Ichnofacies. The lowermost ～50 m of the formation consists of a single deepening upward cycle formed as the basin transitioned from glacioisostatic rebound following the Asselian to early Sakmarian glacial to a regime dominated by regional extensional subsidence without significant glacial influence. The upper ～45 m of the formation can be subdivided into three shallowing upward cycles (parasequences) that formed in the aftermath of rapid, possibly glacioeustatic, rises in relative sea‐level or due to autocyclic progradation patterns. The shift to a parasequence‐dominated architecture and progressive decrease in ice‐rafted debris upwards through the succession records the release from glacioisostatic rebound and amelioration of climate that accompanied the transition to broadly non‐glacial conditions.     

Sediment transport in analogue flume models compared with real‐world sedimentary systems: a new look at scaling evolution of sedimentary systems in a flume

Evolution of sedimentary systems at large temporal and spatial scales cannot be scaled down to laboratory dimensions by conventional hydraulic Froude scaling. Therefore, many researchers question the validity of experiments aiming to simulate this evolution. Yet, it has been shown that laboratory experiments yield stratigraphic responses to allocyclic forcing that are remarkably similar to those in real‐world prototypes, hinting at scale independency with strong dependence on boundary conditions but weak dependence on the actual sediment transport dynamics. This paper addresses the dilemma by contrasting sediment transport rules that apply in the laboratory with those that apply in real‐world geological systems. It is demonstrated that the generation of two‐dimensional stratigraphy in a flume can be simulated numerically by the non‐linear diffusion equation. Sediment transport theory is used to demonstrate that only suspension‐dominated meandering rivers should be simulated with linear diffusion. With increasing grain‐size (coarse sand to gravel) and shallowness of river systems, the prediction of long‐term transport must be simulated by non‐linear, slope‐dependent diffusion to allow for increasing transport rates and thus change in stratigraphic style. To point out these differences in stratigraphic style, three stages in infill of accommodation have been defined here: (i) a start‐up stage, when the system is prograding to base level (e.g. the shelf edge) with no sediment flux beyond the base‐level point; (ii) a fill‐up stage, when the system is further aggrading while progressively more sediment is bypassing base level with the progression of the infill; and (iii) a keep‐up stage, when more than 90% of the input is bypassing the base level and less than 10% is used for filling the accommodation. By plotting the rate of change in flux for various degrees of non‐linearity (varying the exponent in the diffusion equation) it was found that the error between model and real‐world prototype is largest for the suspension‐dominated prototypes, although never more than 30% and only at the beginning of the fill‐up stage. The error reduces to only 10% for the non‐linear sandy‐gravelly and gravelly systems. These results are very encouraging and open up ways to calibrate numerical models of sedimentary system evolution by such experiments.     

Discriminating between pore‐filling load and bed‐structure load: a new porosity‐based method,exemplified for the river Rhine

Sediments contained in the river bed do not necessarily contribute to morphological change. The finest part of the sediment mixture often fills the pores between the larger grains and can be removed without causing a drop in bed level. The discrimination between pore‐filling load and bed‐structure load, therefore, is of practical importance for morphological predictions. In this study, a new method is proposed to estimate the cut‐off grain size that forms the boundary between pore‐filling load and bed‐structure load. The method evaluates the pore structure of the river bed geometrically. Only detailed grain‐size distributions of the river bed are required as input to the method. A preliminary validation shows that the calculated porosity and cut‐off size values agree well with experimental data. Application of the new cut‐off size method to the river Rhine demonstrates that the estimated cut‐off size decreases in a downstream direction from about 2 to 0·05 mm, covariant with the downstream fining of bed sediments. Grain size fractions that are pore‐filling load in the upstream part of the river thus gradually become bed‐structure load in the downstream part. The estimated (mass) percentage of pore‐filling load in the river bed ranges from 0% in areas with a unimodal river bed, to about 22% in reaches with a bimodal sand‐gravel bed. The estimated bed porosity varies between 0·15 and 0·35, which is considerably less than the often‐used standard value of 0·40. The predicted cut‐off size between pore‐filling load and bed‐structure load (D_c,p) is fundamentally different from the cut‐off size between wash‐load and bed‐material load (D_c,w), irrespective of the method used to determine D_c,p or D_c,w. D_c,w values are in the order of 10^?1 mm and mainly dependent on the flow characteristics, whereas D_c,p values are generally much larger (about 10⁰ mm in gravel‐bed rivers) and dependent on the bed composition. Knowledge of D_c,w is important for the prediction of the total sediment transport in a river (including suspended fines that do not interact with the bed), whereas knowledge of D_c,p helps to improve morphological predictions, especially if spatial variations in D_c,p are taken into account. An alternative to using a spatially variable value of D_c,p in morphological models is to use a spatially variable bed porosity, which can also be predicted with the new method. In addition to the morphological benefits, the new method also has sedimentological applications. The possibility to determine quickly whether a sediment mixture is clast‐supported or matrix‐supported may help to better understand downstream fining trends, sediment entrainment thresholds and variations in hydraulic conductivity.     

Sedimentary evolution of a Late Pleistocene temperate red algal reef (Coralligène) on Rhodes,Greece: correlation with global sea‐level fluctuations

Autochthonous red algal structures known as coralligène de plateau occur in the modern warm‐temperate Mediterranean Sea at water depths from 20 to 120 m, but fossil counterparts are not so well‐known. This study describes, from an uplifted coastal section at Plimiri on the island of Rhodes, a 450 m long by 10 m thick Late Pleistocene red algal reef (Coralligène Facies), interpreted as being a coralligène de plateau, and its associated deposits. The Coralligène Facies, constructed mainly by Lithophyllum and Titanoderma, sits unconformably upon the Plio‐Pleistocene Rhodes Formation and is overlain by a Maerl Facies (2 m), a Mixed Siliciclastic‐Carbonate Facies (0·2 m) and an Aeolian Sand Facies (2·5 m). The three calcareous facies, of Heterozoan character, are correlated with established members in the Lindos Acropolis Formation in the north of the island, while the aeolian facies is assigned to the new Plimiri Aeolianite Formation. The palaeoenvironmental and genetic‐stratigraphic interpretations of these mixed siliciclastic‐carbonate temperate water deposits involved consideration of certain characteristics associated with siliciclastic shelf and tropical carbonate shelf models, such as vertical grain‐size trends and the stratigraphic position of zooxanthellate coral growths. Integration of these results with electron spin resonance dates of bivalve shells indicates that the Coralligène Facies was deposited during Marine Isotope Stage 6 to 5e transgressive event (ca 135 to 120 ka), in water depths of 20 to 50 m, and the overlying Maerl Facies was deposited during regression from Marine Isotope Stage 5e to 5d (ca 120 to 110 ka), at water depths of 25 to 40 m. The capping Aeolian Sand Facies, involving dual terrestrial subunits, is interpreted as having formed during each of the glacial intervals Marine Isotope Stages 4 (71 to 59 ka) and 2 (24 to 12 ka), with soil formation during the subsequent interglacial periods of Marine Isotope Stages 3 and 1, respectively. Accumulation rates of about 0·7 mm year^?1 are estimated for the Coralligène Facies and minimum accumulation rates of 0·2 mm year^?1 are estimated for the Maerl Facies. The existence of older red algal reefs in the Plimiri region during at least Marine Isotope Stages 7 (245 to 186 ka) and 9 (339 to 303 ka) is inferred from the occurrence of reworked coralligène‐type lithoclasts in the basal part of the section and from the electron spin resonance ages of transported bivalve shells.