The Architectural Surfaces Characteristics of Sandy Braided River Reservoirs, Case Study in Gudong Oil Field, China

The architecture analysis of the different orders sedimentary bodies is of great significance to the efficient development of oil and gas fields. In order to investigate the effects of the architectural interfaces on reservoir quality and heterogeneity, this study takes the Gudong oil field as a case to investigate the architectural characteristics of hierarchical bounding surfaces using detailed descriptions of core and wireline logs. Architectural models from the 7th-order to the 3rd-order are analyzed, and the developmental characteristics of the 5th-order braided river, 4th-order single sandstone, and 3rd-order accretion are summarized. The interlayer between two braided rivers is floodplain mud deposition, with poor physical properties, stable thickness, and strong blocking capacity. Two models of interlayers are found in the 4th-order deposition. The first interlayer is between the braided filling channel and midchannel bar, which is composed of generally fine-grained sediments with calcium cementation and poor physical properties. The second interlayer is a transformation belt between two midchannel bars and is generally composed of gravel-scoured deposition with penetration capability. The 3rd-order surfaces are defined as the surfaces of accretions within midchannel bars. Two models of interlayers are also found in the 3rd-order surfaces of accretions. The paleocurrent of the sandy braided river is reconstructed by synthesizing the core data, well logging data, and production performance data. A total of 1 fluvial system (7th-order), 2 compound braided rivers (6th-order), 11 braided rivers (5th-order), 41 midchannel bars (4th-order), and 96 accretions (3rd-order) are developed in the study area.

