Improvements to the fuzzy mathematics comprehensive quantitative method for evaluating fault sealing
Improvements to the fuzzy mathematics comprehensive quantitative method for evaluating fault sealing
Da-Wei Dong 0 1 2 3
Ji-Yan Li 0 1 2 3
Yong-Hong Yang 0 1 2 3
Xiao-Lei Wang 0 1 2 3
Jian Liu 0 1 2 3
0 Working Stations for Post Doctors , Dongying 257000, Shandong , China
1 Research Institute of Petroleum Exploration and Development, Shengli Oilfield Company , Dongying 257000, Shandong , China
2 Shengli College of China University of Petroleum , Dongying 257061, Shandong , China
3 Oil Production Plant of Dongxing, Shengli Oilfield Company , Sinopec, Dongying 257061, Shandong , China
Fuzzy mathematics is an important means to quantitatively evaluate the properties of fault sealing in petroleum reservoirs. To accurately study fault sealing, the comprehensive quantitative evaluation method of fuzzy mathematics is improved based on a previous study. First, the single-factor membership degree is determined using the dynamic clustering method, then a single-factor evaluation matrix is constructed using a continuous grading function, and finally, the probability distribution of the evaluation grade in a fuzzy evaluation matrix is analyzed. In this study, taking the F1 fault located in the northeastern Chepaizi Bulge as an example, the sealing properties of faults in different strata are quantitatively evaluated using both an improved and an un-improved comprehensive fuzzy mathematics quantitative evaluation method. Based on current oil and gas distribution, it is found that our evaluation results before and after improvement are significantly different. For faults in ''best'' and ''poorest'' intervals, our evaluation results are consistent with oil and gas distribution. However, for the faults in ''good'' or ''poor'' intervals, our evaluation is not completely consistent with oil and gas distribution. The improved evaluation results reflect the overall and local sealing properties of target zones and embody the nonuniformity of fault sealing, indicating the improved method is more suitable for evaluating fault sealing under complicated conditions
Fault sealing property; Fuzzy mathematics; Dynamic clustering method; Quantitative study
Throughout geological history, faults have played a very
important role in hydrocarbon migration and
accumulation. Large faults, which form the boundaries of oil and
gas fields, can generally control hydrocarbon
accumulation (Agosta et al. 2012; Allan 1989; Balsamo et al.
2010; Bouvier et al. 1989; Braathen et al. 2009; Brogi
and Novellino 2015; Choi et al. 2015; Collettini et al.
2014; Fisher and Jolley 2007). However, small faults
generally separate oil and gas, resulting in an increase in
well spacing that poses challenges in oil and gas
exploration. Thereby, studies of fault sealing have
become a very important focus for oil and gas
exploration and development (Ciftci et al. 2013; Collettini
et al. 2014; Davatzes and Aydin 2005; Fachri et al.
2013a, b). Studies of fault sealing have evolved from
trap theory, which is used to identify hydrocarbon
sealing mechanism(s) and physical parameters of faults
(Hubbert 1953; Smith 1980). Based on different
mathematical principles and modes, the sealing properties of
one or more given faults are quantitatively identified and
predicted by considering various sealing mechanisms and
factors. (Childs et al. 2009; Egholm et al. 2008; Fachri
et al. 2011; Faulkner et al. 2003; Fisher and Knipe 1998;
Wu et al. 2010; Knipe et al. 1997).
With constant improvements in petroleum geology,
scholars in China and abroad have gradually realized that
fault sealing is jointly controlled by multiple factors, rather
than just one or two single factors (L u¨ et al. 2007; Knott
1993; Knipe et al. 1997; Faulkner et al. 2010; Gibson 1994;
Zhang et al. 2013; Pei et al. 2015). With a deepening of
quantitative single-factor studies of fault sealing, many
scholars tried to comprehensively evaluate fault sealing
using mathematical theory methods, including nonlinear
mapping methods, gray relational, logical information,
fault connectivity probabilistic methods and fuzzy
comprehensive evaluation methods (Fu et al. 2005, 2008, 2012;
Lu¨ et al. 1995; Lu¨ and Fu 2002; Zhang et al. 2007; Jiang
et al. 2008; Li 2009; Zhang et al. 2015). However, for all
these methods, some key parameters are difficult to obtain
and have regional limitations, so that human factors are
generally introduced for value assignment or adjustment.
Thus, scholars manage to optimize these parameters using
a variety of mathematical approaches. Among them, the
fuzzy comprehensive evaluation method, which is strongly
systematic and provides definite results, is accepted by
many scholars and has been gradually improved over the
In this paper, based on the principles of fuzzy
comprehensive evaluation, the single-factor membership degree
was established using the dynamic clustering method.
Then, a single-factor evaluation matrix was constructed
using the continuous grading function in order to optimize
the probability distribution of the evaluation grade in the
fuzzy comprehensive evaluation method. Next, we
considered a fault in the Chepaizi Bulge in the northwestern
margin of the Junggar Basin as an example, and the sealing
properties in the vertical and strike directions were
evaluated. Finally, our results are compared with oil and gas
2 Un-improved fuzzy mathematics evaluation of fault sealing
The physical principles that affect fault sealing include
principal stress sealing, lithological allocation sealing,
shale smear sealing, time allocation sealing and occurrence
allocation sealing (Liu 1998). Single factors include fault
properties, fault plane pressure, lithological allocation,
fault dip angle, shale smear and fault active time. Thus,
being affected by numerous factors, fault sealing shows
very strong complexity and randomness. Fortunately, the
fuzzy evaluation method can take various types of single
factors that affect fault sealing into consideration, and
hence realize a comprehensive evaluation using the
principle of fuzzy transformation and maximum membership
degree. Using this mathematical method, fault sealing
properties can be effectively identified. The process is as
follows: First, based on regional geology and fault
development, the single factors that affect fault sealing are
selected and quantified in order to obtain a single-factor
quantization matrix Un 1 (n is the number of single
factors). Then, based on the general division of fault sealing
and regional requirements, the number of evaluation grades
is divided (m). Finally, a single-factor membership degree
matrix is constructed based on the number of evaluation
grades and the oil and gas exploration results of the work
area, V1 m, based on which, Un 1 is evaluated in order to
construct a fuzzy evaluation matrix Rn m ¼ Un 1 V1 m.
Since each single factor has different contributions to fault
sealing, once the weighted matrix of each single factor is
provided, W1 n, the fuzzy evaluation matrix B1 m ¼ W1 n
Rn m can be obtained. Finally, the maximum value of the
fuzzy evaluation matrix B19m is taken as the corresponding
However, there are various methods that can be used to
acquire the parameters mentioned above, both qualitatively
and quantitatively, and which are constantly optimized (Lu¨
and Fu 2002; Russell et al. 2003; Fu et al. 2012; Li et al.
2009). There are many single factors that affect fault
sealing, generally including fault properties, fault plane
pressure, fault lithological allocation and shale smear.
Evaluation grades can be divided into five ranks: good,
better, moderate, fairly poor and poor. The single-factor
membership degree can be obtained via a discrete function
method and/or a continuous function method. The former
applies to qualitative single-factor membership degree,
such as fault properties, while the latter applies to
quantitative single-factor membership degree, such as fault plane
pressure. Weighting coefficients can be obtained via the
expert survey method, analytic hierarchy process (AHP),
Delphi method and/or the weight matrix method (Ding and
Jin 2012; Sun and Wu 1995; Qiu et al. 2007). The larger
the weighting coefficient of a single factor, the greater the
impact of the single factor on fault sealing. The
mathematical model of fuzzy comprehensive evaluation includes
the weighted-average type, the main factor highlight type
and the main factor determination type (Sun et al. 2010).
The weighted-average type can take a variety of single
factors into consideration in order to avoid information
loss. The main factor highlight type and the main factor
determination type emphasize the main controlling factors
and prevent interference factors. In the evaluation process
of fault sealing, many scholars use the weighted-average
type. However, determining the above parameters and
effectively predicting fault sealing must be based on the
actual petroleum geology of the study area and from
suggestions made by local experts in order to reduce
3.2 The establishment of a single-factor evaluation
matrix using the continuous grading function
3 Improved fuzzy mathematics of fault sealing
3.1 The establishment of the single-factor
membership degree using the dynamic
In this study, the membership degree of each single
factor is established by applying the dynamic clustering
method. Based on a large sample dataset,
pre-classification is roughly performed. Then, gradual adjustment is
made until a reasonable classification is obtained. The
maximum and minimum values of each class constitute
the range of membership degrees. For fault sealing, the
single-factor membership degree is established as
follows: Firstly, a large number of samples of one single
factor of the work area and their corresponding oil–gas
show are constructed. Then, taking best, good, poor and
poorest oil and gas show as the standard, the range of
the membership degrees of the fault gouge ratio is
determined. This means that the fault gouge ratio of a
given area is clustered. In the process of clustering, the
maximum and minimum values are the range limits of
membership degree. The so-called dynamic clustering
means the gradual adjustment of a value to a reasonable
range using the dynamic clustering principle. This
method has the advantages of incorporating the oil and
gas geology of the study area, so that it has regional
eigenvalues. It is closely related to the fault sealing
properties of the study area and hence has a certain
degree of predictability. The advantage of this method is
that the regional petroleum geology has a regional
characteristic value which is related to fault sealing, and
After the single-factor membership degree is determined
via the dynamic clustering method, fuzzy evaluation should
be carried out for all single factors in order to establish a
single-factor evaluation matrix. The division of the
singlefactor membership degree SðiÞ only gives evaluation
results of good and poor intervals, and it does not provide
probabilistic evaluation results of each interval. For
example, the range of ‘‘good’’ and ‘‘poor’’ is (1–0.5) and
(0.5–0). If the single factor value is 0.7, it cannot be
regarded simply as ‘‘good,’’ thus implying that a further
probabilistic evaluation is needed. Hence, the definition
‘‘good’’ accounts for 70% good, and ‘‘poor’’ accounts for
30%. eðiÞ is the boundary value of the re-probabilistic
evaluation. The precise method followed is as
follows: First, a range of values for fault sealing is determined
using single factors, SðiÞ. Then, the threshold value of each
adjacent class is determined, namely the classification
representative value eðiÞ, which is determined as per the
following principles (Zhao 2001):
eð1Þ ¼ Sð1Þ >9
eð2Þ ¼ ½Sð1Þ þ Sð2Þ =2 =>
eð4Þ ¼ Sð3Þ
Fault properties, fault gouge ratios and fault sealing are
not always linearly related with each other. However, they
have a piecewise function relation based on single-factor
membership degree. According to Newton’s iteration
principle, it can be assumed that the single-factor
evaluation criteria rðvÞ are a linear piecewise function, which can
be solved by successive approximation. Then, fuzzy
subsets of the fault evaluation criteria can be determined using
the following methods:
Four Trees Sag
Piedmont Fault Zone
First-order tectonic unit
First-order tectonic unit
Industrial Oil flow well
4 Case studies
Fig. 1 The position of the seismic profile study of fault sealing
The fuzzy subset obtained based on N single-factor
evaluation indices constitutes fuzzy sets:
Taking the fault F1 located in the northeastern Chepaizi
Bulge as an example (Fig. 1), the fault sealing properties of
the main target zones were evaluated using both the
improved and un-improved methods. This approach was
used to determine whether the improved comprehensive
evaluation method for fault sealing is more reasonable. The
Chepaizi Bulge is located at the southern end of the
northwestern margin of the Junggar Basin. It is a secondary
structural unit in the western uplift of the Junggar Basin,
adjacent to the Changji Sag and the Zhongguai bulge to the
east, Sikeshu Sag to the south and Zaire Mountain to the
northwest (Fig. 1). Conditions in the Chepaizi Bulge make
this a good place for hydrocarbon accumulation. The main
reservoirs are Cretaceous (K) and Neogene Sha1 Member
(N1s1). The F1 fault in the northeastern Chepaizi Bulge is a
very important oil control fault, and its sealing properties
directly affect hydrocarbon accumulation in this area. The
fault has NNE-trending, EW-trending compression in the
Yanshanian Period and NWW-trending extension in the
Himalayan Period. Vertically, there are reverse faults in the
Cretaceous strata and basement and normal faults in
Cenozoic strata, which represent a negative inverted
structure as a whole.
4.1 Single-factor selection and quantitative
Since the study area has undergone multiple phases of
tectonic evolution, various types of faults with different
occurrences have developed. There are many single factors
that affect fault sealing in response to tectonic evolution, so
that the weighting coefficients of the single factors are
different. However, in this paper, and taking fault F1 as an
example, the sealing properties in the vertical and strike
directions were analyzed. Four single factors with
comparable significance were selected, including fault surface
normal stress (r), fault properties, sand/formation ratio (N/
G) and the shale content of fault zone fillings. In addition,
four seismic profiles vertical to fault F1 were selected (the
location is shown in Fig. 1), which were converted into
four geological sections via a time-depth conversion. Then,
lithological sections were recovered using well logs in
order to obtain the quantitative and qualitative values of
each single factor (Table 1).
Fault surface normal stress, which is an important
parameter used to characterize the fault opening degree, is
crucial to fault sealing (Fu et al. 2005). Typically, this
value is the vector sum of the gravitational force imparted
by the overlying strata, the regional principal compressive
stress, the regional principal stress as well as the fault dip
angle. Fault F1 is a shovel fault, where the fault plane
normal stress increases from the bottom to the top, where
the fault sealing properties become better. However, the
Chepaizi Bulge is located in a slope belt of a foreland
basin, where the buried depth is shallow, so that the fault
surface pressure is small.
In general, the sealing properties of compressional, shear
and compressional-shear faults are good, while the sealing
properties of extensional and extensional-shear faults are
poor. Fault F1 is a negative inverted fault, and it is an
extensional fault that formed during the Neogene. However,
the fault occurrence changes significantly, which is steep in
the upper part and gentle in the lower part. That is, regional
tensile stress can only cause the Neogene ‘‘steep’’ fault
section to extend, but cannot completely make the Cretaceous
‘‘gentle’’ fault section extend. In addition, the pressure
experienced by the Cretaceous fault plane is significantly
larger than that of the Neogene fault plane, which can also
prove this change. After undergoing compressional
deformation, the Cretaceous fault section has better sealing
properties than the Neogene extensional fault section due to its
fault structure and mudstone smearing. Therefore, the
Neogene fault section has been assigned as ‘‘extensional’’, while
the Cretaceous fault section is assigned as ‘‘compressional’’.
Sand shale contraposition is an important process for oil
and gas lateral sealing. In fault dislocation, if sand shale
contraposition occurs in one interval, the fault in this
interval is closed, while if multiple sand layers are
connected in one interval, the fault in this interval is open. The
parameter characterizing possible sand shale contraposition
is the sand/formation ratio, which is the ratio between the
sandstone and layer formation thicknesses, N/G. When
N/G is large, the possibility of a sand–sand connection is
large, and hence the fault sealing ability is poor. When the
N/G ratio is small, the fault sealing ability is in fact good.
However, this single factor does not take the fault throw
generated by fault slipping into consideration, Therefore,
the following single factor is added.
Faults will form fault zones during displacement and
dislocation processes. When a given fault zone is filled
with shale, it can seal oil and gas laterally due to the plastic
flow and compaction of shale. When a fault zone is filled
with sandstone, it can act as a hydrocarbon migration
pathway since sandstone compacts poorly and has good
porosity and permeability. Accordingly, Fu et al. (2012)
improved the fault gouge ratio (SGR) proposed by Yielding
et al. (1997) and proposed to characterize the sealing
property of faults using the ratio between the mudstone
thickness of a fault belt and the sum of the fault throw and
Table 1 Evaluation parameter list of single factors on the sealing ability of the F1 fault
Section line no. Formation
r, MPa Fault surface
Section line no. Formation
r, MPa Fault surface
faulted formation thicknesses, namely the shale content of
the fillings in a given fault belt. Due to the negative
inversion of fault F1, Cretaceous shale smeared the fault
surface repeatedly. Therefore, the shale in the Cretaceous
fault zone fillings was formed in compressional and
extensional stages. Hence, the faults that formed in
compressional and extensional stages were obtained from
structural section restoration.
4.2 Single-factor weighting coefficients
and membership degree
Single-factor weighting coefficients are significantly
different in different areas. The numerous methods that are
used to determine these coefficients are all based on the
expert investigation method. Therefore, in this paper, the
weighting coefficients of four single factors, including the
normal stress of the fault plane (w1), fault properties (w2),
sand/formation ratio (w3) and the shale content of the
fillings in the fault belt (w4), were determined using the expert
investigation method (Table 2). It is worth mentioning that
the formation pressure of the oil and gas discovery wells in
this area is abnormal, while the formation pressure of dry
wells and water wells is normal, meaning pressure is
sensitive to oil and gas accumulation so that the weighting
coefficient of fault plane normal stress is the largest.
The statistics used in this study were applied to the well
logging data of 85 wells around fault F1 and the
corresponding oil and gas shows. Moreover, the values of four
single factors were calculated, which were classified into
four categories using the dynamic clustering method: best,
good, poor and poorest. The maximum and minimum
values of each category were taken as the boundary values
of each membership degree. Based on the advice of oilfield
experts, membership degree was then slightly adjusted.
The final results are shown in Table 3.
4.3 Optimized fuzzy evaluation matrix
The advantage of the dynamic clustering method in
comprehensive fault sealing property evaluation is to establish
a single-factor evaluation matrix. In order to illustrate the
establishment process of the un-improved and improved
single-factor evaluation matrices, the establishment of the
Table 2 Single-factor weight coefficients of fault sealing in the
northeast of Chepaizi Uplift
Table 3 The single-factor membership list in the study area
single-factor matrix in the first member of the Shawan
Formation in profile I was analyzed.
Before improvement, a single value was directly
assigned to a single-factor membership degree evaluation
matrix using the discrete function and continuous function
methods. When the N/G of the first member of the Shawan
Formation was 0.56, (i.e., the grade of the maximum
membership degree is ‘‘poor’’), this grade was directly
assigned a value of 0.5 when using the discrete function
method and 0.3 when using the continuous function
method. The single-factor membership degree evaluation
matrix is RN=G ¼ ð0:5Þ or ð0:3Þ. However, although a
value of 0.56 implies a grade of ‘‘poor,’’ which is close to a
grade of ‘‘good,’’ a value of 0.56 should have some
probability evaluation in the grade ‘‘good.’’ However, the
unimproved single-factor membership degree evaluation
matrix does not embody the probability in this transitional
Instead, a single-factor evaluation matrix that is
constructed by the continuous grading membership function
can properly address the problem mentioned above. Its
calculation process after improvement is as follows: First,
SðiÞ is determined based on N/G, namely S 1
ð Þ ¼ 0:3,
ð Þ ¼ 0:55, Sð3Þ ¼ 0:8. Then, the threshold value of each
adjacent evaluation level is calculated, i.e., grading
representative value eðiÞ:
> eð1Þ = Sð1Þ = 0:3
< eð2Þ ¼ ½Sð1Þ þ Sð2Þ =2 ¼ 0:425
>> eð3Þ ¼ ½Sð2Þ þ Sð3Þ =2 ¼ 0:675
: eð4Þ ¼ Sð3Þ ¼ 0:8
After improvement, the evaluation matrix of the
sand/formation ratio of the first member of the Shawan
Formation is RN=G ¼ ð 0 0:46 0:54 0 Þ.
The comprehensive evaluation matrix of the four factors
of the first member of the Shawan Formation is:
2 0:00 0:51 0:49 0:00 3
6 0:00 0:00 0:80 0:20 7
R ¼ 64 0:00 0:46 0:54 0:00 75
0:00 0:00 0:00 1:00
The fuzzy evaluation matrix of the first member
of the Shawan Formation is B ¼ W R ¼
ð 0 0:27 0:44 0:29 Þ. According to the principle of the
maximum membership degree, the evaluation result of the
first member of the Shawan Formation is poor. However,
the probability distribution of this evaluation grade is
scattered, indicating the fuzziness of the fault sealing
property is large. Conversely, if the probability distribution
of the evaluation grade is concentrated, the evaluation
result of fault sealing property is obvious.
4.4 Fuzzy evaluation results and analysis
Using the improved fuzzy evaluation method described
above, the sealing properties of fault F1 in the vertical and
strike directions were evaluated. The results are shown in
Table 4. Cross sections I, II, III and IV are well-tied cross
sections that span from south to north along the fault strike.
Cretaceous (K), the first member of Shawan Formation
(N1s1), the second member of Shawan Formation (N1s2)
and the third member of Shawan Formation (N1s3)
Table 4 Evaluation form of fault sealing
constitute intervals from top to bottom in the vertical
direction. The evaluation results show that the valuation
result of Cretaceous faults is ‘‘good,’’ but the maximum
probability is only 0.58, indicating oil and gas are sealed
when they migrate to the footwall, but it may be oil and gas
bearing in the hanging wall. The valuation result of the first
member of the Shawan Formation is ‘‘poor’’ or ‘‘poorest.’’
The probability distribution of the evaluation grade is
dispersed, showing the ambiguity of footwalls and hanging
walls on its oil-bearing properties. The evaluation result of
the second member of the Shawan Formation is ‘‘good.’’ A
regional survey showed that this interval is a stably
distributed shale bed, namely regional caprock. The
evaluation result of the third member of the Shawan Formation is
‘‘poor,’’ but mud logging results show this interval is
conglomerate. There were no oil and gas shows discovered
in this interval in the work area.
Next, we compared the variations between the improved
and un-improved fuzzy evaluations and their impact on the
evaluation results. In this paper, a single-factor
membership degree evaluation matrix was directly assigned a value
using the discrete function method. The assigned values of
evaluation grades R ¼ ðbest, good, poor, poorest) = (1, 0
:66; 0:33; 0Þ and the fuzzy evaluation results are shown in
Table 4. Upon inspection, the fuzzy evaluation values of
the first and second members of the Shawan Formation
change significantly, but they do have consistent results.
However, the fuzzy evaluation values and evaluation
Improved evaluation value
Unimproved evaluation value
Oil and gas display
‘‘–’’ in the table means there are no drilling test data or wells are not drilled in this horizon
8 N1t 9 N1s3 10 N1s2 11 N1s1 12
Fig. 2 Analysis of fault sealing and its oil reservoir in sections I–IV 1 Conglomerate, 2 Sandstone, 3 Mudstone, 4 Breccia, 5 Igneous rocks, 6
Hydrocarbon migration direction, 7 Faults open, 8 Faults closed
results of the first member of the Shawan Formation
change significantly before and after improvement, and the
evaluation results before improvement are not well
matched with the oil and gas shows in footwall and hanging
wall. Taking cross section III as an example, the evaluation
results before improvement are poor, but the footwall and
hanging walls are both oil bearing. After improvement, the
probability distribution of the evaluation grades is
scattered, indicating oil and gas may be trapped locally after
having continuously migrated from the footwall to the
hanging wall in the local open interval (Fig. 2). The above
analysis shows that the improved method is more
advantageous for regions with larger degrees of vagueness of
their fault sealing properties, and hence higher exploration
The evaluation of fault sealing properties is an important
part of oil and gas exploration and development. Fuzzy
comprehensive evaluation is a systematic evaluation
method for fault sealing. Since traditional methods are
greatly affected by human factors during the establishment
of the single-factor membership degree, the assignment
method results in fuzzy evaluation values, so that the result
only reflects the target interval. In this paper, the dynamic
clustering method was introduced to determine the
singlefactor membership degree. Then, a single-factor evaluation
matrix was constructed using the continuous grading
function in order to determine the optimum comprehensive
evaluation matrix and make a fuzzy evaluation of the fault
sealing properties. A comparison of the fuzzy evaluation
and its result before and after improvement, combined with
current oil and gas distribution regularity, showed that the
evaluation results before and after improvement are
significantly different. For faults designated as ‘‘best’’ and
‘‘poorest,’’ the evaluation results are consistent with oil and
gas distribution. However, for faults designated as ‘‘good’’
or ‘‘poor,’’ the evaluation results of sealing property are not
completely consistent with the oil and gas distributions.
The improved results reflect the overall and local sealing
properties of target zones and embody the fuzziness of fault
sealing, indicating the improved method is more precise for
evaluating fault sealing properties under complicated
However, the improved method still has its limitations.
First of all, this method still cannot solve the multi-scale
fault sealing problem. As mentioned above, the N1s1
section in profile III displays strong heterogeneity, and
although the evaluation results of large scale are ‘‘poor,’’
the fault sealing properties exist at small scales. Secondly,
this method can only explain the sealing properties of the
faults in profile. When oil and gas migrate in multiple
directions, the evaluation results will not be consistent with
drilling results. For example, the evaluation results of
section K in the III and IV profiles are good, but the fault is
oil bearing in the hanging wall. This phenomenon is related
to a two-way oil supply. Oil source comparison shows that
the oil and gas migration in the study area is from east to
west in the Changji Sag and from south to north in the
Sikeshu Sag (Song et al. 2007). Therefore, this method is
only applicable for studies of single-scale fault sealing
properties. In addition, the direction of oil and gas
migration should be vertical to the fault plane.
Acknowledgements This study is supported by the Science and
Technology Project of Universities and Colleges in Shandong
Province ‘‘Investigation on diagenetic environment and transformation
pattern of red-bed reservoirs in the rift basins’’ (No. J16LH52).
Open Access This article is distributed under the terms of the Creative
Commons Attribution 4.0 International License (http://creative
commons.org/licenses/by/4.0/), which permits unrestricted use,
distribution, and reproduction in any medium, provided you give appropriate
credit to the original author(s) and the source, provide a link to the
Creative Commons license, and indicate if changes were made.
Agosta F , Ruano P , Rustichelli A , et al. Inner structure and deformation mechanisms of normal faults in conglomerates and carbonate grainstones (Granada Basin, Betic Cordillera, Spain): inferences on fault permeability . Struct Geol . 2012 ; 45 : 4 - 20 . doi:10.1016/j. jsg. 2012 .04.003.
Allan US . Model for hydrocarbon migration and entrapment within faulted structures . AAPG Bull . 1989 ; 70 ( 7 ): 803 - 11 . doi:10.1306/ 94885962 - 1704 - 11D7 -8645000102C1865D.
Balsamo F , Storti F , Salvini F , et al. Structural and petrophysical evolution of extensional fault zones in low-porosity, poorly lithified sandstones of the Barreiras Formation , NE Brazil. Struct Geol . 2010 ; 32 ( 11 ): 1806 - 26 . doi:10.1016/j.jsg. 2009 .10.010.
Bouvier JD , Kaars-Sijpesteijn CH , Kluesner DF , et al. Three-dimensional seismic interpretation and fault sealing investigations , Nun River Field, Nigeria. AAPG Bull . 1989 ; 73 ( 11 ): 1397 - 414 .
Braathen A , Tveranger J , Fossen H , et al. Fault facies and its application to sandstone reservoirs . AAPG Bull . 2009 ; 93 ( 7 ): 891 - 917 . doi:10.1306/03230908116.
Brogi A , Novellino R. Low angle normal fault (LANF)-zone architecture and permeability features in bedded carbonate from inner Northern Apennines (RapolanoTerme, Central Italy) . Tectonophysics . 2015 ; 638 ( 638 ): 126 - 46 . doi:10.1016/j.tecto. 2014 .11.005.
Childs C , Manzocchi T , Walsh JJ , et al. A geometric model of fault zone and fault rock thickness variations . Struct Geol . 2009 ; 31 ( 2 ): 117 - 27 . doi:10.1016/j.jsg. 2008 .08.009.
Choi JH , Yang SJ , Han SR , et al. Fault zone evolution during Cenozoic tectonic inversion in SE Korea . Asian Earth Sci . 2015 ; 98 : 167 - 77 . doi:10.1016/j.jseaes. 2014 .11.009.
Ciftci NB , Giger SB , Clennell MB . Three-dimensional structure of experimentally produced clay smears: implications for fault seal analysis . AAPG Bull . 2013 ; 97 ( 5 ): 733 - 57 . doi:10.1306/10161211 192.
Collettini C , Carpenter BM , Viti C , et al. Fault structure and slip localization in carbonate-bearing normal faults: an example from the Northern Apennines of Italy . Struct Geol . 2014 ; 67 : 154 - 66 . doi:10.1016/j.jsg. 2014 .07.017.
Davatzes N , Aydin A. Distribution and nature of fault architecture in a layered sandstone and shale sequence: an example from the Moab fault, Utah . AAPG Memoir . 2005 ; 85 : 153 - 80 . doi:10. 1306/1033722M853134.
Ding WL , Jin WZ . Research and application of the comprehensive evaluation system of multi information to fault sealing . Beijing: Geological Publishing House ; 2012 . p. 70 - 104 (in Chinese).
Egholm DL , Clausen OR , Sandiford M , et al. The mechanics of clay smearing along faults . Geology . 2008 ; 36 ( 10 ): 787 - 90 . doi:10. 1130/G24975A.1.
Fachri M , Rotevatn A , Tveranger J. Fluid flow in relay zones revisited: towards an improved representation of small-scale structural heterogeneities in flow models . Mar Pet Geol . 2013a; 46 ( 3 ): 144 - 64 . doi:10.1016/j.marpetgeo. 2013 .05.016.
Fachri M , Tveranger J , Braathen A , et al. Sensitivity of fluid flow to deformation-band damage zone heterogeneities: a study using fault facies and truncated Gaussian simulation . Struct Geol . 2013b; 52 ( 1 ): 60 - 79 . doi:10.1016/j.jsg. 2013 .04.005.
Fachri M , Tveranger J , Cardozo N , et al. The impact of fault envelope structure on fluid flow: a screening study using fault facies . AAPG Bull . 2011 ; 95 ( 4 ): 619 - 48 . doi:10.1306/09131009132.
Faulkner DR , Jackson CAL , Lunn RJ , et al. A review of recent developments concerning the structure, mechanics and fluid flow properties of fault zones . Struct Geol . 2010 ; 32 ( 11 ): 1557 - 75 . doi:10.1016/j.jsg. 2010 .06.009.
Faulkner DR , Lewis AC , Rutter EH. On the internal structure and mechanics of large strike-slip fault zones: field observations of the Carboneras fault in southeastern Spain . Tectonophysics. 2003 ; 367 ( 3 ): 235 - 51 . doi:10.1016/ S0040-1951(03)00134-3.
Fisher QJ , Jolley SJ . Treatment of faults in production simulation models . Geol Soc Lond Spec Publ . 2007 ; 292 ( 1 ): 219 - 33 . doi:10. 1144/SP292.13.
Fisher QJ , Knipe RJ . Fault sealing processes in siliciclastic sediments . Geol Soc Lond Spec Publ . 1998 ; 147 ( 1 ): 117 - 34 . doi:10.1144/ GSL.SP. 1998 .147.01.08.
Fu G , Liu HX , Du HF . Seal mechanisms of different transporting passways of fault and their research methods . Pet Geol Exp . 2005 ; 27 ( 4 ): 404 - 8 . doi:10.3969/j.issn. 1001 - 6112 . 2005 .04. 016 (in Chinese).
Fu G , Shi JJ , Lu¨ YF. An improvement in quantitatively studying lateral seal of faults . Acta Pet Sin . 2012 ; 33 ( 3 ): 414 -8 (in Chinese).
Fu G , Yin Q , Du Y. A method studying vertical seal of fault with different filling forms and its application . Pet Geol Oilfield Dev Daqing . 2008 ; 27 ( 1 ): 1 - 5 . doi:10.3969/j.issn. 1000 - 3754 . 2008 .01. 001 (in Chinese).
Gibson RG . Fault-zone seals in siliciclastic strata of Columbus Basin, Offshore Trinidad . AAPG Bull . 1994 ; 78 ( 9 ): 1372 - 85 . doi:10. 1306/ A25FECA7-171B-11D7-8645000102C1865D.
Hubbert MK . Entrapment of petroleum under hydrodynamic conditions . AAPG Bull . 1953 ; 37 ( 8 ): 1954 - 2026 . doi:10.1306/ 5CEAD D61-16BB-11D7-8645000102C1865D.
Jiang Z , Dong Y , Li H , et al. Limitation of fault-sealing and its control on hydrocarbon accumulation-an example from the Laoyemiao Oilfield of the Nanpu Sag . Pet Sci . 2008 ; 5 : 295 - 301 . doi:10.1007/s12182- 008 - 0049 -6.
Knipe RJ , Fisher RJ , Jones G , et al. Fault seal analysis: successful methodologies, application and future directions . Nor Pet Soc Spec Publ . 1997 ; 7 ( 97 ): 15 - 38 . doi:10.1016/S0928-8937(97)80004- 5 .
Knott SD . Fault seal analysis in the North Sea . AAPG Bull . 1993 ; 77 ( 5 ): 778 - 92 . doi:10.1306/ BDFF8D58-1718-11D7-86450 00102C1865D.
Li CH , Luo LY , Chen PY , et al. Analysis and application of weight coefficient and attribution function during fuzzy evaluation of fault sealing . Complex Hydrocarb Reserv . 2009 ; 2 ( 1 ): 5 - 13 . doi:10.16181/j.cnki.fzyqc. 2009 .01. 018 (in Chinese).
Li Y. Shale smearing and its quantitative characterization in the perspective of fault rocks-a case study of the Ying 32 fault in the Dongxin field of Jiyang depression . Acta Geol Sin . 2009 ; 83 ( 3 ): 426 - 34 . doi:10.3321/j.issn: 0001 - 5717 . 2009 .03. 010 (in Chinese).
Liu ZR . The mechanism and structural model of fault-block reservoir . Beijing: Petroleum Industry Press ; 1998 (in Chinese).
Lu¨ YF, Chen ZM , Chen FJ . Evaluation of sealing ability of faults using nonlinear mapping analysis . Acta Pet Sin . 1995 ; 16 ( 2 ): 36 - 41 (in Chinese).
Lu¨ YF, Fu G . Study of fault sealing. Beijing: Petroleum Industry Press ; 2002 (in Chinese).
Lu¨ YF, Sha ZX , Fu XF , et al. Quantitative evaluation method for fault vertical sealing ability and its application . Acta Pet Sin . 2007 ; 28 ( 5 ): 34 -8 (in Chinese).
Pei YW , Paton DA , Knipe RJ , et al. A review of fault sealing behaviour and its evaluation in siliciclastic rocks . Earth Sci Rev . 2015 ; 150 : 121 - 38 . doi:10.1016/j.earscirev. 2015 .07.011.
Qiu YB , Zha M , Qu JX . Fuzzy comprehensive evaluation of the sealability of the faults in Chenbao Oilfield . J Xi'an Shiyou Univ (Nat Sci Ed) . 2007 ; 22 ( 3 ): 28 - 32 . doi:10.3969/j.issn. 1673 -064X. 2007.03. 006 (in Chinese).
Russell KD , An LJ , Paul J , et al. Fault-seal analysis south Marsh Island 36 field, Gulf of Mexico . AAPG Bull . 2003 ; 87 : 479 - 91 . doi:10.1306/08010201133.
Smith DA . Sealing and nonsealing faults in Louisiana Gulf coast salt basin . AAPG Bull . 1980 ; 64 ( 2 ): 145 - 72 . doi:10. 1306/2F918946- 16CE-11D7-8645000102C1865D.
Song CC , He LJ , Ma LQ , et al. Characteristic of hydrocarbon accumulation of Chepaizi swell in Junggar Basin . Xinjiang Pet Geol . 2007 ; 28 ( 2 ): 136 -8 (in Chinese).
Sun H , Xu TW , Fan SW , et al. Synthetic fuzzy judgment of fault sealing in MQ area . J Oil Gas Technol . 2010 ; 32 ( 2 ): 200 -3 (in Chinese).
Sun SX , Wu Z. Application of fuzzy synthesized evaluation to the hydrocarbon generating conditions . J Chengdu Univ Technol . 1995 ; 22 ( 3 ): 6 - 10 (in Chinese).
Wu ZP , Wei C , Xue Y. Structural characteristics of faulting zone and its ability in transporting and sealing oil and gas . Acta Geol Sin . 2010 ; 84 ( 4 ): 570 -8 (in Chinese).
Yielding G , Freeman B , Needham DT . Quantitative fault seal prediction . AAPG Bull . 1997 ; 81 ( 6 ): 897 - 917 . doi:10.1016/ S09 28-8937(97)80010-0.
Zhang LK , Luo XR , Liao QJ , et al. Quantitative evaluation of fault sealing property with fault connectivity probabilistic method . Oil Gas Geol . 2007 ; 28 ( 2 ): 181 - 90 (in Chinese).
Zhang LK , Luo XR , Song GQ , et al. Quantitative evaluation of parameters to characterize fault opening and sealing during hydrocarbon migration . Acta Pet Sin . 2013 ; 34 ( 1 ): 92 - 100 . doi:10.7623/syxb201301010 (in Chinese).
Zhang J , Wu ZP , Li W , et al. The fault system of the north section of Liaodong Salient, constraints from seismic data and physical modeling experiment . Acta Geol Sin . 2015 ; 89 (s1): 217 - 27 . doi:10.1111/ 1755 - 6724 .12303_ 28 (in Chinese).
Zhao HP . The comparison of four operators in synthetic fuzzy judgment of environment . Environ Technol Guizhou . 2001 ; 7 ( 3 ): 28 - 35 (in Chinese).