CN104749623B - A kind of imaging of seismic data processing method - Google Patents
A kind of imaging of seismic data processing method Download PDFInfo
- Publication number
- CN104749623B CN104749623B CN201410587602.1A CN201410587602A CN104749623B CN 104749623 B CN104749623 B CN 104749623B CN 201410587602 A CN201410587602 A CN 201410587602A CN 104749623 B CN104749623 B CN 104749623B
- Authority
- CN
- China
- Prior art keywords
- ripples
- gaussian beam
- imaging
- elastodynamics
- formula
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 43
- 238000003672 processing method Methods 0.000 title claims abstract description 6
- 238000013508 migration Methods 0.000 claims description 12
- 230000005012 migration Effects 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 8
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000010287 polarization Effects 0.000 claims description 2
- 230000000630 rising effect Effects 0.000 claims description 2
- 238000000034 method Methods 0.000 description 19
- 238000012937 correction Methods 0.000 description 8
- 238000000926 separation method Methods 0.000 description 6
- 238000012876 topography Methods 0.000 description 5
- 238000007796 conventional method Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011230 binding agent Substances 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000001615 p wave Methods 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004613 tight binding model Methods 0.000 description 1
Abstract
The invention provides a kind of imaging of seismic data processing method, according to the dominant frequency of seismic data and frequency range selection reference frequency and highest frequency, and then vector Gaussian beam original width is determined according to exploratory area P- and S-wave velocity field size, so as to give initial ray parameter space；It is then based on earth's surface obliquity information and elastodynamics Gaussian beam is directly emitted at focus and receiving point, and carry out kinematics ray tracing and kinetics ray-tracing, during asks for the travelling of elastodynamics Gaussian beam and amplitude information；Finally, according to the travelling of elastodynamics Gaussian beam when, amplitude and polarity information, corresponding single track PP ripples and PS ripples are calculated into picture value using cross-correlation image-forming condition, and then single-shot PP ripples and PS ripples are tried to achieve into picture value to per pass circulation, the PP ripples of all big guns are overlapped respectively with PS ripples into picture value and just obtain final PP ripples and PS ripple imaging sections.
Description
Technical field
The invention belongs to seism processing field, specifically a kind of imaging of seismic data processing method.
Background technology
With deepening continuously for exploration and development, the center of gravity of oil-gas exploration, which just turns to, has complicated landform and complicated subsurface geology
Double complex regions of body, geological data also gradually develop into elastic wave data by P wave data.However, near surface elevation and speed
Spend under the MODEL OVER COMPLEX TOPOGRAPHY of acute variation, conventional static correction distorts to caused by wave field can have a strong impact on migration imaging quality.
Because multi-wave and multi-component data comprehensively and more can truly reflect underground medium and formation lithology information, develop double complicated bars
The method that Elastic Wave Migration imaging is directly carried out under part has been increasingly becoming the focus in seismic data process field.
To solve influence of the MODEL OVER COMPLEX TOPOGRAPHY to migration imaging, many scholars have carried out substantial amounts of research.The first kind is
Complicated landform offset imaging method based on ACOUSTIC WAVE EQUATION.Reference plane continuation of the Wiggins (1984) based on geological data, propose
It is applied to the Kirchhoff offset methods of MODEL OVER COMPLEX TOPOGRAPHY, and Alkhalifah and Bagaini is equivalent by setting one
Speed carries out base-level correction；Gray (2005) then proposes a kind of local static correction Gaussian beam suitable for complicated earth surface and offset
Method.Berryhill(1979；1984) method for then proposing wave equation datuming.Then, Yilmaz and Lucas
(1986), Schneider and Pillip etc. (1995), binder course replaces on this basis by Bevc (1997) and Yang etc. (1999)
The method is further developed for thought.Above method is all based on " skew indirectly " imaging method of base-level correction, although
Base-level correction can effectively eliminate relief surface and be distorted to caused by lineups, but the cross directional variations of near-surface velocity can be tight
Ghost image rings the precision of base-level correction.The defects of for base-level correction, Reshef (1991), Beasley and Lynn (1992),
The bright grade (2008) of He Ying etc. (2002), Shragge (2005), Cheng Jiubing (2006), Ye Yue and Lv Bin etc. (2011) successively carry
A variety of " directly skew " imaging methods based on wave equation are gone out to solve the problems, such as complicated earth surface.Xu Yi (2008) then uses triangle
Gridding method realizes sound wave reverse-time migration in the case of relief surface.Liu Hongwei etc. (2011) is to earth's surface free boundary using filtering
Method, relief surface sound wave reverse-time migration is realized with GPU.However, under complex near surface conditionss, ray class method have concurrently flexibly,
The advantages that efficient.Gray etc. (1995) has developed directly carries out the method, Jager etc. of Kirchhoff skews in complicated earth surface
(2003) (2009) and Liu Guofeng etc. such as guarantor's width Kirchhoff offset methods, Dong Chunhui under the conditions of relief surface are proposed
(2010) flow and the application of the direct time migration of relief surface are successively proposed based on Kirchhoff operators.Yue etc. (2010)
Propose a kind of method for directly carrying out protecting the skew of width Gaussian beam in relief surface.Conventional Gaussian beam migration operator described above
All grown up according to local dip superposition theory, Huang Jianping etc. (2014) has then been developed a kind of double multiple based on ACOUSTIC WAVE EQUATION
Non-inclined is superimposed accurate beam offset imaging method under the conditions of miscellaneous, and has obtained preferable imaging effect.Second class is based on elasticity
The complicated landform offset imaging method of wave equation.During traditional multi-component earthquake data processing, ripple is often first passed through
Field separation, is then individually imaged to the P ripples after separation and S ripples.However, remnants during wave field separation often be present
Non- type wave energy, the much noise in imaging results can be caused, have a strong impact on image quality.In recent years, domestic and foreign scholars
Based on horizontal surface conditions, it is proposed that elastic wave reverse-time migration and elastic wave Kirchhoff skews.Yue Yubo (2011) is then proposed
Elastic wave Gaussian beam offset method, this method have taken into account imaging precision and computational efficiency.And the bullet based on MODEL OVER COMPLEX TOPOGRAPHY
Property ripple offset imaging method research be still in starting stage, Huang Jianping etc. (2014) according to local dip superposition theory, propose
A kind of Gaussian beam offset method based on non-separation multi-component earthquake data under the conditions of relief surface.However, due to beam center
There are horizontal range difference and depth displacement between position and receiving point, the method needs to introduce phase correction factor with approximation sign Green
Function, when horizontal range difference is larger, amplitude error caused by this approximation can not be ignored, especially when landform big rise and fall and closely
When earth's surface velocity variations are violent, above-mentioned approximate processing can greatly reduce the effect of migration imaging.
The content of the invention
It is an object of the invention at western part of China and southern mountain front exploratory area multiwave multicomponent earthquake data imaging
The problem of reason, it is proposed that a kind of imaging of seismic data processing method, in the elastic wave height that double complex conditions are handled based on single track
This beam Depth Domain imaging technique, using non-separation multiwave multicomponent earthquake data and p-and s-wave velocity field as input, utilize earth's surface
Obliquity information, the backward extension vector wave field of decoupling is characterized by being emitted elastodynamics Gaussian beam at each receiving point
(longitudinal and transverse wave field), and picture value is asked for into using cross-correlation image-forming condition.
The present invention comprises the following steps：
(1) reference frequency ω is chosen according to the dominant frequency of multiwave multicomponent earthquake data and frequency bandwidthrAnd highest frequency
ωh；
(2) according to following criterion, elastodynamics Gaussian beam original width is calculated：
(formula 1)
In formula：It is the geometrical mean of V-shaped ripple (compressional wave or shear wave) velocity field,
(3) V-shaped ripple initial ray parameter space is：
(formula 2)
(4) earth's surface elevation, inclination angle and p-and s-wave velocity information are based on, and with reference to selected by above-mentioned 1,2,3 steps and is counted
The parameter of calculation, using single track processing mode, elastodynamics Gaussian beam is directly emitted at focus and receiving point respectively, and pass through
Solve kinematics ray tracing equation 3 and kinetics ray-tracing equation 4 and then elastodynamics Gaussian beam is calculated
Travelling when and amplitude information, then try to achieve the backward extension vector wave field decoupled under complex near surface conditionss, as shown in Equation 5；
(formula 3)
(formula 4)
(formula 5)
In formula, wherein, W1 v(xr)、W2 v(xr) it is the weights relevant with p-and s-wave velocity and polarization vector,
For receiving point xrDisplacement caused by place at x, ρ (xr) it is xrLocate Media density, S is the free earth's surface of rising and falling, u1(xr, ω) and u2
(xr, ω) and it is respectively horizontal and vertical component seismic record；
(5) PP ripples and PS ripples based on single track processing under the complex near surface conditionss derived according to cross-correlation image-forming condition into
As formula, using elastodynamics Gaussian beam travelling when, amplitude and polarity information calculate corresponding single track into picture value；
(6) the computing flow of 4,5 steps is repeated, calculate single-shot into picture value, most all big guns are folded into picture value at last
Add to obtain required imaging section, the PP ripples and PS ripple single-shot migration imaging formula such as formula that double complex conditions are handled based on single track
Shown in 6a-6b：
(formula 6a)
(formula 6b)
In formula,It is to correspond to V-shaped ripple in xrThe initial slowness vector at place, φ are that P ripples enter
Firing angle degree, sgn (φ) are sign function；
Single track processing side is utilized during the beneficial effect of the present invention vector wave field that backward extension decouples under MODEL OVER COMPLEX TOPOGRAPHY
Formula, it can effectively avoid amplitude error caused by the common elastic ripple Gaussian beam skew based on local dip superposition theory, divide
Resolution error and noise problem, it is a kind of high-precision migration and imaging techniques based on non-separation multiwave multicomponent earthquake data.
Brief description of the drawings
Fig. 1 a and Fig. 1 b are respectively Zhongyuan Oil Field FAULT MODEL compressional wave and shear wave velocity field.This model has typical small hill
Cheuch complicated earth surface, and the medium and deep of model includes well-developed fault tectonic, is the target area of oil-gas exploration.Model
Parameter is as follows：Model meshes are 436 × 600, and vertically and horizontally sampling interval is respectively 4 meters and 10 meters；
Fig. 2 a and 2b are respectively corresponding single-shot z-component and x-component record.Earthquake big gun data 157 big gun altogether, big gun is at intervals of 20
Rice, Mei Bao 121 are received, and 4000 sampled points of per pass, sample rate is 0.5 millisecond, and road is at intervals of 10 meters；
Fig. 3 a and 3c are respectively the PP ripples and PS ripple imaging sections of conventional method；
Fig. 3 b and 3d are respectively using the PP ripples after the present invention and PS ripple imaging sections.
Embodiment
For make the present invention above and other purpose, feature and advantage become apparent, it is cited below particularly go out preferable implementation
Example, and coordinate institute's accompanying drawings, it is described in detail below：
(1) P- and S-wave velocity model (Fig. 1 a-1b) and corresponding unsegregated multi-wave and multi-component big gun record (Fig. 2 a- are chosen
2b) as input.Rate pattern grid is 436 × 600, and transverse and longitudinal spacing distance is respectively 4 meters and 10 meters.Earthquake big gun data are total to
Have 157 big guns, big gun is at intervals of 20 meters, and Mei Bao 121 is received, 4000 sampled points of per pass, and sample rate is 0.5 millisecond, road at intervals of
10 meters；
(2) suitable reference frequency 10Hz and highest frequency are chosen according to the dominant frequency of earthquake big gun data and frequency bandwidth
50Hz, and then calculate initial beamwidth and initial ray parameter space according to formula 1 and formula 2；
(3) focus and reception are determined according to the observation system described in the information such as earth's surface elevation, inclination angle and above-mentioned steps 1
Point is in the position of complicated earth surface, then direct outgoing elastodynamics Gaussian beam and when calculating its travelling and amplitude information；
(4) when being travelled according to the elastodynamics Gaussian beam asked in step 3, amplitude information and its polarity information, application
Cross-correlation image-forming condition calculates corresponding single track PP ripples and PS ripples into picture value；
(5) repeat step 3 and 4, single-shot PP ripples and PS ripples are calculated into picture value, most at last the PP ripples of all big guns and PS ripples into
Picture value is overlapped to obtain required PP ripples and PS ripple imaging section Fig. 3 b and Fig. 3 d respectively.
It can be seen that from Fig. 3 b and 3d：(1) the either imaging of PP ripples or the imaging of PS ripples clearly reflects rough ground
The well-developed fault tectonic of table and underground；(2) because conversion S wave velocities are smaller than p wave interval velocity thus same in PS ripple imaging results
Phase axle has the longitudinal frame higher than PP imaging results；(3) PS areas imagings are wider than PP areas imaging, and this is due to turn
It is smaller than P wave reflections angle to change the angle of reflection of S ripples, so wave detector can receive the conversion S from the broader scope of subsurface interface
Ripple information.(4) in model deep layer, PS ripples imaging energy ratio PP ripple imaging energy is weak, because incident in deep layer P ripples low-angle
In the case of caused conversion S wave energies it is weaker caused by.
Comparison diagram 3a and 3b and Fig. 3 c and 3d respectively, it can be found that：Relative to conventional relief surface elastic wave Gaussian beam
The imaging of offset method, either PP ripples or the imaging of PS ripples, imaging effect of the invention is preferable, eliminates due to conventional method
The amplitude error caused by relatively large distance between middle beam center position and receiving point be present, the deep layer shown in dashed rectangle especially in figure
Target exploration region shows that the imaging energy of context of methods is significantly stronger than conventional method, has more clearly been imaged well-developed disconnected
Layer construction, imaging results then embody superiority of the present invention compared to conventional method.
Claims (1)
1. a kind of imaging of seismic data processing method, it is characterised in that comprise the following steps：
(1) reference frequency ω is chosen according to the dominant frequency of multiwave multicomponent earthquake data and frequency bandwidthrAnd highest frequency ωh；
(2) elastodynamics Gaussian beam original width is calculated：
In formula：It is the geometrical mean of V-shaped ripple (compressional wave or shear wave) velocity field,
(3) V-shaped ripple initial ray parameter space is：
(4) according to earth's surface elevation, inclination angle and p-and s-wave velocity information, with reference to selected by above-mentioned (1), (2), (3) step and count
The parameter of calculation, using single track processing mode, elastodynamics Gaussian beam is directly emitted at focus and receiving point respectively, and pass through
Solve kinematics ray tracing equation 3 and kinetics ray-tracing equation 4 and then elastodynamics Gaussian beam is calculated
Travelling when and amplitude information, then try to achieve the backward extension vector wave field decoupled under complex near surface conditionss, as shown in Equation 5；
In formula, wherein, W1 v(xr)、It is the weights relevant with p-and s-wave velocity and polarization vector,To receive
Point xrDisplacement caused by place at x, ρ (xr) it is xrLocate Media density, S is the free earth's surface of rising and falling, u1(xr, ω) and u2(xr,
ω) it is respectively horizontal and vertical component seismic record；
(5) the PP ripples under the complex near surface conditionss derived according to cross-correlation image-forming condition based on single track processing and PS ripples are imaged public
Formula, using elastodynamics Gaussian beam travelling when, amplitude and polarity information calculate corresponding single track into picture value；
(6) (4), the computing flow of (5) step are repeated, calculate single-shot into picture value, most all big guns are folded into picture value at last
Add to obtain required imaging section, the PP ripples and PS ripple single-shot migration imaging formula such as formula that double complex conditions are handled based on single track
Shown in 6a-6b：
In formula,It is to correspond to V-shaped ripple in xrThe initial slowness vector at place, φ are P ripple incidence angles
Degree, sgn (φ) is sign function.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410587602.1A CN104749623B (en) | 2014-10-20 | 2014-10-20 | A kind of imaging of seismic data processing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410587602.1A CN104749623B (en) | 2014-10-20 | 2014-10-20 | A kind of imaging of seismic data processing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104749623A CN104749623A (en) | 2015-07-01 |
| CN104749623B true CN104749623B (en) | 2017-12-08 |
Family
ID=53589593
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410587602.1A Active CN104749623B (en) | 2014-10-20 | 2014-10-20 | A kind of imaging of seismic data processing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104749623B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106842308A (en) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | A kind of non-horizontal earth's surface direct continuation beam offset method in seismic prospecting |
| CN106597540B (en) | 2016-12-30 | 2017-10-31 | 中国科学院地质与地球物理研究所 | Gaussian beam offset imaging method and device |
| CN106772585B (en) * | 2017-01-26 | 2018-11-09 | 中国科学院地质与地球物理研究所 | A kind of quasi- analysis method and device of the optimization decoupling equation based on elastic wave |
| CN109655882A (en) * | 2017-10-10 | 2019-04-19 | 中国石油化工股份有限公司 | Seismic forward modeling method and apparatus based on Gaussian beam wave-field simulation |
| CN112068197B (en) * | 2020-08-24 | 2022-03-18 | 西南科技大学 | Fracture cracking scale description method based on seismic wave dynamic parameters |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1797038A (en) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | Method for shifting depth before superposition in seismic data process of undulating the earth's surface |
| CN102103216A (en) * | 2009-12-16 | 2011-06-22 | 中国石油天然气集团公司 | Prestack migration method of two-dimensional Gaussian ray bundle |
| CN102183786A (en) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | Fidelity amplitude gaussian beam pre-stack depth migration method under double complex conditions |
| CN102890290A (en) * | 2012-09-25 | 2013-01-23 | 中国石油天然气股份有限公司 | Prestack depth migration method under condition of undulating surface |
-
2014
- 2014-10-20 CN CN201410587602.1A patent/CN104749623B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1797038A (en) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | Method for shifting depth before superposition in seismic data process of undulating the earth's surface |
| CN102103216A (en) * | 2009-12-16 | 2011-06-22 | 中国石油天然气集团公司 | Prestack migration method of two-dimensional Gaussian ray bundle |
| CN102183786A (en) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | Fidelity amplitude gaussian beam pre-stack depth migration method under double complex conditions |
| CN102890290A (en) * | 2012-09-25 | 2013-01-23 | 中国石油天然气股份有限公司 | Prestack depth migration method under condition of undulating surface |
Non-Patent Citations (3)
| Title |
|---|
| 各向异性介质共炮域高斯束叠前深度偏移;段新意 等;《石油物探》;20140930;第53卷(第5期);579-586 * |
| 复杂地表条件下保幅高斯束偏移;岳玉波 等;《地球物理学报》;20120430;第55卷(第4期);1376-1383 * |
| 格林函数高斯束逆时偏移;黄建平 等;《石油地球物理勘探》;20140228;第49卷(第1期);101-106 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104749623A (en) | 2015-07-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104749623B (en) | A kind of imaging of seismic data processing method | |
| CN102053261B (en) | Method for processing seismic data | |
| CN102841379B (en) | Method for analyzing pre-stack time migration and speed based on common scatter point channel set | |
| CN111158049B (en) | Seismic reverse time migration imaging method based on scattering integration method | |
| CN108196305B (en) | Mountain land static correction method | |
| CN108363101A (en) | A kind of inclined shaft crosshole seismic Gaussian beam pre-stack depth migration imaging method | |
| CN102692648A (en) | Refraction wave residual static correction method based on cannon first arrival | |
| CN104459794A (en) | Common-reflection-point gather time-variable time difference correction method and device | |
| CN104483705A (en) | Three-dimensional residual static correction method | |
| CN102854526B (en) | Multi-component seismic data processing method | |
| CN107561589A (en) | A kind of near surface S-wave interval velocity method for establishing model | |
| CN104533396A (en) | Remote exploration sound wave processing method | |
| CN105093319A (en) | Ground micro-seismic static correction method based on three-dimensional seismic data | |
| CN102053260B (en) | Method for acquiring azimuth velocity of primary wave and method for processing earthquake data | |
| CN104808242B (en) | A kind of imaging of seismic data processing method | |
| CN105866833A (en) | VSP-CDP stacking method and three-dimensional VSP coverage number calculation method | |
| CN104199103A (en) | Static correction method and static correction device for converted wave | |
| CN107656308B (en) | A kind of common scattering point pre-stack time migration imaging method based on time depth scanning | |
| CN103454681A (en) | Method and equipment for evaluating imaging effect of three-dimensional earthquake observing system | |
| CN105510975A (en) | Method and device for improving signal-to-noise ratio of seismic data | |
| CN106338766A (en) | Pre-stack time migration method based on split-step Fourier algorithm | |
| CN107356972A (en) | A kind of imaging method of anisotropic medium | |
| CN106125139A (en) | A kind of D seismic modeling method and system | |
| CN104155690B (en) | The 3D seismic data stack velocity acquiring method deployed based on ellipsoid | |
| CN115144899A (en) | Rugged seabed OBN elastic wave combined deflection imaging method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| DD01 | Delivery of document by public notice |
Addressee: Yuan Maolin Document name: Notification of Passing Preliminary Examination of the Application for Invention |
|
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Qian Fanglin Inventor before: Yuan Maolin Inventor before: Huang Jianping Inventor before: Zhao Shengtian Inventor before: Li Jiguang Inventor before: Yang Jidong Inventor before: Li Zhenchun Inventor before: Duan Xinbiao Inventor before: Guo Shujuan |
|
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20171110 Address after: Hunan Zhongzhou prosperous village 414104 County of Yueyang province Xiang Jiang Jia Village Group No. 10 Applicant after: Qian Fanglin Address before: 266000 Qingdao Changjiang Road, Huangdao District, Shandong, No. 66 Applicant before: Yuan Maolin |
|
| TA01 | Transfer of patent application right |
Effective date of registration: 20171116 Address after: 400056, Chongqing District, Banan City, South Village Red Star Village 12 clubs Applicant after: CHONGQING SHUOAO TECHNOLOGY CO.,LTD. Address before: Hunan Zhongzhou prosperous village 414104 County of Yueyang province Xiang Jiang Jia Village Group No. 10 Applicant before: Qian Fanglin |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180112 Address after: Chongqing District of Banan city bridge 401344 Solitaire Town Street No. 11 6-8 Patentee after: Liu Mengxiang Address before: 400056, Chongqing District, Banan City, South Village Red Star Village 12 clubs Patentee before: CHONGQING SHUOAO TECHNOLOGY CO.,LTD. |
|
| CP02 | Change in the address of a patent holder | ||
| CP02 | Change in the address of a patent holder |
Address after: 274300 Industrial Park, zhangwuzhuang village, Caozhuang Township, Shan county, Heze City, Shandong Province Patentee after: Liu Mengxiang Address before: Chongqing District of Banan city bridge 401344 Solitaire Town Street No. 11 6-8 Patentee before: Liu Mengxiang |
|
| CP02 | Change in the address of a patent holder | ||
| CP02 | Change in the address of a patent holder |
Address after: 401320 6-8, No.11, Jinqiao street, Jielong Town, Banan District, Chongqing Patentee after: Liu Mengxiang Address before: 274300 Industrial Park, zhangwuzhuang village, Caozhuang Township, Shan county, Heze City, Shandong Province Patentee before: Liu Mengxiang |
|
| TR01 | Transfer of patent right |
Effective date of registration: 20240104 Address after: No. 621, 5th Floor, Building C, Lanchi Building, Yaodian Street Office, Qinhan New City, Xixian New District, Xi'an City, Shaanxi Province, 710075 Patentee after: Shaanxi Yibufeng Technology Co.,Ltd. Address before: 401320 6-8, No. 11, Jinqiao street, Jielong Town, Banan District, Chongqing Patentee before: Liu Mengxiang |
|
| TR01 | Transfer of patent right |