Eruption Cycle of Volcanic Rocks Based on Virtual Reality Technology

How to cite item Yang, J., & Wang, X. (2020). Eruption Cycle of Volcanic Rocks Based on Virtual Reality Technology. Earth Sciences Research Journal, 24(3), 277-284. DOI: https://doi.org/10.15446/ esrj.v24n3.89468 The eruption cycle of volcanic rocks based on virtual reality technology is studied to reveal the spatial distribution of volcanic rocks and the genetic relationship between different lithologic assemblages, and effectively guide volcanic exploration. Taking Xujiaweizi fault depression and Yingshan depression in Songliao Basin as the research area, the monocular vision method of virtual reality technology is used to collect data from the study area. Then, the threedimensional dynamic scene of the study area was constructed. Based on this, combined with the research steps and division basis of the volcanic eruption cycle, the eruption cycle of the volcanic region in the study area is studied. According to the characteristics of volcanic activity, volcanic facies sequence, and rock rhythm combination in Xujiaweizi fault depression, the first and three sections of Yingcheng formation are divided into three eruption periods according to the contact relationship between sedimentary strata and volcanic rocks in the eruptive intermittent period. Y1C1, Y1C2, Y1C3, Y3C1, Y3C2, and Y3C3 are respectively from bottom to top. Yingshan volcanic rocks in Yingshan depression are divided into three cycles and six periods. Each stage has significant stages and differences, and the eruption intensity is different. The overall performance is gradually enhanced, and the internal cycle is characterized by the first strong and then weak. ABSTRACT Eruption Cycle of Volcanic Rocks Based on Virtual Reality Technology


Introduction
With the development of oil and gas exploration, volcanic oil and gas reservoirs have become an essential part of oil and gas exploration. Volcanic rocks are complex in rock type, uneven in regional distribution, and fast in lateral lithology. Volcanic oil and gas reservoir exploration technology is a difficult point in petroleum technology research today (Zhu et al., 2017). The volcanic gas reservoirs in the Songliao basin are controlled by space-time matching relationship between source rocks, volcanic traps, and fault transport channels, which provides a basis for finding desirable exploration areas for volcanic rocks. Stages and differences characterize the volcanic eruption, lithology, and structure are complex and diverse, and physical properties are different. Each facies belt has different lithological combinations, forming processes, diagenetic modes and structural structures, and has different pore structure characteristics. Volcanic eruption cycle study can reveal the spatial distribution of volcanic rocks and the genetic relationship among diverse lithologic assemblages, which provides useful guidance for volcanic exploration (Tian et al., 2017).
Chen Jianwen et al. found that oil and gas accumulation is not only related to the lithology of volcanic rocks, but also related to the general characteristics of volcanic rocks, and natural gas is usually enriched in the central part of volcanic structures. By Planke and others, the favorable reservoir is mainly distributed in the favorable lithofacies area of the volcanic crater near fire pass facies belt. In recent years, favorable lithologic and lithofacies zone drilling productivity has been found to be closely related to volcanic type. Tang Jianren et al. have discovered that volcanic reservoirs are controlled by volcanic structures and lithofacies. Different types of volcanic structures have different lithofacies distribution characteristics and hydrocarbon bearing properties. Therefore, studying the eruption cycle of volcanic eruption and defining the distribution characteristics and evolution rules of volcanic apparatus have important guiding and practical significance for volcanic oil and gas exploration.
Volcanic rocks are formed by magma rising from the weakest zone of the strata along the cut through basin basement and the deep fault of the crust. Compared with sedimentary rocks, volcanic rock deposits are not controlled by the accommodation of water bodies. The formation characteristics of this "water and fire" symbiosis make volcanic magma form superposed volcanic bodies in multi-stage and intermittent eruptions (Skelin et al., 2018). The change of volcanic rocks in material composition and lithology is the basis for dividing cycles and periods. In different volcanic eruption cycles or different periods of the same cycle, the material composition is different, resulting in the regularity of volcanic lithology in longitudinal direction, and the corresponding logging curve has the characteristics of prosodic variation (Yang et al., 2018). The seismic reflection characteristics of different lithological characteristics are different in continuity, stratification and amplitude intensity.
Songliao basin is a deep structural depression with large scale and controlled by tensional and torsional faults. Yingshan depression and Xujiaweizi fault depression are the main areas of deep natural gas exploration. The volcanic reservoirs are of strong heterogeneity and lateral variation. The separate volcanic reservoirs can form an independent reservoir forming system, which may appear in a plane with overlapping distribution. Favorable reservoirs are controlled by facies belts (Wendt et al., 2017), the exploration of volcanic gas reservoirs is rather difficult. In order to identify the macroscopic rules of the distribution of volcanic reservoirs in different stages of the study area, in order to guide the exploration of volcanic rocks, this paper divides the eruption characteristics of volcanic rocks in Songliao basin into cycles (Geurden et al., 2017). Due to the less drilling and core data in the study area, the virtual reality technology was used to collect the three-dimensional scene data of Songliao Basin, and the characteristics of the distribution and evolution of each sub volcano were obtained. Based on the study of the step and division of volcanic cycles, the eruption cycle of Songliao Basin was studied in detail, and the subsequent reservoir prediction and meticulous research were carried out. It has important guiding significance for the following reservoir prediction and fine research work.

