Hybrid seismic modeling based on discrete-wave number and finite-difference methods

Any calculation of seismic wave propagation comprising the seismic source, the travel path, and the receiver site in a single finite-difference (FD) model requires a considerable amount of computer time and memory. Moreover, the methods currently available for including point sources in the 2D FD calculations are far-field approximations only. Therefore we have developed a new hybrid method for treating the seismic wave fields at localized 2D near-surface structures embedded in a 1D background medium, and excited by a point source. The source radiation and propagation in the background model is solved by the discrete-wave number (DW) method, while the propagation in the local 2D structure is calculated by the FD method. The coupling between the two sets of calculations is performed on a rectangular excitation box surrounding the local structure. We show the usefulness of the method in ground-motion studies where both near-field source effects and local site effects are important. Technical problems connected with the inconsistency between the 3D source radiation and the 2D FD calculation are minor for the relatively distant in-plane point explosive sources, but are more serious for the in-plane dislocation sources.     

不均匀植被分布对地表面和大气边界层影响的数值试验

研究陆地与大气间相互作用的方法之一是建立联系地表面层与大气间各种过程的数值模式进行模拟。本文是建立一个陆面过程与二维大气边界层相耦合的模式,耦合模式中包含了发生在大气边界层、植被冠层和土壤表层各种动力、热力和水文过程。运用这一模式模拟了荒漠环境中一片绿洲的不均匀地表面形成的局地气候。由于绿洲植被与周围荒漠有着显著不同的水份与能量平衡关系,使绿洲表面与边界层较四周荒漠冷而湿,并形成了相应的局地环流,即所谓“绿洲效应”。试验结果表明,模拟的气候状况与观测现象是一致的。模式可以用于陆气相互作用的研究。     

Pore pressure profiles in fractured and compliant rocks1

Fluid permeability in fractured rocks is sensitive to pore-pressure changes. This dependence can have large effects on the flow of fluids through rocks. We define the permeability compliance γ= 1/k(k/δp_p)_pc, which is the sensitivity of the permeability k to the pore pressure p_p at a constant confining pressure p_c, and solve the specific problems of constant pressure at the boundary of a half-space, a cylindrical cavity and a spherical cavity. The results show that when the magnitude of permeability compliance is large relative to other compliances, diffusion is masked by a piston-like pressure profile. We expect this phenomenon to occur in highly fractured and compliant rock systems where γ may be large. The pressure profile moves rapidly when fluids are pumped into the rock and very slowly when fluids are pumped out. Consequently, fluid pressure, its history and distribution around injection and production wells may be significantly different from pressures predicted by the linear diffusion equation. The propagation speed of the pressure profile, marked by the point where δp_p/δx is a maximum, decreases with time approximately as and the amplitude of the profile also dissipates with time (or distance). The effect of permeability compliance can be important for fluid injection into and withdrawal from reservoirs. For example, excessive drawdown could cause near-wellbore flow suffocation. Also, estimates of the storage capacity of reservoirs may be greatly modified when γ is large. The large near-wellbore pressure gradients caused during withdrawal by large γ can cause sanding and wellbore collapse due to excessive production rates.     

Dynamics of the plasmasphere in the dusk sector as observed in the outer ionosphere

Summary Direct measurements of the thermal plasma parameters in the topside ionosphere reveal variations of the plasmasphere boundary in the dusk sector. The ACTIVE satellite's near-polar orbits at altitudes of 500 – 1800 km around winter solstice 1989 were used to study the bulge region of the plasmasphere during intervals with different levels of geomagnetic agitation. The narrow, sharply defined trough in electron concentration corresponding to the plasmapause under quiet conditions situated at L = 6 – 7 moved to lower L-values with increasing geomagnetic activity. This narrow trough can be found in all main ion constituents. During periods of moderate geomagnetic activity, following the onset of a weak magnetic storm, a portion of the plasmaspheric bulge region was separated from the main plasmaspheric body. This can be seen in the outer ionosphere as an inner narrow trough at lower L-value. Troughs in light ions need no longer coincide with this in electron concentration. He⁺ is the most sensitive constituent reflecting the dusk sector plasmaspheric situation at this altitude.Dedicated to the Memory of Professor Karel P     

Densities of Huanghe River estuary and Bohai Bay waters

