Study on Hydrocarbon Accumulation Periods Based on Fluid Inclusions and Diagenetic Sequence of the Subsalt Carbonate Reservoirs in the Amu Darya Right Bank Block

In order to provide paleo ﬂ uid evidence of hydrocarbon accumulation periods in the Amu Darya Right Bank Block, microexperiments and simulations related to the Middle-Upper Jurassic Callovian-Oxfordian carbonate reservoirs were performed. On the basis of petrographic observation, the diagenetic stages were divided by cathodoluminescence, and the entrapment stages of ﬂ uid inclusions were divided by laser Raman experiment and UV epi ﬂ uorescence. The hydrocarbon generation (expulsion) curve and burial (thermal) history curve of source rocks were simulated by using real drilling data coupled with geochemical parameters of source rocks, such as total organic carbon (TOC) and vitrinite re ﬂ ectance ( R o ). The above results were integrated with microthermometry of ﬂ uid inclusions by inference the timing of hydrocarbon migration into the carbonate reservoirs. The horizon- ﬂ attening technique was used to process the measured seismic pro ﬁ le and restore the structural evolution pro ﬁ le. Four diagenetic periods and three hydrocarbon accumulation periods were identi ﬁ ed. (i) For Syntaxial stage, the ﬂ uid captured by the overgrowing cement around particles is mainly seawater; (ii) for (Early) Mesogenetic burial stage, the calcite cements began to capture hydrocarbon ﬂ uids and show yellow ﬂ uorescence under UV illumination; (iii) for (Late) Mesogenetic burial stage, two sets of cleavage ﬁ ssures developed in massive calcite cements, and oil inclusions with green ﬂ uorescence were entrapped in the crystal; (iv) for Telogenetic burial stage, blue ﬂ uorescent inclusions along with hydrocarbon gas inclusions developed in dully luminescent calcite veins. Based on the accurate division of hydrocarbon migration and charging stages, combined with the structural evolution history of the traps, the hydrocarbon accumulation model was established. Because two of the three sets of source rocks are of marine origin, resulting in the lack of vitrinite in the kerogen of those source rocks, there may be some deviation between the measured value of R o and the real value. Some systematic errors may occur in the thermal history and hydrocarbon generation (expulsion) history of the two sets of source rocks. Due to the limitations of seismic horizon- ﬂ attening technique — such as the inability to accurately recover the inclined strata thickness and horizontal expansion of strata — the ﬁ nal shape of the evolution process of structural pro ﬁ le may also deviate from the real state in geological history. The accumulation model established in this study was based upon the ﬂ uid inclusion experiments, which can e ﬀ ectively characterize the forming process of large condensate gas reservoirs in the Amu Darya Right Bank Block and quantify the timing of hydrocarbon charging. However, the hydrocarbon migration and accumulation model does not take the oil-source correlation into account, but only the relationship between the mature state of source rocks and the timing of hydrocarbon charging into the reservoirs. Subsequent research needs to conduct re ﬁ ned oil-source correlation to reveal the relationship between gas, condensate, source rocks, and recently discovered crude oil and more strictly constrain and modify the accumulation model, so as to ﬁ nally disclose the origin of the crude oil and oil reservoir forming process in the Amu Darya Right Bank Block, evaluate the future exploration potential, and point out the direction of various hydrocarbon resources (condensate gas and crude oil).