Survey of research areas
The study area is located in Songliao Basin. The deep layer of Songliao Basin consists of three first-order tectonic units: eastern fault depression zone, central paleo-uplift zone and western fault depression zone. Take the Xujiaweizi fault depression located in the northern Songliao Basin Eastern fault depression and the main part of Yingshan depression zone in the north of the Songliao Basin as the key research area, the volcanic eruption cycle study is carried out (Chen, 2017).
Xujiaweizi fault depression is close to the central paleo-uplift belt and is a deep and large-scale fault depression. It extends nearly north-south with a length of 90 km. The widest part in the middle is 5 km and its area is 4300 km 2 . The fault depression is a single fault graben fault depression in the west to the east. It can be divided into many structural units, such as the western fault depression belt, the sag depression structural belt, the eastern slope zone, the funyue low uplift and the low uplift. According to the results of the third resource evaluation, the total deep resources of Xujiaweizi fault depression are 6772×108 m 3 . As of the end of 2007, 10 5 wells were drilled in the deep part of the Xujiaweizi fault depression. A total of 5058.8 km 2 of pre-stack time migration processing has been completed in the deep layer, with a total proved natural gas reserve of 2457×108 m 3 handed in, with a proven rate of 36.2% (Bracquart et al., 2017).
Yingshan depression covers an area of 400 km 2 , and its syncline and Xujiaweizi fault depression are divided by Chaoyang Valley anticline belt, which is an asymmetric graben type double fault depression under the control of four station fault and Linjiang fault. Before the deposition of the Doutlouku formation, the structure is controlled by volcanic action, forming two ancient landforms with four concave uplifts, namely, the Paleo uplift, the Qingshan uplift and the Yingshan depression. The north is a low-lying area. Before the deposition of the group, the North South differential subsidence further strengthened the southern structure, and the amplitude of the western Paleo long structure decreased. Up to now, strong tectonic inversion has occurred in the north, slight inversion in the south, and stable subsidence in the middle, resulting in the formation of today's surrounding high, middle low tectonic form. Yingcheng formation is located in the fault depression period of the basin. At present, there are five exploration wells in Yingshan depression. There are few logging and drilling data, and the overall exploration and development degree is weak (Pavey et al., 2017).

Research method for eruption cycle of volcanic rocks
Virtual reality technology is used to reconstruct the volcanic terrain in the study area, and the monocular vision method is applied to collect the data collected in the study area. The image machine is used as the image receiving collector to carry out the data measurement and acquisition in the research area needed for comprehensive three-dimensional modeling, such as image processing, data acquisition and visual capture. According to the three dimensional reconstruction data of dynamic scenes in the study area, combined with the research steps of volcanic eruption cycle and volcanic eruption cycle stages, the lithology, transport characteristics and lithofacies periodic changes of volcanic rocks in the study area were analyzed, and the results of the eruption cycle of volcanic rocks in the study area were obtained.

