Developing a small lake: The case of Palakpakin Lake,San Pablo City,Philippines

This study assesses the development issues of Palakpakin Lake. Premised on the scarcity of development studies on lakes (as the field is dominated by limnology and aquaculture studies) and small-lake studies (as the field is heavily concentrated on big-lake studies), the article argues that the lake has long been wanting of development initiatives from its administrative agencies. It further contends that the Laguna Lake Development Authority and the City Government of San Pablo must now act to promote development of the lake; specifically, by modifying their model-template approach, by formulating a zoning-development plan, and by fostering ecotourism. These basic actions are necessary in managing, developing and conserving the natural resource, as well as in preparing, informing and precipitating other actions on the lake.     

Multi-resolution graph-based clustering analysis for lithofacies identification from well log data: Case study of intraplatform bank gas fields,Amu Darya Basin

In this study, we used the multi-resolution graph-based clustering (MRGC) method for determining the electrofacies (EF) and lithofacies (LF) from well log data obtained from the intraplatform bank gas fields located in the Amu Darya Basin. The MRGC could automatically determine the optimal number of clusters without prior knowledge about the structure or cluster numbers of the analyzed data set and allowed the users to control the level of detail actually needed to define the EF. Based on the LF identification and successful EF calibration using core data, an MRGC EF partition model including five clusters and a quantitative LF interpretation chart were constructed. The EF clusters 1 to 5 were interpreted as lagoon, anhydrite flat, interbank, low-energy bank, and high-energy bank, and the coincidence rate in the cored interval could reach 85%. We concluded that the MRGC could be accurately applied to predict the LF in non-cored but logged wells. Therefore, continuous EF clusters were partitioned and corresponding LF were interpreted, and the distribution and petrophysical characteristics of different LF were analyzed in the framework of sequence stratigraphy.     

Calculation method of prestack FAGVO and its applications

Frequency attenuation occurs when seismic waves propagate through the porous reservoirs containing hydrocarbons. Current researches on the seismic frequency attenuation mainly focus on the post-stack domain instead of the prestack domain. Here we propose the frequency attenuation gradient vs. offset (FAGVO) based on the amplitude variation with offset and frequency attenuation integral equations. We derive the FAGVO equation that equals to zero in a full-elastic medium and is negative in a viscoelastic medium. FAGVO is affected by the viscosity of the medium, the coefficients of reflection, the frequency variation, and high-frequency attenuation. FAGVO uses the differences of partially stacked data to decrease the interference caused by subsurface strata affecting the frequency attenuation, highlights the frequency attenuation gradient anomalies in hydrocarbon-bearing reservoir pores, and finally realizes the hydrocarbon fluid identification. The method was verified using a two-dimensional wave equation forward model and was found to be cost effective. Furthermore, the method does not require well information, which can be applied in the stage of seismic exploration, especially, in the exploration of a none-well project.     

High-Resolution Seismic Characterization of Shallow Gas Accumulations in the Southern Shelf of Marmara Sea,Turkey

High-resolution seismic survey was conducted to investigate acoustic characteristics of gassy sediments along the southern shelf of the Sea of Marmara. The acoustic turbidity zones outlined within the study area are generally below 2–9 m (2-10 ms TWT) the seafloor whilst this vertical distance varies between 9 and 21 m (10–25 ms TWT) for acoustic blanket type reflections. The gassy sediments cover an area of sea floor of about 45, 110, and 75 km² in front of Gönen River, Kocasu River, and Gemlik Bay, respectively. The gassy sediments in the center of Gemlik Bay exhibited an elliptical geometry similar to its basin while the others have deltaic forms in front of the rivers. The sea bottom and near surface sedimentary units are made-up of organic-rich sediments, mostly transported by the southern rivers. The gas observed in sediments is thought to be of biogenic origin, which may be caused by degradation of organic matter in the sediment.     

Comparison of Two Nonstationary Flood Frequency Analysis Methods within the Context of the Variable Regime in the Representative Polish Rivers

Changes in river flow regime resulted in a surge in the number of methods of non-stationary flood frequency analysis. Common assumption is the time-invariant distribution function with time-dependent location and scale parameters while the shape parameters are time-invariant. Here, instead of location and scale parameters of the distribution, the mean and standard deviation are used. We analyse the accuracy of the two methods in respect to estimation of time-dependent first two moments, time-invariant skewness and time-dependent upper quantiles. The method of maximum likelihood (ML) with time covariate is confronted with the Two Stage (TS) one (combining Weighted Least Squares and L-moments techniques). Comparison is made by Monte Carlo simulations. Assuming parent distribution which ensures the asymptotic superiority of ML method, the Generalized Extreme Value distribution with various values of linearly changing in time first two moments, constant skewness, and various time-series lengths are considered. Analysis of results indicates the superiority of TS methods in all analyzed aspects. Moreover, the estimates from TS method are more resistant to probability distribution choice, as demonstrated by Polish rivers’ case studies.     

