聚合物驱油浓度与粘度的研究

本文对聚全物驱油试验获得聚合物溶液的方式和聚合物浓度、粘度的测定进行了研究和比较，得到了辉锐公司自动混配装置中的可编程控制器设置浓度的修正值。此研究对控制聚合物浓度稳定、准确地注入和聚合物驱油效果的评价具有实际价值。     

Analytical models for bubble growth during decompression of high viscosity magmas

Analytical models for decompressional bubble growth in a viscous magma are developed to establish the influence of high magma viscosity on vesiculation and to assess the time-scales on which bubbles respond to decompression. Instantaneous decompression of individual bubbles, analogous to a sudden release of pressure (e.g. sector collapse), is considered for two end-member cases. The infinite melt model considers the growth of an isolated bubble before significant bubble interaction occurs. The shell model considers the growth of a bubble surrounded by a thin shell and is analogous to bubble growth in a highly vesicular magmatic foam. Results from the shell model show that magmas less viscous than 10⁹ Pa s can freely expand without developing strong overpressures. The timescales for pressure re-equilibration are shortened by increased ratios of bubble radius to shell thickness and by larger decompression. Time-scales for isolated bubbles in rhyolitic melts (infinite melt model) are significantly longer, implying that such bubbles could experience internal pressures greater than the ambient pressure for at least a few hours following a sudden release of pressure. The shell model is developed to assess bubble growth during the linear decompression of a magma body of constant viscosity. For the range of decompression rates and viscosities associated with actual volcanic eruptions, bubble growth continues at approximately the equilibrium rate, with no attendant excess of internal pressure. The results imply that viscosity does not have any significant role in preventing the explosive expansion of high viscosity foams. However, for viscosities of >10⁹ Pa s there is the potential for a viscosity quench under the extreme decompression rates of an explosive eruption. It is proposed that the typical vesicularities of pumice of 0.7–0.8 are a consequence of the viscosity of the degassing magmas becoming sufficiently high to inhibit bubble expansion over the characteristic time-scale of eruption. For fully degassed silicic lavas with viscosities in the range 10¹⁰ to 10¹² Pa s time-scales for decompression of isolated bubbles can be hours to many months.     

Stratigraphic architecture of sedimentary basin induced by mantle diapiric upwelling and eustatic event

The stratigraphic architecture of sedimentary basin provides important constraints on the rheological structure of the upper mantle and the crust, eustatic events and tectonic movements. In this study, we examined the convective coupling between the uppermost mantle and ductile lower crust as a formation mechanism of sedimentary basins. In this mechanism, the lower crust is squeezed by the upper mantle flow associated with mantle diapiric upwelling, resulting in the surface subsidence and formation of sedimentary basin. We investigated the stratigraphic architecture by taking into account the horizontal migration of the upwelling in time and spatial domains. The sedimentary basin is characterized by both the thickness and formation period for its gross feature and the sequence stratigraphy. The decay of the upwelling causes the surface uplift, and some parts of sediments deposited during the surface subsidence are consequently eroded. The subsidence area moves as the activity of upwelling horizontally migrates, resulting in the formation of unconformity for the uplifted and eroded area over the previous upwelling. We also incorporated the effects of third order eustasy, with amplitude of 100 m and period of 1 Myr, into convective coupling model. An application of our model to Karatsu-Sasebo coalfield in the Tertiary of the northwest Kyushu, Japan, indicates that the stratigraphic architecture of sedimentary basin including two effects, i.e. convective coupling and eustasy, may provide important information about the viscosity structure of the lower crust and uppermost mantle and spatio-temporal growth and decay histories of the mantle diapiric upwelling.     

Deformation mechanism maps for feldspar rocks