Introduction
The Amu Darya Right Bank Block (the block is denoted by ADRBB in the following text) project is the largest one for natural gas cooperation between China and Turkmenistan. According to the existing data and knowledge, the source rocks of the widely distributed condensate reservoirs in the ADRBB are primarily gas prone. In 2020, some reservoirs that mainly produce normal crude oil were discovered around large condensate gas reservoirs which produce condensate oil (this paper will be the first to report the discovery of crude oil in the ADRRB), indicating that we have not fully understood the regular pattern of hydrocarbon accumulation in the ADRBB. To ascertain the laws regarding the formation and distribution of hydrocarbon reservoirs in the ADRBB, there is an urgent need to restore and finely describe the formation and evolution process of reservoirs which should be regarded as individual units. We conduct various microscopic experiments such as cathodoluminescence of the Middle-Upper Jurassic Callovian-Oxfordian carbonate reservoirs in the ADRRB (the reservoir is represented by J 2 k-J 3 o in the following text) and the microthermometry and Raman spectrum of fluid inclusions of the J 2 k-J 3 o reservoirs. The experimental results are integrated with the burial (thermal) history to determine the reservoir charging period and stages. Then, based on the hydrocarbon generation (expulsion) history of source rocks and the structural evolution history of the strata, the forming process of hydrocarbon reservoirs is restored (Figure 1). The subsequent exploration and development of hydrocarbon reservoirs can be carried out on the basis of the formation process established in this paper.
Hydrocarbon charging and reservoir diagenesis are simultaneous geological processes in the petroleum system, and the constructive and destructive actions of diagenesis on the reservoir determine whether there is enough space to accommodate oil and gas. The various processes of porosity generation, preservation, and reduction "compete" with one another, and the present porosity and permeability characteristics represent the net result of all past processes [1]. Hence, to study the formation period of oil and gas reservoirs, it is necessary first to divide the diagenetic process of reservoirs in a certain area in detail. Constructive diagenesis includes dissolution, dolomitization [2], and tectonic faulting. Destructive diagenesis includes compaction, cementation, stylolitization, and recrystallization. Other diagenesis, such as hydrocarbon generation overpressure in shale, also belongs to positive diagenesis. The generated gas will make the laminated shales open more in the bed-normal direction, which favors incipient vein opening between laminae, even in the absence of deviatoric stress [3]. Veins filled with brittle minerals, such as quartz or calcite veins, are more conducive to the generation of fractured media, which is highly significant for the production of shale gas [4] and the immiscible two-phase flow (in oil zone with a gas cap) displacement [5]. Bacterial sulfate reduction (BSR) and thermochemical sulphate reduction (TSR) have the nature of a double-edged sword for reservoir reconstruction. On the one hand, these two reactions will dissolve sulfate, produce acid gas (H 2 S, CO 2 ), and form secondary pores. On the other hand, solid bitumen may form as a by-product of both BSR and TSR to block pores. Therefore, for blocks with these two redox reactions, the effect of BSR/TSR on reservoir reconstruction needs to be comprehensively evaluated [6]. Cathodoluminescence is a powerful means to determine diagenetic minerals and formation stages, because different minerals precipitated in the same period, or the same mineral precipitated in different periods, show different cathodoluminescence properties [7]. The microthermometry and Raman spectrum of secondary fluid inclusions can provide direct data for the timing of the entrapment, and the structural recovery from real data of traps can show the formation process of hydrocarbon reservoirs specifically and help delimit the migration pathways of oil and gas [8].

