Optimization of grid-drilling using computer simulation

A computer simulation method has been developed to find efficient drilling grids for mineral deposits. A well-known ore deposit is used as a model to develop an efficient pattern for undiscovered ore bodies in the same area or in other prospects where similar geometry is suspected. The model for this study is the Austinville, Virginia deposit, a Mississippi Valley-type deposit composed of 17 ore bodies totaling 34 million short tons (30 million metric tons). The method employs a computer program that simulates drilling the model deposit with different patterns, including various levels of follow-up drilling. Follow-up holes are drilled in fences at one half the original spacing around holes in the grid that show ore-grade mineralization. Each pattern is drilled 100 times from random starting locations to provide a range of outcomes of drilling, including the best, worst, and most likely. For this study, patterns of 100 drill holes were composed of 10 fences spaced 1000–5000 feet (305–1524 m) apart, each with 10 holes spaced 200–1000 feet (61–305 m) apart. In all, 25 grids were used with zero to three levels of follow-up drilling. The 600/2000 grid, with drill holes spaced 600 feet (183 m) apart in fences spaced 2000 feet (610 m) apart, was compared with the 200/5000 grid because they represented contrasting outcomes. The 600/2000 grid penetrated many ore bodies consistently but with few multiple hits to individual ore bodies; whereas the 200/5000 grid inconsistently penetrated few ore bodies with many multiple hits. The 600/2000 grid was more efficient than the 200/5000 grid at hitting large ore bodies of 1,000,000 short tons or greater (900,000 metric tons or greater) and was made more effective by adding one cycle of follow-up drilling. The 600/2000 grid had a 97% chance of hitting one or more large ore bodies with at least one drill hole per ore body, and the 200/5000 grid had a 64% chance. Once hit, there was an 82% chance that the largest ore body would be penetrated by three or more holes when using the 600/2000 grid and an 88% chance using the 200/5000 grid.     

δ34S isotope values of dissolved sulfate (SO4 2−) as a tracer for battery acid (H2SO4) contamination in groundwater

The application of sulfur isotope (³⁴S) values of sulfate in groundwater provided the information necessary to evaluate the source, transport and fate of battery acid and associated contaminants at the Gulf Coast Recycling (GCR) facility. The chemical and isotopic composition of groundwater beneath the (GCR) property, a battery recycling facility in east Tampa, Florida, varies more than expected for an area of comparable size. Sulfate (SO₄^2–) values, for example, range from 1.2 to 11,500 mg/L and oxygen and hydrogen isotopes do not attenuate towards the weighted annual mean. Those samples that are high in sulfate generally have a low pH, which immediately indicates battery acid (H₂SO₄) contamination as a potential source for the sulfate. The low pH and high reactivity of the sulfuric acid groundwater cause the formation of hydrogeological microenvironments due to preferential dissolution of carbonate minerals, which in turn causes enhanced recharge and groundwater flow in certain areas; thus, the extreme scatter in the data set. Because of the difficult hydrogeology it is not straightforward to delineate the point-sources of contamination and up to five potential scenarios have to be evaluated: (1) seawater intrusion, (2) upwelling of high-sulfate groundwater, (3) local dissolution of gypsum, (4) an up-gradient contaminant source to the northeast of the GCR property and (5) battery acid contamination.     

Palaeozoic structures in the Xayacatlán area,Acatlán Complex,southern Mexico: transtensional rift‐ and subduction‐related deformation along the margin of Oaxaquia

The Xayacatlán area (eastern Mixteca terrane, southern Mexico) was previously inferred to preserve the Ordovician‐Silurian thrust contact between vestiges of the Iapetus Ocean and the para‐autochthon bordering Oaxaquia. Detailed remapping indicates that the rocks occur in four vertically‐bounded, NS fault blocks. The latter record the following tectonothermal events that post‐date Iapetus and occurred along the margins of the Rheic (1) and Pacific (2 and 3) oceans: (1) dextral transtension accompanying intrusion of an NS, tholeiitic dike swarm at ～442 Ma; (2) penetrative, greenschist‐facies deformation during the Mississippian related to extrusion of high‐pressure rocks; and (3) subgreenschist‐facies dextral transtension on NS faults during the generation of Middle Permian fabrics.     

Geophysical methods applied to fault characterization and earthquake potential assessment in the Lower Tagus Valley, Portugal

The study region is located in the Lower Tagus Valley, central Portugal, and includes a large portion of the densely populated area of Lisbon. It is characterized by a moderate seismicity with a diffuse pattern, with historical earthquakes causing many casualties, serious damage and economic losses. Occurrence of earthquakes in the area indicates the presence of seismogenic structures at depth that are deficiently known due to a thick Cenozoic sedimentary cover. The hidden character of many of the faults in the Lower Tagus Valley requires the use of indirect methodologies for their study. This paper focuses on the application of high-resolution seismic reflection method for the detection of near-surface faulting on two major tectonic structures that are hidden under the recent alluvial cover of the Tagus Valley, and that have been recognized on deep oil-industry seismic reflection profiles and/or inferred from the surface geology. These are a WNW–ESE-trending fault zone located within the Lower Tagus Cenozoic basin, across the Tagus River estuary (Porto Alto fault), and a NNE–SSW-trending reverse fault zone that borders the Cenozoic Basin at the W (Vila Franca de Xira–Lisbon fault). Vertical electrical soundings were also acquired over the seismic profiles and the refraction interpretation of the reflection data was carried out. According to the interpretation of the collected data, a complex fault pattern disrupts the near surface (first 400 m) at Porto Alto, affecting the Upper Neogene and (at least for one fault) the Quaternary, with a normal offset component. The consistency with the previous oil-industry profiles interpretation supports the location and geometry of this fault zone. Concerning the second structure, two major faults were detected north of Vila Franca de Xira, supporting the extension of the Vila Franca de Xira–Lisbon fault zone northwards. One of these faults presents a reverse geometry apparently displacing Holocene alluvium. Vertical offsets of the Holocene sediments detected in the studied geophysical data of Porto Alto and Vila Franca de Xira–Lisbon faults imply minimum slip rates of 0.15–0.30 mm/year, three times larger than previously inferred for active faults in the Lower Tagus Valley and maximum estimates of average return periods of 2000–5000 years for M 6.5–7 co-seismic ruptures.     