Pareto Joint Inversion of 2D Magnetotelluric and Gravity Data — Towards Practical Applications

In this paper, a Pareto inversion based global optimization approach, to obtain results of joint inversion of two types of geophysical data sets, is formulated. 2D magnetotelluric and gravity data were used for tests, but presented solution is flexible enough to be used for combination of any kind of two or more target functions, as long as misfits can be calculated and forward problems solved. To minimize dimensionality of the solution, space and introduce straightforward regularization Sharp Boundary Interface (SBI) method was applied. As a main optimization engine, Particle Swarm Optimization (PSO) was used. Synthetic examples based on a real geological model were used to test proposed approach and show its usefulness in practical applications.     

Scaling Analysis of Time Distribution between Successive Earthquakes in Aftershock Sequences

The earthquake inter-event time distribution is studied, using catalogs for different recent aftershock sequences. For aftershock sequences following the Modified Omori’s Formula (MOF) it seems clear that the inter-event distribution is a power law. The parameters of this law are defined and they prove to be higher than the calculated value (2–1/p). Based on the analysis of the catalogs, it is determined that the probability densities of the inter-event time distribution collapse into a single master curve when the data is rescaled with instantaneous intensity, R(t; M _th ), defined by MOF. The curve is approximated by a gamma distribution. The collapse of the data provides a clear view of aftershock-occurrence self-similarity.     

Digitized Database of Old Seismograms Recorder in Romania

The aim of this paper is to describe a managing system for a unique Romanian database of historical seismograms and complementary documentation (metadata) and its dissemination and analysis procedure. For this study, 5188 historical seismograms recorded between 1903 and 1957 by the Romanian seismological observatories (Bucharest-Filaret, Foc?ani, Bac?u, Vrincioaia, Câmpulung-Muscel, Ia?i) were used. In order to reconsider the historical instrumental data, the analog seismograms are converted to digital images and digital waveforms (digitization/ vectorialisation). First, we applied a careful scanning procedure of the seismograms and related material (seismic bulletins, station books, etc.). In a next step, the high resolution scanned seismograms will be processed to obtain the digital/numeric waveforms. We used a Colortrac Smartlf Cx40 scanner which provides images in TIFF or JPG format. For digitization the algorithm Teseo2 developed by the National Institute of Geophysics and Volcanology in Rome (Italy), within the framework of the SISMOS Project, will be used.     

The Influence of a Local Fault Zone on High Energy Tremor Occurrence During Longwall Mining of a Coal Seam

Underground mining of coal seams in the Upper Silesian Coal Basin in Poland is accompanied by seismic activity of varying magnitude. The investigations which have been performed for several years distinguished high energy mine tremors connected directly with mining or coupled with geological structures, such as large faults. In mined seams, local fault zones occur. Faults in these zones are usually small, with throws comparable with coal seams thicknesses. Local fault zone may be responsible for the occurrence of high energy tremors as well as large faults, as presented in this article. An analysis of source mechanism of high energy tremors generated during longwall mining of the coal seam No. 510, with presence of a local fault zone, in one of the Polish hard coal mines in the Upper Silesian Coal Basin was performed. For this purpose, the seismic moment tensor inversion method was used. In most of foci, the process of shear predominated. Determined nodal plane parameters were correlated with parameters of faults forming the local fault zone. High energy tremors were generated mostly by dislocations on faults of the local fault zone. Weakening of roof rocks in the neighborhood of local fault zone takes an important role too, and was responsible for share of implosion in the focal mechanism.     

Tectonic Plate Parameters Estimated in the International Terrestrial Reference Frame ITRF2008 Based on SLR Stations

This paper concerns an analysis of the accuracy of estimated parameters Ω (Φ, Λ, ω) which define the tectonic plate motions. The study is based on the velocities of station positions published by ITRF2008 for Satellite Laser Ranging (SLR) technique. The Eurasian, African, North American and Australian plates were used in the analysis. Influence of the number and location of stations on the plate surface on estimation accuracy of the tectonic plate motion parameters was discussed. The results were compared with the APKIM 2005 IGN model. In general, a remarkable concurrence agreement between our solutions and the APKIM 2005 model was found.