Geology Overview
The Amu Darya Basin, located in eastern Turkmenistan and southwestern Uzbekistan, is a large gas-rich foreland basin in Central Asia. It covers an area of 42:70 × 10 4 km 2 and has proven recoverable oil reserves of 1:7 × 10 8 t and proven recoverable natural gas reserves of 22 × 10 12 m 3 [9]. The Amu Darya Basin has been in the northern margin of the Tethyan tectonic domain since the Late Paleozoic. The basin has mainly experienced the Late Paleozoic-Early Mesozoic basement development period, the Permian-Triassic rifting period, the Jurassic-Paleogene Eocene postrifting depression period, and the Neogene Miocene-Quaternary foreland basin period [10]. In terms of the current basement undulation and sedimentary caprock structure state, the Amu Darya Basin is divided into structural units such as Kopet-Dag piedmont depression, Chardzhou terrace, Carlabekaul depression, Bieshikent depression, Murgab depression, and central Karakum uplift. The ADRBB is located in the northeast of Amu Darya basin, covers an area of 1:43 × 10 4 km 2 and traverses 3 secondary structural units such as Chardzhou terrace, Bieshikent depression, and southwest Gissar thrust belt [11,12]. The ADRBB can be further divided into 6 structural units such as Chardzhou uplift, Kennedykiddskurt uplift, Sandykly uplift, Carlabekaul depression, Bieshikent depression, and southwest Gissar piedmont thrust belt (Figure 2(a)).
The ADRBB has experienced the evolution process of continental-marine-lagoon facies. After being uplifted in the Late Triassic, the ADRBB reentered a subsidence stage in the Early-Middle Jurassic. The depositional environment of the ADRBB gradually changed from delta to marine facies. The Middle-Lower Jurassic coal-bearing mudstones with a basin scale deposited in this period in the Amu Darya Basin are the most important set of source rocks in the basin [13]. The ADRRB experienced a large-scale transgression in the Callovian period of the Middle-Late Jurassic, and a Callovian gentle slope carbonate platform developed on the basis of the formed accommodation space in this period. During the Oxfordian period of the Middle-Late Jurassic, the ADRRB suffered another large-scale transgression. The platform was generally submerged, and it retreated, and a rimmed shelf-type carbonate platform formed. The two 2 Geofluids stages of platform facies deposition in the ADRBB have resulted in the development of superimposed intraplatform mound-beach carbonate reservoirs in the western section and the development of longitudinally aggradational gentle-slope reef-beach carbonate reservoirs in the mideastern section [14,15]. The deposited carbonate reservoirs are mostly argillaceous limestones with high organic matter content and are a set of main source rocks in the Jurassic in the ADRBB. As the seawater grew deeper, a deep mudstone layer (Gap layer) with high gamma-ray value developed in the mideastern section of the ADRBB at the end of the Oxfordian period. The thickness of the Gap layer is mostly less than 20 m; it is a set of secondary source rocks of the Jurassic in the ADRBB. That is, working upwards, there are three sets of source rocks in ADRBB: coal-bearing mudstone, argillaceous limestone, and the Gap layer. At the same time of the development of the Gap layer, the ADRBB began to retreat on a large scale, and the reefs disappeared in most regions. In the Late Jurassic Kimmeridgian period, the ADRBB began to enter the gypsum deposition stage, and the carbonate platform gradually degenerated into a gypsum platform. From the Late Jurassic Kimmeridgian period to the Tithonian period, the ADRBB experienced multiple periods of transgression and regression, where thick "three anhydrite units and two salt units" layers were deposited [16]. Three sets of source rocks and J 2 k-J 3 o carbonate reservoirs, together with the ultrathick anhydrite-salt layer, constitute a complete set of source-reservoir-caprock assemblages of continuous sedimentation from bottom to top. The excellent hydrocarbon accumulation assemblage provides the material basis for the formation of large-scale subsalt hydrocarbon reservoirs in the ADRBB [17,18] ( Figure 2(b)). After ten years of exploration and development of the ADRBB, more than 40 hydrocarbon reservoirs have been discovered, which are distributed in the west, middle, and east sections and are mainly natural gas reservoirs. The res-ervoirs are substantially structural traps controlled by anticlines and faults, followed by structural-lithological traps (Figure 2(c)). The western reservoirs are pore-vug type. These are multilayer large-scale intraplatform beach bodies, with secondary pores and residual primary pores as the dominating reservoir space and auxiliary solution cavities. The middle reservoirs are fracture-pore type. These are extensive gentle-slope reef-beach groups, and the reservoir space is dominated by secondary dissolution vugs and fractures. The east reservoirs are fracture-vug type. These are massive fracture-vug carbonate reservoirs related to thrust structures, and the reservoir space is fractures; vugs dissolved along fractures and irregular solution holes. Fractures and faults are pervasive in the subsalt reservoirs and play a significant role in improving the migration pathways and reservoir space. The temperature of the primary pay formation is 100°C-140°C; the salinity of NaCl of formation water is 3-7% [19], and the pressure factor ranges from 0.8 to 2.0; the composition of the well fluids is dominated by CH 4 , and their average volume fraction is below 95%, showing normal/high pressure wet gas reservoirs with a normal temperature system.