Snow cover area change and its relations with climatic variability in Kashmir Himalayas,India

Snowmelt makes an essential component of the hydrological system of Kashmir Himalayas. The present study was carried out to examine the status of Snow Cover Area (SCA) using Moderate Resolution Imaging Spectroradiometer (MODIS) 8-day Snow Cover Product between 2000 and 2016. The intra- and inter-annual variability in SCA and in meteorological parameters was observed and various statistical tests were used to study the interrelationship. Results of statistical analysis indicate decrease in maximum temperature (?0.05 °C/year) and minimum temperatures (?0.02 °C/year) while rise in precipitation (19.13 mm/year). It also showed an increase in annual mean SCA (43.5 sq km) during the study period. The analysis was also carried out on a seasonal basis. The results revealed that in Kashmir Himalayas, climate plays a dominating role in controlling the SCA. The results depict the short-term fluctuations in SCA and show the magnitude of change between two successive values being very large in SCA.     

A review on advances in seismology in Mexico after 30 years from the 1985 earthquake

The 19 September 1985 (M_w8.1) earthquake, located on the Michoacán coast, Mexico, generated great damage in Mexico City, more than 300 km away from the epicentral area. Other important cities near the coast and in central Mexico also suffered severe damage. Thirty years after this important event, the National Autonomous University of Mexico (UNAM), the Centro de Investigación Científica y de Educación Superior de Ensenada, Baja California (CICESE) and other institutions organized a conference to discuss the scientific advances, particularly in seismology, that had taken place in Mexico since then.     

A simplified method for simulation of strong ground motion using finite rupture model of the earthquake source

We present a simplified method to simulate strong ground motion for a realistic representation of a finite earthquake source burried in a layered earth. This method is based on the stochastic simulation method of Boore (Boore, D. M., 1983, Bull. Seism. Soc. Am. 73, 1865–1894) and the Empirical Greens Function (EFG) method of Irikura (Irikura, K., 1986, Proceedings of the 7th Japan Earthquake symposium, pp. 151–156). The rupture responsible for an earthquake is represented by several subfaults. The geometry of subfaults and their number is decided by the similarity relationships. For simulation of ground motion using the stochastic simulation technique we used the shapping window based on the kinetic source model of the rupture plane. The shaping window deepens on the geometry of the earthquake source and the propagation characteristics of the energy released by various subfaults. The division of large fault into small subfaults and the method for accounting their contribution at the surface is identical to the EGF. The shapping window has been modified to take into account the effect of the transmission of energy released form the finite fault at various boundaries of the layered earth model above the source. In the present method we have applied the correction factor to adjust slip time function of small and large earthquakes. The correction factor is used to simulate strong motion records having basic spectral shape of ² source model in broad frequency range. To test this method we have used the strong motion data of the Geiyo earthquake of 24th March 2001, Japan recorded by KiK network. The source of this earthquake is modelled by a simple rectangular rupture of size 24 × 15 km, burried at a depth of 31 km in a multilayered earth model. This rupture plane is divided into 16 rectangular subfaults of size 6.0 × 3.75 km each. Strong motion records at eight selected near-field stations were simulated and compared with the observed records in terms of the acceleration and velocity records and their response spectrum. The comparison confirms the suitability of proposed rupture model responsible for this earthquake and the efficacy of the approach in predicting the strong motion scenario of earthquakes in the subduction zone. Using the same rupture model of the Geiyo earthquake, we compared the simulated records from our and the EGF techniques at one near-field station. The comparison shows that this technique gives records which matches in a wide frequency range and that too from simple and easily accessible parameters of burried rupture.     

APPLICATION OF MAP AND FILE INFORMATIONVISUALIZATION SYSTEM TO COMPREHENSIVEDIVISION OF NATURAL DISASTERS- Takingthe ChangjiangValleg as an Example

In recent years，scholars at home and abroad have method，and that vanous disasters would be analyzed assynthetically studied natural dlsaste‘theoretically and an Integral．SHI Peilun（1991）putforwad a scientificmethodologlcally，as well as its cases analys。s（CND，term—regional disaster system，which Indicates that1987； MA et al，1990； MARBLE，1990； NE et al，the situation of a disaster（calamity loss）results from1999； PATAK et al，1982； SHI，1991；VAN et al，h…     

基于InSAR形变数据分析台湾纵谷断层南段现今运动特征

本文利用2007—2010年花东纵谷南段区域的InSAR形变数据作为约束, 采用分段断层模型和层状介质模型, 反演中国台湾东部纵谷断层南段滑动速率空间分布, 并据此分析断层运动特征。 研究结果表明, 纵谷断层南段整体以逆冲运动为主, 兼具左旋走滑运动。 纵谷断层南段的滑动速率具有空间非均匀性, 在空间上可以细分为深浅两个极值区, 浅部(0~15 km)最大滑动速率为10 cm/a, 位于深度2.5 km左右; 深部(15~30 km)最大滑动速率为21 cm/a, 位于深度25 km左右。 反演结果与用重复地震估算的深部滑动速率基本吻合。