Three-dimensional reconstruction of dynamic scene in research area based on virtual reality technology
By using monocular vision method, the collected image can be either a single image of a single viewpoint or multiple images of several different viewpoints (Lozano-Bilbao et al., 2017). This method can recognize the degree of light and shade, texture lines, distance points, outlines and so on. This method is simple in equipment, simple in structure, and can reconstruct the three-dimensional model of the object using only a single or a few images. These matching points are used to describe and obtain the coordinate information points of the three-dimensional object in space, so as to realize the three-dimensional reconstruction of the research area (Sun, 2017).
(1) Design of photometric stereo vision virtual reality Although conventional brightness methods can reconstruct a threedimensional model from a single image, there is less data information available from a single image. Therefore, the SFS method was improved in Woodham, and the photometric stereo vision method (Englberger & Dörnbrack, 2018) was introduced. The method is based on the brightness equation, as shown in formula (1): dimensional coordinates of the model. After the first display, the collected data will be allocated to a certain area. After the display, the collected data can be calibrated by photometric stereo YYUE. The process is shown in formula (5): x S a a S a L y L b bS b L In the formula, a 0 , a 1 and a 2 are the three-dimensional coordinate valid values of points in the photometric stereo YYUE calibration process, and S is the optimum photometric stereo YYUE calibration parameter value. The point coordinates ML(x, y) and F(x, y) can be obtained intuitively, which greatly improves the occurrence of visual disturbance and ensures the threedimensional dynamic rendering of the dynamic scene in the study area (Rocholl et al., 2017). The effect of three-dimensional reconstruction of the scene in the study area is shown in Figure 2.
In the formula: 1O i ( ) denotes the intensity coefficient of the light source reflected from the surface of the object; g j is the normal vector of the kinetic energy of the object; g i is the refraction vector of the light source, value, which is the value of the response acquisition data to ensure the integrity of the data in the process of data analysis; j N i  ( ) denotes the value range of the data acquisition point. Photometric stereo vision uses multiple sources of light to obtain multiple images of the same object. By combining the conjugate equations of different images, the normal vector g j of the kinetic energy of the object and the refraction vector g i of the light source can be solved. Finally, the three-dimensional shape of the object can be restored. So, the current methods basically use multiple (4-6) light sources for 3D reconstruction, using matrix representation as shown in formula 2.
In the formula, GH matrix represents the mode of data storage calculation and  1 2 represents the number of convertible and simplified conditions of light source. Pixels simplified by transformation can ensure the clarity of the picture in the process of data conversion, and it is more suitable for multi-object acquisition data conversion under different conditions (Xiang et al., 2018). The transformed data need to be distinguished before they can be used, as shown in formula (3): In the formula, exp − β i ( ) represents the high-order attribute distortion vectors after data transformation, which can reflect the validity of the transformation process; X 0 represents the variable parameters that data can effectively retain in the transformation process; X 1 , X 2 , X 3 and X 4 represent the conversion success rates. Through the above process, the design data of photometric stereo vision virtual reality is obtained, and the image of information particle operator distribution is shown in Figure 1.
(2) Three-dimensional reconstruction The principle of 3D reconstruction is mainly to get images from the RGB camera and the far-infrared camera in the same direction, so that the direction of the two cameras can be as parallel as possible. This will help reduce the amount of computation when the camera is corrected, and get the RGB and depth images (Zhang et al., 2019) respectively. In order to solve the problem of image center deviation caused by different positions of two cameras, RGB images are first transformed into coordinates to align with the depth image coordinates. Then the coordinate positions of space points x and y are calculated. Finally, the three-dimensional coordinate point cloud data (x, y, z, R, G, B)is displayed through LOD hierarchical model, as shown in formula (4):

Study on the step of volcanic cycle
According to the 3D reconstruction data obtained from virtual reality technology in the study area, the same volcanic cycle in the volcanic strata in the study area is represented by a set of transitional lithological assemblages, and there are obvious changes in the occurrence of different volcanic cycles. Therefore, the interface of the cycle is manifested as an integrated surface from an unconformity to a sedimentary area in the volcanic area. On the top of the interface, volcanic rocks are obviously overlying. The thickness of a volcanic cycle is generally more than 80m, which can be used for lateral correlation of volcanic cycles with seismic profiles. In the cycle, we can use single well phase contrast  according to the change of lithofacies. The steps of the study of volcanic cycles and phases are: (1) comprehensively using coring and thin slices to identify various data, determine the rock types and the lithofacies types of coring segments in the area, and identify all kinds of geological interfaces (weathering crust, volcanic ash layer and sedimentary interlayer); (2) Make full use of oil field drilling, logging and logging data to establish single well volcanic lithologic columnar section, according to rock assemblage.
(3) Based on the classical volcanic activity facies sequence model, the primary division of volcanic cycles and stages of single well is carried out flexibly; (4) Based on the results of the zircon dating and the division and correlation of volcanic cycles and periods, we should set up separate wells and intervals and push them out according to the principle of "proximity" in order to clarify the eruption cycle and eruption period of the whole area.