Materials and Methods
The core samples were collected from the major gas-bearing horizons of J 2 k-J 3 o in the 7 high yield wells. Wells AWS-5 and AWM-2 are located in the west section of ADRRB, while wells AMP-1, AMO-2, and AMB-1 are located in the middle section. Wells AEA-1 and AEA-2 are located in the east section (Figure 2(a)). Core samples were made into standard double-sided polished slices for inclusion observation.
Petrographic observation was performed using Zeiss AxioCamERc5s microscopes. The vacuum of the cathodoluminescence instrument is 0.3 Pa. The voltage and current of the electron beam are 7.5 kV and 0.63~0.76 mA, respectively. The microthermometry data were obtained using Linkam  3 Geofluids THMS-600 cooling/heating stage, with the temperature measurement range of -196~600°C, and the measurement error of ±0.1°C. Raman spectra were collected by using a JY/Horiba LabRam HR800 Raman system equipped with a frequency doubled Nd: YAG laser (532.06 nm), whose output laser power is 45 mW and a 50× longworking-distance Olympus objective with a 0.5 numerical aperture. Spectra were collected in the range of~100-4000 cm −1 . The acquisition time was approximately 1 s, with two accumulations for each spectrum to maintain a high signal-to-noise ratio. Before Raman testing, the instrument was calibrated with monocrystalline silicon whose Raman shift should be 520.70 cm -1 .
In order to gain a deeper understanding of accumulation period of the ADRRB, we used the PetroMod 2016 software to recover the hydrocarbon generation and expulsion history of source rocks. At the same time, the stratigraphic burial and thermal history of source rocks and structural evolution history were restored. The concepts and rules of restoration are based on the methods discussed by Hantschel and Kauerauf [20]. The maturity data of source rocks are obtained from the measured vitrinite reflectance (R o ). The present strata thickness is obtained from drilling and logging data. The total organic carbon (TOC) and Hydrogen Index (HI) data are obtained from pyrolysis testing. The thermal history curves are restored according to the geothermal gradient of 4.0°C/ 100 m in the uplift area and 2.5°C/100 m in the depression area and the surface temperature of 20°C, and the thermal history curves are constrained and corrected by R o .   5 Geofluids may develop in intraclast and bioclast. However, the fluid inclusions developed in these two kinds of particles are primary inclusions, not secondary inclusions trapped during burial diagenesis. So we do not focus on these three carbonates.

Results and Discussion
According to the luminescent characteristics and relative positions, the sparry calcite cements that can develop secondary fluid inclusions can be divided into four stages.
(i) Syntaxial Cements. It is difficult to identify syntaxial overgrowth cement in plane-polarized light. Syntaxial overgrowth cements keep in optical continuity with their substrate under the condition of cathodoluminescence [21]; that is, they are nonluminescent. The In short, the cement in reservoir pores is calcite which can be determined by cathodoluminescence and petrographic analysis. We have defined four diagenetic periods. The diagenetic periods are (i) Syntaxial stage, (ii) (Early) Mesogenetic burial stage, (iii) (Late) Mesogenetic burial stage, and (iv) Telogenetic burial stage. Figure 4 is a comprehensive sketch of Figure 3 showing all diagenetic stages, in which fluid inclusion assemblages developed in sparry calcite cements are also shown. The stages and corresponding characteristics of fluid inclusions developed at different locations will be discussed in Section 4.2. The three fluorescence colors may indicate that ADRRB experienced at least three stages of oil and gas charging. Additionally, inclusions showing green and blue fluorescence are evidently less abundant than those showing yellow fluorescence, while the hydrocarbon fluid charging intensity is the dominating factor affecting inclusions abundance [25,26], indicating that the hydrocarbon charging scale with yellow fluorescent inclusions formed is the largest The homogenization temperature (T h ) of hydrocarbon inclusions is generally hard to detect due to the thermal instability of the organic matter and T h differences of each   When the formation is not saturated with gas, the T h of inclusions is often lower than the actual trapped temperature and needs to be corrected. When the formation is saturated with gas, the T h of inclusions can represent the trapped temperature. Since the Amu Darya Basin is rich in natural gas, it is assumed that the T h is tested when the formation is saturated with gas.
Using histograms enables effective visual display of the data on the thermal history of multiple phases. The T h histogram of fluid inclusions was plotted (Figures 6(a) and 6(b)) in accordance with the method of FIA and microthermometry data (Table 1). Figure 6(a) is the T h histogram of all aqueous fluid inclusions. Figure 6(b) is the T h histogram of aqueous inclusions which belong to the same FIA with hydrocarbon inclusions showing three fluorescent colors in Well AMO-2. It can be seen from Figure 6(a) that the prevailing T h of the fluid inclusions in the ADRRB is 90~140°C, with continuous distribution characteristics, indicating that on the whole the basin is in a continuous charging state. The main peak temperature is 100~110°C, reflecting the overall large-scale hydrocarbon charging of one stage. The T h of the inclusions in the samples from Well AMO-2 is 90~160°C, and the two peak temperatures are 100~110°C and 130~140°C (Figure 6(b)). The T h data of 100-110°C correspond to yellow fluorescent hydrocarbon inclusions, while T h data of 130-140°C correspond to green and blue fluorescent hydrocarbon inclusions. The T h of some fluid inclusions reaches 150-160°C, which exceeds the highest formation temperature in geological history (the strata have not experienced large-scale uplift, so the current formation temperature is the highest temperature in geological history). The reason may be that the deep hydrothermal fluids were trapped. Or when the fluid inclusions were formed, the fluids were not in a homogeneous phase state but immiscibly trapped, so that the T h measured is far higher than the formation temperature when the fluid inclusions were trapped [28,29].
The salinity of NaCl is determined by Bodnar's formula [30] using freezing point depression data (Figures 6(c) and 6(d)). It can be seen from Figure 6(c) that the T h and the salinity of these secondary inclusions show a certain degree of positive correlation. The T h of fluid inclusions in AMO-2 well in the range of 130-140°C can be divided into two periods by using salinity data (Figure 6(d)). The salinity of the fluid inclusions of the three stages is 2.90~8.41%, 8.14~12.73%, and 4.49~9.86%.

