A computer-implemented method for imaging a three dimensional geological region of interest where gas-bearing sand may be found, relative to a predetermined geological structure comprising at least one of limestone or dolomite, the method comprising the steps of: defining a boundary of the geological structure; defining the region of interest relative to the boundary of the geological structure; selecting from a GUI interface a first attribute representing a first physical property of an acoustic signal and a second attribute representing a second physical property of an acoustic signal, the first attribute and second attribute representing the region of interest; calculating a first attribute volume and a second attribute volume for the region of interest, the first attribute volume and the second attribute volume each comprising a plurality of voxels, each voxel being defined by a set of x,y,z coordinates and a data value associated with the first and second physical property, respectively; selecting a first set of voxels from the first attribute volume that have a data value within a first attribute data value range, the first set of voxels representing at least one of sand or sandstone of the region of interest; selecting a second set of voxels from the second attribute volume that have a data value within a second attribute data value range, the second set of voxels representing the at least one of sand or sandstone; and imaging the first set of voxels and the second set of voxels; wherein the boundary is defined by: selecting from a GUI interface a third attribute representing a third physical property of an acoustic signal, the third attribute representing the geological structure; calculating a third attribute volume (402, 404) for the geological structure, the third attribute volume (402, 404) comprising a plurality of voxels, each voxel being defined by a set of x,y,z coordinates and a data value associated with the third physical property; selecting a third set of voxels from the third attribute volume (402, 404) that have a data value within a third attribute data value range, the third set of voxels representing the boundary of the geological structure; and imaging the third set of voxels. defining a boundary of the geological structure; defining the region of interest relative to the boundary of the geological structure; selecting from a GUI interface a first attribute representing a first physical property of an acoustic signal and a second attribute representing a second physical property of an acoustic signal, the first attribute and second attribute representing the region of interest; calculating a first attribute volume and a second attribute volume for the region of interest, the first attribute volume and the second attribute volume each comprising a plurality of voxels, each voxel being defined by a set of x,y,z coordinates and a data value associated with the first and second physical property, respectively; selecting a first set of voxels from the first attribute volume that have a data value within a first attribute data value range, the first set of voxels representing at least one of sand or sandstone of the region of interest; selecting a second set of voxels from the second attribute volume that have a data value within a second attribute data value range, the second set of voxels representing the at least one of sand or sandstone; and imaging the first set of voxels and the second set of voxels; wherein the boundary is defined by: selecting from a GUI interface a third attribute representing a third physical property of an acoustic signal, the third attribute representing the geological structure; calculating a third attribute volume (402, 404) for the geological structure, the third attribute volume (402, 404) comprising a plurality of voxels, each voxel being defined by a set of x,y,z coordinates and a data value associated with the third physical property; selecting a third set of voxels from the third attribute volume (402, 404) that have a data value within a third attribute data value range, the third set of voxels representing the boundary of the geological structure; and imaging the third set of voxels. selecting from a GUI interface a third attribute representing a third physical property of an acoustic signal, the third attribute representing the geological structure; calculating a third attribute volume (402, 404) for the geological structure, the third attribute volume (402, 404) comprising a plurality of voxels, each voxel being defined by a set of x,y,z coordinates and a data value associated with the third physical property; selecting a third set of voxels from the third attribute volume (402, 404) that have a data value within a third attribute data value range, the third set of voxels representing the boundary of the geological structure; and imaging the third set of voxels.