Classification basis
The division of volcanic eruption periods should be based on the principle of isochronous comparison and hierarchical control, that is, in the same eruption cycle, the lithology, lithofacies, electricity, eruption sedimentary cycle and seismic reflection characteristics can be divided into two parts (Bolgar et al., 2018). The recognition of volcanic eruption stage should pay attention to unconformity, weathering crust, volcanic rock composition, color, structure, sedimentary rock interlayer, and so on, and combine seismology, geology and logging. In addition, the periodic changes of volcanic facies and the cyclical changes of volcanic components, such as volcanic channel facies, explosive phases, intrusive facies, overflow facies, volcanic sedimentary facies, periodic variations of basic, neutral and acidic volcanic rocks, and alkaline and subalkaline periodic changes are also the basis for dividing volcanic eruption stages .

Results
The eruption period of volcanic rocks in the Xujiaweizi fault depression.

Lithology and prosody characteristics of volcanic rocks
In different stages of volcanic activity, volcanic matter composition, different eruption modes, energy differences, and later crystallization and differentiation make each stage of volcanic rock have different lithological types, showing a certain prosody feature in lithological combination and logging response characteristics. The periodic variation of volcanic lithology and the characteristics of rhythmic combination are the reflection of the intensity of volcanic eruption, so it can be used to divide volcanic eruption times. Photographs of volcanic cores at different eruption stages are shown in Figure 3. The lithology and logging response characteristics of different types of volcanic rocks are shown in Figure 4.
The same period of secondary volcanic rocks in Xujiaweizi fault depression usually has the following lithologic variations: (1) Mainly volcaniclastic. Usually the central eruption or central fissure eruption, its eruption energy is strong, the eruption product is mainly volcaniclastic, and has positive rhythmic lithological characteristics. In general, volcanic eruption material is composed of caking, volcanic breccia and tuff at the top. Figure 3 (a) is the third stage volcanic rock of the Sheng Shen 2 well with this lithological feature. The electrical characteristics of this kind of combination, such as Xushen 8 well, are shown in Figure 4 (a). From the electrical characteristics of Xushen 8 well, the block rock has high resistivity, medium and low natural gamma, while the top tuff has low resistivity and high natural gamma. The top and bottom of this type of volcanic lithologic assemblage are usually volcanic sedimentary rocks or tuff, and the electrical characteristics are the decrease of resistivity and the increase of natural gamma, indicating the beginning and end of a volcanic eruption.
(2) Lava is dominant. The energy of fissure eruptive volcanic rocks is weak, and the eruptive products are mainly volcanic lavas, which generally contain only a small amount of pyroclastic rocks. This lithologic assemblage usually occurs in the Shengping-Anda area and the southern Fung Lok area in the northern part of the Xujiaweizi fault depression. The volcanic lavas are mainly dominated by acid rhyolite, while the northern part of the andesite is dominated by basic basalt and neutral andesite (Zong et al., 2018). Figure 3 (b) core image of the Shenshengeng 203 well shows that the top of the lava is usually very porous, reflecting the end of the eruption of the volcanic eruption, and the lava exposing the surface; at the end of the eruption, due to weathering and surface water communication, the pores on the top of lava are very developed. These large numbers of primary stomata are signs of the end of the eruption period. Figure 4 (b) shows the lithology and electrical characteristics of the basic lava in the northern part of the study area. It can be seen that the third stage volcanic rocks in Dashen 6 well are mainly basalts with a small amount of andesite. It has the characteristics of high resistivity and low natural gamma. Figure 4 (c) shows the lithological and electrical characteristics of acid lava in the central and southern part of the study area. It can be seen that the second stage of Shengshen 4 well is mainly fissure eruption, and the rock type is mainly rhyolite. If the lithologic assemblage of these volcanic rocks is thicker, the middle part tends to show higher natural gamma and high resistivity characteristics. The top and bottom are usually volcanic sedimentary rocks or tuff, and the electrical characteristics are the decrease of resistivity and the increase of natural gamma, indicating the beginning and end of the eruption. (3) The interbedding of lava and pyroclastic. The energy of volcanic eruption in the center fissure type is slightly weaker than that in the center. Volcanic eruption products are usually dominated by pyroclastic and lava, indicating that the energy of volcanic eruption from strong to weak or weak to strong. As shown in Figure 4 (d), the deep well of WangShen 903, with low energy at the beginning of the eruption, is dominated by rhyolite, and then gradually increases, the volcanic ejecta is breccia, and then the energy gradually weakens. Overall, the process of weakening one strong and one weak indicates the ending of a volcanic activity.