Characteristics of Laser Raman Spectrum.
Laser Raman spectroscopy can be used for qualitative identification of the components of minerals and fluid inclusions to provide evidence for petrographic analysis. The fluorescence intensity of liquid hydrocarbon is several orders of magnitude higher than that of Raman spectrum intensity. So the laser Raman experiments can only be carried out in gas inclusions and (two-phase) aqueous fluid inclusions.

Geofluids
The four stages of sparry cements classified by petrography were tested separately. The results show that the Raman shift of the highest peak of stage (i) (Figure 7(a)) and stage (ii) (Figure 7(b)) is 1086.01 cm -1 , and that of stage (iii) (Figure 7(c)) and stage (iv) (Figure 7(d)) is 1085.40 cm -1 . That is to say, the minerals of the sparry cements are calcites [31]. Although the contents of Fe 2+ and Mn 2+ can affect the properties of calcite cements under cathodoluminescence, the different contents of Fe 2+ and Mn 2+ in the calcite cement have little effect on the Raman shift of calcite.
The peaks of the host minerals will not be discussed in the subsequent contents. There are five peaks in the gas inclusion (Figure 8(a)). The Raman shifts of the peaks are 1270.66 cm -1 , 1370.89 cm -1 , 2322.85 cm -1 , 2911.93 cm -1 , and 2952.79 cm -1 . The corresponding compounds are CO 2 (1270.66 cm -1 and 1370.89 cm -1 ), N 2 , CH 4 , and C 2 H 6 , respectively.
There are a bump and a peak in the vapor phase of another two-phase aqueous fluid inclusion (Figure 8(c)). The Raman shift of the bump is 1311.24-1640.53 cm -1 , which indicates the vapor phase may contain asphaltine (D band of disordered graphite~1355 cm -1 ; G band of wellordered graphite~1575-1582 cm -1 ) and light hydrocarbon components such as C 6 H 14 (1451 cm -1 ) and C 7 H 16 (1302 cm -1 ) [31]. The Raman shift of the peak is 2912.02 cm -1 , which represents CH 4 .
There are a bump and three peaks in the liquid phase of the two-phase aqueous fluid inclusion (Figure 8(d)). The Raman shift of the bump is 3200-4000 cm -1 , which represents water. The Raman shifts of the three peaks are 1040.74 cm -1 , 1063.41 cm -1 , and 2916.32 cm -1 . The corresponding compounds (ions) are NO 3 -, CO 3 2-, and CH 4 , respectively. Through the identification of hydrocarbon gases, it is shown that the two-phase aqueous inclusions which were used in microthermometry were formed in the same period as hydrocarbon inclusions. The two-phase aqueous inclusions in the microthermometry were not syntaxial inclusions which were formed by the trapping of seawater. That is to say, the petrographic analysis results of FIA are verified, and the validity of microthermometry data is also confirmed.

