Seismic_Attributes_Geologic
12 Pages

Seismic_Attributes_Geologic

Course Number: C 551, Fall 2009

College/University: McGill

Word Count: 8122

Rating:

Document Preview

Seismic Attributes a promising aid for geologic prediction Satinder Chopra+ and Kurt Marfurt* Arcis *Allied Corporation, Calgary, Alberta, Canada Geophysical Laboratories, University of Houston, Houston, Texas, USA expectation is enunciated in that seismic attributes will continue to improve accuracy of interpretations and predictions in hydrocarbon exploration and development. The 1990s saw 3D attribute...

Unformatted Document Excerpt
Coursehero >> Canada >> McGill >> C 551

Course Hero has millions of student submitted documents similar to the one
below including study guides, practice problems, reference materials, practice exams, textbook help and tutor support.

Course Hero has millions of student submitted documents similar to the one below including study guides, practice problems, reference materials, practice exams, textbook help and tutor support.

Attributes Seismic a promising aid for geologic prediction Satinder Chopra+ and Kurt Marfurt* Arcis *Allied Corporation, Calgary, Alberta, Canada Geophysical Laboratories, University of Houston, Houston, Texas, USA expectation is enunciated in that seismic attributes will continue to improve accuracy of interpretations and predictions in hydrocarbon exploration and development. The 1990s saw 3D attribute extractions become commonplace in the interpretation work place. During this time seismic interpreters were making use of dip and azimuth maps, amplitude extractions and seismic sequence attribute mapping concepts were also established. Seismic coherence was introduced, which distinctly revealed buried deltas, river channels, reefs and dewatering features. The remarkable detail with which stratigraphic features show up on coherence displays, without any interpretation bias, and some unidentifiable even with close scrutiny, appealed to the interpreters and this has significantly changed the way geophysicists interpret 3D seismic data. In 1999, Greg Partyka introduced the application of spectral decomposition to stacked seismic data. This allowed interpreters to utilize the discrete frequency components of the seismic bandwidth to interpret and understand the subtle details of subsurface stratigraphy. During the same year Connolly (1999) introduced elastic impedance, which computes conventional acoustic impedance for non-normal angle of incidence. Crossplotting of attributes was introduced to visually display the relationship of two or three variables (White, 1991). Such crossplots have been used since as AVO anomaly indicators (Hilterman and Verm, 1994). Automated pattern recognition on attributes through neural network applications, geostatistics and multivariant regression analysis have helped combine attributes sensitive to relevant geological features. Waveform classification schemes were introduced by Elf Acquitaine in the mid 1990s, where the interpreter defines a zone of interest pegged to an interpreted horizon, and then asks the computer to define a suite of approximately 10-20 waveforms that best express the data. Keeping pace with such emerging technologies were the advancements in visualization. Starting at the seed voxels, a seed tracker will search for connected voxels that satisfy the user-defined search criteria, thereby generating a 3D geobody within the 3D seismic volume. Significant developments: Recent past Spectral decomposition This work continues actively today, with most workers preferring the wavelet transform based a p p roach introduced by Castagna et al. (2003) over the original discrete Fourier transform. Continued on Page 111 Seismic attributes have come a long way since their intro d u ction in the early 1970s and have become an integral part of seismic interpretation projects. Today, they are being used widely for lithological and petrophysical prediction of re s e rvoirs and various methodologies have been developed for their application to broader hydrocarbon exploration and development decision making. Beginning with the digital re c o rding of seismic data in the early 1960s and the ensuing bright spot analysis, the 1970s saw the introduction of complex trace attributes and seismic inversion along with their color displays. This was followed by the development of response attributes, introduction of texture analysis, 2D attributes, horizon and interval attributes and the pervasive use of color. 3D seismic acquisition dominated the 1990s as the most successful exploration technology of several decades and along with that came the seismic sequence attributes. The coherence technology introduced in the mid 1990s significantly changed the way geophysicists interpreted seismic data. This was followed by the introduction of spectral decomposition in the late 1990s and a host of methods for evaluation of a combination of attributes. These included pattern recognition techniques as well as neural network applications. These developments continued into the new millennium, with enhanced visualization and 3D computation and interpretation of texture and curvature attributes coming to the forefront. Of course all this was possible with the power of scientific computing making significant advances during the same period of time. A detailed re c o nstruction of these key historical events that lead to the modern seismic attribute analysis may be found in Chopra and Marfurt (2005). The proliferation of seismic attributes in the last two decades has led to attempts to their classification and to bring some order to their chaotic development. We pick up the threads of this development over the last five to six years that has lead to the present seismic attribute analysis tools and workflows and also discuss some applications of seismic attributes that demonstrate the direction in which the attribute development is headed. Finally, a reasonable Figure 1. A time slice (a) at 2000 ms showing the channel as seen on the 3D seismic (b) 17 Hz amplitude slice from spec tral decomposition showing channel complex and the fault trend. (After Ojo and Sindiku, 2003). 110 CSEG RECORDER 2006 Special Edition Seismic attributes Seismic attributes a promising aid... Continued from Page 110 Contd Figures 1a and b show examples where the application of spectral decomposition on 3D seismic data from Southern Nigeria helped map the geometry of incised channels, one of the primary objectives of the study (Ojo and Sindiku, 2003). Conventional spectral decomposition (using Fourier Transform) was applied to 3D seismic data and viewed in map form. The spectral responses within the zone-of-interest showed the lateral variability of the channel complex clearly and indicated the fault trend as well, features not evident on the traditional time slices. An interesting fact that emerged out of this analysis was that both the time and depth maps generated for the channel bases and the isopach maps of the channel-fill thickness conform to the geometry as obtained from spectral decomposition. Besides imaging the subsurface features, spectral decomposition also helps discriminate the most significant frequency content of a reservoir and so optimizes the seismic image for stratigraphic and/or DHI interpretation (Fahmy et al., 2005). Figure 2 shows a segment of a seismic section from a field in offshore West Africa. The discovery well drilled on this field encountered a 45 m net oil pay and a wet sand downdip by the sidetracked well. Another well was later proposed to prove the presence of a thicker reservoir updip, and evaluate its quality. However, since the seismic amplitude extraction did not support the presence of thicker or better quality sands updip from the first well, the drilling of this well was put on hold. Instantaneous spectral analysis was carried out on several seismic profiles that intersected the proposed well location. Analysis of the resultant frequency gathers showed that the majority of the signals strength for the reservoir reflections was confined to a narrow band of frequencies centered around 11 Hz and below 20 Hz. Consequently, a bandpass filter was designed and applied to the seismic data to preserve the signal bandwidth of the reservoir and attenuate all energy above 20 Hz. The resultant data showed lower resolution, but a pronounced low impedance anomaly corresponding to the reservoir showed up clearly. This significantly helped the interpretation of the stratigraphy from the seismic data. The well in question was later drilled and it encountered 154 m of net oil sands in the target reservoir, validating the pre-drill prognosis of better sand facies than the earlier well. This example demonstrates the application of spectral decomposition for delineating a reservoir specific bandwidth and optimizing the seismic data for DHI and stratigraphic interpretation. Spectral decomposition has been applied not only in the time domain, but in the depth domain as well. Montoya et al., (2005) demonstrate the application of spectral decomposition in the depth domain to an area in the Gulf of Mexico, which helped understand the distribution and classification of deep-water geological elements. Seismic inversion revisited The original recursive or trace-integration seismic inversion technique for acoustic impedance also evolved during the late 1980s and 1990s, with developments in model-based inversion, sparsespike inversion, stratigraphic inversion, and geostatistical inversion, providing accurate results (Chopra and Kuhn, 2001). The earlier techniques used a local optimization method that produced good results when provided with an accurate starting model. Local optimization techniques were followed by global optimization methods that gave reasonable results even with sparse well control. Figure 2. A seismic segment showing comparison of the (a) original seismic with that (b) after optimizing reservoir bandwidth. (After Fahmy et al., 2005) Connolly (1999) introduced elastic impedance, which computes conventional acoustic impedance for non-normal angle of incidence. This was further enhanced by Whitcombe (2002) to reflect different elastic parameters such as Lams parameter , bulk modulus , and shear modulus . Figure 3 shows maps of average extended elastic impedance averaged over a 25 ms time gate measured from top reservoir for the Forties field in the central North Sea. The images are tuned to shear and bulk modulus. These images show the clear imaging of the channel systems and, potentially, remaining hydrocarbons within these channels. Crossplotting of attributes Figure 3. Maps generated using extended elastic impedance for a dataset from North Sea showing the (a) optimized lithology and (b) fluids (After Whitcombe et al., 2002). Crossplotting of attributes was introduced to visually display the relationship between two or three variables Continued on Page 112 2006 Special Edition CSEG RECORDER 111 Seismic attributes Continued from Page 111 Contd gamma ray logs in the area were diagnostic of sands, and there was a fairly even sampling of well data across the field. A non-linear multi-attribute determinant analysis was employed between the derived multiple seismic attribute volumes and the measured gamma ray values at wells. By training a neural network with a statistically representative population of the targeted log responses (gamma ray, sonic, and bulk density) and the multiple seismic attribute volumes available at each well, a non-linear multiattribute transform was computed to produce gamma ray and bulk density inversions across the 3D seismic volume. In Figures 6a and b we show the Lambda-Rho and Mu-Rho sections with the anomaly enclosed in a yellow polygon. Probabilistic neural network solutions were used to derive gamma inversion and are shown in Figure 6c. The data are scaled to API gamma units in Figure 6c and converted to p o rosity in Figure 6d using the standard linear density relationship. From log data, the sand filled channels are interpreted as having gamma values less than 50 API units. This cut off value was used to mask out inverted density values for silts and shales. Analysis of Figures 6c and d shows three distinct sand bearing channels. Cubic B-spline curves (mathematical re p resentation of the a p p ro x i m a t i n g curves in the form of polynomials) have also been used for determination of mathem a t i c a l re l a t i o n s h i p s between pairs of variables from well logs and then using them to invert Seismic attributes a promising aid... (White, 1991). Hilterman and Verm (1994) used crossplots in AVO analysis, which have been used since as AVO anomaly indicators. When appropriate pairs of attributes are cross-plotted, common lithologies and fluid types often cluster together, providing a straightforward interpretation. The off-trend aggregations can then be more elaborately evaluated as potential hydrocarbon indicators, keeping in mind the fact that data that are anomalous statistically are geologically interesting - the essence of successful AVO crossplot analysis. Extension of crossplots to three dimensions is beneficial as data clusters hanging in 3D space are more readily diagnostic, resulting in more accurate and reliable interpretation. In Figure 4 we illustrate the use of modern crossplotting software of three attributes that help identify a gas anomaly P-reflectivity on the x-axis, S-reflectivity on the y-axis and fluid stack on the z-axis. In Figure 4a, we indicate a gas anomaly on a time slice through the P-reflectivity volume by a yellow-red patch. We then draw a yellow polygon on the time slice to select live data points pertaining to the anomaly to be displayed in the crossplot. The red cluster of points in Figure 4b corresponds to the red polygon on the time slice (enclosing the live data points). As the crossplot is rotated toward the left on the vertical axis, the fluid stack shows the expected negative values for the gas sand (Figure 4c). In Figure 5 the AVO attribute Lambda-Rho, Lambda/Mu and Fluid stack are crossplotted. On the Lambda-Rho - Lambda/Mu crossplot the anomaly is seen as yellow cluster showing low values of Lambda-Rho and Lambda/Mu. As the crossplot is turned on the vertical axis, the fluid stack indicates negative values corresponding to the low values of Lambda-Rho and Lambda/Mu, indicating the gas sand clearly. Automated pattern recognition on attributes The attribute proliferation of the 1980s resulted in an explosion in the attribute data available to geophysicists. Besides being overwhelming, the sheer volume of data defied attempts to gauge the information contained within those data using conventional analytical tools, and made their meaningful and timely interpre t ation a challenge. For this reason, one school of geophysicists examined automated pattern recognition techniques, wherein a computer is trained to see the patterns of interest and then made to sift through the available bulk of data seeking those patterns. A second school of geophysicists began combining attributes sensitive to relevant geological features through multiattribute analysis. Neural network application for multiattribute analysis One attempt at automated pattern recognition took the form of neural networks (Russell et al., 1997), wherein a set of input patterns is related to the output by a transformation that is encoded in the network weights. In Figure 6 we show an example on how multivariate statistical analysis can be used in determining whether the derived property volumes are related to gas saturation and lithology (Chopra and Pruden, 2003). For the case study from southern Alberta, it was found that the Figure 4. A time slice through (a) the P-reflectivity (Rp) volume. The suspected gas anomaly is indicated by high (yellow-red) values on the slice (b) Cross-plot of Rp vs. Rs. The red polygon encloses all the live data points on both time slices while the yellow polygon encloses the suspected anomaly. The cross-plot shows the yellow points corresponding to high values of Rp exhibiting negative values of fluid stack as is expected of a gas anomaly. (After Chopra et al., 2003). Continued on Page 113 112 CSEG RECORDER 2006 Special Edition Seismic attributes Seismic attributes a promising aid... Continued from Page 112 Contd Figure 5. Crossplot of Lambda-Rho and Lambda/Mu attributes derived from Rp and Rs shown in Figure 4. Points within the red and yellow polygons (Figure 4) show up as different clusters. The gas anomaly (enclosed by the yellow polygon) shows up with low values of Lambda-Rho and Lambda/Mu and negative values for the fluid stack. (b) 3D crossplot seen from the fluid stack side (After Chopra et al., 2003). attribute volumes into useful inversion volumes like gamma ray and porosity (Chopra et al., 2004). In F i g u re 7 we show spline curve inverted porosity. Notice that the results from the cubic B-spline inversion and the neural network analysis look similar. It may be mentioned that the results could differ significantly if t h e re is a significant variation in geology laterally, in which case the single well used for the spline analysis is no longer representative of the entire area. Enhanced visualization helps attribute interpretation Gradually as geophysicists realized that the additional benefits provided by 3D seismic were beneficial for stratigraphic interpretation of data, seismic interpretation methods also shifted from simple horizon-based to volume-based work. This provided interpreters new insights that were gained by studying objects of different geological origins and their spatial relationships. In Figure 8 we display strat cubes (sub-volumes bounded by two not necessarily parallel horizons), generated from the c o h e rence volumes. The coherence strat cube indicates the N-S channel very clearly, the E-W fault on the right side as well as the down thrown side of the N-S fault on the left. Of course with all this also came the complexity and the magnitude of identification work and the need for faster and more accurate tools. This brought about the significant introduction of techniques for automated identification of seismic objects and stratigraphic features. Geobody tracking is now possible by selecting an amplitude value corresponding to a feature of interest and letting the seed tracker search for connected voxels that satisfy a user- d e f i n e d search criterion. Complete geobodies corresponding to a porous sand accumulation, or a salt body or the channel-fill within a drainage system can be conveniently determined. Figure 9 shows the geobody-tracking of a porous zone in 3D space by selecting a seed amplitude at the point indicated by the small re d sphere. Such geobody-tracking helps Continued on Page 114 2006 Special Edition CSEG RECORDER 113 (a) (b) (c) (d) Figure 6. Time slices through (a) and (b) and volumes. The suspected gas anormaly is indicated by low (blue) values of and high values of and . Neural network inverted gamma ray response is shown in (c). Note the distinct separation of sand from silt and shale. (d) Neural network computed porosity from the inverted density response. The density values have been masked out for gamma ray values representative of silt or shale, giving a rela tive porosity indicator for the sands. (After Chopra and Pruden, 2003). Seismic attributes Continued from Page 113 Contd The network transforms the chosen attributes into a new metaattribute, which indicates the probability of occurrence of the identified feature at different seismic positions. Such highlighted features definitely benefit from the knowledge of shapes and orientations of the features that can be added to the process. Trace shape While spectral decomposition and wavelet analysis compare seismic waveform to precomputed waveforms (typically windowed tapered sines and cosines), an important development was released by Elf Acquitaine in the mid 1990s trace shape classification. In this approach, the interpreter defines a zone of interest pegged to an interpreted horizon, and then asks the computer to define a suite of approximately 10-20 waveforms that best express the data. The most useful of these classifiers is based on Self-Organized Maps (Coleou et al., 2003) that provides maps whose appearance is relatively insensitive to the number of classes. Although the results can be calibrated to well control through forward modeling, and although actual well classes can be inserted, this technology is particularly well suited to a geomorphology driven interpretation, whereby the interpreter identifies depositional and structural patterns from the images and from these infers reservoir properties. Figure 10 shows how waveform identification was used for identification of anomalies. Notice that the anomalies seen in purple are close to what was found on the AVO attribute slices shown in Figure 6. Seismic attributes a promising aid... Figure 7. Spline curve inverted porosity corresponding to the time slices shown in Figure 6 (After Chopra et al., 2004). in assessing the hydrocarbon reserves in place in a given prospect. While one given attribute will be sensitive to a specific geologic feature of interest, a second attribute may be sensitive to a different kind of feature. We can therefore combine multiple attributes to enhance the contrast between features of interest and their surroundings. Diff e rent methodologies have been developed to recognize such features. Meldahl et al. (2001) used neural networks trained on combinations of attributes to recognize features that were first identified in a seed interpretation. (a) (b) (c) (d) (e) (f) Figure 8. Strat cubes from a coherence sub-volume of 160 ms thickness below a horizon of interest at (a) 68 ms (b) 72 ms (c) 88 ms (d) 120 ms (e) 140 ms and (f) 158 ms. Notice the clarity with which different features show up as the strat-cube is sliced through (Data courtesy: Arcis Corporation, Calgary). Continued on Page 115 114 CSEG RECORDER 2006 Special Edition Seismic attributes Seismic attributes a promising aid... Continued from Page 114 Contd Texture attributes More recently the idea of studying seismic textures has been revived. While the term was earlier applied to seismic sections to pick out zones of common signal character (Love and Simaan, 1984), studies are now underway to use statistical measures to classify textures using gray-level co-occurrence matrices (Vinther et al., 1995; Vinther, 1997; Whitehead et al., 1999; West et al., 2002; Gao, 2003, 2004). Some of the statistical measures used are energy (denoting textural homogeneity), entropy (measuring predictability from one texel or voxel to another), contrast (emphasizing the difference in amplitude of neighboring voxels) and homogeneity (highlighting the overall smoothness of the amplitude). Energy, entropy and homogeneity have been found to be the most effective in characterizing seismic data. Figure 11a shows a strat-slice display from a 3D seismic volume at the reservoir level where wells W1 and W2 are producing from the same Lower Cretaceous glauconitic gas sand (Chopra and Alexeev, 2005). However, there is no indication of this on the strat-slice. Texture attribute displays comprising the energy, entropy and the homogeneity are shown in Figure 11b to d respectively. Notice the energy attribute indicates high-energy values for the producing formation for wells W1 and W2. These high values of energy are also associated with low entropy and high homogeneity, a typical combination that is expected for fluvial deposits (Gao, 2003). Curvature With the wide availability of 3D seismic and a renewed interest in fractures, we have seen a rapid acceleration in the use of curvature maps. The structural geology relationship between curvature and fractures is well established (Lisle, 1994) though the exact relationship between open fractures, paleo structure, and present-day stress is not yet clearly understood. Roberts (2001), Hart et al. (2002), Sigismondi and Soldo (2003), Massafero et al. (2003) and others have used seismic measures of reflector curvature to map subtle features and predict fractures. Curvature (a three dimensional property of a quadratic surface that quantifies the degree to which the surface deviates from being planar) attribute analysis of surfaces helps to remove the effects of regional dip and emphasizing small-scale features that might be associated with primary depositional features or small-scale faults. Figure 12a shows a comparison of a time surface from a 3D seismic volume in northwest Alberta. The time surface indicates a pro m inent fault crossing a couple of narrow vertical channels. On the equivalent coherence display (Figure 12b) the main fault and boundaries of some of the channels can be seen clearly. The most negative curvature display (Figure 12c) indicates the base of the narrow channels and how the fault intersects them is evident. The most-positive curvature images the edges of the channels clearly. Notice the crisp definition of one of the channels shown with blue arrows on the most-positive curvature (Figure 12d) and the base of the channel imaged clearly on the most-negative curvature (red arrows in Figure 12c). A similar comparison is shown in Figure 13 where a meandering channel is the only prominent feature seen to the left of the time surface (Figure 13a). The equivalent coherence display (Figure 13b) shows a crisper definition for the channel and a couple of other features indicated with arrows. The most-negative curvature (Figure 13c) shows the base of the channels seen in Figure 13a, but the lower main channel boundary is seen in a different direction (magenta arrow). In fact, there is another channel seen to the southwest of the image. The edges of the channels are seen clearly in Figure 13d. Examples of present-day workflows (1) Attributes used to generate sand probability volumes: When attributes are tied to the available well control, they can be correlated to petrophysical properties, and this helps the interpreter to identify and associate high correlations with specific properties. For Continued on Page 116 2006 Special Edition CSEG RECORDER 115 Figure 9. Geobody-tracking of a porous zone in 3D space by selecting a seed ampli tude at the point indicated by the small red sphere. (Data courtesy: Arcis Corporation, Calgary). Figure 10. Waveform classification used for identification of anomalies for the dataset used in Figure 5, 6 and 7. (After Chopra et al., 2004). Seismic attributes Continued from Page 115 Contd areas of a reservoir, where changes from time-lapse are evident (Galikeev and Davis, 2005). As for example, reservoir condition due to CO2 injection could be detected. They generate attributes that represent reservoir heterogeneity by computing short time window seismic attributes parallel to the reservoir. Such an analysis in short temporal windows ensures that attribute carries an overprint of geology (Partyka et al., 1999). Figure 17 shows the dynamic changes within the Weyburn re s e rvoir (Canada), due to the increased CO2 saturation. This was done by computing the inverted impedance model of the reservoir on the differenced volume of the baseline (2000) and second monitor (2002) surveys. Figure 18 is a computed CO2 saturation map, where the values do not represent absolute CO2 saturation, but rather an estimation of part porosity occupied by CO2 after irreducible water and oil were taken into account. (c) 4-D seismic attributes together with 4-D rock and fluid analysis and incorporation of production engineering information, have been used for pressure-saturation inversion for time-lapse seismic data producing quantitative estimates of reservoir pressure and saturation changes. Application of such an analysis to the Cook reservoir of the Gulfaks field, off s h o re Norway (Lumley et al., 2003), shows that a strong pressure anomaly can be estimated in the vicinity of a horizontal water injector, along with a strong water saturation anomaly drawing towards a nearby producing well (Figure 19). This is in addition to strong evidence of east-west fault block compartmentalization at the time of the seismic survey. Seismic attributes a promising aid... example, Figure 14 shows how attributes from prestack inversion of a high-resolution seismic allowed dataset mapping of sand bodies in a geologically complex area. A key step in the workflow was the petro-elastic analysis of well data that demonstrated that seismic attributes derived from prestack seismic inversion could discriminate between sands and shales. A multi-attribute classification approach, incorporating neural network training techniques, was used to generate sand probability volumes derived from P-wave and S-wave impedances estimated using AVO inversion. The study demonstrated that high-resolution seismic data coupled with targeted inversion can increase confidence and reduce uncertainty. A crucial problem in any multi-attribute analysis is the selection and the number of seismic attributes to be used. Kalkomey (1997) showed that the probability of observing a spurious correlation increases as the number of control points decreases and also if the number of seismic attributes being used increases. A way out of such a situation is to withhold a percentage of the data during the training step and then later to use this hidden data to validate the predictions (Schuelke and Quirein, 1998). (2) Time-lapse analysis (a) Seismic attributes are being used effectively for time-lapse data analysis (4-D). Time-lapse data analysis permits interpretation of fluid saturation and pressure changes, and helps understanding of reservoir dynamics and the performance of existing wells. Figure 15 shows an example from east of Schiehallion field West of the Shetlands (Parr and Marsh, 2000). The pre-production surveys in 1993 in (a) and 1996 in (b) show a high degree of similarity, (a) but the 4-D survey (in (c)) shows large changes around producers and injectors. The poor production rates and low bottom hole-flowing pressures led to the conclusion that well C was located in a compartment that is poorly connected to injection support. The areal extent of this compartment could be picked by the amplitude increase seen on 4-D image and interpreted due to gas liberated from solution. This area is consistent with predictions from material-balance calculations. Figure 15c from the 1999 survey suggested the possibility of a connection (marked by an arrow) between producers C and D. The existence of such a connection was also suspected from the material-balance analysis. Figure 16 shows a coherence display at the required level and depicts the expected connection (marked by a circle). While a plausible explanation for this is not known, it is postulated that the attributes on 4-D seismic provide a clue that a transmissibility barrier may have been bro k e n between the injector and producer. (b) Reservoir based seismic attributes are being used to help delineate anomalous (c) (b) (d) Figure 11. Strat-slices at the reservoir level from (a)seismic (b) energy (c) entropy and (d) homogeneity (texture attrib utes). Figure 11b shows high values of energy corresponding to the producing sandstone. These high values of energy a re associated with low values of entropy and high values of homogeneity. (After Chopra and Alexeev, 2005) Continued on Page 117 116 CSEG RECORDER 2006 Special Edition Seismic attributes Seismic attributes a promising aid... Continued from Page 116 Contd (3) Attributes for detection of gas zones below regional velocity inversion: Evaluation of gas distribution in the fluid system in the Wind River Basin where there are anomalously pressured gas accumulations is a challenge. Using the available logs and seismic data, the regional velocity inversion surface has been mapped in this (a) (b) (c) (d) area (Surdam et al., 2004b), which is the pressure surface separating the anomalously pressured rocks below from the normally pressured rocks above. Seismic attributes have been successfully used to evaluate the distribution of sandstone rich intervals within the prospective reservoir units. The Frenchie Draw gas field in the Wind River Basin is an example of an area where detecting and delineating gas zones below the regional velocity inversion surface is d i fficult. The stratigraphic interval of interest is the Upper Cretaceous-Paleocene Fort Union/Lance lenticular fluvial sandstone formations on a north plunging structural nose. The gas distribution pattern in the formations has been found to be complex, and so the exploitation proved to be risky. Surdam (2004a) has demonstrated that a good corre l a t i o n exists between seismic frequency and gamma ray logs (lithology) in the lower Fort Union/Lance stratigraphic interval. The frequency attribute was used to distinguish sandstone-rich intervals from shalerich intervals. F i g u re 20 shows a frequency attribute section (with seismic data overlaid) covering the Fort Union/Lance stratigraphic interval intersecting the anomalously slow velocity domains (outlined by white dots). In addition to the northplunging structural nose seen in the area, a shale rich sequence (seen in orange) is seen near the upper edge of the gas production. The important thing seen here is the lenticular distribution of the sandstone-rich intervals in blue that stand out against the shale-rich intervals in orange, yellow and g reen. This distributional pattern of lithologies corresponds well with the initial interpretations carried out by geoscientists who discovered the field. The present There are many areas in the attribute world where active development is underway and we mention here some of the prominent ones. 1. Volumetric estimates of curvature: Cracks or small discontinuities are relatively small and fall below seismic re s o l u t i o n . However, the presence of open and closed cracks is closely related to reflector curvature (since tension along a surface increases with increasing curvature and therefore leads to fracture). Until now, such curvature estimates have been limited to the analysis of picked horizons, which previously may be affected by unintentional bias or picking errors introduced during interpretation. Volumetric curvature Continued on Page 118 2006 Special Edition CSEG RECORDER 117 Figure 12. Coherence and curvature computation results shown for a time horizon shown in (a). Figure (b) shows the coherence along the time surface indicating the main fault running NE-SW. The base of the N-S channels are most pronounced on the most-negative curvature (Figure (c ) and the edges of the channels imaged well on the most-positive curvature (Figure (d)) (Data courtesy: Arcis Corporation, Calgary). (a) (b) (c) (d) Figure 13. (a) Time surface (b) coherence (c) most-negative curvature and (d) most-positive curvature. Curvature attributes indicate a better focusing of the base and edges of the channels and other features as compared with coherence. (Data courtesy: Arcis Corporation, Calgary). Seismic attributes Continued from Page 117 Contd decomposition to common offset volumes to better predict changes in lithology and flow barriers. While the interpretation is entirely consistent with AVO analysis, the 3D volumes tend to show discrete geologic features and the limits of hydrocarbon distribution. There is a similar relationship between amplitude vs. azimuth and geometric attributes applied to common azimuth volumes. In this latter case, we often find faults and f r a c t u res better illuminated at those acquisition azimuths perpendicular to structural trends. Gazing into the Future Given the current data explosion both of large regional 3D surveys, but also of new time-lapse and multicomponent surveys, we envision an increasingly rapid evolution of seismic attributes and computer-aided interpretation technology. Some of the clear signs on the horizon are as follows: Seismic attributes a promising aid... computation entails, first, the estimation of volumetric reflection dip and azimuth that represents the best single dip for each sample in the volume, followed by computation of curvature from adjacent measures of dip and azimuth. The result is a full 3D volume of curvature values at different scales of analysis (al Dossary and Marfurt, 2006). 2. 3D classifiers: Supervised 3D classification is being used for integrating several seismic attributes into a volume of seismic facies (Sonneland et al., 1994, Carrillat et al., 2002). Some of the most successful work is in imaging by-passed pay (Xue et al., 2003) using attributes sensitive to amplitude and trace shape. Workers at deGroot-Bril, Paradigm, Rock Solid Images and many others have made progress in using geometric attributes as the basis to automatically pick seismic textures in 3D. 3. Structurally-oriented filtering: Recently, several workers have used the volumetric estimates of dip and azimuth to improve the signal to noise ratio, yet preserve discontinuities such as (a) faults and stratigraphic discontinuities. Hoecker and Fehmers (2002) use an a n i s o t ropic diffusion algorithm that smoothes along dip azimuth only if no discontinuity is detected. Luo et al. (2002) use a multiwindow analysis technique, smoothing in that window containing the (c) analysis point that has the smallest variance. Duncan et al. (2002) built on this latter technique, but instead of a mean filter, applied a principal component filter in the most coherent analysis window. These structurally-oriented filtering (alternatively called edge-preserving smoothing) algorithms improve not only behavior of autotrackers but also the results of coherence and other attributes sensitive to changes in reflector amplitude, (a) waveform, and dip. 4. Volumetric estimation of Q: The most common means of estimating Q is through sample by sample estimates of spectral ratios. Typical workflows use a smoothly varying estimate of Q to increase the spectral bandwidth at depth. Most workers have focused on intrinsic vs. effective Q, the later including effects such as geometric scattering and friendly multiples. The most promising work is done with the aid of VSPs and/or well logs. In these flows, one can interpret not only the spectral amplitude but also the spectral phase compensation necessary to improve seismic resolution. 5. Pre-stack attributes: In addition to AVO, which measures the change in reflection amplitude and phase as a function of offset at a fixed location, we can apply attributes such as coherence and spectral (b) Figure 14. (a) A cross-section from the final processed seismic data volume. The dipping event in the center of the panel is interpreted as a sand injection feature. (b) The same cross-section from the sand probability volume derived from multi-attribute classification. The classifi cation has predicted that the feature is sand being injected from the main sand body seen below. (c) The sand probability volume and amplitude data displayed using 3D visualization. The figure shows seismic ampli tude data in the background; the base reservoir surface is shown in blue and a possible sand injection feature mapped from the inversion results. Note the complexity of the injected sand bodies. (Images courtesy: Steve McHugo, WesternGeco) (b) (c) Figure 15. An example from east of Schiehallion. The net sand maps based on seismic amplitudes on the Pre-production surveys (1993) in (a) and 1996 in (b) show a high degree of similarity. Comparison with the 1999 4-D survey (c) shows large changes around producers and injectors. (After Parr and Marsh, 2000). Continued on Page 119 118 CSEG RECORDER 2006 Special Edition Seismic attributes Seismic attributes a promising aid... Continued from Page 118 Contd Figure 16. A time slice through a coherence volume corresponding to Figure 15 and depicts the expected connection (marked by a circle) between the producers C and D. (After Parr and Marsh, 2000). 1. We expect continued development of texture attributes that can quantify or enhance features used in seismic stratigraphy and seismic geomorphology leading to computer-aided 3D seismic stratigraphy. One of the major challenges is that due to tectonic deformation and sedimentary compaction, such patterns may be arbitrarily rotated from their original position. 2. We expect enhanced emphasis will be placed on time-lapse applications for delineation of flow-barriers, so that reservoir simulation provides more realistic estimates and information regarding the dynamic behavior of reservoirs. 3. We expect continued advances in 3D visualization and multiattribute analysis including clustering, geostatistics, and neural networks to alleviate the problems interpreters face due to an overwhelming number of attributes. Conclusions Figure 17. Position of the time-lapse impedance anomalies in depth relative to CO2 injectors (black) and vertical water injectors (blue). The overall size of the area shown is 9 sq. km. (After Galikeev and Davis, 2005). Figure 18. CO2 saturation map computed from time-lapse (2000-2002 inversion of the difference) impedance values. Shown are areas that responded to CO2 injection wells (white), horizontal injectors (black), unresponded to CO2 wells (blue) and vertical water injectors (yellow). The overall size of the area shown is 9 sq. km. (After Galikeev and Davis, 2005). A seismic attribute is a quantitative measure of a seismic characteristic of interest. Good seismic attributes and attribute analysis tools mimic a good interpreter. Over the past decades, we have witnessed attribute developments track bre a k t h roughs in seismic acquisition and reflector mapping, fault identification, bright spot identification, frequency loss, thin bed tuning, seismic stratigraphy, and geomorphology. More recently, interpreters have used crossplotting to identify clusters of attributes that are associated with either stratigraphic or hydrocarbon anomalies. Once again, the attribute community has worked hard to first duplicate such human-driven clustering through the use of self-organized maps, geostatistics, and neural nets, and then to extend this capability beyond the three dimensions easily visualized by interpreters. Tentative steps have been made towards computer- a s s i s t e d seismic stratigraphy analysis, whereby an interpreter trains the computer on a suite of structural or depositional patterns and asks the computer to find others like them. Pro g ress has been made in automated fault tracking, though currently the technology requires an expert user. In the not too distant future, we can envision an interpreter seeding a channel on a time slice after which the computer paints it in 3D. Although it may take decades, we expect computers will eventually be able to duplicate all the re p e titive processes performed by an interpreter. In contrast we do not Continued on Page 120 2006 Special Edition CSEG RECORDER 119 Seismic attributes Continued from Page 119 Contd thank many individuals for the help we received in the writing of that paper, amongst which the following names would figure prominently: Nigel Anstey, Tury Taner, Jamie Robertson, Alistar Brown, Roy Lindseth, Brian Russell, Roy White, Les Hatton, Bob Sheriff, Sven Treitel, Bee Bednar, Jerry Schuster and Art Barnes. Finally, we thank Ronald Parr for Figures 14 and 15, Tom Davis for Figures 16 and 17, Mark Meadows for Figure 19, and Ron Surdam for Figure 20. We have tried to cover almost all the prominent developments in this vast field of seismic attributes. However, in spite of our best efforts, it is likely we may have missed out on some. Hence any errors or omissions were unintentional and solely our responsibility. One of the authors (S.C.) expresses his appreciation to Arcis Corporation, Calgary for permission to use some of the data examples and publish this paper. R Seismic attributes a promising aid... expect them ever to replicate the creative interpreter imagining depositional environments, structural evolution, diagenetic alteration, and fluid migration. The human interpreter is here to stay, but may be off s h o red in the future. Acknowledgements This paper is an off-shoot of a more comprehensive review paper published in Geophysics (Chopra and Marfurt, 2005). We wish to References Al Dossary, S. and K. J. Marfurt, 2006, 3D volumetric multispectral estimates of reflector curvature and rotation: submitted to Geophysics. Carrillat, A., T. Randen, and L. Sonneland, 2002, Seismic stratigraphic mapping of carbonate mounds using 3D texture attributes: 64th Annual International Meeting, European Association of Geoscientists and Engineers, Z-99. Castagna, J. P., S. Sun. and R. W. Siegfried, 2003, Instantaneous spectral analysis: Detection of low-frequency shadows associated with hydrocarbons: The Leading Edge, 22, 120-127. Figure 19. Probability map on a scale of zero (blue) to 0.6 (white) that water satu ration (left) and pore pressure (right) have increased within the Cook reservoir of the Gulfaks field, Norwegian North Sea. Note that there is a strong pre s s u re anomaly surrounding the B-33 horizontal injector, along with east-west sealing fault compartmentalization. Water saturation change, however, is weak in most of the compartment since well B-33 injects into the water leg. The saturation change is stronger to the southeast of the compartment where water is drawn toward nearby producing well B-1. (After Lumley et al., 2003). Chopra, S. and D. Pruden, 2003, Multiattribute seismic analysis on AVO-derived param eters: The Leading Edge, 22, 998-1002. Chopra, S. and O. Kuhn, 2001, Seismic inversion: CSEG RECORDER: 26, no. 1, 10-14. Chopra, S., V. Alexeev, and Y. Xu, 2003, 3D AVO crossplotting - An effective visualiza tion technique: 73th Annual International Meeting, Society of Exploration Geophysicists Expanded Abstracts, 189-192. Chopra, S., D. Pruden, and V. Alexeev, 2004, Multi-attribute seismic analysis tackling non-linearity: First Break, 22, no.12, p 43-47. Chopra, S. and K. Marfurt, 2005 Seismic attributes- a histor ical perspective, Geophysics, 70, 3SO-28SO. Chopra, S. and V. Alexeev, 2005, Application of texture attrib utes to 3D seismic data, CSEG RECORDER, 30, September, 28-32. Coleou, T., M. Poupon, and K. Azbel, 2003, Unsupervised seismic facies classification: A review and comparison of tech niques and implementation: The Leading Edge, 22, 942-953. Connolly, P., 1999, Elastic impedance: The Leading Edge, 18, 438-452. Duncan, W., P. Constance, and K. J. Marfurt, 2002, Comparison of 3D edge detection seismic attributes to Vinton Dome, Louisiana, 72nd Annual International Meeting Society of Exploration Geophysicists, 577-580. Fahmy, W.A., G. Matteucci, D. Butters, J. Zhang, and J. Castagna, 2005, Successful application of spectral decomposi tion thechnology toward drilling of a key offshore development well: 75th Annual International Meeting, Society of Exploration Geophysicists Expanded Abstracts, 262-264. Galikeev, T., and T. Davis, 2005, Time-lapse seismic attributes and reservoir volumetric calculation: 67th Annual International Meeting, European Association of Geoscientists and Engineers, Z99. Gao, D., 2003, Volume texture extraction for 3D seismic visual ization and interpretation: Geophysics, 68, 1294-1302. Gao, D., 2004, Texture model regression for effective feature discrimination: Application to seismic facies visualization and interpretation: Geophysics, 69, 958-967. Figure 20. Seismic data display superimposed on a frequency attribute section at Frenchie Draw field. Shaded region shows anomalous velocity overlap. (After Surdam et al., 2004) Continued on Page 121 120 CSEG RECORDER 2006 Special Edition Seismic attributes Seismic attributes a promising aid... Continued from Page 120 Hart, B.S., R. Pearson, R. Pearson, and G. C. Rawling, 2002, 3D seismic horizon-based approaches to fracture-swarm sweet spot definition in tight-gas reservoirs: The Leading Edge, 21, 28-35. Hilterman, F. J. and R.W. Verm, 1994, Lithologic color-coded sections by AVO crossplots: 64th Annual International Meeting, Society of Exploration Geophysicists Expanded Abstracts, 1092-1095. Kalkomey, C. T., 1997, Potential risks when using seismic attributes as predictors of reser voir properties: The Leading Edge, 16, 247-251. Lisle, R. J., 1994, Detection of zones of abnormal strains in structures using Gaussian curvature analysis: AAPG Bulletin, 78, 1811-1819. Love, P. L. and M. Simaan, 1984, Segmentation of stacked seismic data by the classifica tion of image texture: 54th Annual International Meeting, Society of Exploration Geophysicists, Session:S7.3. Lumley. D., D. Adams, M. Meadows, S. Cole, and E. Ergas, 2003, 4D seismic pressuresaturation inversion at Gulfaks field, Norway: First Break, 21, 49-56. Luo, Y., W. G. Higgs, and W. S. Kowalik, 1996, Edge detection and stratigraphic analysis using 3D seismic data: 66th Annual International Meeting, Society of Exploration Geophysicists Expanded Abstracts 324-327. Masaferro, J. L., M. Bulnes, J. Poblet, and M. Casson, 2003, Kinematic evolution and fracture prediction of the Valle Morado structure inferred from 3D seismic data, Salta P rovince, northwest Argentina, American Association of Petroleum Geologists Bulletin, 87, 1083-1104. Meldahl, P., R. Heggland, B. Bril, and P. de Groot, 2001, Identifying faults and gas chimneys using multiattributes and neural networks: The Leading Edge, 20, 474-478. Montoya, P., R. Tatham, W. Fishwer, R. Steel, and M. Hudec, 2005, Definition of depo sitional geological elements in deep-water mini basins of the Gulf of Mexico using spectral decomposition in depth domain, 75th Annual International Meeting, Society of Exploration Geophysicists Expanded Abstracts, 481-484. Ojo, C. and I. Sindiku, 2003, Interpretation and analysis of a channel using 3D seismic and well log data a case study, 73th Annual International Meeting, Society of Exploration Geophysicists Expanded Abstracts, 462-465. Parr, R.S. and M. Marsh, 2000, Development of 4-D reservoir management West of Shetland, World Oil, 221, no.9, 39-47. Partyka, G., J. Gridley, and J. Lopez, 1999, Interpretational applications of spectral decomposition in reservoir characterization: The Leading Edge, 18, 353-360. Roberts, A., 2001, Curvature attributes and their application to 3D interpreted horizons: First Break, 19, 85-99. Robertson, J. D., 2005, Private communication. Corresponding author: Satinder Chopra (schopra@arcis.com) Contd Robertson, J.D. and D. A. Fisher, 1988, Complex seismic trace attributes: The Leading Edge, 7, no. 6, 22-26. Robertson, J. D. and H. H. Nogami, 1984, Complex seismic trace analysis of thin beds: Geophysics, 49, 344-352. Russell, B., D. Hampson, J. Schuelke, and J. Quirein,, 1997, Multiattribute seismic analysis: The Leading Edge, 16, 1439-1443. Schuelke, J. S. and J.A. Quirein,, 1998, Validation: A technique for selecting seismic attributes and verifying results: 68th Annual International Meeting, Society of Exploration Geophysicists Expanded Abstracts, 936-939. Sigismondi, M. and J. C. Soldo, 2003, Curvature attributes and seismic interpretation: Case studies from Argentina basins: The Leading Edge, 22, 1122-1126. Sonneland, L., P. Tennebo, T. Gehrmann, and O. Yrke, , 1994, 3D model-based Bayesian classification: 64th Annual International Meeting, Society of Exploration Geophysicists Expanded Abstracts, 510-511. Surdam, R. C., Z. S. Jioa, and Y. Ganshin, 2004a, Reducing risk in low-permeability gas formations: understanding the rock fluid characteristics of Rocky Mountain Laramide Basins: DOE Final Technical Pro g ress Report under contract no. DE-FC2601NT41325, 39. Surdam, R. C., Z. S. Jioa, and Y. Ganshin, 2004b, Reducing the risk of exploring for anomalously pressured gas assets: GasTIPS, Winter 2004, 4. Whitehead, P., J. Fairborn, and R. Wentland, 1999, Identifying stratigraphic units by seismic patterns, 69th Annual International Meeting, Society of Exploration Geophysicists Expanded Abstracts, 942-945. Vinther, R., 1997, Seismic texture classification applied to processed 2D and 3D seismic data, 67th Annual International Meeting, Society of Exploration Geophysicists Expanded Abstracts, 721-724. Vinther, R., K. Mosegaard, K. Kierkegaard, I. Abatzis, C. Andersen, and F. If, 1995, Seismic texture classification: A computer-aided approach to stratigraphic analysis, 65th Annual International Meeting, Society of Exploration Geophysicists Expanded Abstracts, 153-155. West, B., S. May, S., J. E. Eastwood, and C. Rossen, 2002, Interactive seismic facies clas sification using textural and neural networks: The Leading Edge, 21, 1042-1049. Whitcombe, D. N., 2002, Extended elastic impedance for fluid and lithology prediction: Geophysics, 67, 63-67. White, R. E., 1991, Properties of instantaneous seismic attributes: The Leading Edge, 10, no. 07, 26-32. (* Discussion and reply in TLE-11-8-45-46 and TLE-11-10-10-12) Satinder Chopra received M.Sc. and M.Phil. degrees in physics from Himachal Pradesh University, Shimla, India. He joined the Oil and Natural Gas Corporation Limited (ONGC) of India in 1984 and served there till 1997, before joining CTC Pulsonic at Calgary the following year, which later became Scott Pickford and Core Laboratories Reservoir Technologies. Currently, he is working as Manager, Reservoir Services at Arcis Corporation, Calgary. In the last 21 years Satinder has worked in regular seismic processing, and interactive interpretation, but has spent more time in special processing of seismic data involving seismic attributes including coherence and textures, seismic inversion, AVO, VSP processing and frequency enhancement of seismic data. His research interests focus on techniques that are aimed at characterization of reservoirs. He has published 3 books and more than 90 papers and abstracts and likes to make presentations at any beckoning opportunity. He is the Chairman of the SEG Publications Committee, the editor of the CSEG RECORDER and a member of the SEG The Leading Edge Editorial Board. He received several awards at ONGC; he and his colleagues have received the CSEG Best Poster Awards in successive years from 2002 to 2005, and very recently has been awarded the 2005 CSEG Meritorious Services Award. He is a member of SEG, CSEG, EAGE, AAPG, APEGGA (Association of Professional Engineers, Geologists and Geophysicists of Alberta) and TBPG (Texas Board of Professional Geoscientists). Kurt Marfurt began his geophysical career as an Assistant Professor teaching mining geophysics at Columbia Universitys Henry Krumb School of Mines in New York. After five years he joined Amoco at what was then called their Tulsa Research Center. Through successive reorganizations at Amoco, Marfurt has obtained diverse experience in seismic modeling, migration, signal analysis, basin analysis, seismic attribute analysis, reflection tomography, seismic inversion and multicomponent data analysis. Through Amoco he won 5 patents, 2 in seismic coherence technology. He is the author of several dozen articles and coeditor of two books published by the SEG, and has received both the SEG best presentation award for work on seismic modeling (1989) and SEG best paper award for work on seismic coherence (1999). Marfurt joined the University of Houston in 1999 as a Professor in the Department of Geosciences and as Director of the Allied Geophysics Laboratories where his primary emphasis is on the development and calibration of new seismic attribute technology. He is currently a member of the GSH, SEG, EAGE, AAPG and AGU, and has served 16 years as either an Associate or Assistant editor for Geophysics. 2006 Special Edition CSEG RECORDER 121

Textbooks related to the document above:

Find millions of documents on Course Hero - Study Guides, Lecture Notes, Reference Materials, Practice Exams and more. Course Hero has millions of course specific materials providing students with the best way to expand their education.

Below is a small sample set of documents:

Sveriges lantbruksuniversitet - ENSC - 805
1788IEEE TRANSACTIONS ON COMMUNICATIONS, VOL. 48, NO. 11, NOVEMBER 2000Further Analytical Results on the Joint Detection of Cochannel Signals Using Diversity ArraysStephen J. Grant and James K. CaversAbstractWe derive a fully analytical expres
Sveriges lantbruksuniversitet - ENSC - 380
East Los Angeles College - SFOP - 0118
Fictionalism: General IntroductionAntony Eagle Week 1, 18 January 2008the discourse have gotten started in the rst place? If the discourse isn't meant literally, what does it mean? And do possible worlds exist? Without answering any of these quest
East Los Angeles College - SFOP - 0118
Fictionalism: ModalityAntony Eagle Week 5, 15 February 20081Fictionalism about Modality and Fictionalism about Possible WorldsMathematical ctionalists begin with apparent truths like this: (1) There are n Fs i the number of Fs is n.and argue
Sveriges lantbruksuniversitet - ENSC - 805
200IEEE TRANSACTIONS ON COMMUNICATIONS, VOL. 46, NO. 2, FEBRUARY 1998Some New Twists to Problems Involving the Gaussian Probability IntegralMarvin K. Simon, Fellow, IEEE, and Dariush Divsalar, Fellow, IEEEAbstractUsing an alternate form of the
Sveriges lantbruksuniversitet - ENSC - 805
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 805 Techniques of Digital Communications Assignment 2 Due: 5 October, 20061. Differentiation With Respect to a Complex Variable (a) In the function below, z is a complex variable. Calculate
Sveriges lantbruksuniversitet - ENSC - 380
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II One-Sided Laplace Transforms1. Introduction 2. Definition of One-Sided Transform 3. The Laplace Transform on the Complex Plane 3.1 Its Origins 3.2 Inversion by Line
Sveriges lantbruksuniversitet - ENEL - 673
UNIVERSITY OF CANTERBURY Dept. of Electrical and Computer Engineering ENEL 673 Assignment 2 1. Two Seemingly-Different MMSE Solutions If you are not familiar with the singular value decomposition (SVD) of matrices, you can either skip part (b) of thi
Sveriges lantbruksuniversitet - ENSC - 805
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 805 Techniques of Digital Communications Assignment 3 All questions have equal marks. Due: 27 October, 20051. Estimating Decision Reliability Usually we calculate the average BER or SER of
Sveriges lantbruksuniversitet - ENSC - 805
3.3-13.3 Optimal Coherent Detection3.3.1 The Received Signals By linearity, the projection coefficients are the sum of signal and noise contributions, so the equivalent vector channel is r = s m + n .o Coherent reception means that the receiv
Sveriges lantbruksuniversitet - ENSC - 805
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 805 Techniques of Digital Communications Solutions to Assignment 4 1. Timing Recovery and Reality (a) In the interpolator-based tracker, the receivers clock makes its symbol rate Rr and its s
Sveriges lantbruksuniversitet - ENSC - 805
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 805 Techniques of Digital Communications Assignment 4 Due: 13 November, 2007 All questions have equal marks. Subquestion marks are as shown. 1. Linear Multiuser Detection In Section 6.1, we s
Sveriges lantbruksuniversitet - ENSC - 805
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 805 Techniques of Digital Communications Assignment 1 Due: 18 September, 2006 All questions except 1 and 2 have equal marks. Individual work. 1. Bandwidth, Centre Frequency and Noise Realizat
Sveriges lantbruksuniversitet - ENSC - 805
5.4-15.4 Joint Carrier and Timing Recovery In previous sections, when we estimated , we assumed that f o and had already been corrected, and when we estimated , we assumed that f o and had already been corrected. Now its time to look at realit
Sveriges lantbruksuniversitet - ENSC - 805
8.5-18.5 Detection of CPM One of the attractions of CPM is the wide range of ways in which it can be detected: cheap and good enough at high SNR, to Big DSP and good enough at SNR several dB lower. Well look at:o a few simple oneso and the ful
Sveriges lantbruksuniversitet - ENSC - 805
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 805 Techniques of Digital Communications Solutions to Assignment 2 10 October, 20061. Biorthogonal vs Orthogonal Well compare orthogonal and biorthogonal with the same b . Use M o for the
Sveriges lantbruksuniversitet - ENSC - 805
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 805 Techniques of Digital Communications Assignment 2 All questions have equal marks. Due: 10 October, 20061. Biorthogonal vs Orthogonal On page 4.2-30 of the lecture notes, the SER of bior
Sveriges lantbruksuniversitet - ENSC - 805
5.1-15. SYNCHRONIZATION[Proakis Chapt. 6]5.1 Formulation of the Problem We have elegant techniques for signal detection in Sections 3 and 4 but if we don't have accurate timing for sampling or low phase jitter, they degrade quickly. Synchroni
Sveriges lantbruksuniversitet - ENSC - 805
3.2-93.2.3 Channel Filter and Coloured Noise What if the transmitted signals experience filtering in the channel? How should our projection-based front end be modified?o From the receivers perspective, there is a new signal set: sm (t ) = sm (t
Sveriges lantbruksuniversitet - ENSC - 805
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 805 Techniques of Digital Communications Assignment 4 Due: 13 November, 2007 All questions have equal marks. Subquestion marks are as shown. 1. Linear Multiuser Detection In Section 6.1, we s
Sveriges lantbruksuniversitet - ENSC - 805
3.2-13.2 Acquisition of Sufficient Statistics3.2.1 Sufficient Statistics and Irrelevance The receiver picks up a continuous-time waveform r (t ) . How to reduce it to something more manageable (e.g., a scalar or vector) for detection without sac
Sveriges lantbruksuniversitet - ENSC - 805
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 805 Techniques of Digital Communications Assignment 1 All questions have equal marks. 1. Complex Random Process Due: 21 September, 2006The power spectral density Sv ( f ) of a circularly co
Sveriges lantbruksuniversitet - ENSC - 805
7.6-17.6 The Viterbi Equalizer7.6.1 The ISI Problem First, a pictorial view of the problem, to complement the mathematical descriptions of previous sections.o Matched filter output for 16 QAM with no ISI:MF Output, 16QAM5 4 3 2 1 0 1 2 3 4 5
Sveriges lantbruksuniversitet - ENSC - 805
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 805 Techniques of Digital Communications Assignment 3 All questions have equal marks. 1. Union Bound In the lecture notes we had the union bound on SER of biorthogonal signals as Ps ( M - 2)
Sveriges lantbruksuniversitet - ENSC - 805
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 805 Techniques of Digital Communications Assignment 3 Due: 29 October, 2007 All questions have equal marks. Subquestion marks are as shown. 1. Multiuser Interference We are trying to detect a
Sveriges lantbruksuniversitet - ENSC - 805
ENSC 805 Advanced Digital CommunicationsoYou should already have taken a senior undergrad course in digital communications. Preand co-requisites: o ENSC 428 or equivalent (pre-req) o ENSC 802 or equivalent (pre- or co-req) or permission of instru
Sveriges lantbruksuniversitet - ENSC - 805
3.2-73.2.2 Correlators and Matched Filters The correlators used as part of the projection operationcan also be implemented as matched filters. If h(t ) = g * (t ) , like this:then this filter operation: produces y (t ) = r ( ) h(t ) d = r
Sveriges lantbruksuniversitet - ENSC - 805
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 805 Techniques of Digital Communications Assignment 2 All questions have equal marks. Individual work. Due: 5 October, 20061. Differentiation With Respect to a Complex Variable (a) In the f
Sveriges lantbruksuniversitet - ENSC - 805
IEEE TRANSACTIONSON INFORMATIONTHEORY,VOL.IT-l?,NO.2,MARCH1971149using (64) and (66), E(x(O) = + n =z m ) exp ( 2ni $ t - 2 R, [?I F (;), (67) which is a periodic function with period T. It can be shown that x(t) - E(x(t) is expre
Sveriges lantbruksuniversitet - ENSC - 805
Sveriges lantbruksuniversitet - ENSC - 805
IEEE TRANSACTIONS ON COMMUNICATIONS, VOL. 49, NO. 10, OCTOBER 20011727Transactions Papers Convergence Behavior of Iteratively Decoded Parallel Concatenated CodesStephan ten Brink, Member, IEEEAbstractMutual information transfer characteristics
Sveriges lantbruksuniversitet - ENSC - 805
A Unied Approach to the Performance Analysis of Digital Communication over Generalized Fading ChannelsMARVIN K. SIMON, FELLOW, IEEE, AND MOHAMED-SLIM ALOUINI, STUDENT MEMBER, IEEEPresented here is a unied approach to evaluating the errorrate perfo
Sveriges lantbruksuniversitet - ENSC - 805
Sveriges lantbruksuniversitet - ENSC - 805
788IEEE TRANSACTIONS ON COMMUNICATIONS, VOL. 45, NO. 7, JULY 1997Symbol-by-Symbol MAP Demodulation of CPM and PSK Signals on Rayleigh Flat-Fading ChannelsMichael J. Gertsman, Member, IEEE, and John H. Lodge, Senior Member, IEEEAbstract- Demodul
Sveriges lantbruksuniversitet - ENSC - 805
IEEEVEHICULAR NO.TRANSACTIOKS VT-33. TECHNOLOGY, VOL. ON4, NOVEMBER 1984307Differential Detection of Gaussian MSK in a Mobile Radio EnvironmentAbstract-Minimum shift keying with Gaussian shaped transmit pulses is a strong candidate for modula
Sveriges lantbruksuniversitet - ENSC - 805
IEEE TRANSACTIONS ON COMMUNICATIONS, VOL. 49, NO. 2, FEBRUARY 2001213Transactions Letters_FIR Channel-Shortening Equalizers for MIMO ISI ChannelsNaofal Al-DhahirAbstractFinite-length delay-optimized multi-input multi-output (MIMO) equalizers
Sveriges lantbruksuniversitet - ENSC - 805
Sveriges lantbruksuniversitet - ENSC - 805
Sveriges lantbruksuniversitet - ENSC - 805
communication systems, I E E E Trans. Commun. Tcchnol., vol. COM-19, pp. 113-119, Apr. 1971. [I61 A. J. Viterbi, Phase-locked loop dynamics in the presence of noise by Fokker-Planck techniques, Proc. I E E E , vol. 51, pp. 1737-1753, Dec. 1963. [17]
Sveriges lantbruksuniversitet - ENSC - 428
SIMON FRASER UNIVERSITY School of Engineering Science ENSC 428-4 Data CommunicationsAssignment 2 1. Sum of Correlated Random VariablesDue: 2001 02 06Suppose we sum N equispaced samples of a zero-mean stationary process x(t ) with autocorrelatio
Sveriges lantbruksuniversitet - ENEL - 673
UNIVERSITY OF CANTERBURY Dept. of Electrical and Computer Engineering ENEL 673 Assignment 2 1. SNR Calculations Many analyses end up 3 dB out simply through incorrect calculation of signal and noise powers. This question is intended to give you pract
Sveriges lantbruksuniversitet - ENEL - 673
East Los Angeles College - WOLF - 0758
Some elements of the structure of the Xhosa verb1. PreliminariesWe shall propose the following schematic representation (table 1) for the Xhosa complex verbal form, in which the numerals greater than zero are taken to indicate positions or slots
East Los Angeles College - BIOL - 202
BIOL202 CASE STUDY:The Estuary of the River Mersey Introduction to the River Mersey The Estuary of the River Mersey History as revealed by environmental analysis National / International Legislation relating to the Mersey Recent Developments
Allan Hancock College - SVCAB - 2004396
2004 The Parliament of the Commonwealth of Australia THE SENATEPresented and read a first timeSenate Voters' Choice (Preference Allocation) Bill 2004 No. , 2004(Senator Brown)A Bill for an Act to amend the Commonwealth Electoral Act 1918 to e
East Los Angeles College - COMP - 213
COMP 213Advanced Object-oriented ProgrammingLecture 22Concurrency(Picasso )What is Concurrency?Concurrency is when two or more programs execute at the same time, either: sharing the same processor (multi-threading), or on different process
East Los Angeles College - COMP - 553
Tutorial 12: classesDecember 2, 2003Attendance to todays tutorial will be taken by signing the attendance sheet AND submitting (by e-mail to michele@csc.liv.ac.uk) a le containing the two classes mentioned below. You must send your e-mail message b
Allan Hancock College - UTAB - 2007250
2007LEGISLATIVE ASSEMBLY FOR THE AUSTRALIAN CAPITAL TERRITORYUNIT TITLES AMENDMENT BILL 2007EXPLANATORY STATEMENTPresented by authority Mr Andrew Barr MLA Minister for PlanningAuthorised by the ACT Parliamentary Counselalso accessible at ww
East Los Angeles College - COMP - 309
Text Compression through Huffman CodingHuffman codes represent a very effective technique for compressing data; they usually produce savings between 20% and 90%. Preliminary example. We are given a 100,000-character text (maybe a book, or a long rep
East Los Angeles College - COMP - 110
Lecture 8Case Study: Savings Account InterestLoops: recap + example. Files: abstracting from a specic devise. Streams and Tokens. Examples. Material from the second half of Holmes Chapter 4. Problem. You open a savings account with a specied
Allan Hancock College - VEAWWAWWPB - 2007807
2004-2005-2006-2007THE PARLIAMENT OF THE COMMONWEALTH OF AUSTRALIAHOUSE OF REPRESENTATIVESVETERANS ENTITLEMENTS AMENDMENT (DISABILITY, WAR WIDOW AND WAR WIDOWER PENSIONS) BILL 2007EXPLANATORY MEMORANDUM(Circulated by authority of the Minist
East Los Angeles College - DPLB - 0149
Drug delivery systems based on sugar-macromolecule conjugates Benjamin G Davis* & Mark A RobinsonAddress Dyson Perrins Laboratory Department of Chemistry University of Oxford South Parks Road Oxford OX1 3QY UK Email: Ben.Davis@chem.ox.ac.uk *To whom
East Los Angeles College - OXMA - 0006
Performance Based Compensation and Direct Earnings ManagementAntnio Cmara o aVicky HendersonOklahoma State UniversityWarwick Business SchoolJanuary 12, 2008We would like to thank Yener Kandogan, Seyed Mehdian, Keith Moreland, Nandu
Drake - ECON - 001
Principles of Macroeconomics (Econ 001) Drake University, Fall 2000 William M. BoalSignature: Printed name: ID number:MIDTERM EXAMINATION #1: Version A September 18, 2000INSTRUCTIONS: This exam is closed-book, closed-notes, but calculators are p
East Los Angeles College - DPLB - 0149
Chemical Biology, Option IIMechanism Based Proteomic TaggingCase History CH1Proteome Wide Screening of Serine Protease Activity Proc Natl Acad Sci 1999, 97, 14694; Proteomics 2001, 1, 1067; Proc Natl Acad Sci 2002, 99, 10335; Biochemistry 2001,
East Los Angeles College - SANT - 1114
The History of Egyptian CinemaHilary Term 2008 St. Antony's College Hilda Bess Building, The New Room 2-4pm, weeks 1 and 3-7Walter Armbrust (http:/users.ox.ac.uk/%7Esant1114/Armbrust.html)Film Documents Website: http:/users.ox.ac.uk/%7Esant1114/E
East Los Angeles College - DPLB - 0149
Molecular BioSystemsThis article was published as part of the2008 Emerging Investigators IssueHighlighting the work of outstanding young scientists at the chemical- and systems-biology interfaces Please take a look at the full table of contents t
East Los Angeles College - SSFC - 0041
"We, the Peoples of Europe ."Kalypso Nicoladisa more perfect union Political tremors are shaking the old continent. As the European Union's enlargement brings most of the continent under the same banner, Europeans, like their American cousins two c
Sveriges lantbruksuniversitet - LIB - 2240
Is There a January Effect in the Greater China Area?by Chung-Wen (Sally) Hsu M.Sc., Civil Engineering, National Central University, TaiwanPROJECT SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF BUSINESS ADMINISTRAT
Sveriges lantbruksuniversitet - IR - 3552
ANNUAL REPORT20 03C E L E B R AT I N G 3 0 Y E A R S O F S E R V I C EAs of April 1, 2004 the provincial government is reducing its grant to the Native Courtworkers and Counselling Association of BC from $2.04 million to $1.3 millionBy a lett
Sveriges lantbruksuniversitet - LIB - 724
Does the Barro-Gordon Model Explain the Behaviorof Inflation in Canada?Shan ChenBachelor of Business Administration, Nankai University, 2003PROJECT SUBMITTED IN PARTIAL FULFIJLLMENT OF THE REQUIREMENTS FOR THE DEGREE OFMASTER OF ARTSIn the