CN104453877A - Underground deep-buried meander point dam sand body historical reconstruction method - Google Patents

Underground deep-buried meander point dam sand body historical reconstruction method Download PDF

Info

Publication number
CN104453877A
CN104453877A CN201410618134.XA CN201410618134A CN104453877A CN 104453877 A CN104453877 A CN 104453877A CN 201410618134 A CN201410618134 A CN 201410618134A CN 104453877 A CN104453877 A CN 104453877A
Authority
CN
China
Prior art keywords
amassed
river course
point bar
unit
streamline
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.)
Granted
Application number
CN201410618134.XA
Other languages
Chinese (zh)
Other versions
CN104453877B (en
Inventor
单敬福
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Furuibao Energy Technology Co ltd
Original Assignee
Yangtze University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yangtze University filed Critical Yangtze University
Priority to CN201410618134.XA priority Critical patent/CN104453877B/en
Publication of CN104453877A publication Critical patent/CN104453877A/en
Application granted granted Critical
Publication of CN104453877B publication Critical patent/CN104453877B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00Water conservation; Efficient water supply; Efficient water use
    • Y02A20/40Protecting water resources
    • Y02A20/402River restoration

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Architecture (AREA)
  • Mathematical Analysis (AREA)
  • Structural Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Civil Engineering (AREA)
  • Revetment (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses an underground deep-buried meander point dam sand body historical reconstruction method. A preliminary-stage riverway flow line and terminal-stage riverway flow line enveloping method is adopted, and combined with a sand body thickness diagram, a point dam is recognized; through extraction of the point dam lateral deposition body scale, the tendency, the inclination and lateral deposition body spacing parameter in a research region, a point dam building structure recognition key parameter database is built, through a skeleton section, fine dissection is conducted on the meander riverway sand body of a stratigraphic unit, riverway flow lines are obtained, then the riverway flow lines are sequentially recovered and arranged according to the development sequence, and the riverway sedimentary evolution history is represented. The underground deep-buried meander point dam sand body historical reconstruction method has the advantages that the sand body building structure dissection result for riverway historical process rebuilding is made to be more reasonable and closer to the actual underground situation.

Description

A kind of underground buried point bar of meandering sand body history reconstruction method
Technical field
The invention belongs to oil field basic geology technical field, relate to a kind of underground buried point bar of meandering sand body history reconstruction method.
Background technology
Single river course and single bar river course concept, both are distinguishing, are fully following longitudinally by stages, and on plane branch principle basis, time property of longitudinal phase is emphasized in single river course, and single bar river course is then it is emphasised that plane braiding property.Combination channel inside contains some single river courses, river course in sedimentary microfacies map is also all combination channel rank mostly, for the formation of combination channel, mainly river course usually has migration to swing in the plane, the effect such as discarded and cutting causes, thus define the larger combination channel of distribution, broadening generally can reach even a few km of km, and on log, combination channel often shows as that compound is bell, complex case shape and complication clock box.
Fluvial-delta depositional reservoir configuration is the another hot issue of ground over nearly 20 years, and numerous scholar, based on different visual angles, adopts different research means methods, has carried out fine dissection and build geological model based on this to reservoir.But subsurface reservoir is extremely complicated, most sand body is composed sand body, it is not the single origin cause of formation, and have repeatedly the result of the same period and late-stage superposed deformation, if fluvial facies are in lateral transfer swing process, there are erosion, cutting and discarded phenomenon etc., these make sand body complicate all to a certain extent, if can not start with from Genesis Evolution rule, directly deposition results is described, often easily originally complicated reservoir alligatoring, simplification, thus get the wrong sow by the ear and be familiar with, being unfavorable for the adjustment of oil field later stage measure and the enforcement of development plan.Sand body building structure for the meandering channel reservoir origin cause of formation is dissected, and forefathers have done large quantifier elimination, emphasizes that in dry season and flood period water channel, main stream line position and hydrodynamic force there are differences, and points out what point bar was mainly formed in flood period; By to fluvial depositional reservoir configuration, propose level constraint, the subsurface reservoir conformational analysis of model fitting and multidimensional, interactive and modeling method, mainly system overview and summary have been carried out to the thinking of modeling and research method; By methods such as comprehensive outcrop, made ground and the empirical formulas set up thus, establish Gudong field point bar of meandering endo conformation quantitative model, advantage the mathematical methods such as geostatistics is incorporated into Reservoir Architectural dissect, and attempt carrying out quantification sign to reservoir; Lin Chengyan etc. (2011) are by the fine dissection to the building structure of Song-liao basin Daqing Oil Field fluvial depositional reservoir sand body, and the fork proposing river-delta front distributary channel is that the favourable enrichment region of a kind of potential remaining oil is theoretical; Taking into full account that point bar private side amasss mudstone foundation and affects on reservoir permeability spatial distribution basis, propose point bar of meandering 3-d modelling pattern, wherein emphasize from side deposition, face is amassed in side, the important architectural element in 3, abandoned channel is started with, with reference to river course geometry and river work parameter, identify, recover and calculate (13 kinds of extracting methods) point bar side and amass body (comprising inclination angle, tendency, scale) scale, rule is amassed in summary side deposition rock short lap spacial distribution density and side, river course; Zhang Changmin etc. (2013) by identify abandoned channel, latter stage river course, utilizing the change of river course streamline, analyze river course swaying movement regularity, rising to process evolution level first analyzing sand body building structure.Obtaining of above-mentioned achievement in research, for the enforcement of oil field development measure and adjustment provide important reference, has also enriched reservoir architecture theoretical further.
But, in reservoir characterization process, in achievement in research in the past, all too much focus on being described deposition results, for its evolutionary process of recovery, then seldom mention, just domestic indivedual scholar goes to set forth river evolution transition history by the mode of channel line, this method makes river evolution history obtain reconstruction to a certain extent, has had significant progress relatively in the past.In fact, a proposition inherently fuzzy abstract concept of channel line, its core be by channel line initial position and latter stage position estimation, then interpolation main stream line in the middle of, by rebuilding river channel evolution history to the transition of main stream line, although thinking is correct, but the initial position method of discrimination of main stream line, foundation and how the aspect such as interpolation main stream line then discuss less, in fact, the identification of main stream line is always a difficult point, is also key point.The author is with regard to these problems, attempt by reconstruct channel deposit process study thinking and countermeasure, quantize angle duplicate removal from Semi-qualitative semidefinite and build meander river course history of evolution, make the sand body building structure anatomical results of rebuilding based on river course historical process more rationally, closer to underground truth.
It is main study subject that poplar big city, oil recovery factory of Jilin Oil Field Fuyu County seed oil layer is chosen in this research, Fuyu Oil Field construction location is in east, concavity district, Song-liao basin south edge, No. 3, Fuxin uplift belt Fuyu County is textural, be one by the complicated how high some arched roof anticline of tomography, oil reservoir, mainly by orecontrolling factor, belongs to Fractured hyposmosis structure sandstone oil reservoir.The main order stratum of oilfield exploitation is the FY oil-layer of spring head group four sections, and the more stable lithology of reservoir distribution is mainly siltstone and packsand.Previous karyotype studies confirms, poplar big city, oil recovery factory of Jilin Oil Field Fuyu County seed oil group major developmental meandering stream deposit.Fuyu Oil Field, since scale in 1970 is gone into operation, along with going deep into of oil field development, has entered late high water content period mining phase at present, is faced with severe development situation.Although different times carried out meticulous pool description work for research block, but this Research of Reservoir Description based on substratum can not meet the needs that oil field is deeply developed further, and with increase that is static, dynamic data, and awareness is constantly deepened, development problem constantly increases, and existing achievement constantly need adjust improvement.Because remaining oil distribution in this very complicated in period, in the dispersity of height, in the urgent need to carrying out fine dissection by the new research meanses such as reservoir architecture and method to the building structure of reservoir building, to solve the flooding pattern imperfection, dispensing and the effective problem such as the low and flood effectiveness variation of injection rate that cause because geological knowledge is unclear.
Summary of the invention
The object of the present invention is to provide a kind of underground buried point bar of meandering sand body history reconstruction method, solve the inaccurate problem of sand body building structure anatomical results that prior art is rebuild for river course historical process.
The technical solution adopted in the present invention is carried out according to following steps:
Step 1: adopt initial stage river course streamline and latter stage river course streamline enveloping method, in conjunction with sand thickness figure, point bar is identified;
Step 2: the identification of body scale is amassed in point bar side;
(1) extraction of body scale is amassed in side:
The extraction that body scale is amassed in side comprises length and the width that body space spread is amassed in side, and wherein asking for of length is also asking for of point bar length, point bar length (W l) be the function of the width (W) of movable water channel, and the width of movable water channel (W) is the function of the thickness (H) in the single river course of the vertical announcement of individual well, namely meets:
W=6.8H 1.54
W: the width of movable water channel, unit m;
H: the vertical sand thickness in single river course, unit m;
The length of body is amassed in side:
W L=85Ln(W)+250;
W l: body length is amassed in side, unit m;
W: the width of movable water channel, unit m;
And body width W is amassed in side din the asking for of (body width is amassed in side), W dfollowing functional relation is there is with W:
W d=2/3W;
W d: body width is amassed in side, unit m;
W: the width of movable water channel, unit m;
(2) extraction of body tendency is amassed in side:
According to place, crescent abandoned channel arc length mid point tangent position and circle center line connecting, amass the tendency of body as side;
(3) Ji Ti inclination angle in side is extracted:
Formula is extracted at inclination angle:
θ=arctg(Δh/Δl);
θ: Ji Ti inclination angle, side, unit m;
Δ h: with depth displacement between aboveground lower adjacent interlayers, unit m;
Δ l: antithetical phrase well well spacing, unit m;
(4) extraction of body spacing is amassed in side:
It is the interval that side amasss between mudstone foundation nearly river course tip position that body spacing is amassed in side, and its acquiring method is on extraction basis, lamination inclination angle, side, and increase measures well side and amasss mudstone foundation buried depth, is recorded as h2_2, then:
Δz=h2_2-h1
ΔL=Δl+Δz/tgθ
Δ L: body spacing is amassed in side, unit m;
Δ l: antithetical phrase well well spacing, unit m;
Δ z: antithetical phrase aboveground lower two interlayer two well points place drop, unit m;
θ: Ji Ti inclination angle, side;
H2_2: mudstone foundation buried depth is amassed in well 2 side, unit m;
H1: mudstone foundation buried depth is amassed in well 1 side, unit m;
(5) by amassing to point bar side, study area the extraction that body spacing parameter is amassed in body scale, tendency, inclination angle and side, set up point bar building structure identification key parameter database, by skeleton section, the meandering river sandbody of formation unit carries out fine dissection.
Step 3: meander river course evolutionary process and history are rebuild;
Amass that body identifies by step 2 pair point bar side, fine dissection and combination, obtain river course streamline, then restore successively by growth order and arrange, reproduce channel deposit history of evolution.
Further, in described step 1, judging point dam location method is: the first phase that water channel is grown, its streamline connects sand body depocenter with shortest path, and Final Issue water channel streamline path is the longest, the scope of first phase river course streamline and Final Issue river course streamline institute envelope is exactly point bar present position.
Further, in described step 3, the determination of the first later and decadent stage of a school of thought line position in river course, the envelope method based on the thick figure of sand can be adopted, can effectively identify river course initial position and latter stage river location, middle river course streamline amasss with the side same period track that body is combined to form and replaces, and so just can carry out effective restoration and reconstruction to river course sedimentary evolution history of meandering.
The invention has the beneficial effects as follows river course historical process is rebuild sand body building structure anatomical results more rationally, closer to underground truth.
Accompanying drawing explanation
Fig. 1 is just later and decadent stage of a school of thought line envelope method delineation point bar schematic diagram;
Fig. 2 be side amass body tendency ask for schematic diagram;
Fig. 3 is Ji Ti inclination angle, point bar side acquiring method schematic diagram;
Fig. 4 be point bar side amass the body side same period amass body combination with point bar sand body sedimentary evolution and history reconstruction figure;
Fig. 5 is just later and decadent stage of a school of thought line envelope method delineation point bar operational instances figure;
Fig. 6 be point bar side amass body tendency example ask for figure;
Fig. 7 is that point bar side is amassed body spacing and asked for instance graph;
Fig. 8 is that point bar fine structure dissects instance graph;
Fig. 9 is that meander river course sedimentary evolution history recovers instance graph.
Detailed description of the invention
Below in conjunction with the drawings and specific embodiments, the present invention is described in detail.
1. point bar identification
Point bar is sand body framework in meandering river sandbody, corresponding with Miall level Four interface phase with configuration interface.Point bar be formed at unidirectional cross current in river course, asymmetrical cross current causes the imbalance of sediment transport naturally, concave bank is made to suffer sapping, convex bank deposits, define arranged by some side direction, at an angle stacking, rich sand be the point bar sand body of feature, and wherein side alternately to amass between body separate by shale short lap (side deposition rock sandwich shape is formed in the flood intermittent phase).The present invention is by intensive well pattern condition, based on study of micro-sedimentary phase achievement, adopt initial stage river course streamline and latter stage river course streamline enveloping method, in conjunction with sand thickness figure, point bar is predicted, have here a committed step be exactly the initial stage and latter stage river course streamline identification and the determination of envelope.In addition, dense well pattern post-mortem method and abandoned channel method can also be adopted to assist and to carry out comprehensive descision.The result that the continuous side of meandering stream is amassed is exactly that sand body increases along with river camber index and constantly to both sides Migration And Concentration, causes sand body to be the non-linear spread of beading, based on such deposition process.
Point bar identification employing streamline envelope method in first latter stage is in conjunction with the thick delineation point bar of sand, the result that the continuous side of meandering stream is amassed is exactly that sand body is continuous to both sides Migration And Concentration along with the increase of river camber index, sand body is caused to be the non-linear spread of beading, water channel grow the first phase time, its streamline must meet two conditions by way of route simultaneously: 1. water channel camber index is minimum the first phase, according to the shortest principle of point-to-point transmission air line distance, initial flow thread path should be the shortest; 2. ensure have sand body to deposit in single water channel.Therefore, ensure above-mentioned two conditions, the initial streamline in river course must adopt tangential manner to connect Sha Hou center with shortest path.
For the Final Issue streamline in latter stage that water channel is grown, just in time streamline initial with river course is contrary, its two conditions met: 1. Final Issue water channel camber index is maximum, larger according to point-to-point transmission curvature of curve, the longer principle in path, and water channel streamline in latter stage path should be the longest; 2. ensure have sand body to deposit in single water channel, therefore, ensure above-mentioned two conditions, river course streamline in latter stage must adopt tangential manner to connect Sha Hou center with longest path.Therefore, can according to the scope of first phase river course streamline and Final Issue river course streamline institute envelope, just a dam location can be determined, this concrete operation method is as shown in Fig. 1 (a), (b), in figure, 1 represents that body is amassed in first phase side, 2 represent that body is amassed in second phase side, 3 represent that body is amassed in third phase side, 4 represent that body is amassed in fourth phase side, 5 expression activity water channels, 6 represent just streamline envelope method delineation in latter stage point bar model, and 7 represent just streamline envelope method delineation in latter stage point bar result schematic diagram, 8 represent Sha Hou center, local schematic diagram, and 9 represent initial stage river course streamline, the flow direction.10 represent latter stage river course streamlines, flow to schematic diagram.
2. the identification of body scale is amassed in point bar side
Point bar side amasss 3 grades that body research level is equivalent in the configuration theory of Miall interface, and it is higher to the requirement of geologic information.It is most important for sedimentary evolution process analysis procedure analysis that the present invention proposes an important parameter, that is exactly the estimation that body arrangement pitches is amassed in side, the present invention, on the basis of quoting the ripe Theories and methods of forefathers, increases the research that body spacing aspect is amassed in side, completes study area point bar building structure fine dissection.
The recognition methods of body scale is amassed in point bar side of the present invention:
(1) extraction of body scale is amassed in side
The extraction that body scale is amassed in side comprises length and the width that body space spread is amassed in side, and wherein asking for of length is also asking for of point bar length, point bar length (W l) be the function of the width (W) of movable water channel, and the width of movable water channel (W) is the function of the thickness (H) in the single river course of the vertical announcement of individual well, namely meets:
W=6.8H 1.54(1)
W: the width of movable water channel, unit m;
H: the vertical sand thickness in single river course, unit m.
Then ask on basis at movable channel width (W), according to (1985) such as Lorenz, by the sample of investigation worldwide meandering stream, sum up point bar length (W l) and movable channel width (W) existence function relation, the length that body (being also the length of single point bar) is amassed in side can be obtained:
W L=85L n(W)+250 (2)
W l: body length is amassed in side, unit m;
W: the width of movable water channel, unit m.
And body width W is amassed in side dasking for of (body width being amassed in side), discloses according to made ground sample results, general W dfollowing functional relation is there is with W:
W d=2/3W (3)
W d: body width is amassed in side, unit m;
W: the width of movable water channel, unit m.
According to above-mentioned calculation process, the length (W that body is amassed in side finally can be obtained l) and width (W d).
(2) extraction of body tendency is amassed in side
Disclose according to made ground phenomenon, the discarded direction that meander river course is always pointed in direction is amassed on meander Ce Ji side, river course, therefore, can according to place, crescent abandoned channel arc length mid point tangent position and circle center line connecting, the tendency that body is amassed in side can be represented, in fact river course amasss convex bank sensing concave bank direction body tendency as side is always pointed in meander Ce Ji direction, river course, it identifies that detailed process is shown in Fig. 2, in Fig. 2,1 represents point bar of meandering model, 2 represent abandoned channel, 3 represent abandoned channel mid point tangent position, 4 represent circular arc tangential line, 5 represent that model is asked at inclination angle, abandoned channel, 6 represent that body tendency (radius) is amassed in side.
(3) Ji Ti inclination angle in side is extracted
It is the geologic element that adjacent two sides are amassed mudstone foundation and limited that body is amassed in side, observe according to outcrop and made ground, it is generally steep lower slow configuration pattern in upper delaying that mudstone foundation is amassed in side, so, when measuring side and amassing the inclination angle of mudstone foundation, generally whole mudstone foundation inclination angle can be replaced by the inclination angle size at middle part.Amass bulk-mode figure according to Fig. 3 side, by the tolerance of relevant parameter, inclination angle can be summarized as follows and extract formula:
θ=arctg(Δh/Δl) (4)
θ: Ji Ti inclination angle, side, unit m;
Δ h: with depth displacement between aboveground lower adjacent interlayers, unit m;
Δ l: antithetical phrase well well spacing, unit m.
First choose two mouthfuls of close together, the well of body can be amassed in complete announcement two continuous sides, top, the river course same period is evened up, then, measure the depth value that mudstone foundation 1 is amassed in well 1 side, record h1, measure well 2 side and amass mudstone foundation 1 degree of depth, the large I of record h2_1, Δ h=h2_1-h1, Δ l two well spacing directly quotes well data can, so above-mentioned Δ h, Δ l two parameters are substituted into (4) formula, Ji Ti inclination angle, side can be calculated.As shown in Figure 3.
(4) extraction of body spacing is amassed in side
In fact body spacing is amassed in side is exactly the interval that side amasss between mudstone foundation nearly river course tip position, it is a very important parameter in reservoir architecture plan view, its acquiring method is on extraction basis, lamination inclination angle, side, according to Fig. 3, increase measures well 2 side and amasss mudstone foundation 2 buried depth again, be recorded as h2_2, then:
Δz=h2_2-h1 (5)
According to (4) formula,
ΔL=Δl+Δz/tgθ (6)
Δ L: body spacing is amassed in side, unit m;
Δ l: antithetical phrase well well spacing, unit m;
Δ z: antithetical phrase aboveground lower two interlayer two well points place drop, unit m;
θ: Ji Ti inclination angle, side;
H2_2: mudstone foundation buried depth is amassed in well 2 side, unit m;
H1: mudstone foundation buried depth is amassed in well 1 side, unit m.
The process of asking for refers to Fig. 3.
(5) the isoparametric extraction of body spacing is amassed by amassing body scale, tendency, inclination angle and side to point bar side, study area, establish point bar building structure identification key parameter database, and based on this, by skeleton section, the point bar of meandering sand body of formation unit carries out fine dissection.The extraction of above-mentioned key parameter, for the fine dissection of point bar of meandering establishes criterion of identification.Carrying out in fine dissection process, concrete operation method is as follows: first chosen in region to be dissected, and Duo Kou well fixed-direction and is concatenated by this region of mistake, forms skeleton section, the section set up is more, and so the point bar precision of portraying is higher, meticulousr.This work on paper section or computer, associated software platforms can complete, as Coreldraw software etc., by the point bar two ends boundary identified in multiple section, join end to end successively, composition graphs 1 point bar delineation PRELIMINARY RESULTS, finally can draw a circle to approve out point bar spatial scope.On point bar spatial scope delineation basis, body is amassed to the side of point bar inside and carries out meticulous depiction description, concrete operation method is: portraying from section and plane two aspect, is first that section is portrayed, and has just used Ji Ti inclination angle, side (θ), length (W that body spread is amassed in side l) and width (W d) these 3 parameters, the identification of section can adopt two steps: 1. when section is vertical with paleocurrent direction, the width (W that what use is Ji Ti inclination angle, side (θ), body spread is amassed in side d) these two parameters, Ji Tiding bottom surface, single side is obvious in well log curve feature, namely significantly " recessed " phenomenon is had at side Ji Tidingdi circle critical part log, changes in amplitude is unexpected herein to show as log, then with this top bottom interface identified for starting point, follow the trail of to offset well place, in tracking comparison process, amass the width (W of body spread above by the Ji Ti inclination angle, side (θ) estimated and side d) as borders condition, body identification is amassed in the side finally completing vertical paleocurrent direction section; 2. when section is parallel with paleocurrent direction, now Ji Ti inclination angle in side is in being bordering on level, and side is amassed body width and just can not be estimated with formula (1) and formula (3).Certainly, this is an extreme case.In addition, if the skeleton section set up and paleocurrent out of plumb, the inclination angle apparent dip then asked for, but not true dip angle, and less than true dip angle, the width of body is amassed also than proper length slightly length, therefore when portraying section in the side of estimation, choose the section vertical with paleocurrent direction to portray as far as possible, the error of point bar meticulous depiction could be reduced.According to aforesaid operations method, point bar side can be completed and amass identification on body section on transverse and longitudinal section, then the result will identified, with reference to the method for point bar plane combination, body tendency estimation result is amassed as constraints in conjunction with side, carry out plane combination, finally complete the fine dissection of point bar, see point bar meticulous resulting schema schematic diagram 4a.
According to point bar scale result of calculation, seek out width and the length of point bar extension, like this in the identification of section point bar and comparison process, the scope of point bar spatial extent has just had corresponding edge-restraint condition, avoids point bar recognition result bigger than normal or less than normal.
3. river course evolutionary process of meandering and history are rebuild
The same with structural history, river course also exists from being born to final evolution process of withering away, and therefore by analysis and the history reconstruction of its evolution process, contributes to carrying out fine dissection to sand body building structure, and gives rational explanation from the angle of the origin cause of formation.At present not to the method how river course streamline extracts and build.Each issue side is amassed body and is identified by the present invention, solve the problem of river channel evolution, because every first phase side amass body unit represent one time flood period arenaceous sediment event, and side amass between body side deposition rock interlayer representative be rest water phase argillaceous sediment event, therefore, can amass with single side that the extension track in whole piece river course of body is approximate replaces first phase river course streamline, thus just overcome that channel thread line position cannot be estimated and the problem such as river course interpolation main stream line error is larger.
To rebuild channel deposit Evolutionary History, just first must amass that body identifies to point bar side, fine dissection and combination.Obtain river course streamline, then restore successively by growth order and arrange, to reproduce channel deposit history of evolution.The determination of river course initial position, the envelope method based on the thick figure of sand can be adopted, can effectively identify river course initial position and latter stage river location, middle river course streamline amasss with the side same period track that body is combined to form and replaces, and so just can carry out effective restoration and reconstruction to river course sedimentary evolution history of meandering.
On point bar delineation basis, body scale is amassed by step 2 offside, the extraction of the relevant parameters such as the extraction at tendency inclination angle and spacing, implement point bar side and amass body plane growth and Distribution Pattern, as Fig. 4 (a), (b), (c), (d), (e), (f), shown in (g), in Fig. 4,1 amasss body for point bar 1 the 1st phase side, 2 amass body for point bar 1 the 2nd phase side, 3 amass body for point bar 1 the 1st phase side, 4 amass body for point bar 1 the 4th phase side, 1_1 is that body is amassed in point bar 2 the 1st phase side, 2_2 is that body is amassed in point bar 2 the 2nd phase side, 3_3 is that body is amassed in point bar 2 the 3rd phase side, 4_4 is that body is amassed in point bar 2 the 4th phase side, 5 is that body contrast same period connection schematic lines is amassed in the 1st phase side, 6 is that body contrast same period connection schematic lines is amassed in the 2nd phase side, 7 is that body contrast same period connection schematic lines is amassed in the 3rd phase side, 8 is that body contrast same period connection schematic lines is amassed in the 4th phase side.Now choose two adjacent point bars 1 and point bar 2 (see Fig. 4 a), on the basis that step 2 point bar internal construction quantitatively characterizing parameter is set up (see Fig. 4 a), body is amassed to these two the inner sides identified of adjacent point bar and carries out the combination same period (see Fig. 4 b), " adjacent similar " principle is taked to combine, combined result shows, two contiguous sides are amassed body and are meaned the same period, as body 1 is amassed in point bar 1 first phase side, to amass body 1_1 with point bar 2 side be the same period, therefore, restoration result as illustrated in fig. 4 c, by that analogy, respectively body unit combination is amassed in the side same period of follow-up some phases to identify (Fig. 4 d → e), by carrying out accumulating by the phase to above-mentioned each issue, finally complete and the evolutionary process of point bar 1 and point bar 2 entirety and history are rebuild (Fig. 4 g).
Make full use of the data such as core, well logging, dense well pattern, on the basis such as matching and estimation of the identification of abandoned channel and point bar sand body, recovery, river work parameter, in conjunction with former achievements, differentiate that the scale of body, tendency and inclination angle are amassed in point bar side, and based on this, in conjunction with point bar geologic concept model, in single phase river course multiple point bar, body is amassed in the side belonging to the same period identify, then be numbered combination, carry out process analysis procedure analysis according to the rule of river course meander and reproduce its sedimentary evolution process.Such an approach achieves from cause analysis angle to explain point bar interior building structure and Spatial Overlap rule problem thereof, fluvial depositional reservoir structural analysis is made not only to rest on result aspect, but focus on going restoration and reconstruction meander river course Historical Evolution process from the origin cause of formation, process visual angle, make point bar of meandering inner sand body building structure anatomical results will more rationally, more credible, and effectively direct the prediction of oil field middle and later periods remaining oil and tap the latent power
At present, for meandering river facies reservoir, most domestic oil field is all in mid-later development phase, generally there is intensive well pattern condition, therefore, create condition to the fine dissection of point bar of meandering, for the dissection of point bar, the present invention proposes the method for a whole set of Quantitative anatomy point bar, simple easily operation.In addition, reproduce for final the temporal and spatial evolution problem that body is amassed in point bar side, this invention proposes to amass the identification of body with point bar same period side and combined method is rebuild, and scheme is easy-to-understand, easy to operate and effective.
The above is only to better embodiment of the present invention, not any pro forma restriction is done to the present invention, every any simple modification done above embodiment according to technical spirit of the present invention, equivalent variations and modification, all belong in the scope of technical solution of the present invention.
The present invention will be described to enumerate specific embodiment below:
Embodiment 1: for Jilin Oil Field poplar big city seed oil layer group q4-17 individual layer, be example to point bar of meandering, carries out evolutionary process and history is rebuild.
(1) just, river course streamline envelope method identification in latter stage point bar
In the isopach map of sand bodies that single meander river course is main, sand thickness center is many in non-linearity, asymmetric beading spread, instead of the bar-shape spread in anticipation, this may meander continuous with meandering stream, the enrichment of sand body continuous lateral migration is relevant, causes sand body thus along axis at bilateral asymmetry spread.And for the meandering river sandbody of bar-shape spread, it has more in present combination channel isopath map, may be the stacked cumulative result in river course of many phases.According to envelope rule, identify study area q4-17 individual layer point bar, result have identified 7 point bars, in detail as shown in Figure 5, in Fig. 5,1 represents the thick isopleth of sand, and 2 represent Sha Hou center, local signal scope, 3 represent river course streamline in latter stage, and 4 represent initial stage river course streamline, and 5 represent well location.
(2) extraction of body scale is amassed in side
According to the method that step 2 is invented, according to formula (2), in conjunction with known intensive well pattern and phased data, try to achieve abandoned channel (being equivalent to movable water channel in latter stage) width and be about 100m, bring above-mentioned formula into, try to achieve point bar side to amass body length (being also point bar length) and be about 650m, and body width calculation is amassed according to formula (3) in side, can try to achieve side and amass body length W lbe approximately 67m.
(3) extraction of body tendency is amassed in side
According to the recognition methods of the proposition in step 2, body tendency is amassed in the point bar side that can extract study demonstration district, and tendency extraction the results are shown in Figure 6.
(4) Ji Ti inclination angle in side is extracted
Lamination inclination angle, side is extracted, and that generally commonly uses mainly contains 4 kinds of methods, is rock core method, antithetical phrase well technology, discarded surface technology and continous sampling stage extraction technology (Ma Shizhong etc., 2008) respectively.Rock core method and continous sampling stage extraction technology need to get the dark mud rock section that core section has tool horizontal bedding, so that carry out evening up correction, to eliminate the impact of later structural factor.This research is according to the situation of real data, mainly take discarded face and antithetical phrase well technology, because face of amassing, side is not a plane, but present slow-in steep-lower slow concave surface, and in slow be its approximate average, so the Ji Ti inclination angle, replacement side, inclination angle of Middle face is amassed in general side.Forefathers are by outcrop and made ground and case study, draw average 8 ° of top, Ji Ti inclination angle, side, 12 °, middle part and 20 °, bottom (Ma Shizhong etc., 2008), integrating step 2 formula (4) of the present invention, finally asking for its inclination angle average is arctg (3.0/17.0) ≈ 10 °, sees Fig. 7 in detail.
(5) extraction of body spacing is amassed in side
In fact body spacing is amassed in side is exactly the interval that side amasss between mudstone foundation nearly river course tip position, it is a very important parameter in reservoir architecture plan view, according to step 2 formula (5), formula (6), corresponding data is substituted into, side can be tried to achieve amass body spacing and be approximately 31.2m, see Fig. 7.
(6) body fine dissection is amassed in point bar side
The isoparametric extraction of body spacing is amassed by amassing body scale, tendency, inclination angle and side to point bar side, study area, establish point bar building structure identification key parameter database, and based on this, by nearly 20 skeleton sections, fine dissection is carried out to the meandering river sandbody of study area the 17th substratum stratigraphic unit, fine dissection is carried out to 7 that have identified main point bar internal constructions, has referred to table 1 and Fig. 8.In the point bar unit of Z24 wellblock, wherein have identified 7 sides amass body, long-pending to side, Z24-Z2-12-4 well direction, finally sentence at Z4-1 well and gradually abandon abandoned channel termination.In addition, the point bar of growing in Z2-14-9 wellblock, when body unit is amassed in deposition the 6th side, by cutting bending straight effect, this place's point bar is subject to discarded, thus finishes its deposit development history, defines remaining type point bar.
(7) river channel evolution process and history are rebuild
The same with structural history, river course also exists from being born to final evolution process of withering away, and therefore by analysis and the history reconstruction of its evolution process, contributes to carrying out fine dissection to sand body building structure, and gives rational explanation from the angle of the origin cause of formation.
Body structure parameter quantitative characterization experiments data are amassed in table 1 point bar 1 side
This invention of river course streamline think time river course core portion following current not same period by way of route.If body can be amassed in each issue side to identify, the problem of river channel evolution also can be readily solved, because every first phase side amass body unit represent one time flood period arenaceous sediment event, and side amass between body side deposition rock interlayer representative be rest water phase argillaceous sediment event, therefore, can amass with single side that the extension track in whole piece river course of body is approximate replaces first phase river course streamline, thus just overcome that channel thread line position cannot be estimated and the problem such as river course interpolation main stream line error is larger.
For poplar big city, oil recovery factory of Fuyu County seed oil layer group the 17th single layer cell, integrated use is based on the thick envelope method of sand, the same period, body combined method was amassed in side, study area river channel evolution history is recovered, reconstructed results as shown in Figure 9, a material time node by a → f, at this section of depositional stage, the normal meander in river course, and arrived f → g, bottom-right meander river course is the 6th meander ring wherein, flood period may be in, hydrodynamic force strengthens, in addition river course camber index is herein excessive, cause Breach and diversion, the result that neck cuts bending straight effect is that the point bar at this meander ring place is discarded, current after changing its course start new round meander river course history of evolution at side.

Claims (3)

1. the buried point bar of meandering sand body in a underground history reconstruction method, is characterized in that carrying out according to following steps:
Step 1: adopt initial stage river course streamline and latter stage river course streamline enveloping method, in conjunction with sand thickness figure, point bar is identified;
Step 2: the identification of body scale is amassed in point bar side;
(1) extraction of body scale is amassed in side:
The extraction that body scale is amassed in side comprises length and the width that body space spread is amassed in side, and wherein asking for of length is also asking for of point bar length, point bar length (W l) be the function of the width (W) of movable water channel, and the width of movable water channel (W) is the function of the thickness (H) in the single river course of the vertical announcement of individual well, namely meets:
W=6.8H 1.54
W: the width of movable water channel, unit m;
H: the vertical sand thickness in single river course, unit m;
The length of body is amassed in side:
W L=85Ln(W)+250;
W l: body length is amassed in side, unit m;
W: the width of movable water channel, unit m;
And body width W is amassed in side din the asking for of (body width is amassed in side), W dfollowing functional relation is there is with W:
W d=2/3W;
W d: body width is amassed in side, unit m;
W: the width of movable water channel, unit m;
(2) extraction of body tendency is amassed in side:
According to place, crescent abandoned channel arc length mid point tangent position and circle center line connecting, amass the tendency of body as side;
(3) Ji Ti inclination angle in side is extracted:
Formula is extracted at inclination angle:
θ=arctg(Δh/Δl);
θ: Ji Ti inclination angle, side, unit m;
Δ h: with depth displacement between aboveground lower adjacent interlayers, unit m;
Δ l: antithetical phrase well well spacing, unit m;
(4) extraction of body spacing is amassed in side:
It is the interval that side amasss between mudstone foundation nearly river course tip position that body spacing is amassed in side, and its acquiring method is on extraction basis, lamination inclination angle, side, and increase measures well side and amasss mudstone foundation buried depth, is recorded as h2_2, then:
Δz=h2_2-h1
ΔL=Δl+Δz/tgθ
Δ L: body spacing is amassed in side, unit m;
Δ l: antithetical phrase well well spacing, unit m;
Δ z: antithetical phrase aboveground lower two interlayer two well points place drop, unit m;
θ: Ji Ti inclination angle, side;
H2_2: mudstone foundation buried depth is amassed in well 2 side, unit m;
H1: mudstone foundation buried depth is amassed in well 1 side, unit m;
(5) by amassing to point bar side, study area the extraction that body spacing parameter is amassed in body scale, tendency, inclination angle and side, set up point bar building structure identification key parameter database, by skeleton section, the meandering river sandbody of formation unit carries out fine dissection;
Step 3: meander river course evolutionary process and history are rebuild;
Amass that body identifies by step 2 pair point bar side, fine dissection and combination, obtain river course streamline, then restore successively by growth order and arrange, reproduce channel deposit history of evolution.
2. according to underground a kind of described in claim 1 buried point bar of meandering sand body history reconstruction method, it is characterized in that: in described step 1, judging point dam location method is the first phase that water channel is grown, its streamline connects sand body depocenter with shortest path, and Final Issue water channel streamline path is the longest, the scope of first phase river course streamline and Final Issue river course streamline institute envelope is exactly point bar present position.
3. according to underground a kind of described in claim 1 buried point bar of meandering sand body history reconstruction method, it is characterized in that: in described step 3, the determination of river course initial position, the envelope method based on the thick figure of sand can be adopted, can effectively identify river course initial position and latter stage river location, middle river course streamline amasss with the side same period track that body is combined to form and replaces, and so just can carry out effective restoration and reconstruction to river course sedimentary evolution history of meandering.
CN201410618134.XA 2014-11-05 2014-11-05 A kind of buried point bar of meandering sand body history method for reconstructing in underground Active CN104453877B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410618134.XA CN104453877B (en) 2014-11-05 2014-11-05 A kind of buried point bar of meandering sand body history method for reconstructing in underground

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410618134.XA CN104453877B (en) 2014-11-05 2014-11-05 A kind of buried point bar of meandering sand body history method for reconstructing in underground

Publications (2)

Publication Number Publication Date
CN104453877A true CN104453877A (en) 2015-03-25
CN104453877B CN104453877B (en) 2018-05-15

Family

ID=52900549

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410618134.XA Active CN104453877B (en) 2014-11-05 2014-11-05 A kind of buried point bar of meandering sand body history method for reconstructing in underground

Country Status (1)

Country Link
CN (1) CN104453877B (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105550525A (en) * 2015-12-30 2016-05-04 中国科学院遥感与数字地球研究所 Remote sensing technology based ancient wind strength quantification reconstruction method
CN105607146A (en) * 2015-09-10 2016-05-25 中国海洋石油总公司 Quantitative characterization method of meandering river sand body scale
CN108119133A (en) * 2017-11-13 2018-06-05 中国石油天然气股份有限公司 Method and device for determining reservoir deposition process
CN108304627A (en) * 2018-01-18 2018-07-20 中海石油(中国)有限公司上海分公司 A kind of modeling method for portraying river difference cutting mode
CN109025981A (en) * 2018-06-22 2018-12-18 中国石油天然气股份有限公司 Method and system for comparing side accumulated sand bodies of meandering stream
CN110858260A (en) * 2018-08-17 2020-03-03 中国石油化工股份有限公司 Method for describing plane distribution of sand bodies of single-stage river channel of delta river facies
CN112068213A (en) * 2020-09-21 2020-12-11 中国石油天然气股份有限公司 River delta facies sand body configuration prediction method
CN112147714A (en) * 2019-06-28 2020-12-29 中国石油化工股份有限公司 Estuary dam internal structure analysis method
CN112527939A (en) * 2020-12-17 2021-03-19 天津泰勘工程技术咨询有限公司 Urban area ancient river channel identification method
CN112651547A (en) * 2020-12-11 2021-04-13 中国石油天然气股份有限公司 Curve river sandstone reservoir scale prediction method
CN113536520A (en) * 2020-04-22 2021-10-22 中国石油天然气股份有限公司 Method and device for simulating meandering stream river channel migration
CN114706141A (en) * 2022-04-06 2022-07-05 东北石油大学 Underground meandering river zone quantitative characterization method based on sedimentary analogy

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1944956A (en) * 2006-10-31 2007-04-11 大庆油田有限责任公司 Fluvial phase sand building structure identifying method
US20080126048A1 (en) * 2006-11-23 2008-05-29 Total S.A. Method, program and computer system for paleochannel modeling
EP2525242A2 (en) * 2011-05-20 2012-11-21 Baker Hughes Incorporated Multiscale geologic modeling of a clastic meander belt including asymmetry using multi-point statistics
CN103226847A (en) * 2013-04-01 2013-07-31 中国石油天然气股份有限公司 Method and device for generating three-dimensional data volume of point dam of meandering river
CN103699751A (en) * 2013-12-30 2014-04-02 中国石油大学(北京) Sand body reservoir architecture modeling method and system based on space vectors
CN103942842A (en) * 2014-03-19 2014-07-23 中国石油天然气股份有限公司 Embedded meandering stream sand body modeling method
CN104018828A (en) * 2014-05-20 2014-09-03 长江大学 Meander river sand body reservoir building structure analysis method based on evolution process

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1944956A (en) * 2006-10-31 2007-04-11 大庆油田有限责任公司 Fluvial phase sand building structure identifying method
US20080126048A1 (en) * 2006-11-23 2008-05-29 Total S.A. Method, program and computer system for paleochannel modeling
EP2525242A2 (en) * 2011-05-20 2012-11-21 Baker Hughes Incorporated Multiscale geologic modeling of a clastic meander belt including asymmetry using multi-point statistics
CN103226847A (en) * 2013-04-01 2013-07-31 中国石油天然气股份有限公司 Method and device for generating three-dimensional data volume of point dam of meandering river
CN103699751A (en) * 2013-12-30 2014-04-02 中国石油大学(北京) Sand body reservoir architecture modeling method and system based on space vectors
CN103942842A (en) * 2014-03-19 2014-07-23 中国石油天然气股份有限公司 Embedded meandering stream sand body modeling method
CN104018828A (en) * 2014-05-20 2014-09-03 长江大学 Meander river sand body reservoir building structure analysis method based on evolution process

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
岳大力 等: "曲流河古河道储层构型精细解剖——以孤东油田七区西馆陶组为例", 《地学前缘》 *
张昌民 等: "基于过程的分流平原高弯河道砂体储层内部建筑结构分析——以大庆油田萨北地区为例", 《沉积学报》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105607146A (en) * 2015-09-10 2016-05-25 中国海洋石油总公司 Quantitative characterization method of meandering river sand body scale
CN105607146B (en) * 2015-09-10 2017-10-03 中国海洋石油总公司 A kind of quantitatively characterizing method of meandering river sand body scale
CN105550525B (en) * 2015-12-30 2018-09-21 中国科学院遥感与数字地球研究所 A kind of ancient customs power quantification of intensities method for reconstructing based on remote sensing technology
CN105550525A (en) * 2015-12-30 2016-05-04 中国科学院遥感与数字地球研究所 Remote sensing technology based ancient wind strength quantification reconstruction method
CN108119133A (en) * 2017-11-13 2018-06-05 中国石油天然气股份有限公司 Method and device for determining reservoir deposition process
CN108304627B (en) * 2018-01-18 2021-10-15 中海石油(中国)有限公司上海分公司 Modeling method for depicting different cutting modes of river channel
CN108304627A (en) * 2018-01-18 2018-07-20 中海石油(中国)有限公司上海分公司 A kind of modeling method for portraying river difference cutting mode
CN109025981A (en) * 2018-06-22 2018-12-18 中国石油天然气股份有限公司 Method and system for comparing side accumulated sand bodies of meandering stream
CN110858260A (en) * 2018-08-17 2020-03-03 中国石油化工股份有限公司 Method for describing plane distribution of sand bodies of single-stage river channel of delta river facies
CN112147714A (en) * 2019-06-28 2020-12-29 中国石油化工股份有限公司 Estuary dam internal structure analysis method
CN113536520A (en) * 2020-04-22 2021-10-22 中国石油天然气股份有限公司 Method and device for simulating meandering stream river channel migration
CN112068213A (en) * 2020-09-21 2020-12-11 中国石油天然气股份有限公司 River delta facies sand body configuration prediction method
CN112651547A (en) * 2020-12-11 2021-04-13 中国石油天然气股份有限公司 Curve river sandstone reservoir scale prediction method
CN112651547B (en) * 2020-12-11 2024-05-28 中国石油天然气股份有限公司 Method for predicting size of sandstone reservoir of curtreta
CN112527939A (en) * 2020-12-17 2021-03-19 天津泰勘工程技术咨询有限公司 Urban area ancient river channel identification method
CN114706141A (en) * 2022-04-06 2022-07-05 东北石油大学 Underground meandering river zone quantitative characterization method based on sedimentary analogy

Also Published As

Publication number Publication date
CN104453877B (en) 2018-05-15

Similar Documents

Publication Publication Date Title
CN104453877A (en) Underground deep-buried meander point dam sand body historical reconstruction method
Yang et al. Theories and practices of carbonate reservoirs development in China
Federici et al. Glacial history of the Maritime Alps from the last glacial maximum to the little ice age
CN103472484B (en) Horizontal well path optimization method based on RS three-dimensional sensitive earthquake attributive analysis
CN104749638A (en) Determining method of complex mountain-front effective source rock and source rock structural model
CN103454686A (en) Method and system for reservoir prediction based on small-scale sedimentary facies of stratigraphic slice
CN104018828A (en) Meander river sand body reservoir building structure analysis method based on evolution process
CN105319585B (en) A kind of method hidden using thin-layers interference amplitude recovery identification of hydrocarbon
Khomsi et al. Structural control on the deep hydrogeological and geothermal aquifers related to the fractured Campanian-Miocene reservoirs of north-eastern Tunisia foreland constrained by subsurface data
CN105064990A (en) Structure modeling method by means of horizontal well three-dimensional visualization stratigraphic correlation
CN110489892B (en) Three-phase coupling reservoir structure description technical method
CN104533402A (en) Multi-phase superimposed combined braided channel sand phase determining method
CN104453903A (en) Water-preserved mining method of close-distance coal seam group
CN104018827A (en) Delta sand body internal structure analysis method based on overlapped lobes
Taixian et al. Techniques for high-efficient development of offshore fluvial oilfields
CN110428497A (en) Braided stream training image generation method
Liu et al. The Control Theory and Application for Well Pattern Optimization of Heterogeneous Sandstone Reservoirs
Grauch et al. A shifting rift—Geophysical insights into the evolution of Rio Grande rift margins and the Embudo transfer zone near Taos, New Mexico
Pradhan et al. Shale gas play: possibility and challenges in India as an unconventional energy resource
CN109031443B (en) Method for recovering denudation amount based on fixed years of Milnaciaceae Virginia gyrus
CN107991702B (en) The simple accurate determination method of normal fault timing
Chiocchini et al. The Viterbo Hydrothermal System and Its Sustainable Exploitation, Central Italy
CN204882902U (en) Effective hydrocarbon source rock structure model of complicated mountain front
Caddick et al. Global Overview of Fractured Basement Plays.
CN114594518B (en) Fine stratum contrast method for complex fault blocks in later development period based on well-seismic alternation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20220817

Address after: Room 801-82, Unit B, 8th Floor, Building 33, No. 99, Kechuang 14th Street, Beijing Economic and Technological Development Zone, Beijing 100000 (Beijing Pilot Free Trade Zone High-end Industry Zone Yizhuang Group)

Patentee after: Beijing furuibao Energy Technology Co.,Ltd.

Address before: 434023 No. 1 South Ring Road, Hubei, Jingzhou

Patentee before: Yangtze University

TR01 Transfer of patent right