Subluminal values of the speed of gravitational waves (GW) were obtained that reproduce a hypothetical Hulse-Taylor binary pulsar undergoing circular orbit decay. Those values of speed were used to simulate, in the framework of gravitoelectromagnetism (GEM), the in-spiral process of 3 GW events. The calculated results show a significantly better agreement (with the results of the linear theory of relativity) than the ones obtained using the speed of light. A method was proposed to measure the speed of the GW. Constraints of 0.89c < cg < 1.11c and 0.98c < cg < 1.20c were obtained for GW170814 and GW170817 events respectively. Some assumptions for the extremely high constraints on the deviation of the GW speed from the light speed of the almost simultaneous GW170817/GRB170817A events were referenced and alternative scenarios presented. The need for model-independent GW arrival time delays at the detectors is remarked.