Variation based on the volcanic lithofacies belt cycle
In the same period, the cyclical changes of volcanic rock composition, structure and structure make volcanic rocks exhibit certain prosody features. These prosody features overlap or appear repeatedly, reflecting the pattern and output pattern of spatial distribution of volcanic eruption in the same period, and the appearance characteristics of these products, which is the "facies" of volcanic rocks (Cheng et al., 2017). Each stage of volcanic rock is usually composed of one or more phase sequences, showing a continuous or quasi continuous phase sequence.
A complete volcanic facies sequence assemblage in the longitudinal direction is volcanic channel facies -explosive phase overflow facies -intrusion facies -volcanic sedimentary facies. In fact, only 2-3 lithofacies types or single lithofacies types can be found in the facies sequence assemblages, and the phenomenon of duplication or inversion may also occur. The change of phase sequence is actually the manifestation and result of the process of volcanic activity. Therefore, the periodic variation of phase sequence can be used as a basis for dividing volcanic eruption periods.
Volcanic facies in Xujiaweizi fault depression are usually divided into 5 subfacies, 15 facies subfacies, including volcanic channel facies, outburst facies, overflow facies, intrusive facies and volcanic deposits. Each facies zone is subdivided into 3 sub facies. If the eruption energy of the volcano is strong in the early stage and weakened in the later stage, and the eruption time is long, there may be a complete phase sequence combination near the crater. If the energy of the volcanic eruption is weak in the initial stage and strong in the late stage, the combination of the phase sequence will be reversed (Peng et al., 2016). The phase sequence changes of volcanic rocks in Yingcheng formation of Xujiaweizi fault depression mainly begin with explosive facies or overflow facies, and seldom start by volcanic channel facies. Moreover, the volcanic sequence assemblages are usually composed of 2-3 facies belts, and the complete facies assemblages are rare. The characteristics of different eruptive subvolcanic facies in Xujiaweizi fault depression are shown in Table 1.

Division result
From the analysis of lithologic combination and electrical characteristics of these volcanic rocks, we can see that at the beginning and end of the first eruption, the lithology changes are volcanic sedimentary rocks (or tuff) -pyroclastic rocks (or lavas) -volcanic sedimentary rocks (or tuff). The characteristics of electricity are the decrease of density and resistivity, and the increase of natural gamma (Wang et al., 2020). These characteristics are the signs of volcanic intermission. On the basis of lithology and facies combination, according to the regularity of volcanic activity and the particularity of volcanic strata, and combining the characteristics of volcanic lithology, lithologic assemblage, volcanic facies belt, sedimentary rock intercalation and electrical characteristics of the core and volcanic eruption, the first and three sections of Yingcheng formation in Xujiaweizi fault depression are divided into 3 eruption periods, followed by Y1C1, Y1C2 and Y. 1C3, Y3C1, Y3C2, Y3C3. The results are shown in Table 2.

Eruption stages and characteristics of volcanic rocks in Yingshan depression
The volcanic eruption cycle of Yingshan depression in the study area experienced a quiet period of intense volcanic eruption and quiet period of volcanic activity, while sedimentary intercalation is the representative lithology of the quiescence period and the most obvious indicator of the interface of the cycle. The eruption characteristics of volcanic rocks in Yingshan depression are shown in Table 3.
Cycle 1 is divided into stages 1 and 2 from bottom to top. The period 1 is mainly neutral or moderately acidic volcanic rocks. It shows layered like reflection in seismic sections. Sedimentary facies and sedimentary tuff are the main sedimentary facies in stage 2, and the thickness of sedimentary interbeds is about 150m. It can be used as volcanic eruption times alone, indicating volcanic eruption into a long quiet period. The activity of Fire Mountain is from strong to weak, and it shows good stratification and continuity in seismic profile, and is easy to track and compare horizontally.
At the top of the cycle 2, there is a thin sedimentary rock with a thickness of about 50m, showing a continuous strong reflection coaxial axis on the seismic profile, which is the interface between cycle II and cycle III, indicating that the two cycles experienced only a relatively short eruption period. On the whole, the cycle 2 is mainly volcanic rock, and from bottom to top, it is divided into stage 3 and stage 4. Three of them were mainly volcanic lava or volcanic breccia, locally developed fluffy tuff and granularity. Phase 4 was mainly composed of fused tuff and rhyolite tuff, with smaller granularity. It showed that the eruption of multi crater gradually evolved into the eruption of large scale volcanic eruptions. The obvious alteration of granularity and lithology of these volcanic rocks indicates the vertical migration eruption of volcanic rocks at GR. The curves show rhythmic changes, and the reflection of the upper part relative to the lower part of the layer is disorderly on the seismic section.
Cycle 3 is divided into stages 5 and 6 from bottom to top. The stage5 is dominated by rhyolite and rhyolite breccia lava, and the granularity is large. Stage 6 is mainly tuff, with smaller granularity, indicating vertical migration and eruption of volcanic rocks. Volcanic activity gradually ends, and the upper part of the seismic section shows relatively good stratification, while the lower part reflects chaos and vermicular reflection, and the stratification is poor. On the whole, the first stage is controlled by the four station fault, the sub top interfaces are overlapped to the four station fault, and the area is above the first bottom interface.