Deformation mechanism maps for feldspar rocks were constructed based on recently published constitutive laws for dislocation and grain boundary diffusion creep of wet and dry plagioclase aggregates. The maps display constant temperature contours in stress-grain size space for strain rates ranging from 10⁻¹⁶ to 10⁻¹² s⁻¹.Two fields of dominance of grain boundary diffusion-controlled creep and dislocation creep are separated by a strongly grain size-sensitive transition zone. For wet rocks, diffusion-controlled creep dominates below a grain size of about 0.1–1 mm, depending on temperature, stress, strain rate and feldspar composition. Plagioclase aggregates containing up to 0.3 wt.% water as often found in natural feldspars are more than 2 orders of magnitude weaker than dry rocks. The strength of water-bearing feldspar rocks is moderately dependent on composition and water fugacity.For a grain size range of about 10–50 μm commonly observed in natural ultramylonites, the deformation maps predict that diffusion-controlled creep is dominant at greenschist to granulite facies conditions. Low viscosity estimates of 10¹⁸–10¹⁹ Pa·s from modeling postseismic stress relaxation and channel flow of the continental lower crust can only be reconciled with laboratory experiments assuming dislocation creep at high temperatures >900 °C or, at lower temperatures, diffusion creep of fine-grained rocks possibly localized in abundant high strain shear zones. For similar thermodynamic conditions and grain size, lower crustal rocks are predicted to be less than order of magnitude weaker than upper mantle rocks.     

Viscosity measurements of hydrous rhyolitic melts using the fiber elongation method

We have devised a new, simple and easy technique to measure the viscosity of hydrous silicate melts by combining an autoclave for melt hydration and the fiber elongation method for viscosity measurement. Using this, we measured the viscosity of hydrous rhyolitic melts whose water content ranges from 0.02 to 0.58 wt%. We observed a drastic decrease in viscosity against water content: 0.1 wt% water decreases the viscosity about an order of magnitude. Even when the water content is only 0.02 wt%, the viscosity decreased about half an order of magnitude. These results clearly demonstrate that the effect of water on viscosity should not be ignored even when it occurs as a trace constituent. We compared our experimental data with those derived from a non-Arrhenian viscosity model, which is considered to be applicable to calc-alkaline samples. This model succeeded in expressing the viscosity variation against water content but was unable to accurately predict the measured viscosity of liquids.Editorial responsibility: D. Dingwell     

Elastic and inelastic triggering of earthquakes in the North Anatolian Fault zone

Deformation models used to explain the triggering mechanism often assume pure elastic behaviour for the crust and upper mantle. In reality however, the mantle and possibly the lower crust behave viscoelastically, particularly over longer time scales. Consequently, the stress field of an earthquake is in general time-dependent. In addition, if the elastic stress increase were enough to trigger a later earthquake, this triggered event should occur instantaneously and not many years after the triggering event. Hence, it is appropriate to include inelastic behaviour when analysing stress transfer and earthquake interaction.In this work, we analyse a sequence of 10 magnitude M_s > 6.5 events along the North Anatolian Fault between 1939 and 1999 to study the evolution of the regional Coulomb stress field. We investigate the triggering of these events by stress transfer, taking viscoelastic relaxation into account. We evaluate the contribution of elastic stress changes, of post-seismic viscoelastic relaxation in the lower crust and mantle, and of steady tectonic loading to the total Coulomb stress field. We analyse the evolution of stress in the region under study, as well as on the rupture surfaces of the considered events and their epicentres. We study the state of the Coulomb stress field before the 1999 İzmit and Düzce earthquakes, as well as in the Marmara Sea region.In general, the Coulomb stress failure criterion offers a plausible explanation for the location of these events. However, we show that using a purely elastic model disregards an important part of the actual stress increase/decrease. In several cases, post-seismic relaxation effects are important and greater in magnitude than the stress changes due to steady tectonic loading. Consequently, viscoelastic relaxation should be considered in any study dealing with Coulomb stress changes.According to our study, and assuming that an important part of the rupture surface must be stressed for an earthquake to occur, the most likely value for the viscosity of the lower crust or mantle in this region is 5 · 10¹⁷–10¹⁸ Pa · s. Our results cannot rule out the possibility of other time-dependent processes involved in the triggering of the 1999 Düzce event. However, the stress increase due to viscoelastic relaxation brought 22% of the 1999 Düzce rupture area over the threshold value of Δσ_c ≥ 0.01 MPa (0.1 bar), and took the whole surface closer to failure by an average of 0.2 MPa. Finally, we argue that the Marmara Sea region is currently being loaded with positive Coulomb stresses at a much faster rate than would arise exclusively from steady tectonic loading on the North Anatolian Fault.     

岩石古流变性质的构造研究进展

