FR2872296A1 - Procede destine a amiliorer la resolution sismique - Google Patents
Procede destine a amiliorer la resolution sismique Download PDFInfo
- Publication number
- FR2872296A1 FR2872296A1 FR0506155A FR0506155A FR2872296A1 FR 2872296 A1 FR2872296 A1 FR 2872296A1 FR 0506155 A FR0506155 A FR 0506155A FR 0506155 A FR0506155 A FR 0506155A FR 2872296 A1 FR2872296 A1 FR 2872296A1
- Authority
- FR
- France
- Prior art keywords
- seismic
- recording
- close
- wave
- vertical
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000008569 process Effects 0.000 title abstract description 6
- 239000002689 soil Substances 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 4
- 230000003068 static effect Effects 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 3
- 230000002238 attenuated effect Effects 0.000 claims description 2
- 238000011161 development Methods 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 2
- 238000011156 evaluation Methods 0.000 abstract 2
- 238000000926 separation method Methods 0.000 abstract 1
- 239000006185 dispersion Substances 0.000 description 6
- 230000005284 excitation Effects 0.000 description 5
- 238000005070 sampling Methods 0.000 description 5
- 238000001514 detection method Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 239000010409 thin film Substances 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005474 detonation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 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. analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/003—Seismic data acquisition in general, e.g. survey design
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
Abstract
Le procédé comprend essentiellement (1) l'obtention de l'enregistrement sismique vertical xi(t) proche de la surface ; (2) la séparation de l'onde virtuelle réfléchie de x1(t) par traitement de la section sismique verticale ; (3) l'évaluation d'une équation de type à convolution et de l'opérateur de déconvolution proche de la surface ; (4) l'évaluation de l'enregistrement sismique avec une haute résolution grâce à l'équation de type à déconvolution ; (5) l'interprétation géologique de l'enregistrement sismique avec la haute résolution.Procédé destiné à améliorer la résolution sismique dans le domaine de l'exploration pétrolière et du développement pétrolier.
Description
PROCEDE DESTINE A AMELIORER LA RESOLUTION SISMIQUE
L'invention concerne un procédé dans le domaine de l'exploration pétrolière et le développement pétrolier, en particulier un procédé destiné à améliorer la résolution sismique dans le domaine de l'exploration sismique.
L'enregistrement sismique est le résultat enregistré de l'ondelette sismique, qui est produite par l'hypocentre, et le résultat est enregistré après le processus au cours duquel l'ondelette sismique se propage vers le sous-sol et est renvoyée vers le sol lorsqu'elle rencontre la limite de réflexion. La résolution sismique verticale dépend de la fréquence dominante de l'ondelette sismique, qui est affectée par la condition d'excitation et la condition de propagation. Lorsque l'hypocentre de la dynamite sismique est généré, la dispersion d'onde de l'ondelette sismique est améliorée par la quantité de dynamite ou affaiblie par la strate au niveau du point de détonation, donc, la fréquence dominante de l'ondelette sismique sera abaissée. L'atténuation, telle que l'absorption, la diffusion, etc... a lieu au cours du processus de propagation de l'ondelette sismique. L'ondelette sismique continue à s'améliorer avec la distorsion. En particulier, lorsque l'ondelette sismique se propage proche de la surface, Patté- nuation et la dispersion d'onde de l'ondelette sismique sont relativement importantes. La quantité de l'atténuation de chaque mètre proche de la surface est de relativement 500 à 1 000 fois plus élevée que celle de la couche à grande vitesse. L'atténuation de l'onde sismique est différente en raison de la différence d'altitude proche de la surface et de la lithologie de la strate. Il est utile d'améliorer la résolution sismique verticale en utilisant la technologie de traitement, telle que l'acquisition du champ de la présente résolution élevée et de la correction statique, des spectres blancs, une déconvolution d'ondelettes, et un filtrage anti-Q, et analogues. Cependant, la fréquence dominante apparente ne peut pas être suffisante pour les besoins de reconnaissance de la couche mince dans la présente exploration sismique, étant donné qu'elle peut uniquement être améliorée jusqu'aux environs de 60 Hz.
Pour surmonter le défaut de l'art antérieur selon lequel la fréquence dominante apparente de l'onde sismique ne peut pas être suffisante pour les besoins de la reconnaissance de la couche mince dans la présente exploration sismique, étant donné qu'elle peut uniquement être améliorée jusqu'aux environ de 60 Hz, un procédé de la présente invention est fourni pour améliorer la résolution sismique. Ce procédé d'amélioration de la résolution sismique peut accroître la fréquence dominante apparente de l'onde sismique jusqu'à 200 Hz et plus, améliorer la résolution sismique 12 13 revets",23900".23992 doc IO juin 2005 - LS verticale, et résoudre le problème selon lequel il est difficile de reconnaître la couche mince dans la présente exploration sismique.
Le point de départ de la conception conformément à la présente invention consiste à récupérer l'onde sismique reçue par les données sismiques du sol de l'onde sismique produite par la petite dose de dynamite dans la couche à grande vitesse et l'onde sismique reçue dans la couche à grande vitesse qui est sous la couche à faible vitesse, compenser l'influence de l'onde sismique provoquée par l'atténuation et la dispersion d'ondes de l'onde sismique proche de la surface, donc, envisager le sujet de l'amélioration de la résolution sismique verticale.
Selon l'invention, ce procédé comprend essentiellement les étapes consistant à : (1) obtenir l'enregistrement sismique vertical xi(t) proche de la surface en utilisant le procédé de micro-diagraphie des sondages, s'assurer de la profondeur et de la vitesse de la couche à faible vitesse et de la couche à grande vitesse et la caractéristique du spectre de fréquence de l'onde sismique proche de la surface du sol, où xo(t) est la courbe de l'enregistrement sismique vertical proche de la surface du sol lorsque i=0, et indique l'enregistrement sismique des caractéristiques physiques de la surface du sol; et xn(t) est la courbe de l'enregistrement sismique vertical proche de la surface du sol lorsque i=n, n est une certaine profondeur, et indique l'enregistre-ment sismique des caractéristiques physiques de la couche à grande vitesse à une certaine profondeur de niveau de base sous la couche à faible vitesse; (2) séparer l'onde réfléchie virtuelle de xi(t) de la couche à grande vitesse en utilisant un procédé de section sismique vertical, améliorer le rapport signal bruit (SNR) de l'enregistrement sismique de l'onde de trajet direct xn(t) au niveau de la couche à grande vitesse en utilisant un procédé d'empilement de paliers de la section sismique verticale, et fournir la correction statique à l'enregistrement sismique xo(t) ; (3) établir une équation de type à convolution, qui est utilisée lorsque l'onde sismique s'étend proche de la surface, c'est-à-dire xo(t)=xn(t) *s(t), ou Xo(f)= Xn(f) S(f), évaluer l'opérateur de déconvolution proche de la surface: s-'(t)=xn(t)*xo-'(t), ou S-'(f)=Xn(f)*Xo-I(f)' où Xo(f) est la transformée de Fourier de xo(t), Xn(f) est la transformée de Fourier de xn(t), s(t) est l'opérateur de convolution proche de la surface, S(f) est la transformée de Fourier de s(t) ; (4) établir une équation de type à déconvolution de l'enregistrement sismique autour du sol de la micro-diagraphie en utilisant l'opérateur de déconvolution proche de la surface et l'enregistrement sismique du sol, évaluer la haute résolu- tion de l'enregistrement sismique, c'est-à-dire, w(t)=g(t)*s-'(t), ou W(f) -G(f) S-1(f), où g(t) est l'enregistrement sismique du sol, G(f) est la transformée de Fourier de g(t), w(t) est le résultat, qui est compensé lorsque l'onde sismique de la déconvolulion proche de la surface est atténuée, et W(f) est la transformée de Fourier de w(t) ; R Rrecets123900A23992 doc - 10 juin 2005 - 2 (5) fournir l'interprétation géologique et par filtrage passe-bande à l'enre- gistrement sismique avec une haute résolution évaluée.
L'enregistrement sismique vertical proche de la surface obtenue par microdiagraphie des sondages à l'étape (1) ci-dessus peut être obtenu en utilisant la micro- diagraphie à double sondage, et peut également être obtenu en utilisant la micro-diagraphie à sondage unique, où l'enregistrement sismique xo(t), qui répond à la caractéristique physique du sol, peut être remplacé par la première onde de rupture dans l'exploration sismique classique.
De plus, l'enregistrement sismique vertical proche de la surface à l'étape (1) ci- dessus, qui est obtenu grâce à la micro-diagraphie des sondages, peut adopter la même quantité de dynamite au niveau de différents points d'excitation, et peut également adopter la quantité différente de dynamite au niveau du même point d'excitation.
L'onde de trajet directe xr,(t) à l'étape (2) ci-dessus est l'enregistrement d'une certaine profondeur du rapport signal-bruit élevé dans la couche à grande vitesse, où le rapport signal-bruit de l'onde de trajet direct xn(t), qui est amélioré en utilisant le procédé d'empilement de paliers de la section sismique verticale, est celui des valeurs de... xi_2(t), )4_i(t), xn+l(t), xn+2(t)... qui sont corrigées de façon statique par rapport au même temps de l'enregistrement sismique de xn(t) en utilisant la vitesse de la couche à grande vitesse, puis les valeurs de... xn_2(t), xn_i(t), xn+i(t), xn+2(t)... sont empilées pour améliorer le rapport signal-bruit de l'onde de trajet direct xn(t). Dans cette étape, la correction statique pour l'enregistrement sismique dans un puits consiste à enlever la différence du premier temps de rupture entre xo(t) et xn(t).
La circonférence de la micro-diagraphie des sondages dans l'étape (4) cidessus est une plage de plusieurs centaines de mètre à plusieurs dizaines de kilomètres. Si le changement de l'altitude du sol ou de la lithologie est important, cette plage est faible; si le changement de l'aptitude ou de la lithologie est faible, cette plage peut être très importante. L'opérateur de déconvolution proche de la surface entre les microdiagraphies des sondages peut être obtenu par interpolation.
De plus, l'enregistrement sismique du sol g(t) à l'étape (4) peut être l'enregistrement sismique de l'exploration sismique soit à terre soit en mer. Ledit enregistre-ment sismique du sol est représenté comme l'enregistrement sismique de la sonde de démodulation des données sismiques du sol, ou le résultat qui était corrigé de façon statique par l'enregistrement sismique de la sonde de démodulation. Il peut s'agir de l'enregistrement de la section sismique soit avant que le niveau de l'enregistrement sismique du sol ne soit empilé ou ne migre. Les procédés de traitement habituels, tels qu'une correction de l'obliquité, par empilement et une migration, peuvent avoir lieu après que le procédé destiné à améliorer la résolution sismique a été appliqué.
RVBrevets\23900,23992 doc 10 jmn 2005 - 3 A L'expression de l'interprétation géologique à l'étape (5) ci-dessus est l'interprétation structurelle, la stratigraphie à séquence lithologique, la prévision de roche réservoir, la détection sismique de pétrole, et ainsi de suite, qui sont dans le domaine de développement de l'exploration pétrolière.
Cette invention a les effets utiles suivants: étant donné que l'immense atténuation et la dispersion d'ondes de l'onde sismique proche de la surface du sol peuvent être compensées en utilisant ce procédé, la durée de l'onde sismique peut être forte-ment raccourcie. La fréquence dominante apparente des données sismiques peut atteindre 200 Hz et plus, ce qui accroît fortement la résolution sismique verticale, et résout le problème selon lequel il est difficile de reconnaître la couche mince lors de la présente exploration sismique, et donne de l'importance à l'interprétation structurelle, la stratigraphie à séquence lithologique, la prévision de roche réservoir, la détection sismique de pétrole et ainsi de suite, qui sont dans le domaine du développement de l'exploration pétrolière.
D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description détaillée qui suit, donnée à titre d'exemple uniquement et en référence aux dessins.
La figure 1 est le diagramme de l'enregistrement sismique vertical obtenu conformément à la présente invention.
La figure 2 est le diagramme de la section sismique migrée obtenue conformé- ment à la présente invention.
La figure 3 est le diagramme de la section sismique migrée obtenue conformé-ment à l'art antérieur.
Selon la description ci-dessus sur le contexte de la présente invention et l'expli- cation supplémentaire ci-après de l'exemple, l'homme du métier peut envisager l'objet de la présente invention. Cependant, cet exemple n'est qu'une partie de la présente invention, et ne peut pas limiter la plage de protection de la présente invention.
L'objet de l'exemple suivant consiste à améliorer la résolution sismique par ladite résolution sismique verticale de la section sismique, qui est obtenue en adoptant le principe de la présente invention, afin de montrer des améliorations réelles bien supérieures à la section sismique habituelle.
Selon l'étape (1) dudit procédé de la présente invention, l'enregistrement sismique vertical xi(t) est obtenu proche de la surface dans la zone S en utilisant la micro- diagraphie à double diagraphie, qui est représentée sur la figure 1. La distance entre les deux diagraphies est de 4 m, la profondeur du point d'excitation est entre 0,5 m et 34 ni. La distance entre les points d'excitation est de 1 m, et le nombre total de perforateurs est de 33. La ressource excitée est composée d'un détonateur, la profondeur R 'Rreccts 23900'22492 doc - 10juin 2005 - 4'8 de la diagraphie reçue est de 34 m, et la distance d'adoption du temps enregistré est de 0,25 ms. L'enregistrement excité à 0,5 m est appelé xo(t), l'enregistrement à 10,5 m est appelé x10(t). A partir des résultats mesurés dans cette diagraphie, la profondeur de la couche à faible vitesse est entre 0,5 et 6 m, la vitesse au niveau de la couche à faible vitesse est de 360 m/s, la profondeur de la couche à grande vitesse est supérieure à la plage de 6 à 33 m, et la vitesse de la couche à grande vitesse est de 1 700 m/s, tout ceci étant représenté sur la figure 1. La plage d'énergie dominante de xo(t) est entre 10 et 200 Hz, et la plage d'énergie dominante de xlo(t) est entre 30 et 1 400 Hz, lesquelles sont obtenues à partir de l'analyse spectrale du résultat mesuré.
Selon l'étape (2) dudit procédé de la présente invention, l'onde virtuelle réflé- chie est séparée de xl(t) en utilisant le procédé de la section sismique verticale.
Selon l'étape (2) dudit procédé de la présente invention, x8(t), x9(t), x11(t), x12(1) peuvent être corrigés par rapport au même temps que xlo(t) en utilisant le procédé d'empilement de paliers de la section sismique verticale. Le rapport signal-bruit de l'onde de trajet xlo(t) au niveau de la profondeur de 10 m de la diagraphie peut être amélioré par cinq empilements. xo(t)peut être corrigé par rapport au même temps que xlo(t), xo(t)=xo(t-20).
Selon l'étape (3) dudit procédé de la présente invention, une équation de type à convolution: Xo(f)= X10(f) S(f) est établi, laquelle équation est utilisée lorsque l'onde sismique de la micro-diagraphie des sondages se propage proche de la surface. L'opérateur de déconvolution proche de la surface: S-1(f)= X10(f) Xo-1(f) est évalué (rééchantillonnage de xo(t) et de x10(t) avant la transformée de Fourier de X10(f) et de X10(f), l'intervalle d'échantillonnage est de 1 ms, et les points d'échantillonnage totaux sont de 70).
Selon l'étape (4) dudit procédé de la présente invention, une équation de type à déconvolution de l'enregistrement de la section migrée est établie, laquelle équation est utilisée pour le profil sismique sur la ligne SX autour du sol de la micro-diagraphie des sondages, en utilisant l'opérateur de déconvolution proche de la surface et l'enregistrement sismique du sol. La haute résolution de l'enregistrement sismique de l'enregistrement de la section migrée pour le profil à réfraction sismique sur la ligne SX et évaluée, W(f)=G(f) S-1(f). G(f) est la transformée de Fourier du profil à réfraction sismique sur la ligne SX g(f), et l'intervalle d'échantillonnage du temps de g(f) est de 1 ms. Le résultat w(t), qui a été compensé par l'atténuation et la dispersion d'onde de l'onde sismique de déconvolution proche de la surface, est obtenu par la transformée inverse de Fourier. Le résultat, qui reflète différentes propriétés géologi- ques, est obtenu par le filtrage passe-bande.
L'énorme atténuation et la dispersion d'onde de l'onde sismique proche de la surface du sol peuvent être compensées après que ce procédé a été appliqué. La durée R.Brevets.23900A23992. doc - IO juin 2005 - 5"0 de l'ondelette sismique est fortement raccourcie. La fréquence dominante apparente de la section sismique SX classique g(t) peut atteindre 240 Hz, ce qui est représenté sur la figure 2. Cependant, la fréquence dominante apparente de la section sismique SX classique est uniquement de 55 Hz, ce qui est représenté sur la figure 3. Les données représentées sur la figure 2 s'adaptent bien mieux aux données antérieures de l'exploration. Donc, ce procédé peut améliorer fortement la résolution, et résoudre le problème selon lequel il est difficile de reconnaître la couche mince grâce à la présente exploration sismique. Lorsque ce procédé est mis en oeuvre, les caractéristiques structurelles sur la section, telles qu'une phase, un pli, etc., sont évidentes, les phénomènes de la stratigraphie à séquence lithologique, tels qu'un phénomène supérieur, de progradation, d'évent, etc., sont manifestes, et il a une importance sur la prévision de roche de réservoir sismique, la détection pétrolière sismique, etc...
Il est à noter que le présent intervalle d'échantillonnage des données sismiques superficielles est uniquement de 1 ms. Selon le théorème de Nyquist, la plus grande largeur de fréquence, qui peut être enregistrée est 500 Hz. L'effet de la section sismique, dans laquelle la fréquence dominante obtenue dans le présent exemple peut atteindre 240 Hz, est important. S'il est nécessaire que la résolution sismique soit davantage améliorée, il faudra améliorer l'intervalle d'échantillonnage des données sismiques du sol jusqu'à 0,5 ms, 0,25 ms.
R.A6rcvets\23900_23992 doc - 10 juin 2005 - 6B
Claims (1)
1. Procédé destiné à améliorer la résolution sismique, caractérisé en ce qu'il comprend les étapes consistant à : (1) obtenir l'enregistrement sismique vertical xi(t) proche de la surface du sol en utilisant le procédé de micro-diagraphie des sondages, fournissant la profondeur et de la vitesse de la couche à faible vitesse et de la couche à grande vitesse et la caractéristique du spectre de fréquence de l'onde sismique proche de la surface du sol, où xo(t) est la courbe de l'enregistrement sismique vertical proche de la surface du sol lorsque i=0, et indique l'enregistrement sismique des caractéristiques physiques de la surface du sol; et xn(t) est la courbe de l'enregistrement sismique vertical proche de la surface du sol lorsque i=n, n est une certaine profondeur, et indique l'enregistrement sismique des caractéristiques physiques de la couche à grande vitesse à une certaine profondeur de niveau de base sous la couche à faible vitesse; (2) séparer l'onde réfléchie virtuelle de xi(t) de la couche à grande vitesse en utilisant un procédé de section sismique vertical, améliorer le rapport signal bruit (SNR) de l'enregistrement sismique de l'onde de trajet direct Xd(t) au niveau de la couche à grande vitesse en utilisant un procédé d'empilement de paliers de la section sismique verticale, et fournir la correction statique à l'enregistrement sismique xo(t) ; (3) établir une équation de type à convolution, qui est utilisée lorsque l'onde sismique s'étend proche de la surface, c'est-à-dire xo(t)=xn(t) *s(t), ou Xo(f)= Xn(f) S(f), évaluer l'opérateur de déconvolution proche de la surface: s-'(t)=xn(t)*xo-l(t), ou S-l(f)=Xn(f)*Xo-'(f)' où Xo(f) est la transformée de Fourier de xo(t), Xn(f) est la transformée de Fourier de xn(t), s(t) est l'opérateur de convolution 25 proche de la surface, S(f) est la transformée de Fourier de s(t) ; (4) établir une équation de type à déconvolution de l'enregistrement sismique autour du sol de la micro-diagraphie en utilisant l'opérateur de déconvolution proche de la surface et l'enregistrement sismique du sol, évaluer la haute résolution de l'enregistrement sismique, c'est-à-dire, w(t)=g(t)*s-'(t), ou W(f)-G(f) S-'(f), où g(t) est l'enregistrement sismique du sol, G(f) est la transformée de Fourier de g(t), w(t) est le résultat, qui est compensé lorsque l'onde sismique de la déconvolution proche de la surface est atténuée, et W(f) est la transformée de Fourier de w(t) ; (5) fournir l'interprétation géologique et par filtrage passe-bande à l'enregistrement sismique avec une haute résolution évaluée.
R reocr, 23900"_3992 doc - 10 juin 2005 - 7/8
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2004100499961A CN100383557C (zh) | 2004-06-25 | 2004-06-25 | 一种提高地震分辨率的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2872296A1 true FR2872296A1 (fr) | 2005-12-30 |
| FR2872296B1 FR2872296B1 (fr) | 2016-11-04 |
Family
ID=34665883
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR0506155A Expired - Fee Related FR2872296B1 (fr) | 2004-06-25 | 2005-06-17 | Procede destine a ameliorer la resolution sismique |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US7414919B2 (fr) |
| CN (1) | CN100383557C (fr) |
| FR (1) | FR2872296B1 (fr) |
Families Citing this family (48)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100487488C (zh) * | 2006-05-11 | 2009-05-13 | 中国石油集团东方地球物理勘探有限责任公司 | 零井源距垂直地震剖面纵横波数据深度域走廊叠加剖面处理方法 |
| US7453765B2 (en) * | 2006-05-16 | 2008-11-18 | Ikelle Luc T | Scattering diagrams in seismic imaging |
| CN101201407B (zh) * | 2006-12-12 | 2010-05-19 | 中国石油天然气集团公司 | 相对无高频泄漏等效n点平滑谱模拟反褶积方法 |
| US7672193B2 (en) * | 2007-05-21 | 2010-03-02 | Schlumberger Technology Corporation | Methods and systems for processing acoustic waveform data |
| CN101393270B (zh) * | 2007-09-21 | 2011-04-20 | 中国石油集团东方地球物理勘探有限责任公司 | 一种连续介质地区的表层高精度静校正方法 |
| CN102109612B (zh) * | 2009-12-23 | 2012-10-17 | 中国石油天然气集团公司 | 一种地震波吸收衰减补偿方法 |
| CN101915091B (zh) * | 2010-06-22 | 2013-03-27 | 中国石油天然气股份有限公司 | 一种测井频率匹配薄层校正方法与设备 |
| CN102109617A (zh) * | 2010-12-15 | 2011-06-29 | 大庆油田有限责任公司 | 用双井微测井技术测量近地表地层q值的方法 |
| CN102890288B (zh) * | 2011-07-20 | 2015-05-13 | 中国石油天然气集团公司 | 一种地震波层速度反演方法 |
| CN102393533A (zh) * | 2011-07-29 | 2012-03-28 | 中国石油化工股份有限公司 | 一种利用钻井轨迹校正vsp上行转换波的处理方法 |
| CN102411155A (zh) * | 2011-07-29 | 2012-04-11 | 中国石油化工股份有限公司 | 一种利用钻井轨迹校正vsp下行透射转换波的处理方法 |
| CN102914795B (zh) * | 2011-08-05 | 2015-02-25 | 中国石油天然气集团公司 | 一种炮检对照明度确定方法 |
| CN103076631B (zh) * | 2011-10-26 | 2016-03-23 | 中国石油化工股份有限公司 | 一种基于零脉冲反褶积提频技术的煤层气田陷落柱预测方法 |
| CN103513278B (zh) * | 2012-06-19 | 2016-02-10 | 中国石油化工股份有限公司 | 利用地震波组厚度进行储层预测的方法 |
| CN103713317B (zh) * | 2012-09-29 | 2016-09-14 | 北京诺克斯达石油科技有限公司 | 基于对时变子波进行分频段处理的反褶积方法 |
| CN103728659B (zh) * | 2012-10-12 | 2016-07-06 | 中国石油化工股份有限公司 | 一种提高地下岩溶探测精度的方法 |
| CN104142513B (zh) * | 2013-10-10 | 2017-01-11 | 中国石油化工股份有限公司 | 一种提高地震勘探中层析反演精度的方法 |
| CN104977618B (zh) * | 2014-04-09 | 2017-09-22 | 中国石油集团东方地球物理勘探有限责任公司 | 一种评价页岩气储层及寻找甜点区的方法 |
| CN105093266A (zh) * | 2014-05-14 | 2015-11-25 | 中国石油化工股份有限公司 | 地震地质综合描述地层超覆线方法 |
| CN105301658B (zh) * | 2014-06-26 | 2018-05-22 | 中石化石油工程地球物理有限公司胜利分公司 | 利用微测井记录提取近地表反褶积算子的方法 |
| CN105572732A (zh) * | 2014-10-08 | 2016-05-11 | 中国石油化工股份有限公司 | 一种属性变化率逐次晋级的裂缝发育带检测方法 |
| CN105629297A (zh) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | 预测复杂断块油田微断层分布规律的方法 |
| CN104614767A (zh) * | 2014-12-11 | 2015-05-13 | 中国石油大学(华东) | 基于分段延拓的时变地震子波相位校正方法 |
| CN105759310B (zh) * | 2014-12-16 | 2018-05-08 | 中国石油化工股份有限公司 | 复杂非均质储层介质中地震波衰减及速度频散预测方法 |
| CN104597497A (zh) * | 2015-02-26 | 2015-05-06 | 浪潮电子信息产业股份有限公司 | 一种基于叠前瞬时频率属性分析的储层烃类预测方法 |
| CN104965225B (zh) * | 2015-06-09 | 2017-07-28 | 中国石油化工股份有限公司 | 一种基于子波旁瓣信息的油藏识别方法 |
| CN104914466B (zh) * | 2015-06-26 | 2017-12-12 | 中国石油大学(华东) | 一种提高地震资料分辨率的方法 |
| CN106526664A (zh) * | 2015-09-09 | 2017-03-22 | 娄兵 | 一种时空域几何变换水平叠加地震成像方法 |
| CN105954803B (zh) * | 2016-07-08 | 2018-02-02 | 中国石油天然气股份有限公司 | 叠后地震反演方法及装置 |
| CN106443771B (zh) * | 2016-08-05 | 2018-12-11 | 中国石油集团东方地球物理勘探有限责任公司 | 提高转换波地震数据分辨率的方法及其速度反演方法 |
| CN107780400A (zh) * | 2016-08-24 | 2018-03-09 | 北京市水利规划设计研究院 | 一种大口径输水隧洞地基检测系统及方法 |
| CN106291715B (zh) * | 2016-09-24 | 2018-04-03 | 中国地质大学（北京） | 一种基于断裂自相似理论的低级序断层发育规律预测方法 |
| US10895654B2 (en) * | 2017-04-20 | 2021-01-19 | Exxonmobil Upstream Research Company | Method for generating optimized seismic target spectrum |
| RU2690068C1 (ru) * | 2018-05-14 | 2019-05-30 | Акционерное общество "Башнефтегеофизика" | Способ определения границ ВЧР методом прямого МСК в комплексе с методом преломленных волн |
| CN109239783B (zh) * | 2018-10-10 | 2020-02-07 | 东北石油大学 | 一种井震结合中取舍地震信息的方法 |
| CN111077577B (zh) * | 2018-10-22 | 2022-02-01 | 中国石油天然气股份有限公司 | 井地联合储层描述方法及装置 |
| CN111323813B (zh) * | 2018-12-13 | 2022-06-03 | 中国石油天然气集团有限公司 | 一种走廊叠加剖面生成方法及系统 |
| CN109991662A (zh) * | 2019-05-15 | 2019-07-09 | 中油奥博(成都)科技有限公司 | 浅部地层二维或三维弹性参数测量和计算的装置及方法 |
| CN110554435B (zh) * | 2019-07-22 | 2024-04-26 | 中国石油化工股份有限公司 | 一种利用微测井数据构建品质因子体的方法 |
| CN110532615B (zh) * | 2019-07-26 | 2022-11-15 | 自然资源部第二海洋研究所 | 一种浅海复杂地形的逐级分解方法 |
| CN112444845B (zh) * | 2019-08-30 | 2024-04-09 | 中国石油化工股份有限公司 | 一种提高地震记录质量的方法 |
| CN112649847A (zh) * | 2019-10-11 | 2021-04-13 | 中国石油化工股份有限公司 | 一种提高地震垂向分辨率的调谐预测反褶积方法及系统 |
| CN112946755B (zh) * | 2019-12-11 | 2024-04-30 | 中国石油天然气集团有限公司 | 近地表表层的调查方法和装置 |
| CN113534237A (zh) * | 2020-04-18 | 2021-10-22 | 中国石油化工股份有限公司 | 一种基于正演模型的潜山低序级断层分辨力定量分析方法 |
| CN111708081B (zh) * | 2020-05-29 | 2022-04-15 | 成都理工大学 | 考虑衰减频散的深度域地震记录合成方法 |
| CN114114424B (zh) * | 2020-08-28 | 2023-06-09 | 中国石油化工股份有限公司 | 一种关联监视记录的微测井解释及建立成果图的方法 |
| CN112363225A (zh) * | 2020-11-02 | 2021-02-12 | 中国石油天然气集团有限公司 | 横波微测井解释方法及装置 |
| CN114002736B (zh) * | 2021-09-07 | 2023-07-07 | 中国矿业大学 | 一种基于权重反褶积的地震勘探多频数据融合方法 |
Family Cites Families (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4596005A (en) * | 1983-04-20 | 1986-06-17 | Chevron Research Company | Method of seismic collection utilizing multicomponent processing receivers and processing resultant conventional and converted P- or S-wave data |
| US5191557A (en) * | 1986-12-30 | 1993-03-02 | Gas Research Institute | Signal processing to enable utilization of a rig reference sensor with a drill bit seismic source |
| US5197038A (en) * | 1991-06-21 | 1993-03-23 | Schlumberger Technology Corporation | Method and sonic tool apparatus for investigating properties of earth formations transversed by a borehole |
| USH1307H (en) * | 1991-12-11 | 1994-05-03 | Exxon Production Research Company | Method for continuity logging |
| CN1033665C (zh) * | 1991-12-26 | 1996-12-25 | 切夫里昂研究和技术公司 | 改进地质构造的地震分辨率的方法 |
| US5237538A (en) * | 1992-02-20 | 1993-08-17 | Mobil Oil Corporation | Method for removing coherent noise from an array of seismic traces |
| US5583825A (en) * | 1994-09-02 | 1996-12-10 | Exxon Production Research Company | Method for deriving reservoir lithology and fluid content from pre-stack inversion of seismic data |
| JP3542386B2 (ja) * | 1994-10-13 | 2004-07-14 | キヤノン株式会社 | 磁気ヘッドおよび光磁気記録装置 |
| US5764516A (en) * | 1995-12-29 | 1998-06-09 | Atlantic Richfield Company | Method and system for surface-consistent phase and time lag correction of seismic data |
| US6081482A (en) * | 1996-05-09 | 2000-06-27 | 3Dgeo Development, Inc. | Semi-recursive imaging under complex velocity structures |
| US6049759A (en) * | 1998-01-16 | 2000-04-11 | Bp Amoco Corporation | Method of prestack 3-D migration |
| US6131694A (en) * | 1998-09-02 | 2000-10-17 | Ahlliburton Energy Services, Inc. | Vertical seismic profiling in a drilling tool |
| WO2001029581A1 (fr) * | 1999-10-21 | 2001-04-26 | Core Laboratories, Inc. | Procede pour obtenir des fonctions de transfert lors du traitement et de l'exploration de signaux sismiques |
| CN1271420C (zh) * | 2001-01-25 | 2006-08-23 | 维斯特恩格科地震控股有限公司 | 海洋地震数据处理方法和地震测量方法 |
-
2004
- 2004-06-25 CN CNB2004100499961A patent/CN100383557C/zh not_active Expired - Fee Related
-
2005
- 2005-06-17 FR FR0506155A patent/FR2872296B1/fr not_active Expired - Fee Related
- 2005-06-20 US US11/156,787 patent/US7414919B2/en active Active
Non-Patent Citations (1)
| Title |
|---|
| No Search * |
Also Published As
| Publication number | Publication date |
|---|---|
| FR2872296B1 (fr) | 2016-11-04 |
| CN100383557C (zh) | 2008-04-23 |
| US7414919B2 (en) | 2008-08-19 |
| CN1595201A (zh) | 2005-03-16 |
| US20050286344A1 (en) | 2005-12-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| FR2872296A1 (fr) | Procede destine a amiliorer la resolution sismique | |
| CA2733882C (fr) | Procede de surveillance d'un site de stockage geologique de gaz par inversion stratigraphique de donnees sismiques | |
| FR2991463B1 (fr) | Systeme et procede de tomographie a haute definition et a haute resolution pour une utilisation pour la generation de modeles de vitesse et d'images de reflectivite | |
| EP1963886B1 (fr) | Methode d'evaluation quantitative des pressions de fluides et de detection des surpressions d'un milieu souterrain | |
| FR2773618A1 (fr) | Attenuation de reflexions multiples de donnees du fond marin a composantes multiples | |
| FR2955397A1 (fr) | Procede et dispositif d'acquisition de donnees sismiques marines | |
| FR3016970A1 (fr) | ||
| WO2013153311A1 (fr) | Procédé de détermination de paramètres géomécaniques d'un échantillon de roche | |
| CA2704141A1 (fr) | Methode pour imager une zone cible du sous-sol a partir de donnees de type walkaway | |
| FR2844059A1 (fr) | Procede de diagnostic de gisements de petrole et de gaz | |
| FR2845773A1 (fr) | Procede de determination d'un parametre auquel on s'interesse d'une partie d'un reservoir dans des formations terrestres | |
| WO1985005696A1 (fr) | Procede de diagraphie acoustique instantanee dans un puits de forage | |
| FR2708350A1 (fr) | Méthode d'analyse de signaux sismiques. | |
| WO2012038647A1 (fr) | Procede d'estimation de parametres elastiques par inversion de mesures sismiques 4d | |
| FR2888946A1 (fr) | Procede et appareil pour attenuer le bruit du vent dans des donnees sismiques | |
| FR2779532A1 (fr) | Procede d'attenuation des reverberations de colonne d'eau dans un signal sismique de capteurs doubles | |
| EP0246148A1 (fr) | Procédé de prospection sismique à très haute résolution en forages horizontaux | |
| EP0269482B1 (fr) | Procédé de détection des hétérogénéités et/ou de détermination de caractéristiques pétrophysiques de couches géologiques d'un milieu à explorer | |
| EP2315055B1 (fr) | Méthode pour interpréter des enregistrements sismiques répétitives | |
| CA2731474C (fr) | Procede pour caracteriser un panache de co2 dans un aquifere de stockage geologique | |
| WO2009153443A1 (fr) | Methode pour evaluer des pressions de fluides et detecter des surpressions dans un milieu souterrain | |
| EP3172596A1 (fr) | Procédé de filtrage adaptatif de reflexions sismiques multiples | |
| Vesnaver et al. | A workflow for processing mono‐channel Chirp and Boomer surveys | |
| FR2800472A1 (fr) | Procede de prospection sismique mettant en oeuvre un traitement sur les ondes converties | |
| FR3000224A1 (fr) | Method and device to deghost seismic data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PLFP | Fee payment |
Year of fee payment: 12 |
|
| ST | Notification of lapse |
Effective date: 20180228 |