C T& F Vol . 9 Num . 2 D e c emb er 2 01 9 37
SUBMARINE CANYON-FILL RECONSTRUCTION FROM INTEGRATED SEISMIC-STRATIGRAPHIC ANALYSIS – APPLICATION TO BANQUEREAU FORMATION, SCOTIAN BASIN – OFFSHORE CANADA.Chajin, Patricia a*; Illidge, Erick-J.b; Khurama, Sait b
a Grupo de Investigación en Inversión y Modelado Geofísico IMGeo, Ecopetrol – Instituto Colombiano del Petróleo km7 via Piedecuesta –Bucaramanga, C.P 681011, Piedecuesta,
Colombia . b Universidad Industrial de Santander, Carrera 27 calle 9, Bucaramanga, Colombia
*email: [email protected].
ABSTRACT Building geological models (integrating stratigraphic, structural and paleo-environmental 3D models) that allow the interpretation of sand bodies deposited by turbidity currents along submarine canyons or channels, is one of the most useful tools used by geoscientists for the definition of new drilling opportunities in both exploration and development phases. In this context, the integration of methodologies such as sequence stratigraphy and seismic attributes, together with well-log and core information, outline the basis for the interpretation of sand-body lithostratigraphy and chronostratigraphy. Similarly, these models allow the interpreter to reconstruct the depositional environment and deformation history of a sedimentary basin [1]. Based on a series of chronostratigraphic stages, this paper proposes a 3D model for the sedimentation history of the Banquereau Formation.
This model is based on the integration of seismic stratigraphy, seismic attribute interpretation and well-log analysis. Also, a set of system tracts and corresponding transgression and regression phases were identified for the sedimentary interval of interest. The available dataset provided the information to identify the geometry and changes in the sedimentation patterns of the stratigraphic sequences from the Tertiary to the present, thus defining a 3D model of the sedimentological and structural architecture of this interval. Last but not least, the resulting 3D stratigraphic model made possible the identification and description of an amalgamated channel complex filling a submarine canyon associated with a fluvio-deltaic setting. This sort of analysis might be used as an analog for similar reservoirs, providing key insights and vital information for decision making.
KEYWORDS / PALABRAS CLAVE AFFILIATION
RECONSTRUCCIÓN DEL RELLENO DE UN CAÑÓN SUBMARINO A PARTIR DE ANÁLISIS SISMOESTRATIGRÁFICO INTEGRADO – APLICACIÓN A LA FORMACIÓN BANQUERAU, CUENCA ESCOCIA – COSTA AFUERA CANADÁ
RESUMENLa construcción de modelos geológicos que permitan la interpretación de cuerpos de arena depositados por corrientes de turbidez a lo largo de cañones o canales submarinos, es una de las herramientas más utilizadas por los geocientíficos para la definición de nuevas oportunidades de perforación en las fases de exploración y desarrollo. En este contexto, la integración de metodologías como estratigrafía de secuencias y atributos sísmicos, en conjunto con información de registros de pozo y de núcleos, esbozan las bases para la interpretación litoestratigráfica y cronoestratigráfica de cuerpos de arena.De la misma manera, estos modelos permiten al intérprete reconstruir el ambiente de depósito y la historia de deformación de una Cuenca sedimentaria. Basado en una serie de etapas cronoestratigráficas, este artículo propone un modelo 3D para la historia de sedimentación de la Formación Banquerau. Este modelo se basa en la integración de estratigrafía sísmica, la interpretación de atributos sísmicos y el análisis de registros de pozo.También, se identificó un conjunto de system tracts y sus correspondientes fases de transgresión y regresión para el intervalo sedimentario de interés. Los datos disponibles proporcionaron la información para identificar la geometría y los cambios en los patrones de sedimentación de las secuencias estratigráficas desde el Terciario hasta el presente, definiendo así un modelo 3D de la arquitectura sedimentológica y estructural de este intervalo.Por último, pero no menos importante, el modelo estratigráfico 3D resultante hizo posible la identificación y la descripción de un un complejo de canales amalgamados rellenando un cañón submarino asociado con un entorno fluvio-deltáico. Este tipo de análisis podría se usado como un análogo para reservorios similares, proporcionando información clave y vital para la toma de decisiones.
Seismic Attributes | Sequence Stratigraphy | Well Logs | Submarine Canyon.Atributos sismicos | Secuencia Estatigráfica |Well Logs | Cañón Submarino.
A R T I C L E I N F O : Received : March 16, 2017Revised : April 05, 2019Accepted : May 20, 2019CT&F - Ciencia, Tecnologia y Futuro Vol 9, Num 2 December 2019. pages 37 - 52DOI : https://doi.org/10.29047/01225383.119
C T& F Vol . 9 Num . 2 D e c emb er 2 01 9 39
Ec op e t r o lVol . 9 Num . 2 D e c emb er 2 01 9
38 Ec op e t r o l
The study area is located in Offshore East Canada, more specifically in the Nova Scotia Basin (Figure 1). Now, extended works have improved understanding of the Nova Scotia basin since the 1960s [2] - [3], in order to find the most profitable plays. In this context, approximately 180 exploration wells have been drilled on the shelf of the Scotian basin [3], as well an another eight wells on the slope [4]. Most of the discoveries made on the Scotian Shelf relate mainly to gas or condensate [3]. These discoveries have been principally encountered in the Cretaceous and Jurassic Formations in the Scotian Basin [5]. The main target of this paper is the Banquereau Formation, which is a Late Cretaceous to early Quaternary stacked series of prograding sequences that downlap unconformably onto the Wyandot Formation [6].
The geology and stratigraphy of the Scotian Shelf and the Sable Sub-basin have been described in several papers [6]-[9], but due to the lack of prospectivety for hydrocarbons, little attention has been paid to the structure of the Banquereau Formation, its implications for paleogeography and relative sea level [10]. McIver [6] described the main lithologies of the Banquereau Formation as mainly mudstone with few intervals of silt and argillaceous sandstone, which is found in increasing amounts in the upper half of the Formation. In this context, McIver [6] also stated that the mudstone of the Banquereau Formation is commonly burrowed or highly reworked. On the other
hand, the sandstone of the Banquereau Formation is very fine to fine grained and is finely laminated with channel features identifiable in the upper portion of the Formation [6]. McIver [6] also described the lithology of this Formation in detail, but, due to poor well control at the time of publication, this paper does not mention the internal geometry of the Banquereau Formation. Fensome et al. [8] describes the lithology of the Banquereau Formation as sandstone deposited at the topsets of clinoforms, with mudstone fore- and -bottomsets. Fensome et al. [8] also states that the lower interval of the Banquereau Formation is comprised of "prograding deltaic clinoform units that built basinward" [10].
While most of the previous work done on the Banquereau Formation was based on limited information (well-logs separated from seismic data), the scope of this paper is to show how the integration of these sorts of information through different advanced geological and geophysical data interpretation techniques led to the generation of a robust sequence stratigraphy model that accounts for the characterization of the submarine channels and valleys found in the Penobscot Block (Figure 1), more specifically in the Banquereau prograding Tertiary delta complex.
INTRODUCTION1 2. THEORETICAL FRAME All the reported Formations of the Sable Sub-basin are shown on the stratigraphic column in Figure 2. In chronological order, the Formations are Abenaki, Mississauga, Logan Canyon, Dawson Canyon, Wyandot and Banquereau. Abenaki Formation is Jurassic in age, while Mississauga, Logan Canyon, Dawson Canyon, and Wyandot Formations are all from the Cretaceous. Banquereau Formation is Late Cretaceous to early Quaternary in age [13]. Mississauga Formation is mainly sandstone with some minor shale beds. This Formation is commonly broken into an Upper and Lower unit, with the "Base O-Marker". The member Base O-Marker is a regional stratigraphically correlatable unit, which is easy to identify on seismic data and comprises mainly sandstone, minor carbonates and shale [14]. Mississauga Formation can broadly be defined as comprising fluvio-deltaic sediments [15]. Logan Canyon Formation has several members with variable percentages of sandstone and shale, but the main component is a thick shale succession, interpreted as a marine transgression above the Mississauga Formation [15]. Dawson Canyon Formation is mainly shale, while Wyandot Formation is chalk-dominated, highly bioturbated, with minor amounts of marl and shale [13]. In fact, paleontological evidence indicates there is a hiatus associated with the top of the Wyandot Formation [6] above which lies the Banquereau Formation. As mentioned before, Banquereau Formation is mudstone dominated, with minor amounts of siltstone and sandstone in a prograding series of clinoforms [6], [8].
According to Wach and Mosher [16], Cretaceous sediments and the entire Cenozoic sedimentary succession above the Wyandot are designated the Banquereau Formation [6]. Marine shelf mudstones, sandstones and conglomerates of the Banquereau Formation were
Figure 1. Location of the study area. (Modified from Google Earth) Figure 2. Generalized stratigraphic column for the Sable Sub-basin (Modified and updated from [11] and [12]).
influenced throughout the Cenozoic by several major unconformities related to sea level fall [16]. Unconformities are noted during the Paleocene, Oligocene and Miocene intervals where fluvial and deep-water currents eroded largely unconsolidated sediments, subsequently depositing them on the abyssal plain [7],[17]. Winnowing and reworking of deep-water sediment by bottom currents began in the Oligocene [18],[19], providing the earliest evidence of thermohaline circulation. Sediment distribution of Miocene successions were strongly influenced by the Western Boundary Undercurrent with periods of intensified bottom current activity also occurring in the Late Pliocene [19], [20], followed by widespread gully cutting in the Early Pleistocene. During the Quaternary to recent, several hundred meters of glacial and marine sediment were deposited on the outer shelf and slope [7],[17]-[22]. Wach and Mosher [16] also argues that on the Scotian Shelf and Grand Banks, the widespread hiatus eroding either the upper part or all of the Oligocene is marked by a regional unconformity; the nature of which includes canyon formation. Canyon incision at the shelf edge initiated during the Eocene and was extensive by the Oligocene [8]. Pe-Piper and Piper [23] noted that although Oligocene strata are absent on the Scotian Shelf and Grand Banks, they are present both on the Labrador Shelf [24] and on the New Jersey margin where small hiatuses are correlated with the global eustatic 1987 sea-level curve of Haq [25], [26]-[28]. Wade and MacLean [29] attributed the missing Oligocene strata on the eastern Scotian Shelf to a broad southeasterly- trending canyon associated with a sea-level Lowstand, although Grist and Zentilli [30] argued for thermal inversion in the Late Cretaceous to Early Cenozoic, based on results from apatite fission track analysis.
GeologicalTime Scale
GlobalSea Level
(Mille, 2005)
Long-termStratigraphicSequences
ShelfDepositional
settings
Slope and RiseDepositional
settings NW SELithostratigraphyMain
Events Salt Tectonics SeismicHorizons
MajorAnoxic events
Time scale from Ogg et al. (2004).
Eustatic curves: short-term from Miller (2005); intermeditaeterm, mean sea level from SEPM-Haq´80 synthetic.
Stratigraphy adpted from Weston et al. (2012)
Regional Tectonicand
Geological events
C T& F Vol . 9 Num . 2 D e c emb er 2 01 9 41
Ec op e t r o lVol . 9 Num . 2 D e c emb er 2 01 9
40 Ec op e t r o l
3. EXPERIMENTAL DEVELOPMENT In this paper, a stratigraphic sequence model, as defined by Catuneanu, Hancox, Cairncross, and Rubidge [31] (Figure 3), will be used to generate all the sequence stratigraphic analysis and to build a system tracts model.
Figure 3. Timing of system tracts and sequence boundaries for the sequence models currently in use. (Modified from [31]).
Figure 4. 3D seismic survey and correlation wells used in the application.
Figure 5. Workflow applied to transform petrophysical and geophysical information into a 3D geological model.
Figure 6. Well logs available in the interval of analysis for one of the correlation wells.
Sequencemodel
Events
End oftransgression
DepositionalSequence II
DepositionalSequence III
DepositionalSequence IV
GeneticSequence
T-RSequence
HST
TST TST
LSTlate LST(wedge)
early LST(fan)
LST
TST TST TST
RST
RST
HST
HST
FSSTlate HST(fan)
late LST(wedge)
early LST(fan)
early HST(wedge)HST
sequence boundarysystems tract boundarywithin systems tract surface
onset ofbase-level fall
end ofregression
end oftransgression
end ofbase-level fall
Time
HST
early HST HST
End ofregression
End ofbasw-level
fall
onset ofbasw-level
fall
The principal data of this study is a high-quality 3D PSTM from the Penobscot block, located in Nova Scotia Offshore of Canada (Figure 4). The dataset consists of 2 wells with a set of well logs, time to depth tables for the seismic well tie process and the seismic volume mentioned previously (Figure 4). Figure 5 shows the workflow applied to the dataset to obtain the final model. All steps and partial results are illustrated and briefly described in the following sections. It is worth mentioning that this study was conducted in a non-prospective interval. Nevertheless, this methodology may be applied to seismic surveys that image similar geological features. All dataset interpretation was performed in Petrel & Opendtect software packages.
WELL-LOG INTERPRETATION
Following the workflow illustrated in Figure 5, data processing began with well log data interpretation. Figure 6 illustrates the available well logs for one of the correlation wells (Gamma ray [GR], P-Sonic [DTp] and Deep resistivity [ILD]). The first step in the well log interpretation was the definition of the Vshale model, following Asquith and Krygowski [32], using the Gamma ray [GR] and the Equation 1. In this context, Figure 8 illustrates the resulting Vshale model. The next step involved the interpretation of the facies model based on both the Vshale model defined in this study and neural networks as proposed by Bhatt and Helle [33]. By using neural networks, one can associate as many well log properties as there are available, thus finding clusters of points in a crossplot that share similar ranges of the rock properties used as input data. Therefore, one of the clusters defined during the application of neural networks in classification mode would represent a specific lithology with a given range of porosity, shale fraction, density, P-wave etc.
DSMF_PenobscotPSTM
3000
2500
2000
1500
1000
500
-500
-1000
-1500
-2000
-2500-3000
0
System TractsInterpretation
3D GeologicalModel
Boundary SurfacesInterpretation
Data InventoryLithologic
InterpretationElectrofaciesInterpretation
DepositionalEnvironment
1D Modeling
Seismic FaciesInterpretation
Seismic ImageEnhancing Seismic-Well Tie
3D Modeling
MD C1_in
BS_in4.00 14.00
Banquereau
Upper
Litho_Tops GR_API DTp_usec_ftus/ft 240 0.1
ILD_ohm_mohm.m 200.0gAPI10 100
Internediate-2
Internediate-1
Basal
4.00 14.00
in
in
1:13576
600
800
1000
1200
1400
1600
1800
2000
2200
2400
2600
28487
(1)=−−
C T& F Vol . 9 Num . 2 D e c emb er 2 01 9 43
Ec op e t r o lVol . 9 Num . 2 D e c emb er 2 01 9
42 Ec op e t r o l
Additionally, the classification of different types of claystone may also be performed by applying neural networks. Figure 7 shows the resulting crossplot of the facies modeling where 3 different lithologies were interpreted as sandstone, claystone and shale. Figure 7 shows the 1D facies model in the correlation wells along with the Vshale model proposed.
Once the 1D facies and Vshale model were generated the electrofacies were interpreted, enabling the definition of the sedimentary environments per interval, which correlates with the
200
180
160
140
120
100
80 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
80100
120140
160180
200 200
180
160
140
120
100
80 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
80100
120140
160180
200
P - S
onic
,[us/
ft]
P - S
onic
,[us/
ft]
Vshale,[%] (a) (b)
Vshale,[%]
Z-values: Facies_T
SandstoneShale
LimestoneClaystone Marystone
Calcareous_SsShaly SandstoneCalcareous Siltstone
Calcareous_ShaleCalcareous_Claystone
Figure 7. Final crossplot for the facies modeling obtained by applying neural networks and Vshale. Crossplot (a) shows the range of values of Vshale and P-sonic and crossplot (b) illustrates the facies classification for Banquereau Formation.
Figure 8. 1D Facies model for the correlation wells obtained by applying neural networks and Vshale
Figure 9. Synthetic seismogram resulting from the seismic-well tie process on one of the correlation wells.TVD
1:11339
8001000
1:11339TVD Ltho_Tops Facies_T Vsh_T
0.00 1.00%
L-30[TVD]
1200
1400
1600
1800
2000
2200
2400
2600
2800
3000
3131.5
Upper
Upper
Litho_Tops Facies_T Vsh_T Electrofacies3240 m
Irregular
Irregular
Electrofacies
Cylinder
Cylinder
Coarsening UpwardCoarsening Upward
Fining Upward
Fining Upward
B-41[TVD]
0.00 % 1.00
Intermediate-2
Intermediate-2
Intermediate-1Intermediate-1
Basal
Basal
WyanDot
Banquereau
Banquereau
1000
1200
1400
1600
1800
2000
2200
2400
2400
2800
2953
Time (ms) Tops1000 1.5 3 -100 100
Along well path
10021271 1275 1279 1283 1287 1291 1295 1299 1303
1005 1003 1012 1015 1018 1021 1024 102706000*Vp_Final *RHO_Final
m/s g/cc
wave0 : wave0 - x: Time Y: Ar
200
300
Banq
uere
auFo
rmat
ion
WyanDot
DawsonCanyon
400
500
600
700
800
900
1000
Tops P-wave Density DSMF_Penobscot
facies model defined herein (Figure 8). Now, with the sedimentary environments interpreted, the next step was to carry all the information from the 1D models to the 3D scenario. In order to achieve that, seismic-well tie was performed in the correlation wells using the P-Sonic [DTp] and Density [RHOB] logs available. Figure 9 shows the resulting synthetic seismogram from which the depth-time model was obtained to take the 1D stratigraphic information generated (facies, electrofacies and sedimentary environments) to the 3D seismic information
SEISMIC IMAGE ENHANCEMENT & INPUT DATA FOR SEQUENCE STRATIGRAPHIC INTERPRETATION
After tying the correlation wells to the seismic, the interpretation of the 3D seismic data began with enhancement of the seismic image. This process included reflector relative dip estimation and the enhancement of reflector continuity. Based on the reflector relative dip attribute known as Dip Steering (which was calculated using Opendtect) an amplitude filter (Dip Steered Median Filter, DSMF) was applied to the original PSTM image so that the reflector continuity could be enhanced [34] (Figure 10) and thus facilitate the manual interpretation of stratigraphic surfaces. In this context, using the methodology proposed by Odoh, Ilecgukwu and Okoli [34], a fault enhancement filter (FEF) was also applied so that the fault surfaces could be easily interpreted using geometric seismic attributes based on the result of this fault filter (Figure 11). Following the manual interpretation of the main stratigraphic surfaces and faults in the entire seismic volume, a set of about 240 horizons were automatically tracked in the interval of analysis using the main unconformity horizons and faults as constraints. These horizons neatly follow the reflectors and therefore represent timelines that will be used in the following steps. It is worth mentioning that this set of horizons is truncated if they approach each other up to a certain distance defined by the interpreter, therefore making it even easier to interpret stratigraphic features such as sub-aerial unconformities.
SEQUENCE STRATIGRAPHIC INTERPRETATION
Based on the horizons framework, the corresponding reflector terminations and seismic facies interpretation (Figure 13), a complete set of stratigraphic surfaces was defined that was later used to interpret the systems tracts (Figure 14). Using the
This 3D horizons framework, combined with attributes such as thickness variations, make it possible to define time attributes for the bounding surfaces.
Using both the enhanced seismic data and the 3D horizons framework it is possible to interpret the paleo-geomorphology of the entire seismic volume and identify stratigraphic features such as erosional events, non-depositional hiatuses [35], submarine channels, debris flows and most importantly turbidites. Figure 12 illustrates the enhanced seismic image and the corresponding 3D horizons framework for an Inline of the seismic survey. The prograding delta deposited during the late Cretaceous and Tertiary is incised by recent submarine channels and valleys, which is adequately represented by the horizons framework generated in this study.
4. RESULT
C T& F Vol . 9 Num . 2 D e c emb er 2 01 9 45
Ec op e t r o lVol . 9 Num . 2 D e c emb er 2 01 9
44 Ec op e t r o l
Xline 1010ILXL
Random lineRandom line
Xline 1240Xline 1180 Xline 1230
SW NE1220 1220 1220 1220 1220
3000
2000
1000
-1000
-2000
-3000
0
3000
2000
1000
-1000
-2000
-3000
0
1220
-100
-200
-300
-400
-500
-600
-700
-800
-900
-1000
-1100
-1200
-100
-200
-300
-400
-500
-600
-700
-800
-900
-1000
-1100
-1200
[ms]
Original PSTM PSTM + DSMF
1002 1082 1162
0 1000 2000 3000 4000 5000m 0 1000 2000 3000 4000 5000m
1241 1321 1401
Xline 1010ILXL
Random lineRandom line
Xline 1240Xline 1180 Xline 1230
SW NE1220 1220 1220 1220 12201220
1002 1082 1162 1241 1321 1401
[ms]
Xline 1010ILXL
Random lineRandom line
Xline 1240Xline 1180 Xline 1230
SW NE1220 1220 1220 1220 12201220
1002 1082 1162 1241 1321 1401
-100
-200
-300
-400
-500
-600
-700
-800
-900
-1000
-1100
-1200
3000
2000
1000
-1000
-2000
-3000
0
3000
2000
1000
-1000
-2000
-3000
0
-100
-200
-300
-400
-500
-600
-700
-800
-900
-1000
-1100
-1200PSTM + DSMF + FEFPSTM + DSMF
0 1000 2000 3000 4000 5000m 0 1000 2000 3000 4000 5000m
Xline 1010ILXL
Random lineRandom line
Xline 1240Xline 1180 Xline 1230
SW NE1220 1220 1220 1220 12201220
1002 1082 1162 1241 1321 1401
521001 1100
(a) (b)
1200 1300 1400 1478-2000
250[m]
Time (ms) Cross-line
2000
250
500
750
1000
1200
521001 1100 1200 1300 1400 1478
-20002401
Submarine Valley
Time (ms)
2000
250
500
750
1000
1200 Cross-line
Figure 10. Comparison between the original PSTM image (left) and the enhanced PSTM image (right) using the Dip Steered Median Filter (DSMF).
Figure 11. Comparison between the DSMF image (left) and the FEF (right) using the Dip Steered Median Filter (DSMF).
Figure 12. Enhanced seismic image (a) and its corresponding horizons framework (b) defined from the Dip-Steering attribute. Notice the submarine channel eroding the prograding delta.
Figure 13. 1D Facies model, electrofacies and seismic facies defined for one of the correlation wells.
Figure 14. (a) Arbitrary seismic line extracted from the seismic volume showing the correlation wells. (b) Stratigraphic surfaces associated to reflector terminations.
stratigraphic surfaces and the reflector terminations defined in Figure 14, the seismic data was transformed to Wheeler domain where erosion events and non-deposition hiatus were interpreted (Figure 15). The wheeler diagram, together with the well-log interpretation, made it possible to interpret the system tracts associated with each interval and to propose a base level curve for the study area (Figure 16).
Following the interpretation of the base level curve, system tracts were interpreted in the entire 3D seismic survey, thus creating a 3D sequence stratigraphic model calibrated with well log information (Figure 17). In Figure 14 (b), from bottom to top, it can be seen that
MD Litho_Tops Vsh_T Electrofacies Seismic_Facies Seismic_Facies_Interpretation
Semi-parallelreflector withmoderate amplitudeand moderatecontinuity
Semi-parallel reflectorwith moderate to highamplitude and moderateto high continuity
Predominantly claystonewith some occasionalsiltstone and shale
UnconsolidatedSandstone interbeddedwith thin layers ofclaystone
Fluvial sedimentation withpredominance of dandstoneprobably related withchannels or bars
Semi-parallelreflectors with lowamplitude and lowcontinuity
Predominantly shalegrading upwards tosiltstone and sandstone
Delta front sedimentationassociated to a coarseningupwards sequence.
Semi-parallelreflectors with low tomoderate amplitudeand low to moderatecontinuity.
Predominantly shalewith some claystone anssiltstone
Prodelta clay and siltsedimentation associatedto a coastline regression
Prodelta clay and siltsedimentation associatedto a a coastline regression
Lithology_Interpretation Geological_InterpretationDSMF_Penobscot_Banquereau
-3000.00 PSTM 3000.00
L-30 [MD]
Facies_T1:12592
600
800
1000
1200
1400
1600
1800
2000
2200
2400
2600
28002883
Upper
Subm
arin
e_Va
lley
Irregular SF4
UNDEF
SF3
SF2
SF1
Cylinder
Coarsening Upward
Fining Upward
Intermediate-2
Intermediate-1
Basal
Banquereau
ILXL
W
(a) (b)
E10031366 1330 131
105103820 1202 1184
118011651198
11471216
11291233
1111
0 500
-200
-400
-600
803000
2000
Submarine Valley
Banquereau Fm
WyanDot Chalk Fm
Dawson Canyon Fm
Logan Canyon Fm
Maximum Floding Surface [MFS]Subaerial Unconformity [SU]Maximum Regression Surface [MRS]
TruncationOfflap
Downlap
Onlap
[ms]
1000
-1000
-2000
-3000
0
3000
2000
1000
-1000
-2000
-3000
0
70
60
50
40
30
20
-800
-1000
-1200
-200
-400
-600
-800
-1000
-1200
GR_API[gAPI]80
70
60
50
40
30
20
GR_API[gAPI]
1000 1500 2000 2500m
12511193125945 11621021
1348Inline 1349
XInline 1180 XInline 1240XInline 1240Random line Random line line 1178
ILXL
W E10031366 1330 1311
105612751092
12751109
12381127
129310741038
1220 1202 11841180
11651198
11471216
11291233
1111
0 500 1000 1500 2000 2500m
1251109312691145 11621021
1348Inline 1349 Inline 1310 Inline 1240
Inline 1300 XInline 1180 XInline 1240XInline 1240Random line Random line line 1178
the package bounded by the purple (WyanDot Chalk Fm) and the light green surfaces show a downlapping stacking pattern, which can also be seen in the Wheeler domain (Figure 15). It can also be seen that, in the correlation wells, the purple surface represents the top of the WyantDot Chalk Formation and on top of it there are fining upward electrofacies associated with the base of the Banquereau Formation (Figure 13). Therefore, the light green surface represents the end of the transgression and so it is interpreted as a maximum flooding surface [MFS]. This package also corresponds to a set of prograding clinoforms with a gentle dip that downlaps the WyantDot Chalk Formation. Additionally, this package consists of semi-parallel reflectors with low to moderate amplitude and low
C T& F Vol . 9 Num . 2 D e c emb er 2 01 9 47
Ec op e t r o lVol . 9 Num . 2 D e c emb er 2 01 9
46 Ec op e t r o l
(a) (b)
150052
250
500
750
1000
1200 241
200
150
100
50
0
250[m]
-20001 240
TruncationToplapDownlapOnlap
TruncationToplapDownlapOnlap
In-line In-lineRGTTime(ms)
Tim
e(m
s)
Submarine Valley
Submarine Valley
Erosion
Hiatus
Hiatus
Erosion
Erosion
2000
250[m]
-20001
2000240
1400 Inline 1200 1079 1400 1200 1000
Submarine valley
inline
Figure 15. Wheeler domain along with the horizon framework on Xline 1258. The image (a) illustrates the enhanced seismic image with the horizon framework and image (b) shows its corresponding Wheeler diagram interpreted.
Figure 16. Wheeler domain along with the system tract interpretation and the base level change curve.
Figure 17. 3D System tracts model interpreted in the 3D survey. This model is in the time domain and has a vertical exaggeration of 20.
Figure 18. (a) Inline 1220 along with correlation well L-30. (b) Stratigraphic surfaces and submarine channel geometries interpreted are in Inline 1220. This image shows the amalgamated channels deposited inside the submarine canyon (purple
dashed lines on the right).
Wheeler Domain SystemsTracts
StratigraphicSurfaces
Base LevelRise Fall
LSTLST
HST-2HST-2
HST-1HST-1
TSTTST
FSSTFSST
1600 1400 1200 1300 0 MRS1
Top HST2
Top HST1
SB1
MFS1
<Unassigned>
End of Regression
Onset of Base Level fall
Onset of Base Level fall
End of Base Level fall
End of Transgression
64
112
137
182
219
50
100
150
200
RGT
In-line Cross-line
Systems_Tracts_3DSystem Tracts
LST - Submarine Valley
HST-2
HST-1
TST
4892000
740000738000 736000
734000732000 730000 489200
489400
X-asis
Z-asis
Y-asis
6000
0
-200
-400
-600
-800
FSST
ILXL
ILXL
SW SWNE
(a) (b)
11621220
12021220
12411220
12811220
13211220
13611220
14011220
14411220
0 500
-200
-400
-600
[ms] [ms]
803000
2000
Maximum Floding Surface [MFS]Subaerial Unconformity [SU]Maximum Regression Surface [MRS]
TruncationOfflap
Downlap
Onlap
1000
-1000
-2000
-3000
0
3000
2000
1000
-1000
-2000
-3000
0
70
60
50
40
30
20
-800
-1000
-1200
-200
-400
-600
-800
-1000
-1200
GR_API[gAPI]80
70
60
50
40
30
20
GR_API[gAPI]
1000 1500 2000 2500m
XInline 1180 XInline 1250 XInline 1290 XInline 1350 XInline 1453XInline 1365 XInlin
Random line Random line Random line Random line NE
1162122020
222022 1202
122012411220
12811220
13211220
13611220
14011220
14411220
XInline 1180 XInline 1250 XInline 1290 XInline 1350 XInline 1453XInline 1365 XInlin
Random line Random line Random line Random line
to moderate continuity, which was interpreted as seismic facies 1 (SF1, Figure 13). This package is bounded at the bottom by the previous sequence and at the top by the Maximum Flooding Surface [MFS] (Figure 14 & Figure 16). According to the stacking patterns and all the features just mentioned, this package is interpreted as a Transgressive System Tract [TST] (Figure 16).
The overlying package, between the light green and red surfaces, is characterized by downlapping reflectors of a prograding set of clinoforms overlying the MFS. These prograding clinoforms are being truncated at the top by an erosional surface (Figure 14). Moreover, this package consists of semi-parallel reflectors with low amplitude and low continuity, which was interpreted as seismic
facies 2 (SF2) and also exhibits coarsening upwards electrofacies in the correlation wells (Figure 13). This package was interpreted as a Highstand System Tract [HST] and is bounded at the bottom by a Maximum Flooding Surface [MFS] and at the top by a Sub-aerial Unconformity (SU) (Figure 14).
The next package, between the red and dark blue horizons, consists of semi-parallel reflectors with moderate to high amplitude and moderate to high continuity, and was interpreted as seismic facies 3 (SF3, Figure 13). This package consists of a set of offlap and onlap reflectors that exhibit prograding and aggrading stacking patterns, which together result in a forced regression of the coastline. This is an indicator of two important changes in the base level which are the
end of base level rise and the beginning of the base level fall. Based on the offlap and onlap reflector terminations and the prograding stacking patterns, it was inferred that during the early stages of sedimentation of this package the accommodation space is utterly consumed by the sedimentation rate, thus creating a regression of the coastline and progradation of fluvial and nearshore facies, which is common in early stages of base level fall and is supported by the electrofacies interpreted in the correlation wells (Figure 13). This package was interpreted as a Falling Stage System Tract [FSST] (Figure 16) and is bounded at the bottom by a Sub-aerial Unconformity (SU) and at the top by the Maximum Regression Surface [MRS] (Figure 14).
For the following package, between the dark blue and the red surfaces, it can be seen in the Wheeler domain (Figure 15) that the clinoforms start retrograding, which means that the base-level is rising again. Moreover, this package consists of semi-parallel reflectors with moderate amplitude and moderate continuity, which was interpreted as seismic facies 4 (SF4) and also exhibits an irregular electrofacies sequence in the correlation wells (Figure 13). This package was interpreted as a HST based on the characterization mentioned previously.
The youngest package interpreted in this section, between the red and sea floor surfaces, relates to a prograding stacking pattern which is being incised by a prominent submarine valley that eroded down to the last HST. This package consists of semi-parallel reflectors with moderate to high amplitude and moderate to high continuity, which was interpreted as seismic facies 3 (SF3). Additionally, this package is characteristic as it exhibits truncation reflector termination and submarine channel sedimentation (purple dashed lines in Figure
18). Based on the aforementioned description, this package was interpreted as a LST.
Associated to each system tract, a set of seismic facies was interpreted (Figure 13), which helped in the process of defining the gross deposition environment, which is a fluvio-delta prograding from NW to SE of the study area. Each seismic facies has a particular waveform that may be mapped in the entire seismic survey. Now, using Opendtect’s UVQ workflow®, it was possible to map each seismic facies. Figure 19 illustrates the waveform of the seismic facies and a seismic facies distribution map at the MSR. One can clearly see how the seismic facies change as they go deep in the basin.
C T& F Vol . 9 Num . 2 D e c emb er 2 01 9 49
Ec op e t r o lVol . 9 Num . 2 D e c emb er 2 01 9
48 Ec op e t r o l
0
SF3 SF2 SF1 SF4
0
(a)
(b)
5
10
15
20
25
28 Z(ms) Class
In-lineCross-line
1100
1200
1400
1600
1300
1400
NEZ
2 4 5
Figure 19. Wheeler domain along with the system tract interpretation and the base level change curve.
Figure 20. Seismic section illustrating the sand deposits associated with base level fall during the FSST.
Figure 21. Illustration of how sedimentation occurs in an incised valley during the different stages of base level fall and rise (Modified from [36]).
SUBMARINE CANYON FILLRECONSTRUCTION
The study area is located in the continental shelf of the Sable sub-basin. According to Wach and Mosher [16], marine shelf mudstones, sandstones and conglomerates of the Banquereau Formation were influenced throughout the Cenozoic by several major unconformities related to sea level fall. As mentioned before, sediment distribution of Miocene successions was strongly influenced by the Western Boundary Undercurrent with periods of intensified bottom current activity also occurring in the Late Pliocene [19] – [20], followed by widespread Gully cutting in the Early Pleistocene. These events of bottom current activity and Gully cutting in the Early Pleistocene probably marked the beginning of the formation of the canyon being analyzed. Later, during the Quaternary to recent, several hundred meters of glacial and marine sediment were deposited on the outer shelf and slope [17],[21]-[22]. Wach and Mosher [16] also argues that on the Scotian Shelf and Grand Banks, the widespread hiatus eroding either the upper part or all of the Oligocene, is marked by a regional unconformity; the nature of which includes canyon formation. Canyon incision at the shelf edge was initiated during the Eocene and was extensive by the Oligocene [8]. These erosional events are seen in the seismic information where after the base level fall, thick sandstone sequences (200 – 400 [ft]) were deposited on the shelf at the end of the FSST (Figure 20). In this context, the sedimentation within the canyon occurred in stages where the energy of the source of sediments, and therefore the grain size, varies according to the stage.
ILXL
SW
10031366
10381330
11801184
12161147
12511111
0 500
Maximum Flooding Surface [MFS]Subaerial Unconformity [SU]Maximum Regression Surface [MRS]
-200
-300
-400
-500
-600
-700
-800
-900
-1000
-100
[ms]
803000
2000
1000
-1000
-2000
-3000
0
70
60
50
40
30
20
System Tracts
Sand deposits generatedduring base level fall
LST - Submarine Valley
HST-2
FSST
HST-1
TST
GR_API[gAPI]
1000 1500 2000 2500m
Inline 1348 Inline 1XInline 1180 XInline 1240
XInline 1240Random line Random line Random line Inline 1178
TVD Litho_Tops Vsh_T Electrofacies Seismic_Facies System_traes_3DPenobscot_DSMF-3000.00 PSTM 3000.00
L-30 [MD]
Facies_T1:12370
600
800
1000
1200
1400
1600
1800
2000
2200
2400
2800
Upper
Subm
arin
e_Va
lley
Irregular SF4 HST-2
FSST
HST-1
TST
SF3
SF2
SF1
Cylinder
CoarseningUpward
Fining Upward
Intermediate-2
Intermediate-1
Basal
Banquereau
Highstand System Tract: isolated towards amalgamated channel fills
Fluvial aggradation
Fluvial aggradationtowards MFS
Paleosol formation and/orffluvial aggradation
Paleosol formation
Terrace deposits
Sequence boundary
Braided, amalgamated channels Meandering channels
rise risefall
Lagoonal depositsEstuarine deposits
Transgression: tidally influenced, meander channelsand legoonal deposits
Lowstand System Tract: amalgamated channel fills
Fallingstage System Tract: fluvial incision, terrace deposits,paleosols in interfluve areas
Incised valleys form during Fallingstage
ILXL
SW
(b)(a)
NE1002
-200
-400
730000
Morphology of the submarinecanyon interpreted
A
A´
732000 734000 736000 738000 740000 742000
730000 732000 734000 736000 738000 740000 742000
4902000
4901000
4900000
4899000
4898000
4897000
4896000
4895000
4894000
4893000
4892000
4891000
49020000.500.450.400.350.300.250.200.150.100.050.00 Y-axis
X-axis
Y-axis
FEF[Var] 1 [Realizae] 1Seismic - Variance
4901000
4900000
4899000
4898000
4897000
4896000
4895000
4894000
4893000
4892000
48910002500m
1:63114
-600
-800
-1000
-1200
133010421330
10821330
11221330
11621330
12021330
12411330
12811330
13211330
13611330
14011330
14411330
0 500
803000
2000
1000
-1000
-2000
-3000
0
70
60
50
40
30
20
GR_API[gAPI]
1000 1500 2000 2500m
Random line Random line Xline 1180 Xline 1240 Xline 1290 Xline 1350 Xline 1453XlineXline 1365Xline 1240
Figure 22.(a) Plan view of coherence attribute at time slice 552 [ms] showing the submarine canyon interpreted in the 3D seismic. (b) Section A – A’, Inline 1330, highlighting the submarine canyon interpretation.
C T& F Vol . 9 Num . 2 D e c emb er 2 01 9 51
Ec op e t r o lVol . 9 Num . 2 D e c emb er 2 01 9
50 Ec op e t r o l
Figure 21 illustrates how sedimentation occurs in an incised valley during the different stages of base level fall and rise, which explains how, after the incision of the valley, the sedimentation evolved to a set of amalgamated channels across the valley where few or no levees of the channels were preserved (Figure 18). In Figure 21 it can also be seen that the sedimentation started with what seems like an embedded stack of channel deposits in the lowest part of the canyon and, as mentioned before, it evolved to a set of amalgamated channel deposits upwards. Figure 22 and Figure 23 illustrate the morphology of the canyon across the seismic survey and the 3D geological model built from the integration of stratigraphic and structural seismic interpretation. Now, Figure 19 shows the Seismic Facies 3 (SF3) related to the Falling Stage System Tract (FSST) shown in Figure 20, which possesses the highest percentage of sandstones, and therefore the seismic facies distribution map shown in Figure 19 can be used as a guide for well placement in the decision making process.
It is highly recommended that stratigraphic columns generated from the description of core and cuttings samples be used as input data to calibrate all of the interpretations of litho-types made based on neural networks or any other methodology, since well logs represent an indirect measure of a physical property of the rocks while having a rock sample from the interval of interest is priceless.
Since the 3D seismic stratigraphic model proposed by this methodology is mainly based on seismic data, it needs to have a
4. RESULTS ANALYSIS
good S/N (signal/noise) ratio so that the automatic horizon tracking can produce a horizon framework where the resulting horizons will not cross seismic reflectors. Another way to improve the quality of the horizon framework used for sequence stratigraphic analysis is to incorporate the maximum possible number of horizons as constraints, calibrated with stratigraphic analysis performed with core data and well logs. The more guidelines the algorithm has to perform the automatic tracking of the horizon framework, the better the resulting horizon framework will be. Now, the constraints horizons need to be interpreted based on stratigraphic analysis and preferably, sequence-bounding surfaces, thus the interpreter will have control of how the automatic tracking of the horizon framework would turn out.
Referring specifically to system tracts interpretation from the base level curve (Figure 16), it could be argued that right after the Falling Stage System Tract (FSST) an almost completely eroded Transgressive System Tract (TST) might be interpreted following the high negative amplitude above the Maximum Regression Surface (MRS) at approximately 550 [ms] (Figure 18).
Besides just having a basic 2D sequence stratigraphic model, this methodology provides the interpreter with a 3D geological model that incorporates 3D sequence stratigraphy, structural geology and sedimentation history, thus allowing the extraction of geobodies from the seismic facies associated with a specific system tract that has the highest probability of sandstone discovery. Together, these work as a tool that helps the interpreter when it comes to well placement decisions.
Channe_1Elevetion time [ms]
-150-200-250-300-350-400-450-500-550-600-650-700-750
Morphology of thesybmarine canyoninterpreted
Figure 23.Illustration of the 3D geological model built for the study area. This model is in the time domain and has a vertical exaggeration of 15.
CONCLUSIONS The implementation of advanced techniques for 1D facies modeling, such as neural networks calibrated with core description, represents a step forward in the definition of accurate 1D facies models that are used to calibrate the final 3D facies model.
The workflow proposed in this paper for seismic image enhancing guarantees that the seismic data will have the quality necessary for stratigraphic interpretation. Using the resulting enhanced seismic image, a set of four seismic facies were interpreted and properly linked to the electrofacies in the correlation wells, enabling the identification of two main stratigraphic sequences in the Wheeler domain.
Based on the 3D sequence stratigraphic model calibrated with well data, a set of system tracts such as TST, HST, FSST and LST were identified, corresponding to the transgression and regression phases associated with the sedimentation of the Banquereau Formation. This Formation was deposited in a fluvio-
deltaic environment that was eroded down to the last TST sequence by a submarine canyon during the HST-2.
During the description of the stratigraphic packages in the seismic information, the top of the Banquereau Formation is interpreted as the base of a submarine canyon that was generated throughout the prominent erosion event at the end of the prograding HST-2, which afterwards was filled up by an amalgamated channel complex, indicating that the sedimentary environment interpreted matches the description made by Wach (2010).
The integration of advanced techniques for 1D facies modeling along with seismic facies modeling from waveform and 3D seismic stratigraphic analyses allows the interpreter to draw much deeper conclusions about the gross deposition environment and, most importantlysps, have the right information for decision making in the exploration phase, thus reducing the uncertainty in prospective areas.
ACKNOWLEDGEMENTSThis research was supported and carried out under an agreement between Universidad Industrial de Santander and Ecopetrol S.A., named “Acuerdo de Cooperación No. 9 Derivado del Convenio Marco 5222395.” We thank our colleagues from Ecopetrol S.A. and Universidad Industrial de Santander who provided insight and expertise that greatly assisted the research.
REFERENCES[1] Catuneanu, O., (2006). Principles of sequence stratigraphy. Amsterdam: Elsevier.
[2] OETR., (2001). Play Fairway Analysis, Nova Scotia Offshore. Department of Energy. Halifax, Nova Scotia, Canada.
[3] Mukhopadhyay, P. K., Brown, D. E., Kidston, A. G., Bowman, T. D., Faber, J., and Harvey, P. J., (2003). Petroleum Systems of Deepwater Scotian Basin, Eastern Canada: Challenges for Finding Oil versus Gas Provinces. Offshore Technology Conference, Houston, Texas, U.S.A.
[4] Colletta, B., Monnier, F., Rabary, G., Doublet, S., and Letouzey, P., (2011). Offshore Nova Scotia Play Fairway Analysis: 2D Basin Modeling Results. Offshore Technology Conference, Houston Texas, USA.
[5] Cummings, D. I., and Arnott, R. W. (2005). Growth-faulted shelf-margin deltas: a new (but old) play type, offshore Nova Scotia. In Bulletin of Canadian Petroleum Geology, 53 (3), 211-236, Program and Abstracts Volume, pp. 22. Conference on Sequence Stratigraphic and Allostratigraphic Principles and Concepts, Dallas, Texas, USA.
[6] McIver, N. L. (1972). Cenozoic and Mesozoic Stratigraphy of the Nova Scotia Shelf. Canadian Journal of Earth Sciences, 9, 54-70.
[7] Wade, J.A., and MacLean, B.C. (1990). Geology of the continental margin of eastern Canada. In Geological Survey of Canada, Geology of Canada, (2), 190-238, M. J. Keen, and G.L. Williams.
[8] Fensome, R. A., Crux, J. A., Gard, I. G., MacRae, R. A., Williams, G. L., Thomas, F. C., Fiorini, F., and Wach, G. (2008). The last 100 million years on the Scotian Margin, Offshore Eastern Canada. In Atlantic Geology 44. An event-stratigraphic scheme emphasizing biostratigraphic data, 93–126.
[9] Campbell, C. (2012). The Late Cretaceous and Cenozoic Geological History of the Outer Continental Margin of Nova Scotia Canada: Insights into Margin Evolution from a Mature Passive Margin. PhD Thesis, Dalhousie University, Nova Scotia, Canada.
[10] Tobey, D. (2013). Sequence stratigraphy of the Banquereau Formation, Offshore Nova Scotia. B.S. Thesis, Dalhousie University, Nova Scotia, Canada.
[11] PFA, (2011). Play Fairway Analysis (PFA) Offshore, Department of Energy, Nova Scotia, Canada, http://www.oera.ca/offshore-energy-research/geoscience/play-fairway-analysis/.
[12] PFA, (2014). Play fairway analysis (PFA) Offshore, Department of Energy, Nova Scotia, Canada.
[13] Ings, S. J., MacRae, R. A., Shimeld, J. W., and Pe-Piper, G. (2005). Diagenesis and porosity reduction in the Late Cretaceous Wyandot Formation, Offshore Nova Scotia: a comparison with Norwegian North Sea chalks. Bulletin of Canadian Petroleum Geology, 53 (3), 237-249, doi: 10.2113/53.3.237.
[14] CNSOPB, (2008). Regional Geology Overview, Canada-Nova Scotia Offshore Petroleum Board, Call
for Bids 2007-2008, http://www.callforbids.cnsopb.ns.ca/2007/01/regional_geology.html.
[15] Jansa, L. F., and Noguera, V. H. (1990). Geology and Diagenetic History of Overpressured Sandstone Reservoirs, Venture Gas Field, Offshore Nova Scotia, Canada. The American Association of Petroleum Geologist Bulletin, 74 (10), 1640-1658.
[16] Wach, G., and Mosher, D. (2010). Reservoir distribution and characterization: Shelf to slope linked depositional systems. Offshore Energy Technical Research Association.
[17] Kidston, A.G., Brown, D.E., Altheim, B., and Smith, B.M. (2002). Hydrocarbon Potential of the Deep-water Scotian Slope, Canada-Nova Scotia Offshore Petroleum Board (CNSOPB), Halifax, Nova Scotia, Canada.
[18] Gradstein F.M., Jansa L.F., Srivastava S.P., Williamson, M.A., Bonham C. G., and Stam, B. (1990). Aspects of North Atlantic paleoceanography. In Geology of the continental margin off eastern Canada, pp. 351-389. M.J. Keen and G.L. Williams, Geological Survey of Canada, Geology of Canada.
[19] Piper, D. J. W. (2005). Late Cenozoic evolution of the continental margin of eastern Canada, Norwegian Journal of Geology, 85 (4), 305-318.
[20] Myers, R.A., and Piper, D.J.W. (1988). Seismic stratigraphy of late Cenozoic sediments in the northern Labrador Sea: a history of bottom circulation and glaciation, Canadian Journal of Earth Sciences, 25 (12), 2059–2074, doi: 10.1139/e88-191.
Vol . 9 Num . 2 D e c emb er 2 01 9
[21] Piper D. J. W. and Sparkes, R. (1987). Proglacial sediment instability features on the Scotian Slope at 63°W, Marine Geology, 76, 15–31, doi: 10.1016/0025-3227(87)90014-4.
[22] Mosher, D. C., Moran, K., and Hiscott, R.N. (1994). Late Quaternary sediment, sediment mass flow processes and slope stability on the Scotian Slope. Sedimentology, 41, 1039–1061, doi: 10.1111/j.1365-3091.1994.tb01439.x.
[23] Pe-Piper, G. and Piper, D. J. W. (2004). The effects of strike-slip motion along the Cobequid-Chedabucto SW Grand Banks fault system on the Cretaceous-Tertiary evolution of Atlantic Canada. Canadian Journal of Earth Sciences, 41 (7), 799-808, doi: 10.1139/e04-022.
[24] Balkwill, H.R. and McMillan, N.J. (1990). Mesozoic-Cenozoic geology of the Labrador shelf, Baffin Island, and Davis Strait, In Geology of the continental margin of Eastern Canada, 2, pp. 31-85. M.J. Keen and G.L. Williams, Geological Survey of Canada, Geology of Canada.
[25] Haq, B.U., Hardenbol, J. and Vail, P.R. (1987). Chronology of fluctuating sea levels since the Triassic (250 million years ago to present). Science, 235, 1156-1166.
[26] Poag, C.W. (1985). Depositional history and stratigraphic reference section for central Baltimore Canyon trough. In Geologic Evolution of the United States Atlantic Margin, C.W. Poag, pp. 217-263, New York: Van Nostrand Reinhold Co.
[27] Mountain, G. S., and Tucholke, B. E. (1985). Mesozoic and Cenozoic geology of the U.S. Atlantic continental slope and rise. In Geologic evolution of the United States Atlantic Margin, C. W. Poag, pp. 293–341, New York: Van Nostrand Reinhold Co.
[28] Tucholke, B. E., and Mountain, G. S. (1986). Tertiary paleoceanography of the western North Atlantic Ocean. In The Geology of North America, M, The Western North Atlantic Region (38), P.R Vogt and B.E. Tucholke, pp. 631-650, Geological Society of America, Boulder, CO.
[29] Wade, J.A. and McLean, B.C. (1990a). The geology of the southeastern margin of Canada. In Geology of the Continental Margin of Eastern Canada, The Geology of North America, M.J. Keen and G.L. Williams, I (1), pp. 167-238, Geological Survey of Canada, Geology of Canada 2.
[30] Grist, A.M. and Zentilli, M. (2003). Post-Paleocene cooling in the southern Canadian Atlantic region: evidence from apatite fission track models. Canadian Journal of Earth Sciences, 40 (9), 1279-1297, doi: 10.1139/e03-045.
[31] Catuneanu, O., Hancox, P. J., Cairncross, B., and Rubidge, B. S. (2002). Foredeep submarine fans and forebulge deltas: orogenic off-loading in the underfilled Karoo Basin. Journal of African Earth Sciences, 35 (4), pp. 489–502, doi: 10.1016/S0899-5362(02)00154-9.
[32] Asquith, G. and Krygowski, D. (2004). Basic Well Log Analysis: AAPG Methods in Exploration, 2 nd ed. (16), pp. 244. American Association of Petroleum Geologists, Tulsa, OK.
[33] Bhatt, A., and Helle, H. B. (2002). Determination of Facies from Well Logs Using Modular Neural Networks. Petroleum Geoscience, 8 (3), pp. 217–228, doi: 10.1144/petgeo.8.3.217.
[34] Odoh, B. I., Ilechukwu, J. N., and Okoli, N. I. (2014). The Use of Seismic Attributes to Enhance Fault Interpretation of OT Field, Niger Delta. International Journal of Geosciences, (5), pp. 826-834, doi: 10.4236/ijg.2014.58073.
[35] De Bruin, G., Ligtenberg, H., Hemstra, N. and Tingdahl, K. (2006). Synchronized Sequence Stratigraphic Interpretation in the Structural and Chrono-stratigraphic (Wheeler transformed) Domain. EAGE Research Workshop - From Seismic Interpretation to Stratigraphic and Basin Modelling, Present and Future, doi: 10.3997/2214-4609.201403022.
[36] Shanley, K. and McCabe, P. (1993). Alluvial Architecture in a Sequence Stratigraphic Framework: A Case History from the Upper Cretaceous of Southern Utah, USA. In Quantitative Modeling of Clastic Hydrocarbon Reservoirs and Outcrop Analogues, The International Association of Sedimentologists, S. Flint, and I. Bryant. Special Publication, 15, pp. 21-56.
Frequently, biodiesel production operations are limited because of the precipitation of haze in their finished product tanks or those of their customers, which prevents compliance with the cold soak filtration test. It has been widely proved by literature that these deposits arise from the presence of polar chemical compound traces that cannot be eliminated during the conventional transesterification process.
Ecopetrol's Innovation and Technology Center (ICP) has developed the technology known as Ecoteb, which solves this problem. After an exhaustive research, development and upscaling process, two biodiesel producing companies, which represent 40% of Colombia's biodiesel domestic market, have acquired the right to use the technology developed by ICP, and are currently in the engineering and construction phase of ECOTEB plants.
Con frecuencia los productores de biodiesel ven limitadas sus operaciones por la precipitación de Haze que aparece en sus tanques de producto terminado o los de sus clientes e impiden el cumplimiento de la norma de la prueba de filtración en frío (cold soak filtration test). La literatura ha mostrado ampliamente que estos depósitos surgen por la presencia de trazas de compuestos químicos polares que no es posible eliminar durante el proceso convencional de transesterificación.
El ICP desarrolló la tecnología Ecoteb que permite resolver este problema después de seguir un proceso exhaustivo de investigación, desarrollo y escalado. Dos de las empresas productoras de biodiésel, y que representan el 40 % del mercado nacional de biodiésel de Colombia, han adquirido el derecho de uso de la tecnología desarrollada por el ICP y se encuentran en la fase de ingeniería y construcción de plantas de Ecoteb.
ECOTEBECOTEB