Simulation of Hydrocarbon Generation and Expulsion
History. The hydrocarbon generation and expulsion history of three sets of source rocks in the ADRRB were simulated, so as to provide a basis for the judgment of the hydrocarbon accumulation period. The hydrocarbon generation and expulsion history curves of the three sets of source rocks in ADRRB are shown in Figure 9. The hydrocarbon expulsion time lags behind the hydrocarbon generation time because of the need to satisfy the adsorption by the pores of source rocks. The coal-bearing mudstone began to generate hydrocarbons at 170 Ma and to expel hydrocarbons at about 150 Ma, reaching the hydrocarbon expulsion peak in the Paleogene. The argillaceous limestone and Gap layer began to generate hydrocarbons at 103 Ma and to expel hydrocarbons at around 98 Ma. Compared with the coal-bearing mudstones, the argillaceous limestone and Gap layer have a lower level of hydrocarbon generation and expulsion. Before the Himalayan orogeny, the ADRRB was in a stable subsidence stage, and the rise in the strata temperature promoted the hydrocarbon generation and expulsion in the source rocks. The eastern section of the ADRRB was relatively deep and was the main hydrocarbon kitchen in the ADRRB. With the beginning of the Himalayan orogeny, the eastern section of ADRRB was strongly folded and uplifted, and the most important hydrocarbon kitchen in the ADRRB was destroyed, causing the overall hydrocarbon generation and expulsion in the ADRRB to slow down.
In the geological history period, Murgab depression in the southeast of the basin was located in a relative depression area of the basin; that is, the evolution of the Murgab source rocks was ahead of those in the ADRRB. So in the hydrocarbon expulsion stage of the source rocks in the ADRRB, the oil and gas from Murgab depression had already migrated to the ADRRB to charge the traps.

Hydrocarbon
Charging Stages. The studies summarized above show that there are numerous approaches that can be taken to obtain accumulation periods from microscopic analysis of the hydrocarbon reservoirs. The combination method of comprehensive analysis of reservoir charging stages is as follows: (i) Diagenetic stages are divided by cathodoluminescence; (ii) FIAs are identified according to petrographic characteristics; (iii) The entrapment stages of FIAs are divided by UV epifluorescence and microthermometry; (iv) Because the gas components do not emit fluorescence under ultraviolet illumination, the identification of  12 Geofluids hydrocarbon gas components in gas inclusions by laser Raman can provide evidence for hydrocarbon gas charging. The mineral types of cements can also be determined by laser Raman to verify the conclusions of cement identification by petrographic analysis and cathodoluminescence. Comparisons between fluorescence emission spectra of oil inclusions with API gravity and maturity of related oils show that there is a tendency for petroleum fluorescence to change from yellow to blue colors with increasing maturity [32]. By comprehensively analyzing the occurrence location and fluorescence characteristics of hydrocarbon inclusions, it can be confirmed that: (i) The yellow fluorescent FIAs were captured in the (Early) Mesogenetic burial stage. During this period, the reservoirs were gradually divorced from the seawater environment and began to capture the fluids in the formation, and the low mature hydrocarbon and formation water were sealed in the crystal (ii) The green fluorescent FIAs were entrapped in the (Late) Mesogenetic burial stage. Under the increasing formation temperature, the sparry calcite cements continued to crystallize, formed welldeveloped cleavage fissures, and trapped the maturer oil that had migrated into the reservoirs (iii) The blue fluorescent FIAs and gas inclusions containing CH 4 were captured during the Telogenetic burial stage. The tectonic movement made the reservoirs produce fractures which are filled with calcite precipitated from the formation water. The highly mature oil and hydrocarbon gas were sealed in the sparry calcite crystal The timing of hydrocarbon charging can be quantified by projecting the T h of FIAs onto the burial (thermal) history map [32][33][34]. According to the microthermometry data from Well AMO-2 ( Figure 6(b)), the homogenization temperature of yellow fluorescent FIAs is 100-110°C, and that of green and blue fluorescent FIAs is 130-140°C. Although green and blue fluorescent FIAs cannot be distinguished by T h , it can be determined that the entrapment of blue fluorescent inclusions occurred later than that of green fluorescent inclusions by using freezing point depression data, development location, and fluorescence property. The formation of large-scale yellow fluorescent FIAs occurred in 120-105 Ma (Early Cretaceous). At that time, the reservoirs were no longer affected by seawater. Complex physicochemical processes occurred between formation water and reservoirs to precipitate calcite cements from formation water and capture low mature hydrocarbons. The green fluorescent FIAs formed during 75-65 Ma (Late Cretaceous-Early Paleozoic). The rising formation temperature caused the newly precipitated calcite to develop into automorphic crystal with two sets of cleavage fissures. The extension direction of the trapped maturer hydrocarbon FIAs might have cut the cleavage fissures of the calcite. The blue fluorescent FIAs were entrapped in 30-20 Ma (Oligocene-Miocene). The reservoirs experienced the Himalayan tectonic movement, resulting in extensive fractures that provided accommodation for calcite precipitation from formation water. The migrated highly mature oil and hydrocarbon gas inclusions were sealed in calcite veins ( Figure 10).