The densities of 36 water samples from the Huanghe River estuary and Bohai Bay were determinedby a magnetic float densimcter under three temperatures from 15℃ to 25℃.All the measured densities ofsamples were greater than that of the values calculated from the International Equation of State of Seawater.The differences between the measured and calculated densities increased with the decrease of salinities.The dif-ferences appeared exponentially correlated with[Ca~(2+)]/s,[Mg~(2+)]/s and[SO_4~(2-)]/s,and had"s"type curverelationship with the alkalinity in all salinity range.But in the salinity ranging from 25.72 to 31.57,therelationships were all linear.The density difference can be estimated from the equation △ρ（10~3kg·m~(-3)）=（-2.79+236.5([Ca~(2+)]/s)/(-9.7464×10~（-3）+[Ca~(2+)]/s).It was the high alkalinity and[Ca~(2+)]/s that resulted in the measured densi-ties of seawaters being higher than the calculated densities in the Huanghe estuary and Bohai Bay.     

Magnetic properties,self-reversal remanence and thermal alteration products of smythite

Summary This study is a follow up of the investigation of some magnetic properties and metastability of greigite in samples obtained from Miocene claystones in the Kruné hory (Erzgebirge) Piedmont basins (Bohemia). Three different methods of upgrading the smythite were applied; the magnetic properties of the concentrates are compared. The thermal conversion of smythite sets in at 200°C while greigite converts at 250°C. The first intermediate products to be formed are iron sulphides, marcasite clearly dominating over pyrite and pyrrhotite. Apart from a Fe³⁺ sulphate with a composition of Fe₂(SO₄)₃, oxidation of these sulphides gives rise to -Fe₂O₃. The result of the subsequent decomposition of the mentioned sulphate is the formation of -Fe₂O₃, which retains the sulphate structure. The final product of the thermal decomposition at 800°C is -Fe₂O₃. In the smythite concentrate the conversion to Fe³⁺ sulphate and -Fe₂O₃ is about twice as intensive as in greigite. No direct conversion to -Fe₂O₃ was found. During the thermal process self-reversals of remanence were observed, in various samples as many as four reversals in the temperature interval from 340 to 590°C. The occurrences of self-reversals of remanence were only observed at high degrees of thermal demagnetization, of the order of 10^–2 down to 10^–3 in the temperature interval of sulphide origin (below 400°C), and of the order of 10^–4 down to 10^–6 in the temperature interval of Fe-oxides origin (above 400°C).Presented at the 3rd Conference on New Trends in Geomagnetism, Castle of Smolenice, Czechoslovakia, June 22–29, 1992     

Die Rheingrabenrandverwerfung im Wattkopftunnel bei Ettlingen (Nordschwarzwald)

The Wattkopftunnel, near Ettlingen (Nordschwarzwald), drives through the eastern margin of the Rheingraben. The tunnel passes cenozoic and mesozoic sediments. Early quarternary and tertiary beds are situated west of the main thrust of the Rheingraben. Fossil record indicates upper Oligocene age (Chatt) for parts of the tertiary sediments. At the eastern border of the Rheingraben, wedges of jurassic and middle triassic series are squeezed. East of the Rheingraben the tunnel drives in the lower triassic Bausandstein. The eastern margin of the Rheingraben was investigated in detail during tunneling. Faults of the Rheingraben margin are distributed in an 130 meter wide fault zone in the tunnel. Total stratigraphic separation by the normal faults reaches more than 2 000 meters. The cenozoic sequence suffered synsedimentary to early diagenetic deformation, while the mesozoic series are characterized by ruptural deformation. The fault- and joint system is directed in the rheinische Richtung (SSW-NNE). East of the Rheingraben a second direction occur, running parallel to the lower Albtal (W-E).

     

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

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

The effect of an optically thin region on stability of an accretion disk around a supermassive black hole

We analyse conditions of the innermost portion of an accretion disk and establish a set of equations for this region. A stable innermost region may exist, which can probably explain the observed UV and X-ray spectra, avoiding the unstable emission. We then discuss the detailed radial structure of a disk around a black hole for typical AGN parameters and obtain different kinds of- relationships for different regions of a disk. On the basis of this, we discuss the stability. A new type of cycle is present, which we call a double S shaped cycle. In this cycle, the extent of accretion rate variability is much larger than that in dwarf nova cycles. This probably solves the problem of violent variability of AGN. In the meantime, the very high accretion rate at the hottest state in limit cycles in the unstable region may provide continuous injection of matter to the jet and power the relativistic motion of the jet.     

毫秒时间尺度太阳微波辐射新特征

自1990年6月北京天文台太阳射电偏振仪投入实测以来,观测到一百多个“Spike”和短时标精细结构微波幅射,其中包括一些新的微波幅射特征。本文给出了主要的四点特征:(1)窄带辐射;(2)快速偏振逆转;(3)快速频率漂移;(4)不同时标的准周期振荡。