BE1025828A1 - Procédé de calcul du temps de propagation sismique 2-D sur la base d'une construction de front d'onde à sources virtuelles - Google Patents
Procédé de calcul du temps de propagation sismique 2-D sur la base d'une construction de front d'onde à sources virtuelles Download PDFInfo
- Publication number
- BE1025828A1 BE1025828A1 BE20185525A BE201805525A BE1025828A1 BE 1025828 A1 BE1025828 A1 BE 1025828A1 BE 20185525 A BE20185525 A BE 20185525A BE 201805525 A BE201805525 A BE 201805525A BE 1025828 A1 BE1025828 A1 BE 1025828A1
- Authority
- BE
- Belgium
- Prior art keywords
- wavefront
- propagation time
- seismic
- quadrilateral
- grid
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000010276 construction Methods 0.000 title claims abstract description 26
- 238000004364 calculation method Methods 0.000 claims abstract description 27
- 230000004807 localization Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000012821 model calculation 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
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
- G01V1/305—Travel times
-
- 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
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G01V20/00—
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/12—Signal generation
- G01V2210/125—Virtual source
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/57—Trace interpolation or extrapolation, e.g. for virtual receiver; Anti-aliasing for missing receivers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/67—Wave propagation modeling
- G01V2210/671—Raytracing
Abstract
L’invention concerne un procédé de calcul de temps de propagation sismique 2-D basé sur une construction de front d’onde à sources virtuelles. Grâce à l’adoption du procédé à sources virtuelles, lors du calcul du temps de propagation des points de grille dans les quadrilatères de front d’onde de l’invention, la précision de calcul du temps de propagation des points de grille est améliorée et un procédé de calcul sismique basé sur la construction de front d’onde est réalisé.
Description
DOMAINE TECHNIQUE
L'invention concerne le domaine du calcul du temps de parcours sismique, en particulier un procédé de calcul du temps de propagation sismique 2-D basé sur la construction de front d'onde.
CONTEXTE
Le Journal of Jilin University (Earth Science Edition) a divulgué, dans le 2ème numéro de 2008, un article intitulé RayTracing of Wavefront Construction by Bicubic Convolution Interpolation publié par Han Fuxing et al. L'article présente un procédé de calcul de temps de propagation sismique amélioré basé sur la construction de front d'onde, un procédé d'interpolation à convolution bicubique a été appliqué au calcul du temps de propagation à points de grille de manière à améliorer la précision de calcul du temps de propagation sismique. Par ailleurs, l'analyse d'erreurs sur le procédé de construction du front d'onde basé sur un algorithme d'interpolation à convolution bicubique a été réalisée à travers un milieu homogène, et les résultats d'analyse ont obtenu de bons effets.
Le Chinese Journal of Computational Physics a divulgué, dans le numéro 2 de 2008, un article intitulé Application and Comparison of Different Interpolation Algorithms in Ray-Tracing of Wavefront Construction publié par Han Fuxing et al. Dans l'article, les effets d'application d'un procédé d'interpolation de région la plus proche, une procédé d'interpolation bilinéaire, un procédé d'interpolation par morceaux et un procédé d'interpolation par convolution bicubique dans le tracé de rayons de construction de front d'onde ont été comparés. Les résultats du calcul de modèles reflètent que le procédé d'interpolation à convolution bicubique peut donner des
BE2018/5525 trajectoires de rayons plus précises que les trois autres méthodes.
Le Progress in Geophysics a divulgué, dans le numéro 5 de 2009, un article intitulé Study on Gird Point Positioning and Attribute Evaluation with the Method of Wavefront Construction publié par Han Fuxing et al. L'article introduit l'utilisation d'un procédé par produit vectoriel pour évaluer la relation de position relative d'un point de grille rectangulaire et un quadrilatère de front d'onde non régulier, et a fourni un procédé d'interpolation correspondant pour le calcul de l'information d'attribut du point de grille selon la relation de position du point de grille et le quadrilatère de front d'onde non régulier. L procédé fourni a été vérifiée à l'aide d'un modèle homogène et d'exemples de calcul, et de bons résultats de calcul ont été obtenus.
Comme le montrent les exemples ci-dessus, les procédés conventionnelles de calcul du temps de propagation sismique 2D basés sur la construction de front d'onde peuvent améliorer la précision de calcul dans une certaine mesure, mais la règle de propagation des ondes sismiques n'a pas été prise en compte dans le procédé d'interpolation et la précision de calcul améliorée était également limitée.
RÉSUMÉ
L'invention vise à résoudre le problème technique consistant à proposer un procédé de calcul de temps de propagation sismique basée sur la construction de front d'onde de source virtuelle . Compte tenu de la règle de propagation des ondes sismiques dans un milieu, un procédé de calcul de la source virtuelle est utilisée pour remplacer une méthode d'interpolation de convolution bicubique classique dans le procédé de calcul du temps de Voyage de points de grille dans les quadrilatères de front d'onde, et la précision de calcul et la stabilité du
BE2018/5525 procédé de calcul du temps de propagation sismique basé sur la construction du front d'onde sont améliorées.
Afin de résoudre le problème technique, le schéma technique adopté réside dans le fait que un procédé de calcul de temps de propagation sismique de haute précision basé sur la construction de front d'onde comprend les étapes suivantes :
étape 1 : lire des fichiers de paramètres pertinents et un modèle de vitesse, les fichiers de paramètres comprenant le nombre de points de grille, l'espacement de grille, l'emplacement des sources, la longueur de pas de tracé, la plage d'angles des rayons et étape 2 : tracer des rayons, et insérer de nouveaux rayons pendant le traçage pour assurer la couverture des rayons, le but du traçage de rayons étant de résoudre un groupe d'équations de traçage de rayons cinématiques à l'ide d'un procédé de RungeKutta, comme indiqué dans la formule ci-dessous :
dx. -, ~r = v άτ dp. d ( 1 _! dv —- = — = -V dr dxt\v ) dx{ où x± représente la composante de localisation, p± représente la composante de lenteur, τ représente le temps de propagation sismique, et v représente la vitesse sismique ;
étape 3 : diviser un espace de calcul en une pluralité de quadrilatères de front d'onde par calcul de l'information de position spatiale de points discrets sur les rayons acquis ;
étape 4 : évaluer la relation de position du point de grille et du quadrilatère de front d'onde par un procédé à produit vectoriel, et déterminer le point de grille contenu dans le quadrilatère de front d'onde ;
étape 5 : calculer les emplacements de sources virtuelles correspondant aux points via les informations pertinentes des sommets du quadrilatère de front d'onde, et calculer le temps de
BE2018/5525 propagation sismique des points de grille contenus dans le quadrilatère de front d'onde en fonction des localisations des sources virtuelles ; et étape 6: terminer les calculs du temps de propagation de tous les points de la grille et produire les résultats du calcul du temps de propagation final.
Par rapport à l'art antérieur, le procédé de calcul du temps de propagation sismique divulgué par l'invention a pour effets bénéfiques que, compte tenu de la règle de propagation des ondes sismiques dans le milieu, le temps de propagation des points de grille dans les quadrilatères de front d'onde est calculé par un procédé à source virtuelle, la précision de calcul du temps de propagation du point de grille est améliorée et la précision de calcul du procédé à temps de propagation sismique basé sur la construction de front d'onde est amélioré.
BREVE DESCRIPTION DES DESSINS
La figure 1 est un organigramme d'un procédé de calcul de temps de propagation sismique basé sur la construction de front d'onde de sources virtuelles de l'invention.
La figure 2 est un diagramme schématique de la division d'une région de calcul.
La figure 3 est un diagramme schématique du calcul du temps de parcours sismique de sources virtuelles, A, B, C et D étant quatre sommets d'un quadrilatère de front d'onde, OA, OB, Oc et OD étant les emplacements des sources virtuelles correspondant à A, B, C et D, et R étant un point de grille dans le quadrilatère de front d'onde ABCD.
La figure 4 est l'erreur relative du temps de propagation dans un procédé de construction de front d'onde classique dns un milieu homogène.
La figure 5 est l'erreur relative du procédé de calcul de temps de propagation sismique basé sur la construction de front d'onde de sources virtuelles dans le milieu homogène.
BE2018/5525
DESCRIPTION DÉTAILLÉE
L'invention est décrite en détails ci-dessous en référence aux dessins joints et à des modes de réalisation spécifiques.
La figure 1 représente l'organigramme du procédé de calcul du temps de propagation sismique basé sur la construction de front d'onde de sources virtuelles. Le flux de réalisation du procédé de l'invention est représenté spécialement sur la figure comme suit :
1) lire des fichiers de paramètres pertinents et un modèle de vitesse, les fichiers de paramètres comprenant le nombre de points de grille, l'espacement de grille, l'emplacement des sources, la longueur de pas de traçage, la plage d'angles des rayons et l'intervalle d'échantillonnage ;
2) traçage des rayons, et insertion de nouveaux rayons pendant le traçage pour assurer la couverture des rayons, la gamme d'angles d'émission des rayons étant de -80° à + 80°, l'intervalle d'échantillonnage étant de 3° à 6°, la longueur de l'étape de traçage des rayons étant de 2 ms à 6 ms, le but du traçage des rayons étant de résoudre un groupe d'équations de traçage de rayons cinématiques à l'aide d'un procédé de RungeKutta , comme indiqué dans la formule ci-dessous :
dx. -, ~r = v άτ dp. d ( 1 _! dv —- = — = -V dr dxt\v ) dx{ où Xi représente la composante de localisation, p± représente la composante de lenteur, τ représente le temps de propagation sismique, et v représente la vitesse sismique ;
3) diviser un espace de calcul en une pluralité de quadrilatères de front d'onde avec les quatre points sur des rayons adjacents et des surfaces de front d'onde adjacentes comme sommets de chacun des quadrilatères de front d'onde, par
BE2018/5525 calcul des informations de position spatiale des points discrets sur les rayons acquis (voir figure 2) ;
4) d'une part, cribler grossièrement un point de grille qui peut être positionnée dans la zone de couverture de chacun des quadrilatères de front d'onde via les positions des sommets du quadrilatère de front d'onde, puis évaluer la relation de position du point de la grille et du quadrilatère de front d'onde par un procédé à produit vectoriel, et déterminer le point de grille contenu dans le quadrilatère de front d'onde ;
5) calculer les emplacements de sources virtuelles correspondant aux sommets via les informations telles que les directions des rayons des sommets du quadrilatère de front d'onde, le temps de propagation sismique et la vitesse sismique, comme c'est représenté sur la figure 3: supposer que les quatre sommets d'un quadrilatère de front d'onde sont respectivement A, B, C et D, et que R est un point de grille dans le quadrilatère de front d'onde , après les emplacements des sources virtuelles OA, Ob, Oc et OD correspondant à A, B, 0 et D étant calculés via des informations pertinentes, l'expression du temps de propagation sismique au point R étant :
(|OaR| + |ObR|+|OcR| + |OdR|) Tr= 4*Vr |OaR|, |ObR|, |OcR| et |OdR| représentant la distance de OA, OB, Oc et OD au point R, et VR représentant la vitesse sismique au point R ; et
6) terminer les calculs du temps de propagation de tous les points de la grille et produire les résultats du calcul du temps de propagation final.
La précision de calcul du procédé divulgué par l'invention est analysée et vérifiée ci-dessous par un modèle homogène.
La figure 4 et la figure 5 sont des erreurs absolues d'un procédé de construction de front d'onde classique et d'un procédé de construction de front d'onde de sources virtuelles dans un modèle à milieu homogène, respectivement. Le nombre de
BE2018/5525 points de grille horizontaux du modèle à milieu homogène est 761, le nombre de points de grille verticaux est 777, l'espacement de grille horizontal et l'espacement de grille vertical sont 10m, la vitesse est 1000m/s et une source sismique est positionnée à 3800m. Comme on peut le voir sur les figures, par comparaison au procédé de construction de front d'onde classique , la précision de calcul du procédé de calcul de temps de propagation sismique basé sur la construction de front d'onde de sources virtuelles est grandement améliorée.
Le temps de propagation des points de grille dans les quadrilatères de front d'onde de l'invention est calculé par le procédé des sources virtuelles, la précision de calcul du temps de propagation des points de grille dans les quadrilatères du front d'onde est améliorée et un procédé de calcul sismique basé 15 sur la construction de front d'onde est réalisé.
Claims (2)
- Revendications1. Procédé de calcul du temps de propagation sismique 2-D sur la base de la construction de front d'onde à sources virtuelles, caractérisé en ce qu'il comprend les étapes suivantes :étape 1 : lire des fichiers de paramètres pertinents et un modèle de vitesse, les fichiers de paramètres comprenant le nombre de points de grille, l'espacement de grille, l'emplacement des sources, la longueur de pas de tracé, la plage d'angles des rayons et étape 2 : tracer des rayons, et insérer de nouveaux rayons pendant le traçage pour assurer la couverture des rayons, le but du traçage de rayons étant de résoudre un groupe d'équations de traçage de rayons cinématiques à l'ide d'un procédé de RungeKutta, comme indiqué dans la formule ci-dessous :dx. , ~T = V άτ dp. d ( 1 L —- = — = -V dr dxt\v ) dv5χί où Xi représente la composante de localisation, p± représente la composante de lenteur, τ représente le temps de propagation sismique, et v représente la vitesse sismique ;étape 3 : diviser un espace de calcul en une pluralité de quadrilatères de front d'onde par calcul de l'information de position spatiale de points discrets sur les rayons acquis ;étape 4 : évaluer la relation de position du point de grille et du quadrilatère de front d'onde par un procédé à produit vectoriel, et déterminer le point de grille contenu dans le quadrilatère de front d'onde ;étape 5 : calculer les emplacements de sources virtuelles correspondant aux points via les informations pertinentes des sommets du quadrilatère de front d'onde, et calculer le temps de propagation sismique des points de grille contenus dans le quadrilatère de front d'onde en fonction sources virtuelles ; et étape 6 : terminer les calculs du temps les points de la grille et produire lesBE2018/5525 des localrsatrons des de propagation de tous résultats du calcul du temps de propagation final.
- 2. Procédé de calcul de temps de propagation sismique basé sur la construction de front d'onde 2-D à sources virtuelles selon la revendication 1, caractérisé en ce que, à l'étape 5, lorsque le temps de propagation sismique des points de grille contenus dans le quadrilatère de front d'onde est calculé par un procédé à sources virtuelles, les quatre sommets d'un quadrilatère de front d'onde sont supposés être respectivement A, B, C et D, et R est un point de grille dans le quadrilatère de front d'onde, après les emplacements des sources virtuelles OA, OB, Oc et OD correspondant à A, B, C et D sont calculés via des informations pertinentes, l'expression du temps de propagation sismique au point R étant :(|OaR| + |ObR|+|OcR| + |OdR|) Tr= 4*Vr où |OaR|, |ObR|, |OcR| et |OdR| représentent la distance de OA, OB, Oc et OD au point R, et VR représente la vitesse sismique au point R.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810641651.7A CN108957538A (zh) | 2018-06-21 | 2018-06-21 | 一种虚拟震源二维波前构建地震波走时计算方法 |
CN2018106416517 | 2018-06-21 |
Publications (2)
Publication Number | Publication Date |
---|---|
BE1025828A1 true BE1025828A1 (fr) | 2019-07-18 |
BE1025828B1 BE1025828B1 (fr) | 2019-07-25 |
Family
ID=63141839
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
BE2018/5525A BE1025828B1 (fr) | 2018-06-21 | 2018-07-19 | Procédé de calcul du temps de propagation sismique 2-D sur la base d'une construction de front d'onde à sources virtuelles |
Country Status (5)
Country | Link |
---|---|
CN (1) | CN108957538A (fr) |
BE (1) | BE1025828B1 (fr) |
LU (1) | LU100878B1 (fr) |
NL (1) | NL2021354B1 (fr) |
WO (1) | WO2019242045A1 (fr) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108957538A (zh) * | 2018-06-21 | 2018-12-07 | 成都启泰智联信息科技有限公司 | 一种虚拟震源二维波前构建地震波走时计算方法 |
CN110660135A (zh) * | 2019-09-20 | 2020-01-07 | 西南石油大学 | 一种利用三角网格化射线实现波前构建的方法 |
CN110568496B (zh) * | 2019-09-26 | 2021-02-09 | 核工业北京地质研究院 | 一种复杂介质条件下射线追踪方法 |
CN111257939B (zh) * | 2020-03-26 | 2021-06-01 | 中国石油大学(北京) | 一种时移地震虚拟震源双向波场重构方法和系统 |
CN114924312B (zh) * | 2022-05-10 | 2024-03-12 | 吉林大学 | 一种基于波前初始化射线追踪技术的高斯波束偏移方法及装置 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6035256A (en) * | 1997-08-22 | 2000-03-07 | Western Atlas International, Inc. | Method for extrapolating traveltimes across shadow zones |
EP1423735B1 (fr) * | 2001-09-07 | 2010-08-18 | Shell Internationale Research Maatschappij B.V. | Representation en image d'une formation sous la surface au moyen de sources virtuelles |
US20170363759A1 (en) * | 2016-06-17 | 2017-12-21 | Cgg Services Sa | System and method for seismic interferometry optimized data acquisition |
CN105974470B (zh) * | 2016-07-04 | 2017-06-16 | 中国石油大学(华东) | 一种多分量地震资料最小二乘逆时偏移成像方法及系统 |
CN108072897A (zh) * | 2018-01-23 | 2018-05-25 | 西南交通大学 | 一种混合二维地震波走时计算方法 |
CN108957538A (zh) * | 2018-06-21 | 2018-12-07 | 成都启泰智联信息科技有限公司 | 一种虚拟震源二维波前构建地震波走时计算方法 |
-
2018
- 2018-06-21 CN CN201810641651.7A patent/CN108957538A/zh active Pending
- 2018-07-04 WO PCT/CN2018/094459 patent/WO2019242045A1/fr active Application Filing
- 2018-07-04 LU LU100878A patent/LU100878B1/fr active IP Right Grant
- 2018-07-19 BE BE2018/5525A patent/BE1025828B1/fr not_active IP Right Cessation
- 2018-07-19 NL NL2021354A patent/NL2021354B1/en not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
NL2021354A (en) | 2018-08-14 |
WO2019242045A1 (fr) | 2019-12-26 |
CN108957538A (zh) | 2018-12-07 |
LU100878B1 (fr) | 2019-12-30 |
BE1025828B1 (fr) | 2019-07-25 |
NL2021354B1 (en) | 2019-04-26 |
WO2019242045A9 (fr) | 2021-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
BE1025828B1 (fr) | Procédé de calcul du temps de propagation sismique 2-D sur la base d'une construction de front d'onde à sources virtuelles | |
CN104977607B (zh) | 利用变步长网格声波波场模拟的时间域全波形反演方法 | |
CN109100816B (zh) | 一种重磁数据处理方法及系统 | |
CN109116413B (zh) | 成像域立体层析速度反演方法 | |
CN111220222B (zh) | 一种超声波燃气表流量的测定算法 | |
EP4127457A1 (fr) | Système et procédé de mesure de turbulence de flux de vent par lidar dans un terrain complexe | |
CA2883668A1 (fr) | Systeme et procede permettant d'analyser des concepts de releves sismiques | |
LU100751B1 (fr) | Méthode de calcul de temps de propagation sismique mixte en 2D | |
CN108957553A (zh) | 动校正量递推修正的无拉伸畸变动校正方法及装置 | |
CN105572741B (zh) | 一种计算3d高频静校正量的方法 | |
AU739128B2 (en) | A method of seismic processing, and in particular a 3D seismic prospection method implementing seismic data migration | |
Hart | Sphere tracing: Simple robust antialiased rendering of distance-based implicit surfaces | |
JP2007327898A (ja) | ガウシアンレイバンドルモデルにおける音場エネルギーの広がりの標準偏差算出方法及びそのプログラム | |
CN106932821A (zh) | 地震层析反演中的一种目标射线追踪技术 | |
CN113866821A (zh) | 一种基于照明方向约束的被动源干涉偏移成像方法和系统 | |
Paxman et al. | Optimal simulation of volume turbulence with phase screens | |
Zhang et al. | Experimental study of propeller and ducted propeller’s wake inducted scour | |
Ravve et al. | Eigenray tracing in 3D heterogeneous media using spectral element method | |
Kowalski et al. | Unbiased conversion of passive sensor measurements using closest point of approach | |
Ozendi et al. | Stochastic surface mesh reconstruction | |
Elek et al. | Progressive Spectral Ray Differentials. | |
FR2683345A1 (fr) | Procede de determination de la trajectoire d'un rai a travers un milieu geologique tridimensionnel. | |
Hovad et al. | Slope based Grid Creation using Interpolation of LIDAR Data Sets. | |
CN115184986A (zh) | 不依赖震源的全局包络互相关全波形反演方法 | |
Maleika | The Effect of the Density of Measurement Points Collected from a Multibeam Echosounder on the Selection of IDW Interpolation Points in the Process of Creating Seabed Models. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FG | Patent granted |
Effective date: 20190725 |
|
MM | Lapsed because of non-payment of the annual fee |
Effective date: 20200731 |