Restoration of the Hydrocarbon Accumulation Process.
Based on the current seismic profile (the section direction is shown in Figure 2(a)) ( Figure 11), the structural evolution of the Sandykly uplift in the middle section of the ADRRB is restored by horizon-flattening technique. Combined with diagenetic stages, paleofluid properties, hydrocarbon generation (expulsion) simulation, and the division of oil and gas accumulation stages, the formation process of the hydrocarbon reservoirs in the middle section of the ADRRB can be restored ( Figure 12). 14 Geofluids In the Early Cretaceous, the J 1-2 coal-bearing mudstones in the Murgab depression in the southeast of the Amu Darya Basin entered the oil generation window and quickly reached the hydrocarbon expulsion peak stage. In the Albian Age (120~105 Ma), the oil and gas generated in the Murgab depression migrated along the unconformity surface between the Jurassic strata and the Paleozoic strata, as well as the basement faults of Chardzhou terrace, to the reefbeach bodies of the southern slope belt of the Sandykly uplift [15,35]. At this moment, the buried depth of the reservoirs in the middle section of the ADRRB was about 2500 m, and they entered the (Early) Mesogenetic burial diagenesis stage. Microscopically, the calcites precipitated and crystallized around particles and syntaxial cements, capturing the hydrocarbon that had migrated into the reservoirs and forming the yellow fluorescent FIAs. Macroscopically, the reservoirs were charged with oil and gas to form structurallithological hydrocarbon reservoirs (Figure 12(a)). Although the hydrocarbon generation and expulsion simulation shows that the coal-bearing mudstones of the ADRRB have reached the first hydrocarbon expulsion peak (Figure 9), the value of the hydrocarbon generation and expulsion curve   Figure 9: Hydrocarbon generation and expulsion curves of the source rocks in the ADRRB. 15 Geofluids is mainly contributed by the depression in the east of the ADRRB, which has evolved into the Gissar piedmont thrust belt (Figure 2(a)). The hydrocarbon expelled from the east depression accumulated nearby and had little effect on the middle section of the ADRRB. The buried depth of J 1-2 coal-bearing mudstones in the middle section of the ADRRB was less than 3000 m ( Figure 10); they were in a lowmaturity stage with a negligible quantity of hydrocarbon   Figure 10: Burial-thermal history curves of Well AMO-2 in the middle section of the ADRRB.  16 Geofluids generation and expulsion. Therefore, the oil and gas in this period mainly came from the Murgab depression ( Figure 12(a)).
During the Turonian Age of the Late Cretaceous (75-65 Ma), the Amu Darya Basin continued to subside. The buried depth of the reservoirs in the middle section of the ADRRB was about 3500 m, and they entered the (Late) Mesogenetic burial diagenesis stage. The massive calcite cements with two sets of cleavage fissures began to capture the hydrocarbon and formation water in the reservoirs to form green fluorescent FIAs. The source rocks in the Murgab depression expelled vast wet gas. The burial depth of the coal-bearing mudstones in the middle section of the ADRRB exceeded 4000 m, and they began to expel a lot of hydrocarbon. Driven by buoyancy, the hydrocarbon generated by the coal-bearing mudstones in the middle section of the ADRRB migrated upward directly and accumulated in the reservoirs together with the oil and gas that had migrated from the Murgab depression (Figure 12(b)). The argillaceous limestone and Gap layer of the ADRRB  subsequently entered a low-maturity stage in the east section of the ADRRB and expelled a small quantity of hydrocarbons ( Figure 9).
From the Paleogene-Oligocene to the Neogene-Miocene (30~20 Ma), the buried depth of the reservoirs in the middle section of the ADRRB was about 4000 m, and they entered the Telogenetic burial diagenesis stage. In terms of reservoir scale, a large number of fractures were produced in the middle section of the ADRRB affected by the compression of the Miocene Himalayan tectonic movement [18]. The newly formed spaces were filled with calcites precipitated from formation water. Blue fluorescent FIAs and hydrocarbon gas inclusions are sealed in these dully luminescent calcite veins. In terms of trap scale, the closure height of the traps in the Sandykly uplift increased slightly. The trap volume was increased in the Sandykly uplift. In terms of basin scale, the coal-bearing mudstones in Murgab depression entered the gas generation window and expelled dry gas; meanwhile, the argillaceous limestone and Gap layer in the middle section of the ADRRB entered the oil generation window. By now, the hydrocarbon reservoirs in the middle section of the ADRRB had been charged by four sources (Figure 12(c)). The source rocks were (i) J 1-2 coal-bearing mudstones of Murgab depression, (ii) J 1-2 coal-bearing mudstones of ADRRB, (iii) J 2 k-J 3 o argillaceous limestone of ADRRB, and (iv) the Gap layer of the ADRRB. Although the hydrocarbon reservoirs are sealed and protected by ultrathick anhydrite-salt layer [16], under strong gas dissipation [36], continuous charging is also essential for the preservation of the hydrocarbon reservoirs of the Amu Darya Right Bank Block. In the follow-up research, it will be necessary to carry out refined oil-source correlation for ADRRB, so as to determine the contribution of each source rock to the hydrocarbon reservoirs and further constrain and optimize the hydrocarbon accumulation model. (2) The ADRRB is in a continuous charging state on the whole. The main peak temperature of T h histogram is 100~110°C, representing the overall large-scale hydrocarbon charging stage. In the middle section of the ADRRB, the T h peak of 100~110°C corresponds to yellow fluorescent hydrocarbon inclusions, and the T h peak of 130~140°C corresponds to green and blue fluorescent hydrocarbon inclusions. The salinity of the aqueous fluid inclusions coexisting with yellow fluorescent hydrocarbon inclusions is 2.90~8.41%. The salinity of the green fluorescent FIAs is 8.14~12.73% and that of the blue fluorescent FIAs is 4.49~9.86%

Conclusion
(3) In the Albian Age of the Early Cretaceous (120~105 Ma), the reservoirs in the middle section of the ADRRB were charged with oil and gas migrated from the Murgab depression to form structural-lithological hydrocarbon reservoirs. During the Turonian Age of the Late Cretaceous (75-65 Ma), the hydrocarbon generated by the coalbearing mudstones in the middle section of the ADRRB migrated upward directly and accumulated in the reservoirs together with the oil and gas that had migrated from the Murgab depression. From the Paleogene-Oligocene to the Neogene-Miocene (30~20 Ma), the closure height of the traps in the Sandykly uplift increased affected by the compression of the Miocene Himalayan tectonic movement. The oil and gas generated from J 2 k-J 3 o argillaceous limestone and the Gap layer of ADRRB began to charge the local reservoirs. Because of the sealing of an ultrathick anhydrite-salt layer and the continuous charging from multiple sources, the hydrocarbon reservoirs formed in the Early Cretaceous can be well preserved up to now

Data Availability
The microthermometry data used to support the findings of this study are included within the article.

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