Discussions
Volcanic eruption period is a relatively continuous volcanic activity. Its material composition and eruption intensity also have regular changes in the longitudinal direction. A volcanic eruption cycle consists of multiple eruption periods. The lithology and seismic profile characteristics of the same cycle differ from each other, such as volcanic lithology change, volcanic particle size change and seismic reflection characteristics change, etc., which is the basis of the internal division period of the cycle (Bergroth et al., 2018). In this paper, the virtual reality technology is used to study the eruption cycle of volcanic eruption in the Xujiaweizi fault depression and Yingshan depression. The data of the study area are collected by monocular vision based on the 3D reconstruction of the virtual reality technology, which ensures data blocking, avoids the transition of the connected data, and ensures the maximum degree of clarity. The dynamic state of data such as texture and transparency can achieve the best imaging state and realize the dynamic rendering of 3D scene in the study area. On this basis, the eruption cycle times of volcanic rocks in the study area are divided according to the classification steps and division basis of the cycle stages.
During the Mesozoic, the study area experienced the regional tectonic evolution of "North-South strike, east-west subduction, mantle crust, and opening and closing alternation", and 3 major volcanic eruptions occurred in succession, namely, Jurassic flint rock formation and Cretaceous Yingcheng formation. Each volcanic eruption cycle has erupted intermittently. Compared with the eruption cycle, the eruption period is much shorter than that of the first cycle. A volcanic eruption cycle is usually composed of multiple eruption periods (Zhou & Zhang, 2017). The same sub volcanic rocks in Xujiaweizi fault depression usually have three lithologic variation assemblages, mainly volcaniclastic, mainly lava and interbedding of lava and pyroclastic rocks. They reflect the rhythmic characteristics of lithological assemblage and logging response, reflecting the strong and weak eruption of volcanoes. In the same period, volcanic eruption materials are spatially distributed volcanic rocks. Finally, the first and three sections of Yingcheng formation in the Xujiaweizi fault depression are divided into 3 eruption periods, followed by Y1C1, Y1C2, Y1C3, Y3C1, Y3C2 and Y3C3. The Yingshan section of Yingshan depression The reflection is chaotic and wormlike, with poor stratification and large formation thickness.
The formation thickness is thin and the stratification is relatively good.

Stratigraphic Contact Relation
The stratum overlaps the Sizhan fault and partially exceeds T41.
The stratum overlaps the Sizhan fault and partially exceeds T41. showing only a continuous strong reflection axis on the seismic profile, but it is obviously different from the chaotic reflection characteristics of the volcanic rocks. The lateral correlation of volcanic eruption cycles can be achieved. The eruption of volcanic rocks has obvious stages and differences. The intensity of each eruption is different, and the overall performance is gradually enhanced.

Conclusions
Volcanic activity is characterized by polycyclic and multi-stage activities. The magmatic products of a volcanic activity cycle and their volcanic structural traces have certain unity and coordination in time and space distribution. Volcanic rocks of different stages of the same cycle, or volcanic rocks of different stages of the same period, are different in volcanic lithofacies zones due to different energy and eruption modes. Different volcanic lithofacies zones control the development of reservoirs and the enrichment and high production of oil and gas. Therefore, by using virtual reality technology to build three-dimensional dynamic scenes in the study area, it is of great significance to study volcanic eruptions cycles times and effectively divide volcanic eruption stages to identify favorable facies zones of volcanic rocks, analyze the inducing factors of reservoir development, and hydrocarbon distribution and accumulation rules.