CN107526104A - Fracture medium seismic wave field method for numerical simulation based on multimachine multinuclear - Google Patents
Fracture medium seismic wave field method for numerical simulation based on multimachine multinuclear Download PDFInfo
- Publication number
- CN107526104A CN107526104A CN201710749074.9A CN201710749074A CN107526104A CN 107526104 A CN107526104 A CN 107526104A CN 201710749074 A CN201710749074 A CN 201710749074A CN 107526104 A CN107526104 A CN 107526104A
- Authority
- CN
- China
- Prior art keywords
- numerical simulation
- mpi
- seismic wave
- wave field
- parallel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 68
- 238000004088 simulation Methods 0.000 title claims abstract description 32
- 230000008569 process Effects 0.000 claims description 33
- 230000010363 phase shift Effects 0.000 claims description 13
- 238000012937 correction Methods 0.000 claims description 8
- 238000009826 distribution Methods 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 6
- 238000003384 imaging method Methods 0.000 claims description 4
- 238000003825 pressing Methods 0.000 claims 1
- 230000008901 benefit Effects 0.000 abstract description 4
- 238000004422 calculation algorithm Methods 0.000 abstract description 3
- 208000010392 Bone Fractures Diseases 0.000 description 18
- 206010017076 Fracture Diseases 0.000 description 18
- 230000004044 response Effects 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention discloses a kind of fracture medium seismic wave field method for numerical simulation based on multimachine multinuclear, applied to Fractured seismic wave field field, handled by introducing cluster in Wave Equation Numerical, worked by MPI message passing interfaces on multimachine, so that computer cluster is applied in the fracture medium seismic wave field numerical simulation that fine-structure mesh is formatted, it is very big to improve arithmetic speed；And with reference to the application of multithreading, more nulcear properties of individual node are taken full advantage of, take full advantage of the internal memory and slack resources of computer, and a very big step improves the speed-up ratio of algorithm.
Description
Technical field
The invention belongs to Fractured seismic wave field field, more particularly to a kind of fracture-type reservoir geological model.
Background technology
Fractured hydrocarbon reservoir has become Global Oil, the important component of gas reserves and increasing the storage major fields.Profit
It is considered as always by geoscientist and reservoir engineer in the research for having much challenge with ground forecast for seismic data fractured hydrocarbon reservoir
Hold, one of key issue being directed to is that different shape, density, the seismic wave field response characteristic in size yardstick crack are difficult to
Identification.In the environment of current international low oil price, seismic wave field response characteristic caused by crack is clearly recognized to improving crack
Type prediction of oil-gas reserve precision and fractured hydrocarbon reservoir drilling success are most important.Seismic wave field numerical simulation is that people are correct
Propagation law and underground medium construction, the effective means of structure and lithologic character of seismic wave are solved, is had correctly theoretical and real
Border meaning, it is premise and the basis for identifying and describing oil and gas reservoir using method of seismic prospecting.
Conventional seismic wave field method for numerical simulation has ray tracing and Wave Equation Method.Wave Equation Method can be comprehensively
Reflect feature of both the kinematics and dynamics of seismic wave, and use rational algorithm, can obtain high-precisionly
Wavefield forward modeling record is shaken, in particular for fracture medium.It is well known that the crack in stratum has multiple dimensioned property, size
Between several microns to tens meters.The crack of large scale can only be conventionally simulated, it is incompetent for the crack of small yardstick
For power.The numerical simulation of usual fracture medium is all based on EFFECTIVE MEDIUM THEORY and completed, that is, needs according to fracture medium
Property (density, charges etc.) is divided into some Fractured Zones, and assumes all in same Fractured Zone
Containment body medium between crack and crack has identical property.Fracture medium earthquake is simplified although with EFFECTIVE MEDIUM THEORY
The complexity of wave field numerical, but the seismic wave field that numerical simulation result can only be obtained at fracture development geological objects boundary rings
Answer feature, it is impossible to effectively analyze response characteristic of the different shapes and sizes crack to seismic wave field, it is impossible to after effectively instructing
The detailed predicting of continuous fracture reservoir and description.
The content of the invention
In order to solve the above technical problems, present applicant proposes a kind of fracture medium seismic wave field numerical value based on multimachine multinuclear
Analogy method, handled by introducing cluster in Wave Equation Numerical, by MPI message passing interfaces on multimachine
Work so that computer cluster is applied in the fracture medium seismic wave field numerical simulation that fine-structure mesh is formatted, very big to improve
Arithmetic speed.
The technical scheme that the application uses for：Fracture medium seismic wave field method for numerical simulation based on multimachine multinuclear, bag
Include：
S1, according to the direction in crack, shape, volume size, find out minimum Crack Element, built using Mesh Subdivision Technique
Geological model；
S2, combined using phase shift correction and Phase Shift Interpolation seismic wave field numerical simulation is carried out to the geological model of structure；
S3, in wave equation frequency domain MPI and OpenMP parallel forms are used, realize seismic wave field numerical simulation, had
Body is：The multi-process parallel schema of MPI coarseness is used between the calculate node of cluster, is used in calculate node
The OpenMP fine-grained multi-threaded parallel pattern by frequency subdivision.
Further, the multi-process parallel schema of MPI coarseness is used described in step S3 between the calculate node of cluster
Specially：
A1, MPI_INIT_THREAD is called to start parallel environment；
A2, the quantity according to calculate node, MPI_Comm_size order determination process numbers are called, while call MPI_
Comm_rank orders obtain the numbering of each process；
A3, other processes send the scope for asking to need frequency chunks to be processed to host process, after host process receives request
Frequency chunks scope to be processed is needed to the distribution of other processes；
A4, each process perform Wave equation forward modeling numerical simulation parallel computation；
A5, imaging results file is sent to host process, merged by host process；
A6, MPI_Finalize is called, terminate parallel environment.
Further, the OpenMP fine granularity multithreading by frequency subdivision is used described in step S3 in calculate node simultaneously
Row mode, it is specially：
B1, the CPU core number according to current node, call int_omp_set_num_threads function setup present nodes
The number of threads to be started；
B2, frequency cycle outer layer call #pragma omp parallel num_threads () sentence, by frequency wave
Numerical Simulation is set as parallel computation region；
B3, thread private variable application；
B4, in frequency cycle outer layer by calling omp_get thread_num functions, obtain current thread number, be every
Individual thread distribution starting and termination frequency.
Beneficial effects of the present invention：The fracture medium seismic wave field method for numerical simulation based on multimachine multinuclear of the present invention,
Pass through MPI+OpenMP hybrid parallel pattern so that calculator memory and check figure are all fully used to be accelerated.It is logical
Cross with cluster to realize MPI+OpenMP mixed models so that computational efficiency has obtained significant increase, and can reach following has
Beneficial effect：
1) handled by introducing cluster in Wave Equation Numerical, by MPI message passing interfaces on multimachine
Work so that computer cluster is applied in the fracture medium seismic wave field numerical simulation that fine-structure mesh is formatted, very big to improve
Arithmetic speed；
2) apply multithreading in above process again, take full advantage of more nulcear properties of individual node, take full advantage of
The internal memory and slack resources of computer, but a very big step improves the speed-up ratio of algorithm；
3) from physical significance, this method of the present invention it can be considered that different fractuer directions, shape, volume size etc. over the ground
The influence of seismic wave field, obtained numerical simulation result more conform to actual conditions, are that the seismic wave filed character of fracture reservoir is studied
Provide a kind of idea and method.
Brief description of the drawings
Fig. 1 is the protocol procedures figure of the application；
Fig. 2 is the EFFECTIVE MEDIUM THEORY simulation and the effect contrast figure of the application method that the embodiment of the present application provides；
Wherein, Fig. 2 (a) is EFFECTIVE MEDIUM earthquake record；Fig. 2 (b) is the earthquake record of the application method；
Fig. 3 is the migrated section that the embodiment of the present application provides；
Fig. 4 is each platform run time contrast that the embodiment of the present application provides.
Embodiment
For ease of skilled artisan understands that the technology contents of the present invention, enter one to present invention below in conjunction with the accompanying drawings
Step explaination.
It is the protocol procedures figure of the application as shown in Figure 1, the technical scheme of the application is：Crack based on multimachine multinuclear is situated between
Quality seismic wave Numerical Simulation method, including：
S1, according to the direction in crack, shape, volume size, find out minimum Crack Element, built using Mesh Subdivision Technique
Geological model；
S2, the petrophysical property for considering fracture developing zone internal fissure and surrounding rocks, it is horizontal according to speed is suitable for
To big Wave Equation Numerical theory is changed, propose to combine the Geological Model to structure using phase shift correction and Phase Shift Interpolation
Type carries out seismic wave field numerical simulation；
S3, in the implementation process of program, it is contemplated that to carrying out wave equation number using the geological model after grid subdivision
Value simulation amount of calculation is huge, in wave equation frequency domain using MPI and OpenMP parallel forms, realizes seismic wave field numerical value
Simulation, it is specially：The multi-process parallel schema of MPI coarseness is used between the calculate node of cluster, is adopted in calculate node
With the OpenMP fine granularity multi-threaded parallel pattern segmented by frequency.
Phase shift correction and the continuation operator of Phase Shift Interpolation combination are specially in step S2：
From phase shift correction method, the wave equation of three-dimensional zero shot-geophone distance geological data is as follows：
In the case of medium is heterogeneous, slowness field is decomposed into two parts：
In (2) formula, u0(z) it is that and the change of slowness is all included into Δ u (x, y, z) component with reference to slowness.(2) formula
(1) formula of substitution simultaneously makees Fourier transform to time t, obtains：
Or
In formula (3)：4w2(2u0Δu+Δu2) p=s (x, y, the z) wave equation of non-uniform dielectric (be) is referred to as source letter
Number.
At this moment Δ u is ignored2Effect solve equation (4), obtain another e that mutually transplantsiwΔu(x,y,z)Δz, this mutually transplants just
It is the approximate disturbance source item for reflecting horizontal slowness.The wavelength Extrapolation method of phase shift correction up can simply be summarised as (5) formula and retouch
State：
p(x,y,zn+1, w) and=eiwΔuΔz2DDFT-1(2DFFT(P(x,y,znw))eiθΔz) (5)
In (5) formula：
For the laterally violent change of speed-adaptive wave field, the method that interpolation is used to (5) formula, it is located in Δ z layer most
Small speed is v1, maximal rate v2, then have slowness：
U is used in Δ z intervals1And u2Continuation is carried out, u (x, y, z) is the actual slowness in Δ z intervals.
Wherein,For the phase shift correction item in spatial domain.
Actual wave field is obtained by following interpolation formula：
Finally Seismic imaging result is：
Although can realize the seismic wave field numerical simulation of fracture medium according to above formula, this method needs a large amount of
Computing resource could realize.Substantial amounts of computing resource is consumed to be mainly reflected in：(1) seismic wave number of fields is carried out to each frequency
, it is necessary to call Fast Fourier Transform (FFT) repeatedly in time-domain and space when value is simulated；(2) it is horizontal for crack and country rock speed
Quick change obtain the higher wave field of precision, it is necessary to enter row interpolation using multiple speed；(3) fracture medium is directed to, using net
The geological model of the method structure of lattice subdivision also significantly increases amount of calculation.For this, the application employs multimachine multinuclear
Parallel framework is accelerated.
Process level is counted parallel between the node for the Wave Equation Numerical that phase shift correction and Phase Shift Interpolation based on MPI combine
Calculation realizes that step is as follows：
A1, MPI_INIT_THREAD is called to start parallel environment；
A2, the quantity according to calculate node, calling MPI_Comm_size order determination processes number, (node performs
One process), while call the numbering of each process of MPI_Comm_rank orders acquisition；
A3, other processes send the scope for asking to need frequency chunks to be processed to host process, and No. 0 process is host process, is connect
Receiving other backward distributions of request needs frequency chunks scope to be processed；
A4, each process perform Wave equation forward modeling numerical simulation parallel computation；Wave equation expression formula is shown in formula (1)；
A5, imaging results file is sent to No. 0 process, No. 0 process merges；
A6, MPI_Finalize is called, terminate parallel environment.
The Wave Equation Numerical parallel computation realized according to above-mentioned steps, only start one on each node and enter
A CPU core is only called to participate in calculating in journey, namely each node, therefore it can provide sufficiently large internal memory for each process.
Certainly, each node only has a CPU core to participate in computing under the parallel schema, also result in the waste of computing resource, is abundant
Using the computing resource of each node, the present invention is being realized based on Wave Equation Numerical parallel computation between MPI node
On the basis of, in each intra-node again according to the CPU core number of each node, the one process of each node is divided based on OpenMP
Solve and perform a thread for multiple threads, each CPU core, to make full use of computing resource, improve the parallel of Wave equation forward modeling
Computational efficiency.
Wave Equation Numerical Thread-Level Parallelism calculates and realizes that step is as follows in node based on OpenMP：
B1, the CPU core number according to current node, call int_omp_set_num_threads function setup present nodes
The number of threads to be started；
B2, frequency cycle outer layer call #pragma omp parallel num_threads () sentence, by frequency wave
Numerical Simulation is set as parallel computation region；
B3, thread private variable application, such as：#pragma omp parallel for firstprivate(v,
real,imag,vv,re,im) lastprivate(re,im)；In order to ensure thread-safe, because the variable in each thread
What name was just as, when variable changes in thread, easily clash, cause thread dangerous.
B4, frequency cycle outer layer call omp_get thread_num functions, obtain current thread number, be each line
Journey distribution starting and termination frequency.
Shown in the earthquake record such as Fig. 2 (a) simulated according to EFFECTIVE MEDIUM THEORY.Only it is observed that crack is sent out in Fig. 2 (a)
The various diffracted waves with boundary are educated, without response inside Fractured Zone；Fig. 2 (b) is the calculating obtained using this patent method
Earthquake record；Various diffracted waves can be clearly identified in Fig. 2 (b), not only it is observed that the seismic response of slit band boundary
Feature, the reflectance signature inside slit band can be also shown, therefore these fracture developments can be identified on the seismic profile of reality
Area.During due to actual seismic data migration processing, it is impossible to be distinguish between the speed of crack and containment body, we use equivalent Jie
Matter obtains enter line displacement of the speed to Fig. 2 (b), obtains the migrated section shown in Fig. 3.We can clearly identify from Fig. 3
Seismic response features caused by crack, the correctness and validity of proposition method are also confirmed from another side.
Realized using cluster, this patent body of velocity size 4096*4096, opens 6 fortune for the geological model of structure
Row node, multithreading number are 3, it can be seen that the time comparison diagram of each operation platform is as shown in figure 4, be based on MPI+
OpenMP time substantially improves 23 times than unit monokaryon.
One of ordinary skill in the art will be appreciated that embodiment described here is to aid in reader and understands this hair
Bright principle, it should be understood that protection scope of the present invention is not limited to such especially statement and embodiment.For ability
For the technical staff in domain, the present invention can have various modifications and variations.Within the spirit and principles of the invention, made
Any modification, equivalent substitution and improvements etc., should be included within scope of the presently claimed invention.
Claims (3)
1. the fracture medium seismic wave field method for numerical simulation based on multimachine multinuclear, it is characterised in that including：
S1, according to the direction in crack, shape, volume size, find out minimum Crack Element, geology built using Mesh Subdivision Technique
Model；
S2, combined using phase shift correction and Phase Shift Interpolation seismic wave field numerical simulation is carried out to the geological model of structure；
S3, in wave equation frequency domain MPI and OpenMP parallel forms are used, realize seismic wave field numerical simulation, be specially：
The multi-process parallel schema of MPI coarseness is used between the calculate node of cluster, the pressing using OpenMP in calculate node
The fine-grained multi-threaded parallel pattern of frequency subdivision.
2. the fracture medium seismic wave field method for numerical simulation according to claim 1 based on multimachine multinuclear, its feature exist
In the multi-process parallel schema of the coarseness described in step S3 between the calculate node of cluster using MPI is specially：
A1, MPI_INIT_THREAD is called to start parallel environment；
A2, the quantity according to calculate node, MPI_Comm_size order determination process numbers are called, while call MPI_Comm_
Rank orders obtain the numbering of each process；
A3, other processes send request to host process needs the scope of frequency chunks to be processed, host process receive request it is backward its
He needs frequency chunks scope to be processed at process distribution；
A4, each process perform Wave equation forward modeling numerical simulation parallel computation；
A5, imaging results file is sent to host process, merged by host process；
A6, MPI_Finalize is called, terminate parallel environment.
3. the fracture medium seismic wave field method for numerical simulation according to claim 1 based on multimachine multinuclear, its feature exist
In, the fine granularity multi-threaded parallel pattern by frequency subdivision described in step S3 in calculate node using OpenMP, specially：
B1, the CPU core number according to current node, int_omp_set_num_threads function setups present node is called to open
Dynamic number of threads；
B2, frequency cycle outer layer call #pragma omp parallel num_threads () sentence, by frequency wave number of fields
Value simulation is set as parallel computation region；
B3, thread private variable application；
B4, in frequency cycle outer layer by calling omp_get thread_num functions, obtain current thread number, be each line
Journey distribution starting and termination frequency.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710749074.9A CN107526104A (en) | 2017-08-28 | 2017-08-28 | Fracture medium seismic wave field method for numerical simulation based on multimachine multinuclear |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710749074.9A CN107526104A (en) | 2017-08-28 | 2017-08-28 | Fracture medium seismic wave field method for numerical simulation based on multimachine multinuclear |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107526104A true CN107526104A (en) | 2017-12-29 |
Family
ID=60682718
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710749074.9A Pending CN107526104A (en) | 2017-08-28 | 2017-08-28 | Fracture medium seismic wave field method for numerical simulation based on multimachine multinuclear |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107526104A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110542920A (en) * | 2019-09-03 | 2019-12-06 | 北京云庐科技有限公司 | Seismic data processing method and system |
| CN110927784A (en) * | 2019-11-29 | 2020-03-27 | 中国地质大学（武汉） | Frequency domain Gaussian beam Green function calculation method based on OpenMP |
| CN112632005A (en) * | 2019-10-08 | 2021-04-09 | 中国石油化工股份有限公司 | Seismic data calculation method and system based on MPI |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104461466A (en) * | 2013-09-25 | 2015-03-25 | 广州中国科学院软件应用技术研究所 | Method for increasing computing speed through parallel computing based on MPI and OpenMP hybrid programming model |
| CN102239507B (en) * | 2008-12-03 | 2016-02-24 | 雪佛龙美国公司 | For the mess generation system and method for discrete fracture modeling |
-
2017
- 2017-08-28 CN CN201710749074.9A patent/CN107526104A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102239507B (en) * | 2008-12-03 | 2016-02-24 | 雪佛龙美国公司 | For the mess generation system and method for discrete fracture modeling |
| CN104461466A (en) * | 2013-09-25 | 2015-03-25 | 广州中国科学院软件应用技术研究所 | Method for increasing computing speed through parallel computing based on MPI and OpenMP hybrid programming model |
Non-Patent Citations (3)
| Title |
|---|
| 宋鹏 等: "基于MPI+OpenMP的三维声波方程正演模拟", 《中国海洋大学学报》 * |
| 张选朋 等: "相移及相移内插（PSPI）偏移方法研究", 《城市建设理论研究（电子版）》 * |
| 张邦: "三维复杂介质弹性波场有限差分数值模拟方法研究", 《中国优秀硕士学位论文全文数据库·基础科学辑》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110542920A (en) * | 2019-09-03 | 2019-12-06 | 北京云庐科技有限公司 | Seismic data processing method and system |
| CN110542920B (en) * | 2019-09-03 | 2021-06-22 | 北京云庐科技有限公司 | Seismic data processing method and system |
| CN112632005A (en) * | 2019-10-08 | 2021-04-09 | 中国石油化工股份有限公司 | Seismic data calculation method and system based on MPI |
| CN112632005B (en) * | 2019-10-08 | 2024-01-23 | 中国石油化工股份有限公司 | MPI-based seismic data calculation method and system |
| CN110927784A (en) * | 2019-11-29 | 2020-03-27 | 中国地质大学（武汉） | Frequency domain Gaussian beam Green function calculation method based on OpenMP |
| CN110927784B (en) * | 2019-11-29 | 2020-09-29 | 中国地质大学（武汉） | Frequency domain Gaussian beam Green function calculation method based on OpenMP |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ahmed et al. | Control-volume distributed multi-point flux approximation coupled with a lower-dimensional fracture model | |
| Ahmed et al. | Three-dimensional control-volume distributed multi-point flux approximation coupled with a lower-dimensional surface fracture model | |
| US10995592B2 (en) | Method and system for analyzing the uncertainty of subsurface model | |
| US10318663B2 (en) | Method of reservoir compartment analysis using topological structure in 3D earth model | |
| AU2018277657A1 (en) | Constructing structural models of the subsurface | |
| AU2011356658A1 (en) | Method of reservoir compartment analysis using topological structure in 3D earth model | |
| Wheeler et al. | A parallel multiblock/multidomain approach for reservoir simulation | |
| Zhang et al. | Determination of RVE with consideration of the spatial effect | |
| CN107526104A (en) | Fracture medium seismic wave field method for numerical simulation based on multimachine multinuclear | |
| Xu et al. | Building a three dimensional sealed geological model to use in numerical stress analysis software: a case study for a dam site | |
| Chen et al. | Scalable generation of large-scale unstructured meshes by a novel domain decomposition approach | |
| Wang et al. | Method for generating high-quality tetrahedral meshes of geological models by utilizing CGAL | |
| Wang et al. | Accurate modeling and simulation of seepage in 3D heterogeneous fractured porous media with complex structures | |
| CN114895351A (en) | Method and device for medium modeling for simulating seismic wave propagation at any discontinuous interface | |
| Qian et al. | GAN inversion method of an initial in situ stress field based on the lateral stress coefficient | |
| Mohajerani et al. | An efficient algorithm for generation of conforming mesh for three-dimensional discrete fracture networks | |
| Liu et al. | Flat voxel-based modelling, assessment and visualization of urban underground space resource quality | |
| Zhang et al. | Conditional stochastic simulation of fluvial reservoirs using multi-scale concurrent generative adversarial networks | |
| Gao | Inverse back analysis based on evolutionary neural networks for underground engineering | |
| Tran et al. | Characterizing and modelling of fractured reservoirs with object-oriented global optimization | |
| Hu et al. | An efficient computational framework for height-contained growing and intersecting hydraulic fracturing simulation via SGBEM–FEM | |
| Käser et al. | SeisSol–a software for seismic wave propagation simulations | |
| CN110162804A (en) | The wavefield forward modeling optimization method accelerated based on CPU | |
| Liang et al. | A novel workflow for fracture reconstruction and uncertainty analysis for unconventional reservoir development | |
| Gable et al. | FracMan/dfnWorks: Translating Geological Fracture Characterizations for Massively Parallel Subsurface Flow and Transport Simulations |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171229 |