change properties in faults

[Denisiuskley]

Hi guys. You have created an amazing program for synthetic seismic generation. However, I would like to add property changes in the faults. In which part of the code would it be best to implement this, and which data array should be used to determine the fault location?

[Denisiuskley]

I mean physical properties, of course (vp, vs, rho).

[tpmerrifield]

Hi Denis,

Thanks!

I assume you want to specify some vp, vs, rho values in the fault zones, is that right?
My suggestion:

• The 3D fault plane volume is stored as a binary segmentation as fault_planes.

• You may want to perform some morphological dilation to widen the fault zone, for example the _dilate_faults method in Closures.py.

• A very simple method would be to do something similar to how we overwrite the property values for salt. Check for those voxels which are in a fault zone and enter the values you want to use.

• If you wanted to do something more complicated, which accounted for changes with depth / per layer, then you could add a new block to apply_faulting_to_geomodels_and_depth_maps, checking whether voxels in the layer are in a faulted zone, similar to how the shale and sand properties are created.

  • I'd do this after all of the other lithologies have been generated, so that they don't get overwritten. Also take care that they don't get modified by the scaling with the random scaling factors. It might be simplest to modify apply_scaling_factors to ignore the voxels in the faulted zones.

Best,
Tom

[Denisiuskley]

Hi, Tom. Thanks for your reply. My idea is simply to multiply each property by a random coefficient, varying, for example, from 0.5 to 1. Your idea about modifying properties, as it's done for salt bodies, is brilliant. Thank you for your hard work! You're making the world a better!