Seismic inversion holds significant importance across various domains of geoscience and engineering, including the characterization of energy resource reservoirs, the assessment of polluted sites, and CO2 storage. It is a process of estimating rock properties from seismic data that is inherently uncertain, nonlinear, non‐unique, and highly challenging. Using multiple seismic attributes increases the size of the data, requiring considerable processing resources and time. However, deep learning can accurately fit quantities of nonlinear variables, making it an excellent method for predicting spatially distributed subsurface properties. We trained some multi‐output regression neural networks to carry out porosity inversion from seismic data. We initially computed a series of seismic attributes and generated the corresponding porosity using interpreted horizons, well logs, and seismic data. Subsequently, we proposed a technique to identify the most relevant seismic attributes for porosity inversion. Because our networks work as stochastic modeling entities, we created a weight‐averaging ensemble approach to build a strong model with the highest level of accuracy. We combined realizations from baseline entities, considering their respective performance levels. Using the statistics between these realizations and the robust model, we determined the degree of uncertainty associated with the outcome. We found an R2 of 0.993 and an MAE of 0.00112 in the F3 block offshore the Netherlands, proving the method's effectiveness. The mean porosity was 0.175193, compared to 0.175626 from a reference model, and the mean uncertainty was ±0.0008998.