Many scholars have carried out a series of researches on braided rivers through the observation of modern deposits and outcrops. Yu et al. [14] studied the braided river outcrop section and systematically described the sedimentary characteristics, lithofacies association, hierarchical interface, and sedimentary model of braided river reservoir. Jin et al. [16] divided the channel sand bodies into three types, which include superimposed, lateral, and isolated. Rucsandr (2001) applied GPR technology, core, and outcrop studies to dissect fluvial sand bodies. The radar profile is used to divide the reservoir architecture units.
The architectural analysis mainly includes two aspects: (1) the spatial contact relation and internal structure of architectural elements and (2) the distribution law and development characteristics of architectural bounding surfaces [17][18][19]. The architectural element analysis of reservoirs and sedimentary architectural patterns of sandy braided rivers have currently achieved certain successes [2,[20][21][22][23]. However, the sedimentary structural characteristics, architectural pattern, and hierarchical interfaces characteristics of braided river reservoirs have not been systematically studied due to the complexity and diversity of cutting and superposing relationships for underground reservoirs [14,24,25].
Likewise, the paleocurrent reconstruction of underground reservoirs in dense well-network pattern areas has not been attempted ( [26,27] [28]). Thus, the sedimentary characteristics of architectural elements and hierarchical bounding surfaces are systematically analyzed, and the channel and barform architecture of a sandy braided river is reconstructed using Miall's classification method for the hierarchical architecture of fluvial channels combined with the core data, well logging data, and dynamic data of dense drilling well areas. The architectural characteristics of underground reservoirs should be investigated due to their important theoretical significance in exploration and exploitation.
In this paper, sedimentary environment and sedimentary characteristics were studied by core description and analysis. The architecture interface and architecture units were identified using wireline logs, and the characteristics of different architecture units are further analyzed. Finally, the direction of the paleocurrent was determined.

Geological Setting
The Gudong oil field is a near-north-south-trend anticline structural reservoir in the northeastern part of the Jiyang Depression and the northern edge of the Yellow River estuary in eastern China ( Figure 1). The reservoir is mainly composed of loose sandstones from the Paleogene Guantao Formation. At the end of the Paleogene, the Jiyang Depression was lifted by the Dongying movement, thereby stripping and flattening the top of the Dongying Formation. During the Neogene, the Guantao Formation began to develop extensive fluvial deposits due to the height difference, abundant provenance, and braided river development. Then, the terrain height difference became small, and the sedimentary environment gradually shifted from braided to the meandering river when the upper section of the Guantao Formation began to deposit. The main oil-bearing formation is the upper section of the Guantao Formation. In this study, the 6 3+4 layer of the Guantao Formation from the Neogene is used as an example; this reservoir has high porosity and permeability, and its buried depth is 1220-1340 m ( Figure 2).
The dense well areas are in the east wing of the Gudong oil field. The north, west, and south sides are bounded by faults. The terrain slope is gentle, and the dip angle is generally less than 3°. The study area has 355 wells, which mainly use the antinine point well pattern and have approximate distances between wells of 30-150 m. Furthermore, 6 sealed coring wells and 3 tracer monitoring test well groups covering 31 injection and production wells are situated in the study area.

Methods
A total of 355 drilling wells and 6 coring wells are found in the dense well pattern area. All drilling wells have conventional logging curves with spontaneous-potential (SP), Gamma (GR), caliper log (CAL), sonic (AC), neutron (CNL), density (DEN), true resistivity (RT), and induced conductivity (COND). From the coring wells, we use 100 samples for performing several tests, such as granularity analysis, porosity and permeability analysis, rock thin section, scanning electron microscope, and conventional mercury injection test.
Various scaled architectural elements and the hierarchical bounding surfaces are identified and classified refinedly by observing the cored-well, single-well recognition, and interwell prediction of dense well pattern area (Li et al., 2011; [16,29]). The 7th-order to 3rd-order architectural models are analyzed, and the developmental characteristics of the 5th-order braided river, the 4th-order single sandstone (single midchannel bar and single-braided filling channel), and 3rd-order accretion are summarized by using Miall's classification method of architectural elements and hierarchical bounding surfaces.
The lithofacies variation of various scaled architectural elements and hierarchy of bounding surfaces are summarized, and the genetic mechanism of the hierarchical bounding surfaces is analyzed based on the cored-well data. The penetration capabilities of various hierarchical bounding surfaces are compared according to the physical properties of reservoirs obtained in an analytical laboratory test.
The well logging and production performance data of the dense well pattern area are synthesized to reconstruct the paleocurrent of the sandy braided river in the west 7th block of the Gudong oil field. In the 29-J254 coring example, the sample points 6 31 , 6 32 , 6 41 , and 6 42 are considered the cumulative probability curve and C-M diagram ( Figure 4). The particle size probability curve and C-M diagram show that the cumulative curve of grain size probability comprises two steps. The overall suspension of the upper and lower parts of the overall jump passes the overall slope of line 70-80, thus reflecting good sorting; on the C-M diagram, the sample points are mainly concentrated in the QR section, thereby showing positive rhythm characteristics of the river and onshore tractive current. In summary, the Guan 6 3+4 position has an irregular binary structure, and the underlying sand deposits are thick.

Geofluids
The study area has the same characteristics as a typical sandy braided river and is subdivided into braided filling channel, midchannel bar, and floodplain microfacies. Guan 6 3 is mainly composed of the braided filling channel and floodplain deposits, and Guan 6 4 mainly has the midchannel bar sedimentary microfacies.

Architectural elements.
The sedimentary environment of the study area is a typical sandy braided river. The braided river is divided into 5th-order braided river, 4th-order single sandstone (single midchannel bar and single-braided filling channel sand bodies), and 3rd-order accretion, as categorized by Miall's classification of fluvial architecture. This study focuses on the 4th-order architectural elements. The main 4th-order architectural elements are the midchannel bar and braided filling channel. Based on core description and logging data, seven kinds of midchannel bar deposits and four kinds of braided filling channel deposits are identified ( Figure 5).
The midchannel bar mainly consists of box, compound box, gearbox, high-amplitude bell type, compound bell, com-pound round, and funnel types. Fine sandstones dominate the box type, and the SP and COND curves show a single box-like shape with a high amplitude, a thickness of 4~8 m, and a ratio of 41.8% for the seven types. In the compound box type, the interlayer of the silty mudstone can be found in fine-sand sandstones. The curve shows many sets of box features; the amplitude is high, the thickness is larger than the box shape, which is 4~13 m, and the proportion is 22.9%. The gearbox profile is in the shape of a box. However, the SP has a tooth-like shape and has a noticeable return, thus reflecting the apparent distribution of the silty mudstone interlayer, which accounts for approximately 5.5% of all types. In the high-amplitude bell type, where the lower part is fine sandstone and the upper part is siltstone (silty    5 Geofluids mudstone), the SP and COND curves show a bell-like shape with a thickness of 3-6 m, the curve value ranging from the bottom to the top is weak, and the proportion is approximately 20%. In the compound bell type, a high overall amplitude curve and a large thickness are shown by a group of bellshaped combinations with a thickness of 7~12 m, which accounts for approximately 5.3%. In the compound round type, fine sandstone-siltstone and pebbly sandstone can be found in the top and bottom parts. The overall amplitude of the curve is high, the thickness of the sand body is large, the upper part is bell-shaped, the lower part is funnel-shaped, thereby corresponding to the compound rhythm, and the thickness is 6~11 m, which accounts for approximately 2.8% of all types. Compared with the bell-shaped type, the funnel type has a larger upward curve; furthermore, the upper part is fine-grained sandstone, and the bottom is pebbly sandstone. The curve shows a funnel-like shape and a thickness of 4~7 m, which accounts for 1.7% of all diagram types.
The braided filling channel comes in four types: sand, sand mixed with thin mud, sand at the bottom and mud at the top, and mud filling. In the sand filling type, the SP and COND curves show bell-shaped features; the logging range is high but weaker than that of the midchannel bar. The thickness is 2~5 m. This type accounts for 38.2% of the total channel filling. In the sand mixed with thin mud filling type, the curve is toothed bell-shaped, the thickness is approxi-mately 3~5 m, and the proportion is 14.5%. In the sand at the bottom and mud at the top filling type, the curve is low bell-shaped, and the top amplitude is low. In this type, the thickness is 2~5 m, which accounts for 19.5% of the total. The curve of mud filling is not an apparent finger or lowtooth profile, the amplitude is low, the thickness is 1~3 m, and the proportion is 27.8% of the filling channels.

Hierarchical Bounding Surface Identification.
A hierarchical architectural interface has a certain sedimentologic significance that reflects the heterogeneity of reservoirs in different scales and has several identification marks in core and logging data. Miall's classification method of architectural elements and hierarchical bounding surfaces is used for the sandy braided river. The 5th-order surfaces are defined as the surfaces of the braided river. The 4thorder surfaces are defined as the surfaces of braided filling channels or midchannel bars. The 3rd-order surfaces are defined as the surfaces of accretions in midchannel bars ( Figure 6).
A one-to-one correspondence is observed between architectural bounding surfaces and architectural elements, and the 5th-order bounding surface is the boundary of the single-braided river. Stabilized floodplain mudstone deposits generally exist in the braided river at different periods. Thus, the floodplain is an important symbol for judging the 5thorder bounding surface, reflected by the logging curve. The   Figure 6: Lithological characteristics of hierarchical bounding surfaces. FP: flood plain; CH: channel; DC: channel dam. 6 Geofluids SP curve shows a clear baseline segment, and the ML logging curve also has noticeable returns. The 4th-order architectural bounding surface is the boundary of the single midchannel bar or braided filling channel sediments. Midchannel bars are usually deposited in the middle of the river because of the special bidirectional circulation of braided rivers, with the hydrodynamic flow weakening in the later period, the gullies on two sides of the midchannel bar eventually filled with sandy or muddy sediments, which are called braided filling channel sediments. The symbols for judging the 4th-order architectural bounding surface are the connected interfaces of midchannel bars and braided filling channel sediments called bar-channel transformation interfaces. However, midchannel bars may migrate laterally or the braided filling channel sediments may be completely scoured by a subsequent river. A midchannel bar may directly be adjacent to the latter midchannel bar. This kind of 4th-order architectural bounding surface is called a bar-bar transformation interface. Fine-grained argillaceous siltstones usually exist at the bar-channel interface, with calcium cementation. The SP curve eventually shows a noticeable return. Meanwhile, the ML logging curve has noticeably high values due to calcium cementation deposition. The SP and ML wireline logs are reflected as relatively slight returns because the bar-bar interface generally has gravel-scoured erosion sediments deposited at the bottom of the latter midchannel bars.
The 3rd-order architectural bounding surface is the boundary of the single accretion. The end of one accretion and the start of the next generally reflect the direction and strength changes of the hydrodynamic circumstance. Two types of interfaces are found between the two accretions. The first interface is a thin muddy siltstone deposition, with an increased clay content weakened hydrodynamic. The second interface comprises sedimentary bedding structures that reflect changes in hydrodynamic conditions. The 3rd-order architectural bounding surface has a small effect on reservoir heterogeneity, and the amplitude that responds on logging carves is also relatively weak. The responding difference on the SP logging is too weak to identify the 3rd-order architectural bounding surface and is frequently accompanied by a slight return of the GR and ML curves.
In general, senior interfaces limit the deposition of junior interfaces. Thus, the architectural interface analysis is hierarchical. The 4th-order interface is analyzed in the 5th-order interface, and the 3th-order interface is analyzed in the 4thorder interface. However, the interface usually records an erosion event for fluvial deposits. That is, the sediments of a new period cut through and erode the sediments of an old period. Thus, the analysis of the hierarchical interface of underground braided river reservoirs is based on the characteristics of cutting and erosional channels. For example, if the top surface of a midchannel bar is a 4th-order interface, then the bottom surface is generally recognized as a 5th-order interface located in a braided river.

Geofluids
interfaces have different effects on reservoir quality and heterogeneity. An interface is considered a physical boundary. We select three parameters for investigating their effects on reservoir heterogeneity (Figure 7). The first parameter, which is called k0, is the interlayer's porosity and permeability between two architectural elements. The second parameter, which is called k, is the porosity and permeability of architectural elements. The third parameter, which is called Δk, is the variation of permeability on the two sides of the interface and represents the rate of change in permeability, which is expressed as Δk = ðk − k0Þ/k.
The floodplain mudstones among the 5th-order interfaces are generally deposited in the intermittent period of the river. The hydrodynamics is weak, and the mudstones develop stably with a large thickness and a wide range [30]. The maximum, minimum, and average thicknesses are 5.5, 0.7, and 1.69 m, respectively. Floodplain mudstones with extremely small porosity and permeability cannot effectively store and transport fluid; the maximum, minimum, and average k0 is 736:6, 0, and 67:1 × 10 −3 μm 2 , respectively, and Δk is comparatively large.
An interlayer that is different from the architectural elements exists among the 4th-order hierarchical interfaces. The interlayer among the bar-channel interfaces is in the middle of the midchannel bar and the braided filling channel sediments. The interlayer is usually a fine-grained interlayer or argillaceous siltstone with poor petrophysical properties and calcium cementation. At the end of the midchannel bar 25 Figure 9: Identification profiles of single midchannel bar. 9 Geofluids deposition, the hydrodynamics on the two sides of the bar is weakened, thus, gradually filling the fine-grained sediments. The maximum, minimum, and average thicknesses of the interlayer among the bar-channel interfaces are 0.6, 0.4, and 0.52 m, respectively.
Meanwhile, fine-grained sediments have a certain penetrability; the maximum, minimum, and average k0 is 1212:5 , 196:7, and 693:2 × 10 −3 μm 2 , respectively, and Δk is large but smaller than that of the 5th-order interfaces. Second, the interlayer among the bar-bar interfaces is in the middle of two midchannel bars. However, several differences are found in sedimentation due to the differences in hydrodynamic circumstances between the two bars. The interlayer among the bar-bar interfaces frequently has gravel-eroded sediments that affect the petrophysical properties and heterogeneity of the reservoir. The maximum, minimum, and average thicknesses of the interlayer are 0.7, 0.4, and 0.54 m, respectively. The maximum, minimum, and average k0 is 1374:7, 1271:2, and 1323:4 × 10 −3 μm 2 , respectively, and Δk is comparatively small.
The interlayer among the 3rd-order hierarchical interfaces (or between two stages of accretion) is commonly recognized as the fall silt seam or trench. The interlayer is thin-layered muddy stripes, thin, and generally poor continuity. The maximum, minimum, and average thicknesses of the interlayer are 0.6, 0.37, and 0.43 m, respectively. The maximum, minimum, and average k0 is 2678:8, 2413:1, and 2546:7 × 10 −3 μm 2 , respectively, and Δk is comparatively small. A 3rd-order interlayer with cross-bedding variation that reflects a slight change in hydrodynamics without muddy strips is found between two accretions. The thickness of the interlayer among the interfaces is relatively thin, thus indicating the short intermittent period of the two stages for hydrodynamic changes. The lithology of the sediments is mainly fine-grained muddy sandstones. The maximum, minimum, and average thicknesses of the interlayer are 0.5, 0.21, and 0.33 m, respectively. The interlayer does not affect the petrophysical properties and heterogeneity. The maximum, minimum, and average k0 is 2784:6, 2236:5, and 2610:1 × 10 −3 μm 2 , respectively, and the value of Δk is invariant.

Compound Bar
Interpretation. The paleocurrent of the sandy braided river in the west 7th block of the Gudong oil field is reconstructed by synthesizing the core, logging, and production performance data of the dense well pattern areas, using Miall's classification method of architectural elements and hierarchical bounding surfaces. The developmental characteristics of the braided river, single midchannel bar, and single accretion are analyzed quantitatively. The width and  thickness of the braided river are ranging from 408 to 1296 m, and from 2.9 to 11.7 m, respectively. The width-tothickness ratio is from 110 to 167. The width and thickness of the single midchannel bar are from 196 to 337 m and from 4.0 to 7.3 m, respectively. The width-to-thickness ratio is 38-51. The thickness of the single accretion is 1.0-1.5 m. The paleocurrent is reconstructed by dividing the compound midchannel bar into single midchannel bars. Sand body deposits with multiperiod single midchannel bars are then superimposed to form a compound bar.
In the thick Sixth Member sandstone of the Guantao Formation in the west 7th block in Gudong oil field, for example, according to the braided filling channel sand bodies among the different bars, the different thicknesses of adjacent bars, different changes in sedimentary rhythm, unmatched interlayers in different bars, and different logging curve characteristics and dynamic data reveal disconnected signs and identify single bars. We identified four single midchannel bars (Figures 8 and 9) and divided them into two different depositional types, simultaneously with different position types and the same layers at different times. In Figure 8, the 1st-, 2nd-, and 3rd-time depositions are bars 1, 2, 3, and 4, respectively. Bars 1 and 2 belong to the depositional type at the same time but in different positions, while bars 2 and 3 belong to the depositional type with the same layer but at different times. A total of 1 fluvial system (7th-order), 2 compound braided rivers (6th-order), 11 braided rivers (5thorder), 41 midchannel bars (4th-order), and 96 accretions (3rd-order) are developed in the study area ( Figure 10).

Conclusions
In this paper, the architecture interface and architecture units were identified using wireline logs, based on the understanding of the sedimentary environment and sedimentary characteristics. And the characteristics of different architecture units are further analyzed combined with wireline logs' shape and core photos. The following conclusions were obtained: (1) The Guan 6 3+4 layer of the study area is mainly typical sandy braided river deposits. The braided river deposits can be subdivided into the braided filling channel, midchannel bar, and floodplain microfacies (2) Sand body deposits with multiperiod single midchannel bars are superimposed to form a compound bar. Four single midchannel bars were identified and then divided into two different depositional types, simultaneously with different position types and the same layers at different times (3) The interlayer between two braided rivers is floodplain mud deposition, with poor physical properties, stable thickness, and strong blocking capacity. The interlayer between the braided filling channel and the midchannel bar is generally fine-grained sediments with calcium cementation and poor physical properties. The interlayer, between two midchannel bars, is a transformation belt and generally composed of gravel-scoured deposition with penetration capability.

Data Availability
All data can be obtained from the corresponding authors.

Conflicts of Interest
The authors declare that they have no conflicts of interest.