Seismic geometric attribute analysis for fracture characterization: New methodologies and applications

In 3D subsurface exploration, detection of faults and fractures from 3D seismic data is vital to robust structural and stratigraphic analysis in the subsurface, and great efforts have been made in the development and application of various seismic attributes (e.g. coherence, semblance, curvature, and flexure). However, the existing algorithms and workflows are not accurate and efficient enough for robust fracture detection, especially in naturally fractured reservoirs with complicated structural geometry and fracture network. My Ph.D. research is proposing the following scopes of work to enhance our capability and to help improve the resolution on fracture characterization and prediction.;For discontinuity attribute, previous methods have difficulty highlighting subtle discontinuities from seismic data in cases where the local amplitude variation is non-zero mean. This study proposes implementing a gray-level transformation and the Canny edge detector for improved imaging of discontinuities. Specifically, the new process transforms seismic signals to be zero mean and helps amplify subtle discontinuities, leading to an enhanced visualization for structural and stratigraphic details. Applications to various 3D seismic datasets demonstrate that the new algorithm is superior to previous discontinuity-detection methods. Integrating both discontinuity magnitude and discontinuity azimuth helps better define channels, faults and fractures, than the traditional similarity, amplitude gradient and semblance attributes.;For flexure attribute, the existing algorithm is computationally intensive and limited by the lateral resolution for steeply-dipping formations. This study proposes a new and robust volume-based algorithm that evaluate flexure attribute more accurately and effectively. The algorithms first volumetrically fit a cubic surface by using a diamond 13-node grid cell to seismic data, and then compute flexure using the spatial derivatives of the built surface. To avoid introducing interpreter bias, this study introduces a new workflow for automatically building surfaces that best represent the geometry of seismic reflections. A dip-steering approach based on 3D complex seismic trace analysis is implemented to enhance the accuracy of surface construction and to reduce computational time. Applications to two 3D seismic surveys demonstrate the accuracy and efficiency of the new flexure algorithm for characterizing faults and fractures in fractured reservoirs.;For robust fracture detection, this study presents a new methodology to compute both magnitude and directions of most extreme flexure attribute. The new method first computes azimuthal flexure; and then implements a discrete azimuth-scanning approach to finding the magnitude and azimuth of most extreme flexure. Specially, a set of flexure values is estimated and compared by substituting all possible azimuths between 0 degree (Inline) and 180 degree (Crossline) into the newly-developed equation for computing azimuthal flexure. The added value of the new algorithm is demonstrated through applications to the seismic data set from Teapot Dome of Wyoming. The results indicate that most extreme flexure and its associated azimuthal directions help reveal structural complexities that are not discernible from conventional coherence or geometric attributes.;Given that the azimuth-scanning approach for computing maximum/minimum flexure is time-consuming, this study proposes fracture detection using most positive/negative flexures; since for gently-dipping structures, most positive is similar to maximum flexure while most negative flexure to minimum flexure. After setting the first reflection derivatives (or apparent dips) to be zero, the localized reflection is rotated to be horizontal and thereby the equation for computing azimuthal flexure is significantly simplified, from which a new analytical approach is proposed for computing most positive/negative flexures. Comparisons demonstrate that positive/negative flexures can provide quantitative fracture characterization similar to most extreme flexure, but the computation is 8 times faster than the azimuth-scanning approach.;Due to the overestimate by using most positive/negative flexure for fracture characterization, 3D surface rotation is then introduced for flexure extraction in the presence of structural dip, which makes it possible for solving the problem in an analytical manner. The improved computational efficiency and accuracy is demonstrated by both synthetic testing and applications to real 3D seismic datasets, compared to the existing discrete azimuth-scanning approach.;Last but not the least, strain analysis is also important for understanding structural deformation, predicting natural fracture system, and planning well bores. Physically, open fractures are most likely to develop in extensional domains whereas closed fractures in compressional ones. The beam model has been proposed for describing the strain distribution within a geological formation with a certain thickness, in which, however, the extensional zone cannot be distinguished from the compression one with the aid of traditional geometric attributes, including discontinuity, dip, and curvature. To resolve this problem, this study proposes a new algorithm for strain reconstruction using apparent dips at each sample location within a seismic cube

Model driven segmentation and the detection of bone fractures

Bibliography: leaves 83-90.The introduction of lower dosage image acquisition devices and the increase in computational power means that there is an increased focus on producing diagnostic aids for the medical trauma environment. The focus of this research is to explore whether geometric criteria can be used to detect bone fractures from Computed Tomography data. Conventional image processing of CT data is aimed at the production of simple iso-surfaces for surgical planning or diagnosis - such methods are not suitable for the automated detection of fractures. Our hypothesis is that through a model-based technique a triangulated surface representing the bone can be speedily and accurately produced. And, that there is sufficient structural information present that by examining the geometric structure of this representation we can accurately detect bone fractures. In this dissertation we describe the algorithms and framework that we built to facilitate the detection of bone fractures and evaluate the validity of our approach

Computer-aided diagnosis of complications of total hip replacement X-ray images

