|
5 | 5 | <head> |
6 | 6 | <meta charset="utf-8" /> |
7 | 7 | <meta name="viewport" content="width=device-width, initial-scale=1.0" /> |
8 | | - <title>fmu.tools.domainconversion.dconvert — fmu-tools 1.32.1.dev1+g79dcf8305 documentation</title> |
| 8 | + <title>fmu.tools.domainconversion.dconvert — fmu-tools 1.32.2.dev1+g858860831 documentation</title> |
9 | 9 | <link rel="stylesheet" type="text/css" href="../../../../_static/pygments.css?v=03e43079" /> |
10 | 10 | <link rel="stylesheet" type="text/css" href="../../../../_static/css/theme.css?v=9edc463e" /> |
11 | 11 |
|
12 | 12 |
|
13 | 13 | <script src="../../../../_static/jquery.js?v=5d32c60e"></script> |
14 | 14 | <script src="../../../../_static/_sphinx_javascript_frameworks_compat.js?v=2cd50e6c"></script> |
15 | | - <script src="../../../../_static/documentation_options.js?v=ac42e60e"></script> |
| 15 | + <script src="../../../../_static/documentation_options.js?v=aa3c8e2d"></script> |
16 | 16 | <script src="../../../../_static/doctools.js?v=9bcbadda"></script> |
17 | 17 | <script src="../../../../_static/sphinx_highlight.js?v=dc90522c"></script> |
18 | 18 | <script src="../../../../_static/js/theme.js"></script> |
@@ -154,6 +154,7 @@ <h1>Source code for fmu.tools.domainconversion.dconvert</h1><div class="highligh |
154 | 154 | <span class="n">_scube_d</span><span class="p">:</span> <span class="n">xtgeo</span><span class="o">.</span><span class="n">Cube</span> <span class="o">=</span> <span class="n">field</span><span class="p">(</span><span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span> <span class="n">init</span><span class="o">=</span><span class="kc">False</span><span class="p">)</span> <span class="c1"># avg slowness cube in depth</span> |
155 | 155 |
|
156 | 156 | <span class="n">_v_surfaces</span><span class="p">:</span> <span class="nb">list</span><span class="p">[</span><span class="n">xtgeo</span><span class="o">.</span><span class="n">RegularSurface</span><span class="p">]</span> <span class="o">=</span> <span class="n">field</span><span class="p">(</span><span class="n">default_factory</span><span class="o">=</span><span class="nb">list</span><span class="p">,</span> <span class="n">init</span><span class="o">=</span><span class="kc">False</span><span class="p">)</span> |
| 157 | + <span class="n">_vi_surfaces</span><span class="p">:</span> <span class="nb">list</span><span class="p">[</span><span class="n">xtgeo</span><span class="o">.</span><span class="n">RegularSurface</span><span class="p">]</span> <span class="o">=</span> <span class="n">field</span><span class="p">(</span><span class="n">default_factory</span><span class="o">=</span><span class="nb">list</span><span class="p">,</span> <span class="n">init</span><span class="o">=</span><span class="kc">False</span><span class="p">)</span> |
157 | 158 | <span class="n">_s_surfaces</span><span class="p">:</span> <span class="nb">list</span><span class="p">[</span><span class="n">xtgeo</span><span class="o">.</span><span class="n">RegularSurface</span><span class="p">]</span> <span class="o">=</span> <span class="n">field</span><span class="p">(</span><span class="n">default_factory</span><span class="o">=</span><span class="nb">list</span><span class="p">,</span> <span class="n">init</span><span class="o">=</span><span class="kc">False</span><span class="p">)</span> |
158 | 159 | <span class="n">_d_surfaces</span><span class="p">:</span> <span class="nb">list</span><span class="p">[</span><span class="n">xtgeo</span><span class="o">.</span><span class="n">RegularSurface</span><span class="p">]</span> <span class="o">=</span> <span class="n">field</span><span class="p">(</span><span class="n">default_factory</span><span class="o">=</span><span class="nb">list</span><span class="p">,</span> <span class="n">init</span><span class="o">=</span><span class="kc">False</span><span class="p">)</span> |
159 | 160 | <span class="n">_t_surfaces</span><span class="p">:</span> <span class="nb">list</span><span class="p">[</span><span class="n">xtgeo</span><span class="o">.</span><span class="n">RegularSurface</span><span class="p">]</span> <span class="o">=</span> <span class="n">field</span><span class="p">(</span><span class="n">default_factory</span><span class="o">=</span><span class="nb">list</span><span class="p">,</span> <span class="n">init</span><span class="o">=</span><span class="kc">False</span><span class="p">)</span> |
@@ -215,6 +216,7 @@ <h1>Source code for fmu.tools.domainconversion.dconvert</h1><div class="highligh |
215 | 216 | <span class="bp">self</span><span class="o">.</span><span class="n">time_surfaces</span><span class="p">,</span> <span class="n">template_surf</span><span class="p">,</span> <span class="n">fill</span><span class="o">=</span><span class="kc">True</span><span class="p">,</span> <span class="n">ensure_consistency</span><span class="o">=</span><span class="kc">True</span> |
216 | 217 | <span class="p">)</span> |
217 | 218 | <span class="bp">self</span><span class="o">.</span><span class="n">_ensure_surfaces_has_msl</span><span class="p">()</span> |
| 219 | + <span class="bp">self</span><span class="o">.</span><span class="n">_calculate_interval_velocity_surfaces</span><span class="p">()</span> |
218 | 220 | <span class="n">_logger</span><span class="o">.</span><span class="n">debug</span><span class="p">(</span><span class="s2">"Check and fix surfaces... DONE"</span><span class="p">)</span> |
219 | 221 |
|
220 | 222 | <span class="nd">@staticmethod</span> |
@@ -321,6 +323,32 @@ <h1>Source code for fmu.tools.domainconversion.dconvert</h1><div class="highligh |
321 | 323 | <span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">_names</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">!=</span> <span class="s2">"__MSL"</span><span class="p">:</span> |
322 | 324 | <span class="bp">self</span><span class="o">.</span><span class="n">_names</span><span class="o">.</span><span class="n">insert</span><span class="p">(</span><span class="mi">0</span><span class="p">,</span> <span class="s2">"__MSL"</span><span class="p">)</span> |
323 | 325 |
|
| 326 | + <span class="k">def</span><span class="w"> </span><span class="nf">_calculate_interval_velocity_surfaces</span><span class="p">(</span><span class="bp">self</span><span class="p">)</span> <span class="o">-></span> <span class="kc">None</span><span class="p">:</span> |
| 327 | +<span class="w"> </span><span class="sd">"""Calculate interval velocity surfaces."""</span> |
| 328 | + <span class="n">_logger</span><span class="o">.</span><span class="n">debug</span><span class="p">(</span><span class="s2">"Calculate interval velocity surfaces"</span><span class="p">)</span> |
| 329 | + |
| 330 | + <span class="n">vel</span> <span class="o">=</span> <span class="p">[]</span> |
| 331 | + <span class="k">for</span> <span class="n">no</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span> <span class="nb">len</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">_d_surfaces</span><span class="p">)):</span> |
| 332 | + <span class="n">t0</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">_t_surfaces</span><span class="p">[</span><span class="n">no</span> <span class="o">-</span> <span class="mi">1</span><span class="p">]</span> |
| 333 | + <span class="n">t1</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">_t_surfaces</span><span class="p">[</span><span class="n">no</span><span class="p">]</span> |
| 334 | + <span class="n">d0</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">_d_surfaces</span><span class="p">[</span><span class="n">no</span> <span class="o">-</span> <span class="mi">1</span><span class="p">]</span> |
| 335 | + <span class="n">d1</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">_d_surfaces</span><span class="p">[</span><span class="n">no</span><span class="p">]</span> |
| 336 | + |
| 337 | + <span class="n">vint</span> <span class="o">=</span> <span class="n">d1</span><span class="o">.</span><span class="n">copy</span><span class="p">()</span> |
| 338 | + <span class="n">tdiff</span> <span class="o">=</span> <span class="n">t1</span><span class="o">.</span><span class="n">values</span> <span class="o">-</span> <span class="n">t0</span><span class="o">.</span><span class="n">values</span> |
| 339 | + <span class="n">tdiff</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">where</span><span class="p">(</span><span class="n">tdiff</span> <span class="o">==</span> <span class="mf">0.0</span><span class="p">,</span> <span class="mf">1e-06</span><span class="p">,</span> <span class="n">tdiff</span><span class="p">)</span> |
| 340 | + <span class="n">vint</span><span class="o">.</span><span class="n">values</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">divide</span><span class="p">((</span><span class="n">d1</span><span class="o">.</span><span class="n">values</span> <span class="o">-</span> <span class="n">d0</span><span class="o">.</span><span class="n">values</span><span class="p">),</span> <span class="n">tdiff</span><span class="p">)</span> |
| 341 | + <span class="n">vint</span><span class="o">.</span><span class="n">values</span> <span class="o">*=</span> <span class="mi">2000</span> |
| 342 | + <span class="n">vel</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">vint</span><span class="p">)</span> |
| 343 | + |
| 344 | + <span class="n">d1_values</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">_d_surfaces</span><span class="p">[</span><span class="mi">1</span><span class="p">]</span><span class="o">.</span><span class="n">values</span> |
| 345 | + <span class="n">t1_values</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">_t_surfaces</span><span class="p">[</span><span class="mi">1</span><span class="p">]</span><span class="o">.</span><span class="n">values</span> |
| 346 | + <span class="k">if</span> <span class="n">np</span><span class="o">.</span><span class="n">allclose</span><span class="p">(</span><span class="n">d1_values</span><span class="p">,</span> <span class="mf">0.0</span><span class="p">)</span> <span class="ow">and</span> <span class="n">np</span><span class="o">.</span><span class="n">allclose</span><span class="p">(</span><span class="n">t1_values</span><span class="p">,</span> <span class="mf">0.0</span><span class="p">):</span> |
| 347 | + <span class="n">vel</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">=</span> <span class="n">vel</span><span class="p">[</span><span class="mi">1</span><span class="p">]</span> |
| 348 | + |
| 349 | + <span class="n">vel</span><span class="o">.</span><span class="n">insert</span><span class="p">(</span><span class="mi">0</span><span class="p">,</span> <span class="n">vel</span><span class="p">[</span><span class="mi">0</span><span class="p">])</span> |
| 350 | + <span class="bp">self</span><span class="o">.</span><span class="n">_vi_surfaces</span> <span class="o">=</span> <span class="n">vel</span> |
| 351 | + |
324 | 352 | <span class="k">def</span><span class="w"> </span><span class="nf">_velo_maps_average</span><span class="p">(</span><span class="bp">self</span><span class="p">)</span> <span class="o">-></span> <span class="kc">None</span><span class="p">:</span> |
325 | 353 | <span class="w"> </span><span class="sd">"""Create average velocities from MSL to surface N"""</span> |
326 | 354 |
|
@@ -812,6 +840,14 @@ <h1>Source code for fmu.tools.domainconversion.dconvert</h1><div class="highligh |
812 | 840 | <span class="k">yield</span> <span class="n">surf</span></div> |
813 | 841 |
|
814 | 842 |
|
| 843 | +<div class="viewcode-block" id="DomainConversion.vint_surfaces"> |
| 844 | +<a class="viewcode-back" href="../../../../fmu.tools.domainconversion.html#fmu.tools.DomainConversion.vint_surfaces">[docs]</a> |
| 845 | + <span class="k">def</span><span class="w"> </span><span class="nf">vint_surfaces</span><span class="p">(</span><span class="bp">self</span><span class="p">)</span> <span class="o">-></span> <span class="n">Generator</span><span class="p">[</span><span class="n">xtgeo</span><span class="o">.</span><span class="n">RegularSurface</span><span class="p">]:</span> |
| 846 | +<span class="w"> </span><span class="sd">"""Return a generator of interval velocity surfaces."""</span> |
| 847 | + <span class="k">for</span> <span class="n">surf</span> <span class="ow">in</span> <span class="bp">self</span><span class="o">.</span><span class="n">_vi_surfaces</span><span class="p">:</span> |
| 848 | + <span class="k">yield</span> <span class="n">surf</span></div> |
| 849 | + |
| 850 | + |
815 | 851 | <span class="nd">@property</span> |
816 | 852 | <span class="k">def</span><span class="w"> </span><span class="nf">average_velocity_cube_in_time</span><span class="p">(</span><span class="bp">self</span><span class="p">)</span> <span class="o">-></span> <span class="n">xtgeo</span><span class="o">.</span><span class="n">Cube</span> <span class="o">|</span> <span class="kc">None</span><span class="p">:</span> |
817 | 853 | <span class="k">return</span> <span class="bp">self</span><span class="o">.</span><span class="n">_vcube_t</span> |
|
0 commit comments