岩石古流变性质研究是大陆动力学研究的一个重要方面。岩石古流变性质研究的一个发展方向是通过比例模型,建立流变参数与构造变形几何参数之间的定量关系,再利用天然岩石构造变形几何参数的定量测定,反推在形成这些变形时所处物理化学条件下的岩石流变参数。对利用露头尺度地质构造的定量测量研究岩石古流变性质的方法与进展做了重点介绍。     

地壳、上地幔变形属性的判别

基于Ｍａｘｗｅｌｌ和Ｋｅｌｖｉｎ流变模型，以延续时间与松弛期的比值，即松弛数，作为衡量指标，对地壳、上地幔不同层次的变形属性进行了判别。结果表明，除软流圈地幔主要为粘性流动外，其上各层的情况是：在漫长的地质过程中，不仅岩石圈地幔和下地壳属于粘性，上地壳也有可能处于粘性流动状态；在延续数十年的情况下，上地壳以及岩石圈地幔都可能呈弹性状态；延续时间介乎二者，通常上地壳为弹性，下地壳和岩石圈地幔为粘性或滞弹粘性。此外，还结合变形属性的判别，就岩石圈地幔的流动和板内的驱动力远程传递等问题进行了讨论     

Late Pleistocene crustal uplift and gravity anomaly in the eastern part of Kyushu, Japan, and its geophysical implications

The Miyazaki Plain, eastern part of Kyushu, Japan, is characterized by both significant negative gravity anomalies and aseismic crustal uplifting (1 mm/year) in the Late Pleistocene and Holocene. We examine the relationship between these two phenomena, which may provide important constraints on the interaction between the collision and/or subduction of the Kyushu-Palau Ridge and the forearc. We estimate the mass deficiency below 11-km depth by using the gravity anomalies and P-wave velocity structure of the upper crust. The onset of the load accumulation, 0.5–0.4 Ma, is inferred from the movement of the fluvial terraces considering the tephrochronology. The loading history is assumed to be a linear function of time. We evaluate the crustal rebound by assuming a viscoelastic plate deformation with an underplating load existing at 20- or 30-km depth. The predicted crustal movement for models with a lithospheric (crustal) viscosity of 10²³–10²⁴ Pa s can explain the observed altitudes of the shoreline of the marine terraces formed at the Last Interglacial of about 125 kyr BP and the middle Holocene of 5–6 kyr BP. Although we cannot restrict the origin of the buoyant body, the subduction of the Kyushu-Palau Ridge, remnant arc associated with back-arc opening of the Shikoku Basin, may be related to the buoyancy for the uplifting region examined here. On the other hand, the buoyant body off the Miyazaki Plain probably plays an important role in the interaction between the subducting oceanic slab and the overriding forearc crust. Thus, the observed lateral variation of the interplate coupling on the convergent boundary along the Nankai Trough may be attributed to the existence of the buoyant body.     

震后GPS观测数据揭示的日本M_W9.0地震周边地区地幔黏滞性结构垂向变化

本文利用大范围的震后GPS数据和黏弹性球形地球位错理论,定量研究了日本M_W9.0地震周边地区地幔黏滞性结构的垂向变化.首先结合陆地和海底的GPS观测数据,以及基于球形地球位错理论格林函数和贝叶斯反演方法,反演了该地震的同震滑动分布,发现其最大错动量高达59m.然后在均一地幔黏滞性结构的假设前提下,确定了震源周边地区地幔黏滞因子的最优解,发现依据该地幔黏滞因子获得的理论远场震后位移和GPS观测结果之间的均方根误差高达0.81cm,不能解释远场观测结果.为解决上述问题,本文对震中周边地区地幔黏滞性结构沿垂向方向进行分层,建立了一个随深度变化的地幔黏滞性构造模型,然后综合利用远近场的GPS数据对该地区地幔黏滞因子进行反演研究,结果表明,震源周边地区岩石圈弹性层厚度最优解为40km,40～220km深度的地幔黏滞因子最优解为6×10~(18)Pa·s,220～670km深度之间的地幔黏滞因子最优解为1.5×10~(19)Pa·s.上述地幔黏滞性构造使远场的均方根误差降为0.12cm,仅为利用均一地幔黏滞性构造所得均方根误差值的15%,大大提高了远场模拟结果的准确性.最后,观测值和模拟值之间的均方根误差分析表明,近场震后形变数据主要约束浅层的地幔黏滞性结构,而远场震后形变数据主要约束深部的地幔黏滞性结构.