Hip replacement surgery has experienced a dramatic evolution in recent years supported by the latest developments in many areas of technology and surgical procedures. Unfortunately complications that follow hip replacement surgery remains the most challenging dilemma faced both by the patients and medical experts. The thesis presents a novel approach to segment the prosthesis of a THR surgical process by using an Active Contour Model (ACM) that is initiated via an automatically detected seed point within the enarthrosis region of the prosthesis. The circular area is detected via the use of a Fast, Randomized Circle Detection Algorithm. Experimental results are provided to compare the performance of the proposed ACM based approach to popular thresholding based approaches. Further an approach to automatically detect the Obturator Foramen using an ACM approach is also presented. Based on analysis of how medical experts carry out the detection of loosening and subsidence of a prosthesis and the presence of infections around the prosthesis area, this thesis presents novel computational analysis concepts to identify the key feature points of the prosthesis that are required to detect all of the above three types of complications. Initially key points along the prosthesis boundary are determined by measuring the curvature on the surface of the prosthesis. By traversing the edge pixels, starting from one end of the boundary of a detected prosthesis, the curvature values are determined and effectively used to determine key points of the prosthesis surface and their relative positioning. After the key-points are detected, pixel value gradients across the boundary of the prosthesis are determined along the boundary of the prosthesis to determine the presence of subsidence, loosening and infections. Experimental results and analysis are presented to show that the presence of subsidence is determined by the identification of dark pixels around the convex bend closest to the stem area of the prosthesis and away from it. The presence of loosening is determined by the additional presence of dark regions just outside the two straight line edges of the stem area of the prosthesis. The presence of infections is represented by the determination of dark areas around the tip of the stem of the prosthesis. All three complications are thus determined by a single process where the detailed analysis defer. The experimental results presented show the effectiveness of all proposed approaches which are also compared and validated against the ground truth recorded manually with expert user input

Shape/image registration for medical imaging : novel algorithms and applications.

This dissertation looks at two different categories of the registration approaches: Shape registration, and Image registration. It also considers the applications of these approaches into the medical imaging field. Shape registration is an important problem in computer vision, computer graphics and medical imaging. It has been handled in different manners in many applications like shapebased segmentation, shape recognition, and tracking. Image registration is the process of overlaying two or more images of the same scene taken at different times, from different viewpoints, and/or by different sensors. Many image processing applications like remote sensing, fusion of medical images, and computer-aided surgery need image registration. This study deals with two different applications in the field of medical image analysis. The first one is related to shape-based segmentation of the human vertebral bodies (VBs). The vertebra consists of the VB, spinous, and other anatomical regions. Spinous pedicles, and ribs should not be included in the bone mineral density (BMD) measurements. The VB segmentation is not an easy task since the ribs have similar gray level information. This dissertation investigates two different segmentation approaches. Both of them are obeying the variational shape-based segmentation frameworks. The first approach deals with two dimensional (2D) case. This segmentation approach starts with obtaining the initial segmentation using the intensity/spatial interaction models. Then, shape model is registered to the image domain. Finally, the optimal segmentation is obtained using the optimization of an energy functional which integrating the shape model with the intensity information. The second one is a 3D simultaneous segmentation and registration approach. The information of the intensity is handled by embedding a Willmore flow into the level set segmentation framework. Then the shape variations are estimated using a new distance probabilistic model. The experimental results show that the segmentation accuracy of the framework are much higher than other alternatives. Applications on BMD measurements of vertebral body are given to illustrate the accuracy of the proposed segmentation approach. The second application is related to the field of computer-aided surgery, specifically on ankle fusion surgery. The long-term goal of this work is to apply this technique to ankle fusion surgery to determine the proper size and orientation of the screws that are used for fusing the bones together. In addition, we try to localize the best bone region to fix these screws. To achieve these goals, the 2D-3D registration is introduced. The role of 2D-3D registration is to enhance the quality of the surgical procedure in terms of time and accuracy, and would greatly reduce the need for repeated surgeries; thus, saving the patients time, expense, and trauma

FracDetect: A novel algorithm for 3D fracture detection in digital fractured rocks

Fractures have a governing effect on the physical properties of fractured rocks, such as permeability. Accurate representation of 3D fractures is, therefore, required for precise analysis of digital fractured rocks. However, conventional segmentation methods fail to detect and label the fractures with aperture sizes near or below the resolution of 3D micro-computed tomographic (micro-CT) images, which are visible in the greyscale images, and where greyscale intensity convolution between different phases exists. In addition, conventional methods are highly subjective to user interpretation. Herein, a novel algorithm for the automatic detection of fractures from greyscale 3D micro-CT images is proposed. The algorithm involves a low-level early vision stage, which identifies potential fractures, followed by a high-level interpretative stage, which enforces planar continuity to reject false positives and more reliably extract planar fractures from digital rock images. A manually segmented fractured shale sample was used as the groundtruth, with which the efficacy of the algorithm in 3D fracture detection was validated. Following this, the proposed and conventional methods were applied to detect fractures in digital fractured coal and shale samples. Based on these analyses, the impact of fracture detection accuracy on the analysis of fractured rocks' physical properties was inferred

Image analysis in medical imaging: recent advances in selected examples

Medical imaging has developed into one of the most important fields within scientific imaging due to the rapid and continuing progress in computerised medical image visualisation and advances in analysis methods and computer-aided diagnosis. Several research applications are selected to illustrate the advances in image analysis algorithms and visualisation. Recent results, including previously unpublished data, are presented to illustrate the challenges and ongoing developments

Surface-Only Simulation of Fluids

Surface-only simulation methods for fluid dynamics are those that perform computation only on a surface representation, without relying on any volumetric discretization. Such methods have superior asymptotic complexity in time and memory than the traditional volumetric discretization approaches, and thus are more tractable for simulation of complex fluid phenomena. Although for most computer graphics applications and many engineering applications, the interior flow inside the fluid phases is typically not of interest, the vast majority of existing numerical techniques still rely on discretization of the volumetric domain. My research first tackles the mesh-based surface tracking problem in the multimaterial setting, and then proposes surface-only simulation solutions for two scenarios: the soap-films and bubbles, and the general 3D liquids. Throughout these simulation approaches, all computation takes place on the surface, and volumetric discretization is entirely eliminated